結果

問題 No.1688 Veterinarian
ユーザー lllllll88938494lllllll88938494
提出日時 2023-04-28 03:56:50
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 295 ms / 3,000 ms
コード長 1,107 bytes
コンパイル時間 1,023 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 155,264 KB
最終ジャッジ日時 2024-04-28 14:12:14
合計ジャッジ時間 2,846 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 46 ms
52,480 KB
testcase_01 AC 38 ms
52,992 KB
testcase_02 AC 37 ms
53,248 KB
testcase_03 AC 33 ms
52,224 KB
testcase_04 AC 40 ms
52,736 KB
testcase_05 AC 37 ms
52,224 KB
testcase_06 AC 54 ms
65,792 KB
testcase_07 AC 62 ms
73,856 KB
testcase_08 AC 295 ms
155,264 KB
testcase_09 AC 283 ms
153,116 KB
testcase_10 AC 114 ms
83,584 KB
testcase_11 AC 107 ms
84,992 KB
testcase_12 AC 55 ms
65,024 KB
testcase_13 AC 96 ms
77,952 KB
testcase_14 AC 67 ms
69,504 KB
testcase_15 AC 117 ms
81,864 KB
testcase_16 AC 55 ms
65,536 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

a,b,c,n=map(int,input().split())

dp=[[[[0]*(c+2) for i in range(b+2)]for i in range(a+2)]for i in range(n+2)]
dp[0][0][0][0] = 1

for x in range(n):
    for i in range(a+1):
        for j in range(b+1):
            for k in range(c+1):
                if dp[x][i][j][k] == 0:
                    continue
                ni=a-i
                nj=b-j
                nk=c-k
                nc2=(ni+nj+nk)*(ni+nj+nk-1)//2
                nic2=ni*(ni-1)//2
                njc2=nj*(nj-1)//2
                nkc2=nk*(nk-1)//2
                
                dp[x+1][i+1][j][k] += nic2 / nc2 * dp[x][i][j][k]
                dp[x+1][i][j+1][k] += njc2 / nc2 * dp[x][i][j][k]
                dp[x+1][i][j][k+1] += nkc2 / nc2 * dp[x][i][j][k]
                dp[x+1][i][j][k] += (1 - (nic2+njc2+nkc2) / nc2) * dp[x][i][j][k]
                
                
ii = 0
jj = 0
kk = 0
for i in range(a+1):
    for j in range(b+1):
        for k in range(c+1):
            ii +=  i * dp[n][i][j][k]
            jj +=  j * dp[n][i][j][k]
            kk +=  k * dp[n][i][j][k]
            
            
print(ii,jj,kk)
0