結果
問題 | No.1688 Veterinarian |
ユーザー | ああいい |
提出日時 | 2022-03-04 10:21:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 706 ms / 3,000 ms |
コード長 | 1,117 bytes |
コンパイル時間 | 1,373 ms |
コンパイル使用メモリ | 86,700 KB |
実行使用メモリ | 131,128 KB |
最終ジャッジ日時 | 2023-09-25 09:18:38 |
合計ジャッジ時間 | 6,188 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 73 ms
71,188 KB |
testcase_01 | AC | 72 ms
71,176 KB |
testcase_02 | AC | 75 ms
71,192 KB |
testcase_03 | AC | 71 ms
71,304 KB |
testcase_04 | AC | 71 ms
71,340 KB |
testcase_05 | AC | 77 ms
71,124 KB |
testcase_06 | AC | 93 ms
76,312 KB |
testcase_07 | AC | 134 ms
79,972 KB |
testcase_08 | AC | 706 ms
131,128 KB |
testcase_09 | AC | 657 ms
122,220 KB |
testcase_10 | AC | 196 ms
79,424 KB |
testcase_11 | AC | 202 ms
79,700 KB |
testcase_12 | AC | 121 ms
78,008 KB |
testcase_13 | AC | 127 ms
78,060 KB |
testcase_14 | AC | 150 ms
78,820 KB |
testcase_15 | AC | 172 ms
78,940 KB |
testcase_16 | AC | 115 ms
77,656 KB |
ソースコード
A,B,C,N = map(int,input().split()) dp = [[[0] * (C + 1) for _ in range(B + 1)] for _ in range(A + 1)] dp[-1][-1][-1] = 1 for _ in range(N): nx = [[[0] * (C +1) for _ in range(B + 1)] for _ in range(A + 1)] for a in range(A + 1): for b in range(B + 1): for c in range(C + 1): d = a + b + c u = dp[a][b][c] if d <= 2:continue if a >= 2: nx[a-1][b][c] += u * a * (a-1) / (d * (d-1)) if b >= 2: nx[a][b-1][c] += u * b * (b-1) / (d * (d-1)) if c >= 2: nx[a][b][c-1] += u * c * (c - 1) / (d * (d - 1)) tmp = a * (a - 1) + b * (b - 1) + c * (c - 1) tmp = 1 - tmp / (d * (d - 1)) nx[a][b][c] += u * tmp dp = nx white = 0 black = 0 brown = 0 for a in range(A + 1): for b in range(B + 1): for c in range(C + 1): u = dp[a][b][c] white += (A - a) * u black += (B - b) * u brown += (C - c) * u print(white,black,brown)