結果
問題 | 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;}