結果

問題 No.174 カードゲーム（Hard）
コンテスト
ユーザー t98slider
提出日時 2022-07-15 05:13:28
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 141 ms / 2,000 ms
コード長 1,661 bytes
コンパイル時間 1,750 ms
コンパイル使用メモリ 177,164 KB
実行使用メモリ 19,956 KB
最終ジャッジ日時 2024-06-26 23:42:12
合計ジャッジ時間 3,723 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

int main(){
    int n;
    double pa, pb;
    scanf("%d %lf %lf", &n, &pa, &pb);
    vector<int> a(n), b(n);
    vector<vector<double>> A(n, vector<double>(n)), B(n, vector<double>(n));
    for(auto &&v: a)cin >> v;
    for(auto &&v: b)cin >> v;
    sort(a.begin(), a.end());
    sort(b.begin(), b.end());
    vector<double> 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);
}
0