結果
問題 | No.1688 Veterinarian |
ユーザー | chineristAC |
提出日時 | 2021-09-24 21:34:03 |
言語 | PyPy3 (7.3.13) |
結果 |
AC
|
実行時間 | 309 ms / 3,000 ms |
コード長 | 1,223 bytes |
コンパイル時間 | 275 ms |
コンパイル使用メモリ | 86,996 KB |
実行使用メモリ | 101,232 KB |
最終ジャッジ日時 | 2023-09-18 20:49:45 |
合計ジャッジ時間 | 3,527 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge12 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 93 ms
72,204 KB |
testcase_01 | AC | 93 ms
71,952 KB |
testcase_02 | AC | 101 ms
72,096 KB |
testcase_03 | AC | 93 ms
72,088 KB |
testcase_04 | AC | 93 ms
72,324 KB |
testcase_05 | AC | 94 ms
72,004 KB |
testcase_06 | AC | 113 ms
78,368 KB |
testcase_07 | AC | 123 ms
78,968 KB |
testcase_08 | AC | 309 ms
101,232 KB |
testcase_09 | AC | 293 ms
97,156 KB |
testcase_10 | AC | 162 ms
79,380 KB |
testcase_11 | AC | 163 ms
79,360 KB |
testcase_12 | AC | 112 ms
78,148 KB |
testcase_13 | AC | 133 ms
79,076 KB |
testcase_14 | AC | 127 ms
79,072 KB |
testcase_15 | AC | 165 ms
79,008 KB |
testcase_16 | AC | 118 ms
78,472 KB |
ソースコード
import sys,random input = lambda :sys.stdin.readline().rstrip() mi = lambda :map(int,input().split()) li = lambda :list(mi()) A,B,C,N = mi() dp = [[[0 for k in range(C+1)] for j in range(B+1)] for i in range(A+1)] dp[A][B][C] = 1 for _ in range(N): ndp = [[[0 for k in range(C+1)] for j in range(B+1)] for i in range(A+1)] for i in range(A+1): for j in range(B+1): for k in range(C+1): if not dp[i][j][k]: continue S = i + j + k if S==0: ndp[i][j][k] += dp[i][j][k] continue base = dp[i][j][k]/(S*(S-1)//2) if i: ndp[i-1][j][k] += i*(i-1)//2 *base if j: ndp[i][j-1][k] += j*(j-1)//2 *base if k: ndp[i][j][k-1] += k*(k-1)//2 *base ndp[i][j][k] += (S*(S-1)//2-i*(i-1)//2-j*(j-1)//2-k*(k-1)//2) * base dp = ndp res = [0,0,0] for i in range(A+1): for j in range(B+1): for k in range(C+1): p = dp[i][j][k] res[0] += (A-i) * p res[1] += (B-j) * p res[2] += (C-k) * p print(*res)