結果

問題 No.1688 Veterinarian
ユーザー kohei2019kohei2019
提出日時 2023-06-04 21:39:10
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 142 ms / 3,000 ms
コード長 931 bytes
コンパイル時間 256 ms
コンパイル使用メモリ 82,280 KB
実行使用メモリ 78,024 KB
最終ジャッジ日時 2024-06-09 04:08:32
合計ジャッジ時間 2,055 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
53,104 KB
testcase_01 AC 32 ms
52,792 KB
testcase_02 AC 33 ms
52,880 KB
testcase_03 AC 31 ms
53,316 KB
testcase_04 AC 32 ms
53,372 KB
testcase_05 AC 32 ms
54,216 KB
testcase_06 AC 40 ms
63,384 KB
testcase_07 AC 47 ms
66,492 KB
testcase_08 AC 142 ms
78,024 KB
testcase_09 AC 134 ms
77,700 KB
testcase_10 AC 83 ms
76,960 KB
testcase_11 AC 79 ms
77,056 KB
testcase_12 AC 43 ms
63,536 KB
testcase_13 AC 57 ms
72,060 KB
testcase_14 AC 53 ms
69,960 KB
testcase_15 AC 72 ms
75,156 KB
testcase_16 AC 41 ms
62,336 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

A,B,C,N = map(int,input().split())
dp = [[[0]*(C+1) for i in range(B+1)] for j in range(A+1)]
dp[A][B][C] = 1
# a,b,c枚あるとき、全パターン(a+b+c)C2,a2枚はaC2,b2枚はbC2,c2枚はcC2

for i in range(N):
    for j in range(1,A+1):
        for k in range(1,B+1):
            for l in range(1,C+1):
                if dp[j][k][l]==0:
                    continue
                alln = (j+k+l)*(j+k+l-1)//2
                dp[j-1][k][l] += dp[j][k][l]*(j)*(j-1)/2/alln
                dp[j][k-1][l] += dp[j][k][l]*(k)*(k-1)/2/alln
                dp[j][k][l-1] += dp[j][k][l]*(l)*(l-1)/2/alln
                dp[j][k][l] -= dp[j][k][l]*((j)*(j-1)/2/alln+(k)*(k-1)/2/alln+(l)*(l-1)/2/alln)

Av = 0
Bv = 0
Cv = 0


for j in range(1,A+1):
    for k in range(1,B+1):
        for l in range(1,C+1):
            Av += (A-j)*dp[j][k][l]
            Bv += (B-k)*dp[j][k][l]
            Cv += (C-l)*dp[j][k][l]

print(*[Av,Bv,Cv])
0