結果
問題 |
No.3174 勝ち残りじゃんけん
|
ユーザー |
|
提出日時 | 2025-06-07 14:05:52 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 336 ms / 2,000 ms |
コード長 | 1,631 bytes |
コンパイル時間 | 1,878 ms |
コンパイル使用メモリ | 192,444 KB |
実行使用メモリ | 66,560 KB |
最終ジャッジ日時 | 2025-06-07 14:05:58 |
合計ジャッジ時間 | 5,804 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef pair<int, int> pii; typedef long long ll; const int N = 2000008, MOD = 998244353, INF = 0x3f3f3f3f; ll res; int n, m, cnt, w[N]; ll f[5086], c[5086], fac[N], inv[N], p2[N], p3[N], f2[5086], c2[5086]; inline ll qmi(ll a, ll b, ll c) { ll res = 1; while (b) { if (b & 1) res = res * a % c; a = a * a % c; b >>= 1; } return res; } void init() { p3[0] = p2[0] = fac[0] = inv[0] = 1; for (int i = 1; i < N; i++) fac[i] = fac[i - 1] * i % MOD, p3[i] = p3[i - 1] * 3 % MOD, p2[i] = p2[i - 1] * 2 % MOD; inv[N - 1] = qmi(fac[N - 1], MOD - 2, MOD); for (int i = N - 2; i; i--) inv[i] = inv[i + 1] * (i + 1) % MOD; } ll C(int a, int b) { if (a < b) return 0; return fac[a] * inv[b] % MOD * inv[a - b] % MOD; } int main() { init(); cin >> n; f2[n] = f[n] = 0, c2[n] = c[n] = 1; for (int i = n; i > 1; i--) { ll sum = p3[i], v = (p2[i] * 3 % MOD - 6 + MOD) % MOD; ll len = sum * qmi(v, MOD - 2, MOD) % MOD; ll pre = 0; ll abc = qmi(v, MOD - 2, MOD), abd = qmi(c[i], MOD - 2, MOD); for (int j = 1; j < i; j++) { ll t = C(i, j) * 3 % MOD; t = t * abc % MOD; c[j] = (c[j] + t * c[i]) % MOD; c2[j] = (c2[j] + (t + pre) % MOD * c[i]) % MOD; ll g = (f[i] * abd + len) % MOD; f[j] = (f[j] + g * c[i] % MOD * t % MOD) % MOD; f2[j] = (f2[j] + g * c[i] % MOD * ((t + pre) % MOD)) % MOD; pre = t + pre; } } for (int i = 1; i < n + 1; i++) printf("%lld ", f2[i] * qmi(c2[i], MOD - 2, MOD) % MOD); return 0; }