結果
問題 |
No.2369 Some Products
|
ユーザー |
👑 ![]() |
提出日時 | 2023-06-30 21:26:51 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 457 ms / 2,500 ms |
コード長 | 1,106 bytes |
コンパイル時間 | 832 ms |
コンパイル使用メモリ | 82,464 KB |
最終ジャッジ日時 | 2025-02-15 03:29:28 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 14 |
ソースコード
#include <iostream> #include <string> #include <vector> #include <algorithm> #include <atcoder/modint> using namespace std; using i32 = int; using u32 = unsigned int; using i64 = long long; using u64 = unsigned long long; #define rep(i,n) for(int i=0; i<(int)(n); i++) const i64 INF = 1001001001001001001; using Modint = atcoder::static_modint<998244353>; int main(){ int N; cin >> N; vector<Modint> P(N); rep(i,N){ int a; cin >> a; P[i] = a; } vector<vector<Modint>> inv(N+1); inv[0].resize(N+1); inv[0][0] = 1; rep(i,N){ auto x = - P[i]; inv[i+1] = inv[i]; rep(j,N) inv[i+1][j+1] += inv[i+1][j] * x; } vector<vector<Modint>> prod(N+1); prod[0].resize(N+1); prod[0][0] = 1; rep(i,N){ auto x = P[i]; prod[i+1] = prod[i]; for(int j=N-1; j>=0; j--) prod[i+1][j+1] += prod[i+1][j] * x; } int Q; cin >> Q; rep(q,Q){ int l,r,k; cin >> l >> r >> k; l--; Modint ans = 0; for(int j=0; j<=k; j++) ans += prod[r][j] * inv[l][k-j]; cout << ans.val() << '\n'; } return 0; }