結果
問題 |
No.3045 反復重み付き累積和
|
ユーザー |
![]() |
提出日時 | 2025-03-02 15:20:07 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,420 bytes |
コンパイル時間 | 6,453 ms |
コンパイル使用メモリ | 332,756 KB |
実行使用メモリ | 10,276 KB |
最終ジャッジ日時 | 2025-03-02 15:20:20 |
合計ジャッジ時間 | 13,185 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 4 TLE * 1 -- * 36 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; //using mint = modint1000000007; //const int mod = 1000000007; using mint = modint998244353; const int mod = 998244353; //const int INF = 1e9; //const long long LINF = 1e18; #define rep(i, n) for (int i = 0; i < (n); ++i) #define rep2(i,l,r)for(int i=(l);i<(r);++i) #define rrep(i, n) for (int i = (n) - 1; i >= 0; --i) #define rrep2(i,l,r)for(int i=(r) - 1;i>=(l);--i) #define all(x) (x).begin(),(x).end() #define allR(x) (x).rbegin(),(x).rend() #define P pair<int,int> template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; } template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; } mint binom(int p, int q) { mint x = 1, y = 1; rep(i, q) { x *= p - i; y *= i + 1; } return x / y; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n, q; cin >> n >> q; vector<mint>a(n); rep(i, n) { int x; cin >> x; a[i] = x; } while (q--) { int t; cin >> t; if (t == 1) { long long k, t; cin >> k >> t; // (1/(1-kx))^t vector<mint>b(n); mint pw = 1; rep(i, n) { b[i] = binom(i + t - 1, i) * pw; pw *= k; } a = convolution(a, b); a.resize(n); } else { int x; cin >> x; x--; cout << a[x].val() << endl; } } return 0; }