题目
求
f(n)
模
998244353(7×17×223+1)
看到模数
这么奇怪的模数,很容易想到ntt啦
正解
首先我们看第二类斯特林数的意义
S(n,m)
:将n个不同的元素拆分成m个集合的方案数。
那么乘上
m!
即
(S(n,m)×m!)
的意义不过就是让每个集合有了区别而已。
然后那么根据意义,我们可以用容斥原理得到:
S(n,m)×m!=∑mi=0(−1)k×(mk)×(m−k)n
详情见百度百科【第二类斯特林数的展开式】,链接:http://baike.baidu.com/link?url=Y-Wz-m8oEmw44OP1F5qi-blDPQwf8rWllPqnOqD1mdmkDsByyZQQZFR7DkkAoqOwsWJcWtz_3enoGtJmXU9dy-5YiL9b7dxlQsSKMGPaH_tVywqvNq48w_qXOhliVWFUH05Yt1ts_eJlB4kZllFukp5qnoq3SVDelpnD1t234saZPH-1mB08G-T0p1ZoktzSpqdc6UEYJi9FzK1_US9ru_
再看原来我们要求的式子:
∑ni=0∑ij=0S(i,j)×2j×(j!)
∵
当
j>i
时
S(i,j)=0
∴∑ni=0∑ij=0S(i,j)×2j×(j!)
=∑ni=0∑nj=0S(i,j)×2j×(j!)
=∑ni=0∑nj=02j×∑jk=0(−1)j−k×(jk)×ki
=∑nj=02j×∑jk=0(−1)j−k×j!k!(j−k)!×∑ni=0ki
=∑nj=02j×j!×∑jk=0(−1)j−k(j−k)!×∑ni=0kik!
很明显,后面的部分
∑jk=0(−1)j−k(j−k)!×∑ni=0kik!
即为一个卷积形式。
设
ai=(−1)ii!,bi=∑nk=0iki!,ci=∑ik=0ak×bi−k
求出c之后,那么
ANS=∑nj=02j×j!×cj
完美解决,时间复杂度
O(nlog2n)
贴个代码:
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long LL;
typedef double db;
int get(){
char ch;
int s=0;
bool pd=0;
while(ch=getchar(),(ch<'0'||ch>'9')&&ch!='-');
if (ch=='-')pd=1;
else s=ch-'0';
while(ch=getchar(),ch>='0'&&ch<='9')s=s*10+ch-'0';
if (pd)return -s;
return s;
}
const int N = 270000;
const int mo = 998244353;
const int G = 3;
int n,flen;
LL js[N],inv[N],be[N],a[N],b[N],c[N],ans,rev[N],w[N],v[N];
LL quickmi(LL x,LL tim){
LL ans=1;
while(tim){
if (tim%2)ans=ans*x%mo;
x=x*x%mo;
tim/=2;
}
return ans;
}
void set(){
fo(i,0,flen-1){
int cur=0;
for(int j=1,x=i;j<flen;j<<=1,x>>=1)cur=(cur<<1)+(x&1);
rev[i]=cur;
}
w[0]=1;w[1]=quickmi(G,(mo-1)/flen);
fo(i,2,flen)w[i]=w[i-1]*w[1]%mo;
}
void dft(LL *a,int sig){
fo(i,0,flen-1)v[rev[i]]=a[i];
for(int i=2,tot=flen,hal,cur;i<=flen;i<<=1){
tot>>=1,cur=sig?flen:0,hal=i>>1;
for(int j=0;j<hal;j++){
for(int k=j,wi=w[cur];k<flen;k+=i){
LL u=v[k],r=v[k+hal]*wi%mo;
v[k]=(u+r)%mo;
v[k+hal]=(u+mo-r)%mo;
}
sig?cur-=tot:cur+=tot;
}
}
fo(i,0,flen-1)a[i]=v[i];
}
void ntt(){
set();
LL invl=quickmi(flen,mo-2);
dft(a,0);dft(b,0);
fo(i,0,flen-1)c[i]=a[i]*b[i]%mo;
dft(c,1);
fo(i,0,flen-1)c[i]=c[i]*invl%mo;
}
int main(){
n=get();
flen=1;
while(flen<=n*2)flen<<=1;
js[0]=1;
fo(i,1,flen)js[i]=js[i-1]*i%mo;
inv[flen]=quickmi(js[flen],mo-2);
fd(i,flen-1,0)inv[i]=quickmi(js[i],mo-2);
b[0]=1;
b[1]=n+1;
fo(k,2,n)b[k]=(quickmi(k,n+1)-1)*quickmi(k-1,mo-2)%mo*inv[k]%mo;
fo(i,0,n)a[i]=i%2?(-inv[i]+mo)%mo:inv[i];
ans=0;
LL mi=1,jc=1;
ntt();
ans=0;
mi=1,jc=1;
fo(i,0,n){
ans=(ans+mi*jc%mo*c[i]%mo)%mo;
mi=mi*2%mo;
jc=jc*(i+1)%mo;
}
printf("%lld\n",ans);
return 0;
}