結果

問題 No.3174 勝ち残りじゃんけん
ユーザー The Forsaking
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

#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;
}
0