#include <algorithm>
#include<bits/stdc++.h>

using namespace std;

using db = long double;
constexpr db eps = 1e-50;

void solve() {
	int n, a, b; cin >> n >> a >> b;
	vector<int>vec(n);
	for(int i = 0 ; i < n ; i ++) cin >> vec[i];

	map<int, int>cnt;
	for(int i = 0 ; i < n ; i ++) cnt[vec[i]]++;

	vector<vector<db>>dp(n+1, vector<db>(b+1));
	auto ndp = dp;

	dp[0][0] = 1;
	for(int i = 0 ; i < n ; i ++) {	
		for(auto &v: ndp) fill(begin(v), end(v), 0);
		for(int x = 0 ; x <= i ; x ++) {
			for(int y = 0 ; y <= b ; y ++) {
				ndp[x][y] += dp[x][y];
			}
			for(int y = 0 ; y <= b-vec[i] ; y ++) {
				ndp[x+1][y + vec[i]] += dp[x][y] * (x+1) / (n-x);
			}
		}
		dp = ndp;
	}

	db ans = 0;

	for(auto &[v, t]: cnt) {
		if(t<=100) {
			ndp = dp;
			for(int x = 1 ; x <= n ; x ++) {
				for(int y = v ; y <= b ; y ++) {
					auto &val = ndp[x-1][y-v];
					if(val<eps) val = 0;
					if(val>1-eps) val = 1;
					ndp[x][y] -= val * x / (n-x+1);
				}
			}

			db temp = 0;
			for(int x = 0 ; x < n ; x ++) {
				db t1 = 0;
				for(int y = max(0, a-v+1) ; y <= min(a, b-v) ; y ++) {
					t1 += ndp[x][y];
				}
				temp += t1/(n-x);
			}
			ans += temp * t;
		}
		else{
			vector<vector<db>>dp(n+1, vector<db>(b+1));
			dp[0][0] = 1;
			bool flag = false;
			for(int i = 0 ; i < n ; i ++) {	
				if(!flag and vec[i]==v) {
					flag = true;
					continue;
				}
				ndp = vector<vector<db>>(n+1, vector<db>(b+1, 0));
				for(int x = 0 ; x <= i ; x ++) {
					for(int y = 0 ; y <= b ; y ++) {
						ndp[x][y] += dp[x][y];
					}
					for(int y = 0 ; y <= b-vec[i] ; y ++) {
						ndp[x+1][y + vec[i]] += dp[x][y] * (x+1) / (n-x);
					}
				}
				dp = ndp;
			}
			db temp = 0;
			for(int x = 0 ; x < n ; x ++) {
				db t1 = 0;
				for(int y = max(0, a-v+1) ; y <= min(a, b-v) ; y ++) {
					t1 += dp[x][y];
				}
				temp += t1/(n-x);
			}
			ans += temp * t;
		}
	}

	cout << (long double)ans << endl;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cout << fixed << setprecision(20);

	solve();

	return 0;
}