#include #define rep(i,a,b) for(int i=a;i=b;i--) #define fore(i,a) for(auto &i:a) #define all(x) (x).begin(),(x).end() #pragma GCC optimize ("-O3") using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); } typedef long long ll; const int inf = INT_MAX / 2; const ll infl = 1LL << 60; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b> N; rep(j, 0, 2) cin >> P[j]; rep(j, 0, 2) { rep(i, 0, N) cin >> A[j][i]; sort(A[j], A[j] + N); } rep(j, 0, 2) { vector dp(1 << N); dp[0] = 1; rep(msk, 0, 1 << N) { int cnt = 0, f = 1; rep(i, 0, N) if (!(msk & (1 << i))) cnt++; rep(i, 0, N) if (!(msk & (1 << i))) { double pp = 0; if (cnt == 1) { pp = dp[msk]; } else { if (f) pp = dp[msk] * P[j], f = 0; else pp = dp[msk] * (1.0 - P[j]) / (cnt - 1); } v[j][N - cnt][i] += pp; dp[msk + (1 << i)] += pp; } } } double ans = 0; rep(i, 0, N) { rep(a, 0, N) rep(b, 0, N) if (A[0][a] > A[1][b]) ans += v[0][i][a] * v[1][i][b] * (A[0][a] + A[1][b]); } printf("%.10f\n", ans); }