#include using namespace std; #define REP(i,a,b) for(int i=a;i<(int)b;i++) #define rep(i,n) REP(i,0,n) int N, T; double pa,pb; int A[22], B[22]; int win, sum; vector getrandvec(double p, default_random_engine& gene, uniform_real_distribution& dist) { int x; vector v,r; rep(i, N) { v.push_back(i); } rep(i, N-1) { double d = dist(gene); if(d < p) x = 0; else x = 1+(N-1-i)*(d-p) / (1-p); r.push_back(v[x]); v.erase(v.begin()+x); } r.push_back(v[0]); return r; } int main() { cin >> N; double pa, pb; cin>>pa>>pb; rep(i, N) cin >> A[i], T += A[i]; rep(i, N) cin >> B[i], T += B[i]; sort(A, A+N); sort(B, B+N); default_random_engine gene; uniform_real_distribution dist(0.,1.); int x; rep(i, 300000) { sum ++; vector va, vb; va = getrandvec(pa, gene, dist); vb = getrandvec(pb, gene, dist); x = 0; rep(y, N) if(A[va[y]] > B[vb[y]]) { x += A[va[y]] + B[vb[y]]; } if(x > T/2) { win ++; } } printf("%.10f\n", (double)win/sum); return 0; }