結果

問題 No.1688 Veterinarian
ユーザー 已经死了已经死了
提出日時 2024-06-29 22:25:37
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 634 ms / 3,000 ms
コード長 2,954 bytes
コンパイル時間 274 ms
コンパイル使用メモリ 81,908 KB
実行使用メモリ 116,952 KB
最終ジャッジ日時 2024-06-29 22:25:44
合計ジャッジ時間 6,267 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 133 ms
92,244 KB
testcase_01 AC 135 ms
92,280 KB
testcase_02 AC 137 ms
92,408 KB
testcase_03 AC 136 ms
92,348 KB
testcase_04 AC 135 ms
92,212 KB
testcase_05 AC 140 ms
92,388 KB
testcase_06 AC 158 ms
93,268 KB
testcase_07 AC 217 ms
95,412 KB
testcase_08 AC 634 ms
116,952 KB
testcase_09 AC 629 ms
112,012 KB
testcase_10 AC 245 ms
95,220 KB
testcase_11 AC 259 ms
94,936 KB
testcase_12 AC 185 ms
94,092 KB
testcase_13 AC 194 ms
93,784 KB
testcase_14 AC 210 ms
94,444 KB
testcase_15 AC 223 ms
93,920 KB
testcase_16 AC 170 ms
93,396 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import os,sys,random,threading
#sys.exit() 退出程序
from random import randint,choice,shuffle
#randint(a,b)从[a,b]范围随机选择一个数
#choice(seq)seq可以是一个列表,元组或字符串,从seq中随机选取一个元素
#shuffle(x)将一个可变的序列x中的元素打乱
from copy import deepcopy
from io import BytesIO,IOBase
from types import GeneratorType
from functools import lru_cache,reduce
#reduce(op,迭代对象)
from bisect import bisect_left,bisect_right
#bisect_left(x) 大于等于x的第一个下标
#bisect_right(x) 大于x的第一个下标
from collections import Counter,defaultdict,deque
from itertools import accumulate,combinations,permutations
#accumulate(a)用a序列生成一个累积迭代器，一般list化前面放个[0]做前缀和用
#combinations(a,k)a序列选k个 组合迭代器
#permutations(a,k)a序列选k个 排列迭代器
from heapq import  heapify,heappop,heappush
#heapify将列表转为堆
from typing import Generic,Iterable,Iterator,TypeVar,Union,List
from string import ascii_lowercase,ascii_uppercase,digits
#小写字母，大写字母，十进制数字
from math import ceil,floor,sqrt,pi,factorial,gcd,log,log10,log2,inf
#ceil向上取整，floor向下取整 ，sqrt开方 ，factorial阶乘
from decimal import Decimal,getcontext
#Decimal(s) 实例化Decimal对象,一般使用字符串
#getcontext().prec=100 修改精度
from sys import stdin, stdout, setrecursionlimit
input = lambda: sys.stdin.readline().rstrip("\r\n")
MI = lambda :map(int,input().split())
li = lambda :list(MI())
ii = lambda :int(input())
mod = int(1e9 + 7) #998244353
inf = int(1e20)
py = lambda :print("YES")
pn = lambda :print("NO")
DIRS = [(0, 1), (1, 0), (0, -1), (-1, 0)]  # 右下左上
DIRS8 = [(0, 1), (1, 1), (1, 0), (1, -1), (0, -1), (-1, -1), (-1, 0),(-1, 1)]  # →↘↓↙←↖↑↗


aa,bb,cc,n=li()

g=[[[0]*(cc+1) for _ in range(bb+1)] for _ in range(aa+1)]
g[aa][bb][cc]=1



for _ in range(n):
    f=[[[0]*(cc+1) for _ in range(bb+1)] for _ in range(aa+1)]
    for a in range(aa+1):
        for b in range(bb+1):
            for c in range(cc+1):
                if a+b+c<2:
                    f[a][b][c]+=g[a][b][c]
                else:
                    tot=(a+b+c)*(a+b+c-1)//2
                    t=tot
                    if a>=2:
                        f[a-1][b][c]+=g[a][b][c]*((a*(a-1)//2)/tot)
                        t-=a*(a-1)//2
                    if b>=2:
                        f[a][b-1][c]+=g[a][b][c]*((b*(b-1)//2)/tot)
                        t-=b*(b-1)//2
                    if c>=2:
                        f[a][b][c-1]+=g[a][b][c]*((c*(c-1)//2)/tot)
                        t-=c*(c-1)//2
                    f[a][b][c]+=g[a][b][c]*(t/tot)
    g=f


ta=tb=tc=0

for a in range(aa+1):
    for b in range(bb+1):
        for c in range(cc+1):
            ta+=(aa-a)*g[a][b][c]
            tb+=(bb-b)*g[a][b][c]
            tc+=(cc-c)*g[a][b][c]

print(ta,tb,tc)
0