bzoj P5093 [Lydsy1711月赛]图的价值

背景：

貌似$\text{luogu}$没有这道题，于是来巨慢的$\text{bzoj}$交了。

题目传送门：

https://www.lydsy.com/JudgeOnline/problem.php?id=5093

题意：

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

思路：

大水题。
发现每一个点的贡献是互相独立的，因此我们考虑先算一个点的贡献，最后$\ast n$即可。
因为每一个点是互相独立的，因此剩余的$n-1$个点的连边与我无关，因此贡献为$2^{\frac{(n-1)(n-2)}{2}}$。
考虑每一个点与其他点的连边数$i$。
从$n-1$个点中选出$i$个点（度数就为$i$）的方案为$C_{n-1}^i$，乘上度数的贡献$i^k$即可。
所以：
$$ans=n\ast 2^{\frac{(n-1)(n-2)}{2}}\sum _{i=0}^{n-1}{C}_{n-1}^i{i}^k$$

$$ans=n\ast 2^{\frac{(n-1)(n-2)}{2}}\sum _{i=1}^{n-1}{C}_{n-1}^i{i}^k$$
我们发现后面的式子跟$\text{CF932E Team Work}$是类似的 ，因此直接抄即可 。
当然这道题要$\text{NTT}$处理第二类斯特林数。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
const int mod=998244353,G=3,inv_G=332748118;
using namespace std;
	int r[1000010],a[1000010],b[1000010],f[1000010],g[1000010],inv[1000010],Inv[1000010];
	int n,k,limit,l,ans=0;
int dg(int x,LL k)
{
	if(!k) return 1;
	int op=dg(x,k>>1);
	if(k&1) return (LL)op*op%mod*x%mod; else return (LL)op*op%mod;
}
int get_inv(int x)
{
	return dg(x,mod-2);
}
void init(int n)
{
	limit=1,l=0;
	while(limit<(n<<1))
		limit<<=1,l++;
	for(int i=1;i<limit;i++)
		r[i]=((r[i>>1]>>1)|((i&1)<<(l-1)));
}
void NTT(int *now,int limit,int op)
{
	for(int i=0;i<limit;i++)
		if(i<r[i]) swap(now[i],now[r[i]]);
	for(int mid=1;mid<limit;mid<<=1)
	{
		int wn=dg(op==1?G:inv_G,(mod-1)/(mid<<1));
		for(int j=0;j<limit;j+=(mid<<1))
		{
			int w=1;
			for(int k=0;k<mid;k++,w=((LL)w*wn)%mod)
			{
				int x=now[j+k],y=(LL)w*now[j+k+mid]%mod;
				now[j+k]=(x+y)%mod;
				now[j+k+mid]=(x-y+mod)%mod;
			}
		}
	}
}
void dft(int *f,int n,int limit)
{
	NTT(f,limit,-1);
	int INV=get_inv(limit);
	for(int i=0;i<n;i++)
		f[i]=(LL)f[i]*INV%mod;
}
void Init()
{
	inv[0]=inv[1]=1;
	for(int i=2;i<=k;i++)
		inv[i]=((LL)mod-mod/i)*inv[mod%i]%mod;

	Inv[0]=Inv[1]=1;
	for(int i=2;i<=k;i++)
		Inv[i]=(LL)Inv[i-1]*inv[i]%mod;
}
int main()
{
	scanf("%d %d",&n,&k);
	Init();
	for(int i=0;i<=k;i++)
	{
		f[i]=(((i&1)?-1ll:1ll)*Inv[i]+mod)%mod;
		g[i]=(LL)dg(i,k)*Inv[i]%mod;
	}
	init(k+1);
	NTT(f,limit,1),NTT(g,limit,1);
	for(int i=0;i<limit;i++)
		f[i]=(LL)f[i]*g[i]%mod;
	dft(f,k+1,limit);
	int tmp=n-1;
	for(int i=1;i<=min(n-1,k);i++)
	{
		ans=((LL)ans+(LL)f[i]*tmp%mod*dg(2,n-1-i)%mod)%mod;
		tmp=(LL)tmp*(n-1-i)%mod;
	}
	printf("%d",(LL)n*dg(2,(LL)(n-1)*(n-2)/2)%mod*ans%mod);
}