結果
問題 | No.1688 Veterinarian |
ユーザー | tcltk |
提出日時 | 2021-09-30 05:29:55 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 856 ms / 3,000 ms |
コード長 | 1,343 bytes |
コンパイル時間 | 1,603 ms |
コンパイル使用メモリ | 87,072 KB |
実行使用メモリ | 149,416 KB |
最終ジャッジ日時 | 2023-09-24 08:00:47 |
合計ジャッジ時間 | 9,241 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge12 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 214 ms
80,912 KB |
testcase_01 | AC | 216 ms
80,872 KB |
testcase_02 | AC | 209 ms
80,628 KB |
testcase_03 | AC | 207 ms
80,792 KB |
testcase_04 | AC | 206 ms
80,916 KB |
testcase_05 | AC | 207 ms
80,984 KB |
testcase_06 | AC | 230 ms
84,316 KB |
testcase_07 | AC | 261 ms
85,828 KB |
testcase_08 | AC | 856 ms
149,416 KB |
testcase_09 | AC | 815 ms
144,724 KB |
testcase_10 | AC | 322 ms
90,948 KB |
testcase_11 | AC | 329 ms
92,312 KB |
testcase_12 | AC | 250 ms
85,364 KB |
testcase_13 | AC | 263 ms
85,596 KB |
testcase_14 | AC | 277 ms
85,784 KB |
testcase_15 | AC | 298 ms
88,264 KB |
testcase_16 | AC | 244 ms
85,100 KB |
ソースコード
#!/usr/bin/env python3 # from typing import * import sys import io import math import collections import decimal import itertools import bisect import heapq def input(): return sys.stdin.readline()[:-1] # sys.setrecursionlimit(1000000) # _INPUT = """4 3 2 13 # """ # sys.stdin = io.StringIO(_INPUT) INF = 10**10 A, B, C, N = map(int, input().split()) dp = [[[[0.0] * (C+1) for _ in range(B+1)] for _ in range(A+1)] for _ in range(N+1)] dp[0][A][B][C] = 1.0 for i in range(N): for a in range(A+1): for b in range(B+1): for c in range(C+1): n = a+b+c if n <= 1: continue dp[i+1][a][b][c] += dp[i][a][b][c] * (a*(b+c) + b*(a+c) + c*(a+b)) / (n*(n-1)) if a >= 2: dp[i+1][a-1][b][c] += dp[i][a][b][c] * (a*(a-1))/(n*(n-1)) if b >= 2: dp[i+1][a][b-1][c] += dp[i][a][b][c] * (b*(b-1))/(n*(n-1)) if c >= 2: dp[i+1][a][b][c-1] += dp[i][a][b][c] * (c*(c-1))/(n*(n-1)) ans_a, ans_b, ans_c = 0.0, 0.0, 0.0 for a in range(A+1): for b in range(B+1): for c in range(C+1): ans_a += (A-a) * dp[N][a][b][c] ans_b += (B-b) * dp[N][a][b][c] ans_c += (C-c) * dp[N][a][b][c] print(ans_a, ans_b, ans_c)