結果

問題 No.2685 Cell Proliferation (Easy)
ユーザー kakel-san
提出日時 2024-03-21 23:43:11
言語 C#
(.NET 8.0.404)
結果
AC  
実行時間 100 ms / 2,000 ms
コード長 1,627 bytes
コンパイル時間 11,768 ms
コンパイル使用メモリ 168,300 KB
実行使用メモリ 188,268 KB
最終ジャッジ日時 2024-09-30 10:24:41
合計ジャッジ時間 15,199 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (109 ms)。
MSBuild のバージョン 17.9.6+a4ecab324 (.NET)
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

diff # 

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;

class Program
{
    static int NN => int.Parse(ReadLine());
    static int[] NList => ReadLine().Split().Select(int.Parse).ToArray();
    public static void Main()
    {
        Solve();
    }
    static void Solve()
    {
        var c = NList;
        var (p1, p2, q1, q2, t) = (c[0], c[1], c[2], c[3], c[4]);
        var mod = 998_244_353;
        var p2rev = Exp(p2, mod - 2, mod);
        var q2rev = Exp(q2, mod - 2, mod);
        var prob = new long[t + 1];
        prob[0] = 1;
        prob[1] = q1 * q2rev % mod;
        for (var i = 2; i <= t; ++i)
        {
            prob[i] = prob[i - 1] * prob[1] % mod;
        }
        var dp = new long[t + 1];
        dp[0] = 1;
        for (var i = 0; i < t; ++i)
        {
            var ndp = new long[t + 1];
            for (var j = 0; j < t; ++j)
            {
                ndp[j + 1] = dp[j] * prob[j + 1] % mod;
                ndp[0] = (ndp[0] + dp[j] * p1 % mod * p2rev % mod) % mod;
            }
            dp = ndp;
        }
        var ans = 0L;
        for (var i = 0; i < dp.Length; ++i) ans = (ans + dp[i]) % mod;
        WriteLine(ans);
    }
    static long Exp(long n, long p, int mod)
    {
        long _n = n % mod;
        var _p = p;
        var result = 1L;
        if ((_p & 1) == 1) result *= _n;
        while (_p > 0)
        {
            _n = _n * _n % mod;
            _p >>= 1;
            if ((_p & 1) == 1) result = result * _n % mod;
        }
        return result;
    }
}
0