結果
| 問題 | 
                            No.66 輝け☆全国たこやき杯
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2024-07-07 08:40:21 | 
| 言語 | Python3  (3.13.1 + numpy 2.2.1 + scipy 1.14.1)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 507 ms / 5,000 ms | 
| コード長 | 682 bytes | 
| コンパイル時間 | 107 ms | 
| コンパイル使用メモリ | 12,544 KB | 
| 実行使用メモリ | 11,008 KB | 
| 最終ジャッジ日時 | 2024-07-07 08:40:23 | 
| 合計ジャッジ時間 | 2,058 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge4 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 10 | 
ソースコード
import math
def dfs(s, n):
    ret = [0.0] * n
    if n == 1:
        ret[0] = 1
        return ret
    
    a = dfs(s, n // 2)
    b = dfs(s + n // 2, n // 2)
    for i in range(len(a)):
        for j in range(len(b)):
            p = a[i] * b[j]
            ap = s + i
            bp = s + n // 2 + j
            ret[i] += p * params[ap] * params[ap] / (params[ap] * params[ap] + params[bp] * params[bp])
            ret[n // 2 + j] += p * params[bp] * params[bp] / (params[ap] * params[ap] + params[bp] * params[bp])
    return ret
M = int(input())
mem = int(math.pow(2, M))
params = []
for _ in range(mem):
    params.append(int(input()))
    
ans = dfs(0, mem)
print(ans[0])