#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const ll CYCLE_PER_SEC = 2400000000; double TIME_LIMIT = 2.8; unsigned long long int get_cycle() { unsigned int low, high; __asm__ volatile ("rdtsc" : "=a" (low), "=d" (high)); return ((unsigned long long int) low) | ((unsigned long long int) high << 32); } double get_time(unsigned long long int begin_cycle) { return (double) (get_cycle() - begin_cycle) / CYCLE_PER_SEC; } 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()); ll start_cycle = get_cycle(); int play_cnt = 0; int ret[3] = {0, 0, 0}; while (get_time(start_cycle) < TIME_LIMIT) { ++play_cnt; int s = simulate(PA * 1000, PB * 1000); ret[s]++; } cout << fixed << setprecision(10) << ret[1] * 1.0 / play_cnt << endl; return 0; }