結果

問題 No.1688 Veterinarian
コンテスト
ユーザー ayaoni
提出日時 2021-09-24 22:18:19
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 2,372 ms / 3,000 ms
コード長 1,046 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 225 ms
コンパイル使用メモリ 85,376 KB
実行使用メモリ 320,640 KB
最終ジャッジ日時 2026-03-26 17:42:26
合計ジャッジ時間 6,714 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge3_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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