#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; unsigned long long xor128() { static unsigned long long rx = 123456789, ry = 362436069, rz = 521288629, rw = 88675123; unsigned long long rt = (rx ^ (rx << 11)); rx = ry; ry = rz; rz = rw; return (rw = (rw ^ (rw >> 19)) ^ (rt ^ (rt >> 8))); } int N; double PA, PB; vector A; vector B; int simulate(int select_min_rate_a, int select_min_rate_b) { bool used_a[N]; bool used_b[N]; memset(used_a, false, sizeof(used_a)); memset(used_b, false, sizeof(used_b)); int scores[2] = {0, 0}; int m_idx_a = 0; int m_idx_b = 0; for (int i = 0; i < N; ++i) { int idx_a; int idx_b; int pa = xor128() % 1000; while (used_a[m_idx_a]) { ++m_idx_a; } if (pa < select_min_rate_a || i == N - 1) { idx_a = m_idx_a; } else { do { idx_a = xor128() % N; } while (used_a[idx_a] || m_idx_a == idx_a); } int a = A[idx_a]; used_a[idx_a] = true; int pb = xor128() % 1000; while (used_b[m_idx_b]) { ++m_idx_b; } if (pb < select_min_rate_b || i == N - 1) { idx_b = m_idx_b; } else { do { idx_b = xor128() % N; } while (used_b[idx_b] || m_idx_b == idx_b); } int b = B[idx_b]; used_b[idx_b] = true; if (a > b) { scores[0] += a + b; } else { scores[1] += a + b; } } if (scores[0] > scores[1]) { return 1; } else if (scores[0] < scores[1]) { return 2; } else { return 0; } } int main() { cin >> N >> PA >> PB; A.resize(N); B.resize(N); for (int i = 0; i < N; ++i) { cin >> A[i]; } for (int i = 0; i < N; ++i) { cin >> B[i]; } sort(A.begin(), A.end()); sort(B.begin(), B.end()); int play_cnt = 5000000; int ret[3] = {0, 0, 0}; for (int i = 0; i < play_cnt; ++i) { int s = simulate(PA * 1000, PB * 1000); ret[s]++; } cout << fixed << setprecision(10) << ret[1] * 1.0 / play_cnt << endl; return 0; }