結果
問題 | No.2871 Universal Serial Bus |
ユーザー |
|
提出日時 | 2024-09-06 21:46:08 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
WA
|
実行時間 | - |
コード長 | 778 bytes |
コンパイル時間 | 280 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 11,392 KB |
最終ジャッジ日時 | 2024-09-06 21:46:26 |
合計ジャッジ時間 | 18,063 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 14 WA * 4 |
ソースコード
from decimal import DecimalH, W = map(int, input().split())S = [list(input()) for _ in range(H)]T = [list(input()) for _ in range(H)]S = [[1 if S[i][j] == '#' else 0 for j in range(W)] for i in range(H)]T = [[1 if T[i][j] == '#' else 0 for j in range(W)] for i in range(H)]a, b = 1, 1for i in range(H):for j in range(W):a &= (S[i][j] + T[i][j] == 1)b &= (S[H - 1 - i][W - 1 - j] + T[i][j] == 1)a = Decimal(a)b = Decimal(b)ans = Decimal(0.0)sp = Decimal(0.0)coef = Decimal(1.0)for k in range(1, 202020):if k % 2:prob = (Decimal(1.0) - sp) * a * (Decimal(1.0) - coef)else:prob = (Decimal(1.0) - sp) * b * (Decimal(1.0) - coef)ans += prob * Decimal(k)sp += probcoef /= Decimal(2.0)print(ans)