結果

問題 No.1688 Veterinarian
ユーザー tassei903
提出日時 2021-09-24 23:35:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 868 ms / 3,000 ms
コード長 1,111 bytes
コンパイル時間 272 ms
コンパイル使用メモリ 81,664 KB
実行使用メモリ 220,016 KB
最終ジャッジ日時 2024-07-05 11:35:57
合計ジャッジ時間 3,943 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = lambda :sys.stdin.readline()[:-1]
ni = lambda :int(input())
na = lambda :list(map(int,input().split()))
sys.setrecursionlimit(10**7)
yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")
no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")
#######################################################################
a,b,c,n = na()

import functools
@functools.lru_cache(maxsize = None)
def saiki(n, a, b, c):
    if n==0:
        return 0,0,0
    s = a+b+c
    pa = (a-1)*a/s/(s-1)
    pb = (b-1)*b/s/(s-1)
    pc = (c-1)*c/s/(s-1)
    po = 1-pa-pb-pc
    ra,rb,rc = 0,0,0
    if a>1:
        za,zb,zc = saiki(n-1,a-1,b,c)
        za+=1
        ra+=za*pa
        rb+=zb*pa
        rc+=zc*pa
    if b>1:
        za,zb,zc = saiki(n-1,a,b-1,c)
        zb+=1
        ra+=za*pb
        rb+=zb*pb
        rc+=zc*pb
    if c>1:
        za,zb,zc = saiki(n-1,a,b,c-1)
        zc+=1
        ra+=za*pc
        rb+=zb*pc
        rc+=zc*pc
    za,zb,zc = saiki(n-1,a,b,c)
    ra+=za*po
    rb+=zb*po
    rc+=zc*po
    return ra,rb,rc

print(*saiki(n,a,b,c))
0