結果

問題 No.174 カードゲーム（Hard）
ユーザー yaoshimaxyaoshimax
提出日時 2015-05-06 18:04:57
言語 PyPy2
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,478 bytes
コンパイル時間 1,160 ms
コンパイル使用メモリ 77,636 KB
実行使用メモリ 96,416 KB
最終ジャッジ日時 2023-09-19 08:37:27
合計ジャッジ時間 13,634 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 85 ms
75,500 KB
testcase_01 AC 81 ms
76,052 KB
testcase_02 TLE -
testcase_03 TLE -
testcase_04 TLE -
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 TLE -
testcase_10 AC 83 ms
76,536 KB
testcase_11 AC 84 ms
76,740 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def bitcount(N):
    ans=0
    while N!=0:
        ans+=1
        N&=N-1
    return ans

N,Pa,Pb=raw_input().split()
N=int(N)
Pa,Pb=float(Pa),float(Pb)
A=map(int,raw_input().split())
B=map(int,raw_input().split())
A.sort()
B.sort()
dpb=[[0.0 for i in range(N+1)] for j in range(N+1)]
dpmaskb=[0.0 for i in range(1<<N)]
dpmaskb[(1<<N)-1]=1.0
dpa=[[0.0 for i in range(N+1)] for j in range(N+1)]
dpmaska=[0.0 for i in range(1<<N)]
dpmaska[(1<<N)-1]=1.0

for i in range((1<<N)-1,-1,-1):
#    print "dpmaska[",i,"]=",dpmaska[i],
    for j in range(N):
        n=bitcount(i)
        if (i&(1<<j))!=0:
            k=i-(1<<j)
            if k==(i&(i-1)):
                if n==1:
                    dpmaska[k]+=dpmaska[i]
                    dpmaskb[k]+=dpmaskb[i]
                    dpa[N-n][j]+=dpmaska[i]
                    dpb[N-n][j]+=dpmaskb[i]
                else: 
                    dpmaska[k]+=dpmaska[i]*Pa
                    dpmaskb[k]+=dpmaskb[i]*Pb
                    dpa[N-n][j]+=dpmaska[i]*Pa
                    dpb[N-n][j]+=dpmaskb[i]*Pb
            else:
                dpmaska[k]+=dpmaska[i]*(1.0-Pa)/(n-1)
                dpmaskb[k]+=dpmaskb[i]*(1.0-Pb)/(n-1)
                dpa[N-n][j]+=dpmaska[i]*(1.0-Pa)/(n-1)
                dpb[N-n][j]+=dpmaskb[i]*(1.0-Pb)/(n-1)
ans=0.0
for ai in range(N):
    for bi in range(N):
        for ti in range(N):
            if A[ai] > B[bi]:
                ans+=(A[ai]+B[bi])*dpa[ti][ai]*dpb[ti][bi]
print "%.11f"%ans
0