結果
問題 | No.66 輝け☆全国たこやき杯 |
ユーザー | rpy3cpp |
提出日時 | 2015-07-22 01:06:25 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 578 bytes |
コンパイル時間 | 98 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-07-08 11:44:56 |
合計ジャッジ時間 | 1,000 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 25 ms
10,752 KB |
testcase_01 | AC | 24 ms
10,496 KB |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
ソースコード
def read_data(): M = int(input()) S = [int(input()) ** 2 for i in range(2**M)] return M, S def solve(M, S): N = 2**M p = [1.0] * N for m in range(0, M-1): np = [0.0] * N for i in range(N): pi = p[i] Si = S[i] start = i ^ (1 << m) for j in range(start, start + (1 << m)): np[i] += pi * p[j] * Si / (Si + S[j]) p = np p0 = p[0] S0 = S[0] ans = p0 * sum(p[j] * S0 /(S0 + S[j]) for j in range(N//2, N)) return ans M, S = read_data() print(solve(M, S))