結果
問題 | No.1688 Veterinarian |
ユーザー | lloyz |
提出日時 | 2021-11-23 17:55:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,689 ms / 3,000 ms |
コード長 | 1,721 bytes |
コンパイル時間 | 422 ms |
コンパイル使用メモリ | 87,188 KB |
実行使用メモリ | 266,096 KB |
最終ジャッジ日時 | 2023-09-07 21:00:05 |
合計ジャッジ時間 | 7,375 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 65 ms
71,296 KB |
testcase_01 | AC | 65 ms
71,392 KB |
testcase_02 | AC | 70 ms
71,392 KB |
testcase_03 | AC | 67 ms
71,336 KB |
testcase_04 | AC | 66 ms
71,372 KB |
testcase_05 | AC | 68 ms
71,240 KB |
testcase_06 | AC | 95 ms
78,380 KB |
testcase_07 | AC | 151 ms
85,524 KB |
testcase_08 | AC | 1,689 ms
266,096 KB |
testcase_09 | AC | 1,639 ms
251,916 KB |
testcase_10 | AC | 286 ms
94,172 KB |
testcase_11 | AC | 321 ms
96,644 KB |
testcase_12 | AC | 145 ms
80,328 KB |
testcase_13 | AC | 157 ms
80,764 KB |
testcase_14 | AC | 167 ms
81,144 KB |
testcase_15 | AC | 218 ms
88,280 KB |
testcase_16 | AC | 110 ms
78,580 KB |
ソースコード
a, b, c, n = map(int, input().split()) dpa = [[[[0.0 for _ in range(c)] for _ in range(b)] for _ in range(a)] for _ in range(n + 1)] dpb = [[[[0.0 for _ in range(c)] for _ in range(b)] for _ in range(a)] for _ in range(n + 1)] dpc = [[[[0.0 for _ in range(c)] for _ in range(b)] for _ in range(a)] for _ in range(n + 1)] for i in range(1, n + 1): for aa in range(a): for bb in range(b): for cc in range(c): all_num = aa + bb + cc + 3 prob_a = (aa + 1) * aa / all_num / (all_num - 1) prob_b = (bb + 1) * bb / all_num / (all_num - 1) prob_c = (cc + 1) * cc / all_num / (all_num - 1) prob_res = 1 - prob_a - prob_b - prob_c dpa[i][aa][bb][cc] = prob_res * dpa[i - 1][aa][bb][cc] + \ prob_a * (dpa[i - 1][aa - 1][bb][cc] + 1) + \ prob_b * dpa[i - 1][aa][bb - 1][cc] + \ prob_c * dpa[i - 1][aa][bb][cc - 1] dpb[i][aa][bb][cc] = prob_res * dpb[i - 1][aa][bb][cc] + \ prob_a * dpb[i - 1][aa - 1][bb][cc] + \ prob_b * (dpb[i - 1][aa][bb - 1][cc] + 1) + \ prob_c * dpb[i - 1][aa][bb][cc - 1] dpc[i][aa][bb][cc] = prob_res * dpc[i - 1][aa][bb][cc] + \ prob_a * dpc[i - 1][aa - 1][bb][cc] + \ prob_b * dpc[i - 1][aa][bb - 1][cc] + \ prob_c * (dpc[i - 1][aa][bb][cc - 1] + 1) print(dpa[-1][-1][-1][-1], dpb[-1][-1][-1][-1], dpc[-1][-1][-1][-1])