結果
問題 |
No.2685 Cell Proliferation (Easy)
|
ユーザー |
![]() |
提出日時 | 2024-03-20 21:58:35 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 146 ms / 2,000 ms |
コード長 | 1,399 bytes |
コンパイル時間 | 230 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 79,360 KB |
最終ジャッジ日時 | 2024-09-30 07:48:34 |
合計ジャッジ時間 | 3,925 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
mod=998244353 table_size=2*10**5 fac=[1]*(table_size+1) finv=[1]*(table_size+1) for i in range(2,table_size+1): fac[i]=fac[i-1]*i%mod finv[table_size]=pow(fac[table_size],mod-2,mod) for i in range(table_size-1,-1,-1): finv[i]=finv[i+1]*(i+1)%mod def rebuild(n): global table_size,fac,finv fac+=[0]*(n-table_size) fac+=[0]*(n-table_size) finv+=[0]*(n-table_size) for i in range(table_size+1,n+1): fac[i]=fac[i-1]*i%mod finv[n]=inv(fac[n]) for i in range(n-1,table_size,-1): finv[i]=finv[i+1]*(i+1)%mod table_size=n def binom(n,k): if n<0 or k<0: return 0 if k>n: return 0 if n>table_size: rebuild(n+10**4) return (fac[n]*finv[k]%mod)*finv[n-k]%mod def fpow(x,k): res=1 while k: if k&1: res=res*x%mod x=x*x%mod k>>=1 return res def inv(a): if a<table_size: return fac[a-1]*finv[a]%mod return fpow(a,mod-2) def modint_to_frac(a,mod): a%=mod if a==0: return '0/1' for X in range(1,10000): Y=a*X%mod if 0<Y<10000: return str(Y)+'/'+str(X) if mod-10000<Y<mod: return '-'+str(mod-Y)+'/'+str(X) return 'inexpressible' P1,P2,Q1,Q2,T=map(int,input().split()) P=P1*inv(P2)%mod Q=Q1*inv(Q2)%mod dp=[0]*(T+1) dp[0]=1 for t in range(1,T+1): for t1 in range(t): dp[t]+=dp[t1]*P dp[t]%=mod for t2 in range(t): dp[t2]=dp[t2]*fpow(Q,t-t2) dp[t2]%=mod print(sum(dp)%mod)