結果
問題 |
No.2685 Cell Proliferation (Easy)
|
ユーザー |
![]() |
提出日時 | 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; }