結果
問題 | No.2685 Cell Proliferation (Easy) |
ユーザー |
![]() |
提出日時 | 2024-03-22 08:23:48 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 15 ms / 2,000 ms |
コード長 | 1,026 bytes |
コンパイル時間 | 2,919 ms |
コンパイル使用メモリ | 248,384 KB |
実行使用メモリ | 10,368 KB |
最終ジャッジ日時 | 2024-09-30 10:30:12 |
合計ジャッジ時間 | 4,011 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#define rep(i,n) for(int i=0;i<(int)(n);i++) #define ALL(v) v.begin(),v.end() typedef long long ll; #include<bits/stdc++.h> using namespace std; const int MOD=998244353; ll modpow(ll x,ll n){ x%=MOD; ll ans=1; while(n){ if(n&1) ans=ans*x%MOD; x=x*x%MOD; n/=2; } return ans; } ll dp[1010][1010]; int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); ll p1,p2,q1,q2,t; cin>>p1>>p2>>q1>>q2>>t; ll p=p1*modpow(p2,MOD-2)%MOD; ll q=q1*modpow(q2,MOD-2)%MOD; vector<ll> Q(t+2); Q[0]=1; rep(i,t+1) Q[i+1]=Q[i]*q%MOD; for(int i=1;i<=t;i++) Q[i]=Q[i]*Q[i-1]%MOD; dp[0][0]=1; for(int i=0;i<=t;i++){ for(int j=0;j<=t;j++){ if(dp[i][j]==0) continue; dp[i+1][j+1]=(dp[i+1][j+1]+dp[i][j])%MOD; } if(i==t) break; for(int j=1;j<=t;j++){ if(dp[i+1][j]==0) continue; dp[i+1][0]=(dp[i+1][0]+dp[i+1][j]*p%MOD*Q[j-1]%MOD)%MOD; } } ll ans=0; for(int j=0;j<=t+1;j++) ans=(ans+dp[t][j]*Q[j]%MOD)%MOD; cout<<ans<<endl; return 0; }