結果
問題 | No.2876 Infection |
ユーザー | SnowBeenDiding |
提出日時 | 2024-09-06 23:48:36 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 165 ms / 2,000 ms |
コード長 | 1,767 bytes |
コンパイル時間 | 5,529 ms |
コンパイル使用メモリ | 317,680 KB |
実行使用メモリ | 113,152 KB |
最終ジャッジ日時 | 2024-09-06 23:48:48 |
合計ジャッジ時間 | 10,836 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 114 ms
112,984 KB |
testcase_01 | AC | 112 ms
112,916 KB |
testcase_02 | AC | 114 ms
112,816 KB |
testcase_03 | AC | 112 ms
113,036 KB |
testcase_04 | AC | 111 ms
112,820 KB |
testcase_05 | AC | 163 ms
112,816 KB |
testcase_06 | AC | 163 ms
112,816 KB |
testcase_07 | AC | 163 ms
112,928 KB |
testcase_08 | AC | 164 ms
112,996 KB |
testcase_09 | AC | 165 ms
113,016 KB |
testcase_10 | AC | 118 ms
112,976 KB |
testcase_11 | AC | 118 ms
113,020 KB |
testcase_12 | AC | 115 ms
113,152 KB |
testcase_13 | AC | 124 ms
112,936 KB |
testcase_14 | AC | 154 ms
112,928 KB |
testcase_15 | AC | 138 ms
112,932 KB |
testcase_16 | AC | 114 ms
112,884 KB |
testcase_17 | AC | 112 ms
112,980 KB |
testcase_18 | AC | 122 ms
112,924 KB |
testcase_19 | AC | 113 ms
113,016 KB |
testcase_20 | AC | 128 ms
112,848 KB |
testcase_21 | AC | 111 ms
112,932 KB |
testcase_22 | AC | 131 ms
112,960 KB |
testcase_23 | AC | 125 ms
112,876 KB |
testcase_24 | AC | 118 ms
112,816 KB |
testcase_25 | AC | 132 ms
112,948 KB |
testcase_26 | AC | 112 ms
112,880 KB |
testcase_27 | AC | 122 ms
112,968 KB |
testcase_28 | AC | 113 ms
112,884 KB |
testcase_29 | AC | 148 ms
113,032 KB |
ソースコード
#include <atcoder/all> #include <bits/stdc++.h> #define rep(i, a, b) for (ll i = (ll)(a); i < (ll)(b); i++) using namespace atcoder; using namespace std; typedef long long ll; using mint = modint998244353; struct Comb { vector<mint> fact, ifact; int MAX_COM; Comb() {} Comb(int n) { MAX_COM = n; // ここでMAX入力を調整 init(998244353, MAX_COM); } void init(long long MOD, long long MAX_COM) { int n = MAX_COM; assert(n < MOD); fact = vector<mint>(n + 1); ifact = vector<mint>(n + 1); fact[0] = 1; for (int i = 1; i <= n; ++i) fact[i] = fact[i - 1] * i; ifact[n] = fact[n].inv(); for (int i = n; i >= 1; --i) ifact[i - 1] = ifact[i] * i; } mint operator()(long long n, long long k) { if (k < 0 || k > n) return 0; return fact[n] * ifact[k] * ifact[n - k]; } }; Comb comb(5000010); int main() { int N, x; cin >> N >> x; vector<int> alc(3000); vector<mint> c(3000); vector als(3000, vector<int>(3000)); vector s(3000, vector<mint>(3000)); mint ans = 0; mint p = (mint)x / 100; mint q = 1 - p; function<mint(int)> fc; function<mint(int, int)> fs; fc = [&](int n) { if (alc[n]) return c[n]; mint ret = 1; rep(i, 1, n) { ret -= fs(n, i); } alc[n] = 1; return c[n] = ret; }; fs = [&](int n, int k) { if (als[n][k]) return s[n][k]; mint ret = comb(n - 1, k - 1); ret *= fc(k); ret *= q.pow(k * (n - k)); als[n][k] = 1; return s[n][k] = ret; }; rep(k, 1, N + 1) { ans += k * fs(N, k); } cout << ans.val() << endl; }