結果
問題 | No.1688 Veterinarian |
ユーザー |
👑 ![]() |
提出日時 | 2021-09-24 23:43:13 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 648 ms / 3,000 ms |
コード長 | 820 bytes |
コンパイル時間 | 174 ms |
コンパイル使用メモリ | 82,688 KB |
実行使用メモリ | 207,312 KB |
最終ジャッジ日時 | 2024-07-05 11:37:37 |
合計ジャッジ時間 | 3,178 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 14 |
ソースコード
A,B,C,N=map(int,input().split())DP=[[[[0]*(C+2) for ___ in range(B+2)] for _ in range(A+2)] for _ in range(N+1)]DP[0][A][B][C]=1for n in range(1,N+1):for a in range(A+1):for b in range(B+1):for c in range(C+1):K=a+b+cif K<=1:continues=(K*(K-1)-(a*(a-1)+b*(b-1)+c*(c-1)))/(K*(K-1))t=((a+1)*a)/((K+1)*K)u=((b+1)*b)/((K+1)*K)v=((c+1)*c)/((K+1)*K)DP[n][a][b][c]=s*DP[n-1][a][b][c]+t*DP[n-1][a+1][b][c]+u*DP[n-1][a][b+1][c]+v*DP[n-1][a][b][c+1]S=T=U=0for a in range(A+1):for b in range(B+1):for c in range(C+1):S+=(A-a)*DP[N][a][b][c]T+=(B-b)*DP[N][a][b][c]U+=(C-c)*DP[N][a][b][c]print(S,T,U)