結果
問題 | No.1688 Veterinarian |
ユーザー | 👑 tamato |
提出日時 | 2021-09-24 21:42:24 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 864 ms / 3,000 ms |
コード長 | 1,393 bytes |
コンパイル時間 | 907 ms |
コンパイル使用メモリ | 87,532 KB |
実行使用メモリ | 166,468 KB |
最終ジャッジ日時 | 2023-09-18 21:00:55 |
合計ジャッジ時間 | 5,143 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 71 ms
71,784 KB |
testcase_01 | AC | 73 ms
71,764 KB |
testcase_02 | AC | 75 ms
72,008 KB |
testcase_03 | AC | 74 ms
71,900 KB |
testcase_04 | AC | 75 ms
71,928 KB |
testcase_05 | AC | 75 ms
71,624 KB |
testcase_06 | AC | 101 ms
77,516 KB |
testcase_07 | AC | 130 ms
80,484 KB |
testcase_08 | AC | 864 ms
166,468 KB |
testcase_09 | AC | 806 ms
155,472 KB |
testcase_10 | AC | 190 ms
79,768 KB |
testcase_11 | AC | 200 ms
80,680 KB |
testcase_12 | AC | 118 ms
78,364 KB |
testcase_13 | AC | 133 ms
78,140 KB |
testcase_14 | AC | 132 ms
78,472 KB |
testcase_15 | AC | 165 ms
78,672 KB |
testcase_16 | AC | 114 ms
78,164 KB |
ソースコード
mod = 1000000007 eps = 10**-9 def main(): import sys input = sys.stdin.readline A, B, C, N = map(int, input().split()) dp = [[[0] * (C+1) for _ in range(B+1)] for _ in range(A+1)] dp[A][B][C] = 1 for _ in range(N): dp_new = [[[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): abc = a + b + c if abc <= 2: continue lower = abc * (abc - 1) pp = dp[a][b][c] pa = a * (a-1) / lower dp_new[a-1][b][c] = dp_new[a-1][b][c] + pp * pa pb = b * (b-1) / lower dp_new[a][b-1][c] = dp_new[a][b-1][c] + pp * pb pc = c * (c-1) / lower dp_new[a][b][c-1] = dp_new[a][b][c-1] + pp * pc dp_new[a][b][c] = dp_new[a][b][c] + pp * (1 - pa - pb - pc) dp = dp_new ans_a = 0 ans_b = 0 ans_c = 0 for a in range(A+1): for b in range(B+1): for c in range(C+1): ans_a = ans_a + dp[a][b][c] * (A - a) ans_b = ans_b + dp[a][b][c] * (B - b) ans_c = ans_c + dp[a][b][c] * (C - c) print(ans_a, ans_b, ans_c) if __name__ == '__main__': main()