結果
問題 |
No.3174 勝ち残りじゃんけん
|
ユーザー |
![]() |
提出日時 | 2025-06-06 21:59:53 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 134 ms / 2,000 ms |
コード長 | 2,873 bytes |
コンパイル時間 | 3,705 ms |
コンパイル使用メモリ | 303,040 KB |
実行使用メモリ | 11,344 KB |
最終ジャッジ日時 | 2025-06-06 21:59:59 |
合計ジャッジ時間 | 5,868 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
コンパイルメッセージ
main.cpp:49:6: warning: ‘{anonymous}::mint {anonymous}::perm(int, int)’ defined but not used [-Wunused-function] 49 | mint perm(int n, int k) { | ^~~~ main.cpp:48:6: warning: ‘{anonymous}::mint {anonymous}::fact(int)’ defined but not used [-Wunused-function] 48 | mint fact(int n) {return Fact[n];} | ^~~~ main.cpp:44:6: warning: ‘{anonymous}::mint {anonymous}::icomb(int, int)’ defined but not used [-Wunused-function] 44 | mint icomb(int n, int k) { | ^~~~~
ソースコード
#include<bits/stdc++.h> namespace { #pragma GCC diagnostic ignored "-Wunused-function" #include<atcoder/all> #pragma GCC diagnostic warning "-Wunused-function" using namespace std; using namespace atcoder; #define rep(i,n) for(int i = 0; i < (int)(n); i++) #define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--) #define all(x) begin(x), end(x) #define rall(x) rbegin(x), rend(x) template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; } template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; } using ll = long long; using P = pair<int,int>; using VI = vector<int>; using VVI = vector<VI>; using VL = vector<ll>; using VVL = vector<VL>; using mint = modint998244353; constexpr int FACT_SIZE = 1000000; mint Fact[FACT_SIZE + 1]; mint iFact[FACT_SIZE + 1]; const auto fact_init = [] { Fact[0] = mint::raw(1); for(int i = 1; i <= FACT_SIZE; ++i) { Fact[i] = Fact[i-1] * i; } iFact[FACT_SIZE] = Fact[FACT_SIZE].inv(); for(int i = FACT_SIZE; i; --i) { iFact[i-1] = iFact[i] * i; } return false; }(); mint comb(int n, int k) { if (k == 0) return mint::raw(1); assert(n >= 0 && k >= 0); if (k > n) return mint::raw(0); return Fact[n] * iFact[n - k] * iFact[k]; } mint icomb(int n, int k) { return iFact[n] * Fact[n - k] * Fact[k]; } mint fact(int n) {return Fact[n];} mint perm(int n, int k) { assert(0 <= n); return Fact[n] * iFact[n - k]; } template <class T> struct cntsum { T a, b; // a+bx mod x^2 friend cntsum operator+(const cntsum& lhs, const cntsum& rhs) { return {lhs.a + rhs.a, lhs.b + rhs.b}; } friend cntsum operator-(const cntsum& lhs, const cntsum& rhs) { return {lhs.a - rhs.a, lhs.b - rhs.b}; } friend cntsum operator*(const cntsum& lhs, const cntsum& rhs) { return {lhs.a * rhs.a, lhs.a * rhs.b + lhs.b * rhs.a}; } friend cntsum operator*(const T& c, const cntsum& x) { return {x.a, x.b + c * x.a}; } cntsum& operator+=(const cntsum& x) { return *this = *this + x; } cntsum& operator-=(const cntsum& x) { return *this = *this - x; } cntsum& operator*=(const cntsum& x) { return *this = *this * x; } }; using CS = cntsum<mint>; } int main() { ios::sync_with_stdio(false); cin.tie(0); int n; cin >> n; vector<CS> dp(n + 1); dp[n] = {1, 0}; vector<mint> ans(n + 1); mint esm; for (int i = n; i >= 1; i--) { ans[i] = esm; esm -= dp[i].b; if (i == 1) break; mint p_suc = 3 * (mint(2).pow(i) - 2) / mint(3).pow(i); auto v = p_suc.inv() * dp[i]; esm += v.b; mint denom = (mint(2).pow(i) - 2).inv(); for (int j = 1; j < i; j++) { mint num = comb(i, j); auto w = v; w.a *= num * denom; w.b *= num * denom; dp[j] += w; } } for (int i = 1; i <= n; i++) cout << ans[i].val() << " \n"[i == n]; }