結果

問題 No.2685 Cell Proliferation (Easy)
コンテスト
ユーザー norioc
提出日時 2025-10-20 22:19:45
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,809 ms / 2,000 ms
コード長 1,057 bytes
コンパイル時間 370 ms
コンパイル使用メモリ 82,904 KB
実行使用メモリ 280,456 KB
最終ジャッジ日時 2025-10-20 22:20:06
合計ジャッジ時間 15,094 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections.abc import Iterable


def mod_inv(n: int, mod: int) -> int:
    return pow(n, mod-2, mod)


def accum_dp1(xs: Iterable, f, op, e, size: int, init: Iterable, *, is_reset=True):
    dp = [e] * size
    for i, v in init:
        dp[i] = v

    for x in xs:
        pp = [e] * size if is_reset else dp.copy()
        dp, pp = pp, dp
        for i in range(size):
            for p, v in f(i, pp[i], x):
                if not (0 <= p < size): continue
                dp[p] = op(dp[p], v)

    return dp


memo = {}
def modpow(a, b, mod):
    if (a, b) in memo: return memo[a, b]

    res = pow(a, b, mod)
    memo[a, b] = res
    return res



MOD = 998244353
P1, P2, Q1, Q2, T = map(int, input().split())

p = P1 * mod_inv(P2, MOD) % MOD
q = Q1 * mod_inv(Q2, MOD) % MOD


def f(i, v, x):
    res = []
    # 新たに誕生
    res.append((0, p * v))
    # 生存
    res.append((i+1, modpow(q, i+1, MOD) * v))

    return res


def op(a, b):
    return a + b


dp = accum_dp1(range(T), f, op, 0, T+1, [(0, 1)])
ans = sum(dp) % MOD
print(ans)
0