結果

問題 No.2369 Some Products
ユーザー vjudge1
提出日時 2025-05-30 17:25:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 933 ms / 2,500 ms
コード長 970 bytes
コンパイル時間 1,549 ms
コンパイル使用メモリ 163,856 KB
実行使用メモリ 394,856 KB
最終ジャッジ日時 2025-05-30 17:25:49
合計ジャッジ時間 8,425 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define ll long long
using namespace std;
constexpr ll maxn=5e3+5,modd=998244353;
int mul[maxn][maxn],inv[maxn][maxn*3];
ll n,q,p[maxn];
ll qpow(ll a,ll b){
	ll ret=1;
	while(b){
		if(b&1)ret=(ret*a)%modd;
		a=(a*a)%modd;
		b>>=1;
	}
	return ret;
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)cin>>p[i];
	mul[1][0]=inv[1][n]=inv[0][n]=1;
	mul[1][1]=p[1];
	for(int i=n+1;i<=3*n;i++)inv[1][i]=(1ll*inv[1][i-1]*(-p[1]))%modd;
	for(int i=2;i<=n;i++){
		mul[i][0]=1;
		ll invp=qpow(p[i],modd-2);
		for(int j=1;j<=n;j++)mul[i][j]=(1ll*mul[i-1][j]+1ll*mul[i-1][j-1]*p[i])%modd;
		for(int j=3*n;j>=0;j--)inv[i][j]=(1ll*inv[i-1][j]*invp+1ll*inv[i][j+1]*(-invp))%modd;
		for(int j=0;j<=3*n;j++)inv[i][j]=inv[i][j+1];
	}
	cin>>q;
	while(q--){
		ll l,r,k,ans=0;
		cin>>l>>r>>k;
		for(int i=-n+k;i<=k;i++)ans=(ans+1ll*mul[r][k-i]*inv[l-1][i+n])%modd;
		cout<<(ans%modd+modd)%modd<<"\n";
	}
	return 0;
}
0