結果

問題 No.3045 反復重み付き累積和
コンテスト
ユーザー kwm_t
提出日時 2025-03-02 15:27:12
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 385 ms / 5,000 ms
コード長 1,594 bytes
コンパイル時間 6,734 ms
コンパイル使用メモリ 333,512 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2025-03-02 15:27:24
合計ジャッジ時間 11,546 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
	}
	vector<mint>ifrac(n + 1, 1);
	rep2(i, 1, n + 1)ifrac[i] = ifrac[i - 1] * mint(i).inv();
	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;
			mint tmp = 1;
			rep(i, n) {
				// b[i] = binom(i + t - 1, i) *  pw;
				if (i != 0)tmp *= i + t - 1;
				b[i] = tmp * ifrac[i] * pw;
				pw *= k;
			}
			a = convolution(a, b);
			a.resize(n);
		}
		else {
			int x; cin >> x;
			x--;
			cout << a[x].val() << endl;
		}
	}
	return 0;
}
0