#include using namespace std; int main(){ int n; double pa, pb; scanf("%d %lf %lf", &n, &pa, &pb); vector a(n), b(n); vector> A(n, vector(n)), B(n, vector(n)); for(auto &&v: a)cin >> v; for(auto &&v: b)cin >> v; sort(a.begin(), a.end()); sort(b.begin(), b.end()); vector dp1(1 << n), dp2(1 << n); dp1[0] = dp2[0] = 1; for(int i = 0; i < (1 << n); i++){ int cnt = __builtin_popcount(i); for(int j = 0, flag = 1; j < n; j++){ if(~i >> j & 1){ if(flag){ if(cnt + 1 == n){ A[j][cnt] += dp1[i]; B[j][cnt] += dp2[i]; continue; } A[j][cnt] += dp1[i] * pa; dp1[i | (1 << j)] += dp1[i] * pa; B[j][cnt] += dp2[i] * pb; dp2[i | (1 << j)] += dp2[i] * pb; flag = 0; }else{ A[j][cnt] += dp1[i] * (1 - pa) / (n - cnt - 1); dp1[i | (1 << j)] += dp1[i] * (1 - pa) / (n - cnt - 1); B[j][cnt] += dp2[i] * (1 - pb) / (n - cnt - 1); dp2[i | (1 << j)] += dp2[i] * (1 - pb) / (n - cnt - 1); } } } } double ans = 0; for(int i = 0; i < n; i++){ for(int j = 0; j < n; j++){ if(a[i] < b[j])continue; for(int k = 0; k < n; k++){ ans += (a[i] + b[j]) * A[i][k] * B[j][k]; } } } printf("%.15lf\n", ans); }