結果
| 問題 |
No.1688 Veterinarian
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-01-25 17:09:01 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 790 ms / 3,000 ms |
| コード長 | 1,025 bytes |
| コンパイル時間 | 399 ms |
| コンパイル使用メモリ | 82,304 KB |
| 実行使用メモリ | 181,608 KB |
| 最終ジャッジ日時 | 2024-12-16 08:49:19 |
| 合計ジャッジ時間 | 4,285 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 14 |
ソースコード
memo = {}
def dfs(a, b, c, n):
if (a, b, c, n) in memo:
return memo[(a, b, c, n)]
tot = a + b + c
if n == 0 or tot <= 1:
return 0, 0, 0
ret = [0, 0, 0]
div = tot * (tot - 1) // 2
dif = div
if a >= 2:
aa = a * (a - 1) // 2
dif -= aa
tmp = dfs(a - 1, b, c, n - 1)
ret[0] += aa / div
for i in range(3):
ret[i] += tmp[i] * aa / div
if b >= 2:
bb = b * (b - 1) // 2
dif -= bb
tmp = dfs(a, b - 1, c, n - 1)
ret[1] += bb / div
for i in range(3):
ret[i] += tmp[i] * bb / div
if c >= 2:
cc = c * (c - 1) // 2
dif -= cc
tmp = dfs(a, b, c - 1, n - 1)
ret[2] += cc / div
for i in range(3):
ret[i] += tmp[i] * cc / div
tmp = dfs(a, b, c, n - 1)
for i in range(3):
ret[i] += tmp[i] * dif / div
memo[(a, b, c, n)] = ret
return ret
a, b, c, n = map(int, input().split())
print(*dfs(a, b, c, n))