【PKUWC2019模拟2019.1.15】Permutation

最新推荐文章于 2019-01-31 15:27:00 发布

YiPeng_Deng

最新推荐文章于 2019-01-31 15:27:00 发布

阅读量209

点赞数

分类专栏：题解计数 FFT和NTT 组合数学文章标签：计数 FFT和NTT 组合数学

本文链接：https://blog.csdn.net/qq_43649416/article/details/86516063

版权

题解同时被 3 个专栏收录

133 篇文章

订阅专栏

计数

32 篇文章

订阅专栏

FFT和NTT

8 篇文章

订阅专栏

Description

在这里插入图片描述

Data Constraint

在这里插入图片描述

一个简单的计数。

在这里插入图片描述

稍微化简一下发现我们需要求
(i-1)/2*(i-2)!/(i-y)!
考场上的时候想到了这里，也想到了NTT，就是不会转换！！！？？
我们可以这样卷：
- 设f[i]=(i-1)/2*(i-2)!,g[i]=1/(n-i)!
- 那么那个式子就是
- sigma(f[i]*g[j])[i+j==n+y]
- 妙啊妙…
但这只是nlogn，还有O（n）的方法，然鹅nlogn也能勉强水过这题。。。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define LL long long
#define maxn 2000005
#define mo 998244353
using namespace std;

int n,q,i,j,k,x,y,jc[maxn],ny[maxn];
int limit,a[maxn*2],b[maxn*2],bt[maxn*2];

LL ksm(LL x,LL y){
	LL s=1;
	for(;y;y/=2,x=x*x%mo) if (y&1) 
		s=s*x%mo;
	return s;
}

void ntt(int *a,int n,int sig){
	for(int i=0;i<n;i++) if (i<bt[i]) 
		swap(a[i],a[bt[i]]);
	for(int mid=1;mid<n;mid<<=1){
		int gn=ksm(3,(mo-1)/(mid<<1));
		if (sig<0) gn=ksm(gn,mo-2);
		for(int j=0;j<n;j+=(mid<<1)){
			int g=1;
			for(int k=0;k<mid;k++,g=(LL)g*gn%mo){
				int x=a[j+k],y=(LL)g*a[j+k+mid]%mo;
				a[j+k]=(x+y)%mo;
				a[j+k+mid]=(x-y+mo)%mo;
			}
		}
	}
}

int main(){
	freopen("permutation.in","r",stdin);
	freopen("permutation.out","w",stdout);
	scanf("%d%d",&n,&q);
	jc[0]=ny[0]=1;
	for(i=1;i<=n;i++) jc[i]=(LL)jc[i-1]*i%mo,ny[i]=(LL)ny[i-1]*ksm(i,mo-2)%mo;
	
	for(i=2;i<=n;i++) a[i]=(LL)(i-1)/2*jc[i-2]%mo,b[i]=ny[n-i];
	for(limit=1;limit<=n*2;limit<<=1);
	for(i=1;i<limit;i++) bt[i]=(bt[i>>1]>>1)|((i&1)?(limit>>1):0);
	ntt(a,limit,1),ntt(b,limit,1);
	for(i=0;i<limit;i++) a[i]=(LL)a[i]*b[i]%mo;
	ntt(a,limit,-1);
	int inv=ksm(limit,mo-2);
	for(i=0;i<limit;i++) a[i]=(LL)a[i]*inv%mo;
	
	for(i=1;i<=q;i++){
		scanf("%d%d",&x,&y);
		printf("%lld\n",(LL)jc[n-y]*a[n+y]%mo);
	}
}