[ZROI322]C酱的表格 组合数学

[ZROI322]

  • 考虑贡献即可

 

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define per(i,a,b) for(int i=(a);i>=(b);i--)
using namespace std;
const int mod=998244353;
const int N=1e7+2e5;
int f[N],fact[N],ifact[N],inva[N],ans=0,n,m,a,b,p,q,x,y;
inline int read(){
	int num=0;char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))num=num*10+ch-'0',ch=getchar();
	return num;
}
inline void Print(int x){if(x>9)Print(x/10);putchar(x%10+'0');}
int power(int a,int b)
{int ans=1;for(;b;b>>=1,a=1LL*a*a%mod)if(b&1)ans=1LL*ans*a%mod;return ans;}

inline int C(int n,int m){
	if(n<0||m<0||n<m)return 0;
	return 1LL*fact[n]*ifact[m]%mod*ifact[n-m]%mod;
} 
void subtask2(int x,int y){
	int ans=0;
	rep(i,max(1,y-x+p),y)ans=(ans+1LL*C(x-p,y-i)*f[i]%mod*power(b,y-i)%mod*power(a,x-p-y+i)%mod)%mod;
	Print(ans);puts("");
}
void subtask1(int x,int y){
	int ans=0;
	rep(i,1,y)ans=(ans+1LL*f[i]*inva[p-x+y-i]%mod*power(-b,y-i)%mod*C(y-i+p-x-1,p-x-1)%mod+mod)%mod;
	Print(ans);puts("");
}
int main()
{

	m=read(),n=read(),a=read(),b=read(),p=read(),q=read();
	fact[0]=1;rep(i,1,n+m)fact[i]=1LL*i*fact[i-1]%mod;
	ifact[0]=ifact[1]=1;rep(i,2,n+m)ifact[i]=(-1LL*mod/i*ifact[mod%i]%mod+mod)%mod;
	rep(i,1,n+m)ifact[i]=1LL*ifact[i]*ifact[i-1]%mod;
	int fuck=power(a,mod-2);inva[0]=1;rep(i,1,n+m)inva[i]=1LL*fuck*inva[i-1]%mod;
	rep(i,1,n)f[i]=read();
	while(q--){
		x=read(),y=read();
		if(x>=p)subtask2(x,y);
		else subtask1(x,y);
	}
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值