#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef pair Pii; #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define FOR(i, a, b) for (int i = (int)a; i <= (int)b; i++) template void checkmin(T &a, T b) { if (b < a) a = b; } template void checkmax(T &a, T b) { if (b > a) a = b; } const int MaxN = 20; int N; double PA, PB; int A[MaxN], B[MaxN]; double bitdp[1 << MaxN][MaxN]; double AtA[MaxN][MaxN], AtB[MaxN][MaxN]; // who, at int bitcount(int x) { int c = 0; while (x != 0) { c++; x &= x - 1; } return c; } void checkProb(double selectMinProb, double res[MaxN][MaxN]) { memset(bitdp, 0, sizeof(bitdp)); bitdp[0][0] = 1; rep(mask, 1 << N) { int num = bitcount(mask); int rem = N - num; int min = -1; double currProb = 0; rep(i, N) { currProb += bitdp[mask][i]; if (min == -1 && (mask >> i & 1) == 0) min = i; } rep(next, N) if ((mask >> next & 1) == 0) { double selectProb = (rem == 1) ? 1 : (next == min) ? selectMinProb : (1 - selectMinProb) / (rem - 1); double p = currProb * selectProb; bitdp[mask | 1 << next][next] += p; res[next][num] += p; } } } void solve() { sort(A, A + N); sort(B, B + N); checkProb(PA, AtA); checkProb(PB, AtB); double ans = 0; rep(ai, N) rep(bi, N) if (A[ai] > B[bi]) rep(pos, N) { ans += (A[ai] + B[bi]) * AtA[ai][pos] * AtB[bi][pos]; } printf("%.12f\n", ans); } int main() { cin >> N >> PA >> PB; rep(i, N) cin >> A[i]; rep(i, N) cin >> B[i]; solve(); return 0; }