結果

問題 No.1688 Veterinarian
ユーザー ayaoni
提出日時 2021-09-24 22:18:19
言語 PyPy3
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,046 bytes
コンパイル時間 149 ms
コンパイル使用メモリ 81,896 KB
実行使用メモリ 276,960 KB
最終ジャッジ日時 2024-07-05 10:49:28
合計ジャッジ時間 9,342 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 12 TLE * 2
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
input = sys.stdin.readline
sys.setrecursionlimit(10**7)

A,B,C,N = map(int,input().split())

dp = [[[(0,0,0)]*(C+1) for _ in range(B+1)] for _ in range(A+1)]
for _ in range(N):
    new_dp = [[[(0,0,0)]*(C+1) for _ in range(B+1)] for _ in range(A+1)]
    for i in range(1,A+1):
        for j in range(1,B+1):
            for k in range(1,C+1):
                s = i+j+k
                ss = s*(s-1)
                ww = i*(i-1)/ss
                bb = j*(j-1)/ss
                cc = k*(k-1)/ss
                q = 1-ww-bb-cc
                prob = [ww,bb,cc,q]
                white,black,brown = ww,bb,cc
                X = [(1,0,0),(0,1,0),(0,0,1),(0,0,0)]
                for l in range(4):
                    p = prob[l]
                    di,dj,dk = X[l]
                    white += p*dp[i-di][j-dj][k-dk][0]
                    black += p*dp[i-di][j-dj][k-dk][1]
                    brown += p*dp[i-di][j-dj][k-dk][2]
                new_dp[i][j][k] = (white,black,brown)
    dp = new_dp

w,b,c = dp[-1][-1][-1]
print(w,b,c)
0