#include using namespace std; #ifdef LOCAL #include "debug.hpp" #else #define debug(...) 1 #endif int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n; double pa, pb; cin >> n >> pa >> pb; vector a(n), b(n); for (int i = 0; i < n; i++) { cin >> a[i]; } for (int i = 0; i < n; i++) { cin >> b[i]; } auto solve = [&] (double p) -> vector> { vector dp(1 << n, 0.0); vector> card(n + 1, vector(n, 0.0)); dp[0] = 1.0; for (int bit = 0; bit < (1 << n); bit++) { bool mi = true; int pc = __builtin_popcount(bit); for (int i = 0; i < n; i++) { if ((bit >> i) & 1 ^ 1) { int nbit = bit | 1 << i; if (mi) { mi = false; dp[nbit] += dp[bit] * (pc == n - 1 ? 1.0 : p); card[pc + 1][i] += dp[bit] * (pc == n - 1 ? 1.0 : p); } else { dp[nbit] += dp[bit] * (1.0 - p) / (n - pc - 1); card[pc + 1][i] += dp[bit] * (1.0 - p) / (n - pc - 1); } } } } return card; }; sort(a.begin(), a.end()); sort(b.begin(), b.end()); auto cardA = solve(pa), cardB = solve(pb); debug(cardA); double ans = 0; for (int i = 1; i <= n; i++) { for (int j = 0; j < n; j++) { for (int k = 0; k < n; k++) { if (a[j] > b[k]) { ans += (double) (a[j] + b[k]) * cardA[i][j] * cardB[i][k]; } } } } cout << fixed << setprecision(15) << ans << '\n'; }