Tjoi2016&Heoi2016 求和

题目

求$f(n)$模$998244353(7×17×2^{23}+1)$

看到模数

这么奇怪的模数，很容易想到ntt啦

正解

首先我们看第二类斯特林数的意义$S(n,m)$：将n个不同的元素拆分成m个集合的方案数。
那么乘上$m!$即$(S(n,m)×m!)$的意义不过就是让每个集合有了区别而已。
然后那么根据意义，我们可以用容斥原理得到：$S(n,m)×m!=\sum_{i=0}^m(-1)^k×(^m_k)×(m-k)^n$
详情见百度百科【第二类斯特林数的展开式】，链接：http://baike.baidu.com/link?url=Y-Wz-m8oEmw44OP1F5qi-blDPQwf8rWllPqnOqD1mdmkDsByyZQQZFR7DkkAoqOwsWJcWtz_3enoGtJmXU9dy-5YiL9b7dxlQsSKMGPaH_tVywqvNq48w_qXOhliVWFUH05Yt1ts_eJlB4kZllFukp5qnoq3SVDelpnD1t234saZPH-1mB08G-T0p1ZoktzSpqdc6UEYJi9FzK1_US9ru_
再看原来我们要求的式子：
$\sum_{i=0}^n\sum_{j=0}^iS(i,j)×2^j×(j!)$
$\because$当$j>i$时$S(i,j)=0$
$\therefore\sum_{i=0}^n\sum_{j=0}^iS(i,j)×2^j×(j!)$
$=\sum_{i=0}^n\sum_{j=0}^nS(i,j)×2^j×(j!)$
$=\sum_{i=0}^n\sum_{j=0}^n2^j×\sum_{k=0}^j(-1)^{j-k}×(^j_k)×k^i$
$=\sum_{j=0}^n2^j×\sum_{k=0}^j(-1)^{j-k}×\frac{j!}{k!(j-k)!}\times \sum_{i=0}^nk^i$
$=\sum_{j=0}^n2^j\times j!\times \sum_{k=0}^j\frac{(-1)^{j-k}}{(j-k)!}\times \frac{\sum_{i=0}^nk^i}{k!}$
很明显，后面的部分$\sum_{k=0}^j\frac{(-1)^{j-k}}{(j-k)!}\times \frac{\sum_{i=0}^nk^i}{k!}$即为一个卷积形式。
设$a_i=\frac{(-1)^i}{i!}，b_i=\frac{\sum_{k=0}^ni^k}{i!}，c_i=\sum_{k=0}^ia_k\times b_{i-k}$
求出c之后，那么$ANS=\sum_{j=0}^n2^j\times j!\times c_j$
完美解决，时间复杂度$O(n log_2n)$

贴个代码:

#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;
}