結果

問題 No.2876 Infection
ユーザー 寝癖寝癖
提出日時 2024-07-06 00:14:46
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 251 ms / 2,000 ms
コード長 859 bytes
コンパイル時間 2,559 ms
コンパイル使用メモリ 247,588 KB
実行使用メモリ 105,984 KB
最終ジャッジ日時 2024-07-06 01:08:39
合計ジャッジ時間 5,890 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/modint>
using namespace std;
using namespace atcoder;

using mint = modint998244353;

mint c[5005], s[5005][5005], p, fact[5005];
bool c_used[5005], s_used[5005][5005];

mint C(int n);
mint S(int n, int k);

mint C(int n) {
    if (c_used[n]) return c[n];
    c_used[n] = true;
    mint res = 1;
    for (int i = 1; i <= n-1; i++) res -= S(n, i);
    return c[n] = res;
}

mint S(int n, int k) {
    if (s_used[n][k]) return s[n][k];
    s_used[n][k] = true;
    mint res = fact[n-1] / fact[k-1] / fact[n-k] * C(k) * (1-p).pow(k*(n-k));
    return s[n][k] = res;
}

int main() {
    int N, x;
    cin >> N >> x;

    p = (mint)x / 100;
    fact[0] = 1;
    for (int i = 1; i <= N; i++) fact[i] = fact[i-1] * i;

    mint ans = 0;
    for (int k = 1; k <= N; k++) ans += k * S(N, k);
    cout << ans.val() << endl;
}
0