結果
| 問題 |
No.2685 Cell Proliferation (Easy)
|
| コンテスト | |
| ユーザー |
SSRS
|
| 提出日時 | 2024-03-20 21:33:02 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 15 ms / 2,000 ms |
| コード長 | 942 bytes |
| コンパイル時間 | 2,016 ms |
| コンパイル使用メモリ | 195,916 KB |
| 最終ジャッジ日時 | 2025-02-20 08:49:18 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 26 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
long long modpow(long long a, long long b){
long long ans = 1;
while (b > 0){
if (b % 2 == 1){
ans *= a;
ans %= MOD;
}
a *= a;
a %= MOD;
b /= 2;
}
return ans;
}
long long modinv(long long a){
return modpow(a, MOD - 2);
}
int main(){
long long P1, P2, Q1, Q2;
int T;
cin >> P1 >> P2 >> Q1 >> Q2 >> T;
long long P = P1 * modinv(P2) % MOD;
long long Q = Q1 * modinv(Q2) % MOD;
vector<long long> POW(T + 1);
POW[0] = 1;
for (int i = 1; i <= T; i++){
POW[i] = POW[i - 1] * Q % MOD;
}
vector<long long> dp(T + 1, 0);
dp[0] = 1;
long long ans = 0;
for (int i = 0; i <= T; i++){
long long tmp = dp[i];
for (int j = 1; j <= T - i; j++){
dp[i + j] += tmp * P;
dp[i + j] %= MOD;
tmp *= POW[j];
tmp %= MOD;
}
ans += tmp;
}
ans %= MOD;
cout << ans << endl;
}
SSRS