結果

問題 No.2685 Cell Proliferation (Easy)
ユーザー umezoumezo
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0