結果

問題 No.174 カードゲーム（Hard）
コンテスト
ユーザー EmKjp
提出日時 2015-03-28 10:51:43
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 311 ms / 2,000 ms
コード長 2,117 bytes
コンパイル時間 684 ms
コンパイル使用メモリ 91,656 KB
実行使用メモリ 11,972 KB
最終ジャッジ日時 2024-07-06 21:37:46
合計ジャッジ時間 3,540 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<sstream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include<cmath>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<numeric>
#include<functional>
#include<complex>

using namespace std;
#define BET(a,b,c) ((a)<=(b)&&(b)<(c))
#define FOR(i,n) for(int i=0,i##_end=(int(n));i<i##_end;i++)
#define SZ(x) (int)(x.size())
#define ALL(x) (x).begin(),(x).end()
#define MP make_pair
#define FOR_EACH(it,v) for(__typeof(v.begin()) it=v.begin(),it_end=v.end() ; it != it_end ; it++)
typedef vector<int> VI;
typedef vector<VI> VVI;

vector<vector<double> > getTable(VI A, double p){
    int n = SZ(A);
    vector<double> prob(1<<n);
    vector<vector<double> > useP(n, vector<double>(n)); 
    prob[0] = 1.0;
    for(int b=0;b<(1<<n);b++){
        if(prob[b] == 0) continue;
        int minp = 1<<28;
        int minIndex = -1;
        int rem = 0;
        FOR(i,n){
            if(b & (1<<i)) continue;
            rem++;
            if(minp > A[i]) {
                minp = A[i];
                minIndex = i;
            }
        }
        if(rem == 1){
            prob[b | (1<<minIndex)] += prob[b];
            useP[n - rem][minIndex] += prob[b];
        }else{
            FOR(i,n){
                if(b & (1<<i)) continue;
                if(i == minIndex){
                    prob[b | (1<<i)] += prob[b] * p;
                    useP[n - rem][i] += prob[b] * p;
                }else {
                    prob[b | (1<<i)] += prob[b] * (1.0 - p) / (rem - 1);
                    useP[n - rem][i] += prob[b] * (1.0 - p) / (rem - 1);
                }
            }
        }
    }
    return useP;
}

int main()
{
    int N;
    double pa,pb;
    cin>>N>>pa>>pb;
    VI A(N); FOR(i,N) cin>>A[i];
    VI B(N); FOR(i,N) cin>>B[i];

    auto t1 = getTable(A, pa);
    auto t2 = getTable(B, pb);
    double ans = 0 ;

    FOR(t,N){
        FOR(i,N) FOR(j,N){
            if(A[i] > B[j]){
                ans += t1[t][i] * t2[t][j] * (A[i] + B[j]);
            }
        }
    }

    printf("%.10f\n", ans);
    
    return 0;
}
0