結果
問題 | No.1547 [Cherry 2nd Tune *] 偶然の勝利の確率 |
ユーザー | googol_S0 |
提出日時 | 2021-06-11 22:50:14 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 553 ms / 2,000 ms |
コード長 | 1,025 bytes |
コンパイル時間 | 277 ms |
コンパイル使用メモリ | 86,992 KB |
実行使用メモリ | 77,968 KB |
最終ジャッジ日時 | 2023-08-21 13:37:52 |
合計ジャッジ時間 | 10,613 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 72 ms
71,308 KB |
testcase_01 | AC | 81 ms
76,404 KB |
testcase_02 | AC | 100 ms
76,448 KB |
testcase_03 | AC | 82 ms
76,168 KB |
testcase_04 | AC | 79 ms
76,356 KB |
testcase_05 | AC | 78 ms
76,176 KB |
testcase_06 | AC | 77 ms
76,180 KB |
testcase_07 | AC | 78 ms
76,180 KB |
testcase_08 | AC | 79 ms
76,100 KB |
testcase_09 | AC | 78 ms
76,236 KB |
testcase_10 | AC | 77 ms
76,348 KB |
testcase_11 | AC | 78 ms
76,288 KB |
testcase_12 | AC | 78 ms
76,236 KB |
testcase_13 | AC | 191 ms
77,300 KB |
testcase_14 | AC | 136 ms
76,424 KB |
testcase_15 | AC | 102 ms
76,216 KB |
testcase_16 | AC | 118 ms
76,496 KB |
testcase_17 | AC | 118 ms
76,416 KB |
testcase_18 | AC | 117 ms
76,416 KB |
testcase_19 | AC | 93 ms
76,560 KB |
testcase_20 | AC | 346 ms
77,468 KB |
testcase_21 | AC | 105 ms
76,388 KB |
testcase_22 | AC | 80 ms
76,412 KB |
testcase_23 | AC | 546 ms
77,968 KB |
testcase_24 | AC | 548 ms
77,860 KB |
testcase_25 | AC | 550 ms
77,896 KB |
testcase_26 | AC | 553 ms
77,724 KB |
testcase_27 | AC | 549 ms
77,740 KB |
testcase_28 | AC | 548 ms
77,592 KB |
testcase_29 | AC | 548 ms
77,772 KB |
testcase_30 | AC | 551 ms
77,796 KB |
testcase_31 | AC | 549 ms
77,836 KB |
testcase_32 | AC | 548 ms
77,584 KB |
testcase_33 | AC | 547 ms
77,792 KB |
testcase_34 | AC | 129 ms
76,484 KB |
testcase_35 | AC | 128 ms
76,580 KB |
ソースコード
mod=998244353 def m(a,b): r=[[0]*len(b[0]) for i in range(len(a))] for i in range(len(a)): for k in range(len(b)): for j in range(len(b[0])): r[i][j]=(r[i][j]+a[i][k]*b[k][j])%mod return r def p(a,n): r=[[0]*len(a) for i in range(len(a))] b=[] for i in range(len(a)): r[i][i]=1 b.append(a[i][:]) l=n while l>0: if l&1: r=m(b,r) b=m(b,b) l>>=1 return r MA,NA,S=map(int,input().split()) MB,NB,T=map(int,input().split()) X,Y=MA*pow(NA,mod-2,mod)%mod,MB*pow(NB,mod-2,mod)%mod K=int(input()) L=S+T+1 A=[[0]*L for i in range(L)] B=[[0]*L for i in range(L)] A[0][0]=1 B[0][0]=1 A[-1][-1]=1 B[-1][-1]=1 for i in range(-T+1,S): x=i y=1 while x<=S: if x==S: A[x+T][i+T]=y break A[x+T][i+T]=(y*(1-X))%mod y=y*X%mod x+=1 x=i y=1 while x>=-T: if x==-T: B[x+T][i+T]=y break B[x+T][i+T]=(y*(1-Y))%mod y=y*Y%mod x-=1 C=[[0]*L for i in range(L)] C[T][-1]=1 Z=m(p(m(B,A),K),C) print(Z[-1][-1],Z[0][-1],sep='\n')