結果
| 問題 | 
                            No.66 輝け☆全国たこやき杯
                             | 
                    
| コンテスト | |
| ユーザー | 
                             convexineq
                         | 
                    
| 提出日時 | 2020-12-13 10:23:52 | 
| 言語 | PyPy3  (7.3.15)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 69 ms / 5,000 ms | 
| コード長 | 501 bytes | 
| コンパイル時間 | 135 ms | 
| コンパイル使用メモリ | 81,904 KB | 
| 実行使用メモリ | 76,336 KB | 
| 最終ジャッジ日時 | 2024-09-19 23:24:42 | 
| 合計ジャッジ時間 | 1,101 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge5 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 10 | 
ソースコード
m,*a = map(int,open(0).read().split())
M = 1<<m
for i in range(M): a[i] *= a[i]
dp = [1.0]*M
for i in range(m):
    ndp = dp[:]
    L = 1<<i
    for j in range(1<<(m-i-1)):
        for k in range(L):
            nk = j*L*2+k
            v = w = 0.0
            for r in range(L):
                v += dp[j*L*2+L+r]*(a[nk]/(a[nk]+a[j*L*2+L+r]))
                w += dp[j*L*2+r]*(a[nk+L]/(a[nk+L]+a[j*L*2+r]))                
            ndp[nk] *= v
            ndp[nk+L] *= w
    dp = ndp
print(dp[0])
            
            
            
        
            
convexineq