码迷,mamicode.com
首页 > 其他好文 > 详细

[HEOI2016/TJOI2016]求和

时间:2018-03-24 20:41:13      阅读:140      评论:0      收藏:0      [点我收藏+]

标签:.com   时间   img   姐姐   tput   lin   inpu   printf   整数   

Discription

在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心。

现在他想计算这样一个函数的值:
技术分享图片
S(i, j)表示第二类斯特林数,递推公式为:
S(i, j) = j ? S(i ? 1, j) + S(i ? 1, j ? 1), 1 <= j <= i ? 1。
边界条件为:S(i, i) = 1(0 <= i), S(i, 0) = 0(1 <= i)
你能帮帮他吗?
Input

输入只有一个正整数

Output

 输出f(n)。

由于结果会很大,输出f(n)对998244353(7 × 17 × 223 + 1)取模的结果即可。

1 ≤ n ≤ 100000

Sample Input
3

Sample Output

87

 

 

   我们知道第二类斯特林数和排列数(下降幂)组合在一起可以表示n^k,又因为排列等于组合乘上一个阶乘,于是我们就可以开开心心的二项式反演,得到一个某一行(其实也可以很多行,鉴于这个式子的特殊性质,我们可以把不同行的同一列合并)某一列的斯特林数的表达式。

    具体的说,S(k,n) = Σ (i^k / i!) * ((-1)^(n-i) / (n-i)!)     [具体推导就不写了,就是一个二项式反演]。

这个式子的特殊性质太多了,首先它是一个卷积的形式,所以我们可以直接用NTT 在 N log N 的时间求出某一行的所有第二类斯特林数分别是多少;

并且只有 i^k 项和行数有关,所以同一列很好合并,于是这个题就做完了2333。

 

 

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=330005;
const int ha=998244353;
const int root=3,inv=ha/3+1;
int a[maxn],b[maxn],jc[maxn];
int r[maxn],N,M,n,INV,l;

inline int add(int x,int y){
	x+=y;
	return x>=ha?x-ha:x;
}

inline int ksm(int x,int y){
	int an=1;
	for(;y;y>>=1,x=x*(ll)x%ha) if(y&1) an=an*(ll)x%ha;
	return an;
}

inline void NTT(int *c,const int f){
	for(int i=0;i<N;i++) if(i<r[i]) swap(c[i],c[r[i]]);
	
	for(int i=1;i<N;i<<=1){
		int omega=ksm(f==1?root:inv,(ha-1)/(i<<1));
		for(int p=i<<1,j=0;j<N;j+=p){
			int now=1;
			for(int k=0;k<i;k++,now=now*(ll)omega%ha){
				int x=c[j+k],y=c[j+k+i]*(ll)now%ha;
				c[j+k]=add(x,y);
				c[j+k+i]=add(x,ha-y);
			}
		}
	}
	
	if(f==-1) for(int i=0;i<N;i++) c[i]=c[i]*(ll)INV%ha;
}

inline void init(){
	jc[0]=1;
	for(int i=1;i<=n;i++) jc[i]=jc[i-1]*(ll)i%ha;
	for(int i=0;i<=n;i++){
		if(!i) a[i]=1;
		else if(i==1) a[i]=n+1;
		else a[i]=add(ksm(i,n+1),ha-1)*(ll)ksm(add(i,ha-1)*(ll)jc[i]%ha,ha-2)%ha;
		if(i&1) b[i]=ha-ksm(jc[i],ha-2);
		else b[i]=ksm(jc[i],ha-2);
	}
	
	M=n<<1;
	for(N=1;N<=M;N<<=1) l++;
	for(int i=0;i<N;i++) r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
}

inline void solve(){
	NTT(a,1),NTT(b,1);
	for(int i=0;i<N;i++) a[i]=a[i]*(ll)b[i]%ha;
	INV=ksm(N,ha-2),NTT(a,-1);
}

inline void output(){
	int ans=0,base=1;
	for(int i=0;i<=n;i++,base=add(base,base)) ans=add(ans,a[i]*(ll)base%ha*(ll)jc[i]%ha);
	printf("%d\n",ans);
}

int main(){
	scanf("%d",&n);
	init();
	solve();
	output();
	return 0;
}

  

[HEOI2016/TJOI2016]求和

标签:.com   时间   img   姐姐   tput   lin   inpu   printf   整数   

原文地址:https://www.cnblogs.com/JYYHH/p/8640895.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!