結果
問題 |
No.3174 勝ち残りじゃんけん
|
ユーザー |
![]() |
提出日時 | 2025-06-07 17:50:59 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 199 ms / 2,000 ms |
コード長 | 2,517 bytes |
コンパイル時間 | 720 ms |
コンパイル使用メモリ | 57,184 KB |
実行使用メモリ | 6,272 KB |
最終ジャッジ日時 | 2025-06-07 17:51:02 |
合計ジャッジ時間 | 2,917 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
コンパイルメッセージ
main.cpp: In function ‘int main()’: main.cpp:88:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 88 | scanf("%d", &n); | ~~~~~^~~~~~~~~~
ソースコード
/* -*- coding: utf-8 -*- * * 3174.cc: No.3174 蜍昴■谿九j縺倥c繧薙¢繧・- yukicoder */ #include<cstdio> #include<vector> #include<algorithm> using namespace std; /* constant */ const int MAX_N = 5000; const int MOD = 998244353; /* typedef */ template<const int MOD> struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator-() const { return MI(MOD - v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(int n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; using mi = MI<MOD>; using vmi = vector<mi>; /* global variables */ vmi fs, invfs; mi e2s[MAX_N + 1], e3s[MAX_N + 1], invds[MAX_N + 1]; mi ps[MAX_N + 1], xs[MAX_N + 1]; /* subroutines */ inline mi nck(int n, int k) { // nCk % MOD if (n < k || k < 0) return 0; return fs[n] * invfs[n - k] * invfs[k]; } inline mi nhk(int n, int k) { // nHk % MOD return nck(n + k - 1, k); } void prepare_fracs(int n) { fs.resize(n + 1), invfs.resize(n + 1); fs[0] = invfs[0] = 1; for (int i = 1; i <= n; i++) { fs[i] = fs[i - 1] * i; invfs[i] = fs[i].inv(); } } /* main */ int main() { int n; scanf("%d", &n); prepare_fracs(n); e2s[0] = e3s[0] = 1; for (int i = 1; i <= n; i++) e2s[i] = e2s[i - 1] * 2, e3s[i] = e3s[i - 1] * 3; for (int i = 2; i <= n; i++) invds[i] = (e2s[i] - 2).inv(); ps[n] = 1; for (int i = n - 1; i >= 1; i--) for (int j = i + 1; j <= n; j++) ps[i] += ps[j] * nck(j, i) * invds[j]; mi s = 0; for (int i = n; i >= 1; i--) { xs[i] = s; s += ps[i] * e3s[i - 1] * invds[i]; } for (int i = 1; i <= n; i++) printf("%d%c", (int)xs[i], (i < n) ? ' ' : '\n'); return 0; }