結果

問題 No.174 カードゲーム（Hard）
ユーザー yaoshimaxyaoshimax
提出日時 2015-05-06 18:04:57
言語 PyPy2
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,478 bytes
コンパイル時間 2,282 ms
コンパイル使用メモリ 76,408 KB
実行使用メモリ 95,336 KB
最終ジャッジ日時 2024-12-28 12:03:35
合計ジャッジ時間 21,362 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 TLE * 8
権限があれば一括ダウンロードができます

ソースコード

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