結果

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

ソースコード

diff # 

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