結果

問題 No.1688 Veterinarian
ユーザー ayaoniayaoni
提出日時 2021-09-24 22:18:19
言語 PyPy3
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,046 bytes
コンパイル時間 306 ms
コンパイル使用メモリ 87,144 KB
実行使用メモリ 276,564 KB
最終ジャッジ日時 2023-09-18 21:37:15
合計ジャッジ時間 10,214 ms
ジャッジサーバーID
(参考情報)
judge12 / judge14
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 73 ms
71,216 KB
testcase_01 AC 73 ms
71,292 KB
testcase_02 AC 82 ms
76,532 KB
testcase_03 AC 73 ms
71,044 KB
testcase_04 AC 73 ms
71,404 KB
testcase_05 AC 75 ms
71,256 KB
testcase_06 AC 89 ms
76,668 KB
testcase_07 AC 164 ms
92,220 KB
testcase_08 TLE -
testcase_09 TLE -
testcase_10 AC 345 ms
116,104 KB
testcase_11 AC 442 ms
139,720 KB
testcase_12 AC 124 ms
79,204 KB
testcase_13 AC 119 ms
78,244 KB
testcase_14 AC 157 ms
83,480 KB
testcase_15 AC 255 ms
90,956 KB
testcase_16 AC 107 ms
78,244 KB
権限があれば一括ダウンロードができます

ソースコード

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