結果

問題 No.174 カードゲーム(Hard)
ユーザー yaoshimaxyaoshimax
提出日時 2015-05-06 18:07:04
言語 PyPy2
(7.3.13)
結果
AC  
実行時間 870 ms / 2,000 ms
コード長 1,508 bytes
コンパイル時間 1,639 ms
コンパイル使用メモリ 77,488 KB
実行使用メモリ 104,472 KB
最終ジャッジ日時 2023-09-21 02:58:54
合計ジャッジ時間 9,543 ms
ジャッジサーバーID
(参考情報)
judge12 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 76 ms
76,468 KB
testcase_01 AC 77 ms
76,228 KB
testcase_02 AC 870 ms
104,396 KB
testcase_03 AC 857 ms
104,336 KB
testcase_04 AC 857 ms
104,184 KB
testcase_05 AC 855 ms
104,348 KB
testcase_06 AC 857 ms
104,472 KB
testcase_07 AC 853 ms
104,368 KB
testcase_08 AC 853 ms
104,284 KB
testcase_09 AC 861 ms
104,304 KB
testcase_10 AC 76 ms
76,208 KB
testcase_11 AC 74 ms
76,620 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
b=[bitcount(i) for i in range(1<<N)]

for i in range((1<<N)-1,-1,-1):
#    print "dpmaska[",i,"]=",dpmaska[i],
    for j in range(N):
        n=b[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):
        if A[ai] > B[bi]:
            for ti in range(N):
                ans+=(A[ai]+B[bi])*dpa[ti][ai]*dpb[ti][bi]
print "%.11f"%ans
0