結果
問題 | No.1688 Veterinarian |
ユーザー | None |
提出日時 | 2021-10-06 05:24:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 510 ms / 3,000 ms |
コード長 | 852 bytes |
コンパイル時間 | 285 ms |
コンパイル使用メモリ | 87,056 KB |
実行使用メモリ | 80,916 KB |
最終ジャッジ日時 | 2023-09-30 08:37:42 |
合計ジャッジ時間 | 4,400 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 69 ms
71,480 KB |
testcase_01 | AC | 71 ms
71,388 KB |
testcase_02 | AC | 70 ms
71,152 KB |
testcase_03 | AC | 68 ms
71,540 KB |
testcase_04 | AC | 73 ms
71,564 KB |
testcase_05 | AC | 70 ms
71,420 KB |
testcase_06 | AC | 89 ms
76,844 KB |
testcase_07 | AC | 117 ms
79,188 KB |
testcase_08 | AC | 510 ms
80,724 KB |
testcase_09 | AC | 469 ms
80,916 KB |
testcase_10 | AC | 156 ms
78,776 KB |
testcase_11 | AC | 168 ms
78,820 KB |
testcase_12 | AC | 110 ms
77,888 KB |
testcase_13 | AC | 120 ms
78,056 KB |
testcase_14 | AC | 125 ms
78,172 KB |
testcase_15 | AC | 139 ms
77,956 KB |
testcase_16 | AC | 115 ms
77,896 KB |
ソースコード
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): for a in range(A+1): for b in range(B+1): for c in range(C+1): if a+b+c<=1: continue if a>=1: dp[a-1][b][c]+=dp[a][b][c]*a*(a-1)/(a+b+c)/(a+b+c-1) if b>=1: dp[a][b-1][c]+=dp[a][b][c]*b*(b-1)/(a+b+c)/(a+b+c-1) if c>=1: dp[a][b][c-1]+=dp[a][b][c]*c*(c-1)/(a+b+c)/(a+b+c-1) dp[a][b][c]*=(1-(a*(a-1)+b*(b-1)+c*(c-1))/(a+b+c)/(a+b+c-1)) res=[0,0,0] for a in range(A+1): for b in range(B+1): for c in range(C+1): res[0]+=(A-a)*dp[a][b][c] res[1]+=(B-b)*dp[a][b][c] res[2]+=(C-c)*dp[a][b][c] print(*res)