結果

問題 No.2369 Some Products
ユーザー vjudge1
提出日時 2025-05-30 16:56:49
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,001 bytes
コンパイル時間 1,665 ms
コンパイル使用メモリ 163,932 KB
実行使用メモリ 813,272 KB
最終ジャッジ日時 2025-05-30 16:56:58
合計ジャッジ時間 5,336 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other WA * 1 MLE * 1 -- * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define ll long long
using namespace std;
constexpr ll maxn=5e3+5,modd=998244353;
ll mul[maxn][maxn<<2],inv[maxn][maxn<<2],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*2]=inv[0][n*2]=1;
	mul[1][1]=p[1];
	for(int i=n*2+1;i<=n*4;i++)inv[1][i]=(inv[1][i-1]*(-p[1]));
	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]=(mul[i-1][j]+mul[i-1][j-1]*p[i])%modd;
		for(int j=0;j<=4*n;j++)inv[i][j]=(inv[i-1][j]*invp)%modd;
		for(int j=4*n;j>=0;j--)inv[i][j]=(inv[i][j]+inv[i][j+1]*(-invp))%modd;
		for(int j=0;j<=4*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=-2*n;i<=k;i++)ans=(ans+mul[r][k-i]*inv[l-1][i+2*n])%modd;
		cout<<(ans%modd+modd)%modd<<"\n";
	}
	return 0;
}
0