bzoj5093图的价值：多项式，斯特林数（二项式反演）

Description

“简单无向图”是指无重边、无自环的无向图（不一定连通）。 一个带标号的图的价值定义为每个点度数的k次方的和。 给定n和k，请计算所有n个点的带标号的简单无向图的价值之和。

因为答案很大，请对998244353取模输出。$n le 10^9,k le 200000$

化学学考时含义推式子+手动打表找规律得到了一个$O(nlogn)$的式子开心的很我以为我要AC了回来看数据范围就升天了。

问NC大神这题用到了什么：斯特林数/伯努利数。然后就自闭了学了一天的知识点还去做了点伯努利数的别的题。

然而最后没想到伯努利数的做法，是用第二类斯特林数解决的。

先列出来答案的式子，分别考虑每个点的贡献：

$ans=2^{frac{(n-1)(n-2)}{2}} imes n imes sumlimits_{i=0}^{n} C_{n-1}^i imes i^k$

就是首先在n个点里选出你当前考虑贡献答案的是哪个点，然后在剩余的n-1个点的子图里联通情况随意就是$2^{frac{(n-1)(n-2)}{2}}$

然后从n-1个点里选出i个与当前贡献答案的点连边，选法是组合数，贡献的价值是$i^k$。

可以发现，这里有个求自然数幂和乘个组合数的形式，而自然数幂和的范围很大不能直接求解。

动用学过的（和刚学的）所有数学知识点看有什么能把幂拆成可以接受的形式。

1，伯努利数：$sumlimits_{i=1}^{n} i^m = frac{sumlimits_{i=1}^{m} B_i imes C_{m+1}^{i} imes n^{m+1-i} }{m+1}$

这不妨是个好主意，但是这道题里要求的是自然数幂乘组合数后再求和，并不能直接适用，所以至少我现在还不会做。

等会上网学习一下再来fix这里。

学习完了，网上并没有这样的做法，但是伯努利数和第二类斯特林数可以互相转化。。。

对这题并没有什么帮助。。。

2，第二类斯特林数：$x^n=sumlimits_{i=0}^{n} left { ^n_k ight } imes frac{x!}{(x-i)!}$

这样的话虽然不能直接求和，但是好歹可以把难以解决的高次幂转化为经常处理的阶乘什么的了。

首先说说关于斯特林数怎么求：先弄明白含义，第二类斯特林数$left { ^n_k ight }$表示把n个元素放进k个无标号非空集合的方案数。

然后我们再弄几个类似的定义：设f(n,k)表示把n个元素放进k个有编号非空集合的方案数，g(n,k)表示n个元素放进k个有标号可空集合的方案数。

对于是否有编号的问题，其实就是多了个内部顺序，所以$left { ^n_k ight }=frac{f(n,k)}{k!}$

然后g其实是可以直接求解的：$g(n,k)=k^n$

然后至于f和g之间的关系，我们可以通过枚举它具体有几个集合是空的，得到$g(n,k)=sumlimits_{i=0}^{k} C_k^i imes f(n,i)$

就是选出到底哪几个不是空的，在这些非空的里面进行分配。

然后发现这个形式就是二项式反演的基本形式，所以反演它：$f(n,k)=sumlimits_{i=0}^{k} (-1)^{k-i} imes C_k^i imes g(n,i)$

这样我们就能得到f，进而求解出第二类斯特林数。

总的式子就是$left { ^n_k ight }=sumlimits_{i=0}^{k} frac{(-1)^{k-i} imes i^k}{(k-i)! imes i!}$

这样的话单点求第二类斯特林数就是$O(n)$的了。而上式是个卷积形式，所以求一行第二类斯特林数就可以用NTT在$O(nlogn)$之内解决。

回到这道题来，求解斯特林数的复杂度是$O(klogk)$可以接受，那么现在答案的式子是什么？

$ans=2^{frac{(n-1)(n-2)}{2}} imes n imes sumlimits_{i=0}^{n} C_{n-1}^{i} sumlimits_{j=0}^{k} left { ^k_j ight } imes frac{i!}{(i-j)!}$

把各种组合数拆一拆放在一起消一消，得到：$ans=2^{frac{(n-1)(n-2)}{2}} imes n! imes sumlimits_{i=0}^{n} sumlimits_{j=0}^{k} left { ^k_j ight } imes frac{1}{(i-j)! imes (n-i-1)!}$

换一下ij的枚举顺序，发现最后那个分母上两个阶乘的底之和是一个与i无关的定值，上下都乘一个$(n-j-1)!$合成一个组合数。它等于$frac{C_{n-j-1}^{i-j}}{(n-j-1)!}$

然后再套上外面那个求和。。。一行组合数求和得到的是2的次幂，所以这个和式就可以扔掉了得到下面形式：

$ans=2^{frac{(n-1)(n-2)}{2}} imes n! imes sumlimits_{j=0}^{k} frac{left { ^k_j ight } imes 2^{n-j-1}}{(n-j-1)!} $

现在的瓶颈就在于式子中的两个阶乘不能在合理复杂度内求出。

然而它们一个在分子上，一个在分母上，而且它们底数的差值在k范围内。。。

抵消啊！就剩下一个下降幂了！

所以最后可做的式子就是$2^{frac{(n-1)(n-2)}{2}} sumlimits_{j=0}^{k} left { ^k_j ight } imes 2^{n-j-1} imes n^{underline{j+1}}$

#include<cstdio>
#define int long long
#define mod 998244353
#define S 555555
int n,k,a[S],s[S],len=1,fac[S],inv[S],rev[S],dw[S],ans;
int pow(int b,int t,int a=1){for(;t;t>>=1,b=b*b%mod)if(t&1)a=a*b%mod;return a;}
void NTT(int *a,int opt){
    for(int i=1;i<len;++i)if(rev[i]>i)a[i]^=a[rev[i]]^=a[i]^=a[rev[i]];
    for(int mid=1;mid<len;mid<<=1)
        for(int i=0,t=pow(3,(mod-1)/2/mid*opt+mod-1);i<len;i+=mid<<1)
            for(int j=i,x,y,w=1;j<mid+i;++j,w=w*t%mod)
                x=a[j],y=a[j+mid]*w%mod,a[j]=(x+y)%mod,a[j+mid]=(x-y+mod)%mod;
    if(opt==-1)for(int i=0,iv=pow(len,mod-2);i<len;++i)a[i]=a[i]*iv%mod;
}
signed main(){
    scanf("%lld%lld",&n,&k);fac[0]=dw[0]=1;
    while(len<=k<<1)len<<=1;
    for(int i=1;i<len;++i)rev[i]=rev[i>>1]>>1|(i&1?len>>1:0);
    for(int i=1;i<=k;++i)fac[i]=fac[i-1]*i%mod;
    inv[k]=pow(fac[k],mod-2);
    for(int i=k-1;~i;--i)inv[i]=inv[i+1]*(i+1)%mod;
    for(int i=0;i<=k;++i)a[i]=pow(i,k)*inv[i]%mod,s[i]=(i&1?mod-1:1)*inv[i]%mod;
    NTT(a,1);NTT(s,1);for(int i=0;i<len;++i)s[i]=s[i]*a[i]%mod;NTT(s,-1);
    for(int i=1;i<=k+1;++i)dw[i]=dw[i-1]*(n+1-i)%mod;
    for(int i=0;i<=k;++i)ans=(ans+dw[i+1]*pow(2,n-i-1)%mod*s[i])%mod;
    printf("%lld
",ans*pow(2,(n-1)*(n-2)/2)%mod);
}