【省选模拟】20/04/09 黑暗（生成函数）（多项式）

第一次考场自己搞生成函数，算是进步吧
好像可以用斯特林推，不过有一种无脑的方法

• 题意：$n$ 个点的所有无向图，每个图的贡献是连通块大小的 $m$ 次方，多次询问 $n\le 3e4,m\le 15$

• 先预处理 $F(x)$ 表示无向联通图的 $EGF$
  设 $G_m(x)$ 为答案的 $EGF$，那么有
  $$G_m(x)=\sum _{i\ge 0}\frac{F(x{)}^i\ast i^m}{i!}$$
  考虑递推，稍微推导发现
  $$G_m(x)=F(x)\frac{\text{d}{G}_{m-1}(x)}{\text{d}F(x)},G_0=\exp F(x)$$
  现在知道 $\frac{\text{d}G}{\text{d}x}$，那么我们可以预处理 $\frac{\text{d}x}{\text{d}F}$，然后就做完了，$O(nm\log n)$，好像不用推式子，询问可以 $O(1)$

• 题解的做法是考虑 $n^m={\sum }_{i=1}^m{S}_{m,i}\left(\genfrac{}{}{0ex}{}{n}{i}\right)i!$，然后让一个图贡献 $\left(\genfrac{}{}{0ex}{}{n}{i}\right)i!$，设这个是 $f_{i,n}$，那么有 $f_{i,n}={\sum }_j{g}_j\left(\genfrac{}{}{0ex}{}{n}{j}\right)f_{i-1,n-j}$，$g$ 是无向联通图个数，转移系数是因为 $i$ 个联通块是有序的
  不过这个询问要 $O(m)$

#include<bits/stdc++.h>
#define cs const
#define poly vector<int>
#define pb push_back
using namespace std;
int read(){
	int cnt=0, f=1; char ch=0;
	while(!isdigit(ch)){ ch=getchar(); if(ch=='-') f=-1; }
	while(isdigit(ch)) cnt=cnt*10+(ch-'0'), ch=getchar();
	return cnt*f;
}
typedef long long ll;
cs int Mod = 998244353;
int add(int a, int b){ return a + b >= Mod ? a + b - Mod : a + b; }
int dec(int a, int b){ return a - b < 0 ? a - b + Mod : a - b; }
int mul(int a, int b){ return 1ll * a * b % Mod; }
void Add(int &a, int b){ a = add(a,b); }
void Dec(int &a, int b){ a = dec(a,b); }
void Mul(int &a, int b){ a = mul(a,b); }
int ksm(int a, int b){ int as=1; for(;b;b>>=1,a=mul(a,a)) if(b&1) as=mul(as,a); return as; }
cs int K = 17, M = 1 << K | 5;
int fc[M], ifc[M], iv[M], T, n, m;
int up, bit; poly rev;
poly w[K+1];
poly F[16];
void NTT_init(){
	for(int i=1; i<=K; i++) w[i].resize(1<<(i-1)); 
	int wn=ksm(3,(Mod-1)/(1<<K)); w[K][0]=1;
	for(int i=1; i<(1<<(K-1)); i++) w[K][i]=mul(w[K][i-1],wn);
	for(int i=K-1;i;i--) for(int j=0; j<(1<<(i-1)); j++) w[i][j]=w[i+1][j<<1];
	iv[0]=iv[1]=1; for(int i=2; i<=(1<<K); i++) iv[i]=mul(Mod-Mod/i,iv[Mod%i]);
}
void Fac_init(int n){
	fc[0]=fc[1]=ifc[0]=ifc[1]=1;
	for(int i=2; i<=n; i++) fc[i]=mul(fc[i-1],i);
	for(int i=2; i<=n; i++) ifc[i]=mul(ifc[i-1],iv[i]);
}
void init(int deg){
	up=1; bit=0; while(up<deg) up<<=1,++bit; rev.resize(up);
	for(int i=0; i<up; i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(bit-1));
}
void NTT(poly &a, int typ=1){
	for(int i=0; i<up; i++) if(i<rev[i]) swap(a[i],a[rev[i]]);
	for(int i=1,l=1; i<up; i<<=1,++l)
	for(int j=0; j<up; j+=(i<<1))
	for(int k=0; k<i; k++){
		int x=a[k+j], y=mul(w[l][k],a[k+j+i]);
		a[k+j]=add(x,y); a[k+j+i]=dec(x,y);
	}
	if(typ==-1){
		reverse(a.begin()+1,a.end());
		for(int i=0; i<up; i++) Mul(a[i],iv[up]);
	}
}
poly operator * (poly a, poly b){
	int deg=a.size()+b.size()-1; 
	if(a.size()<=16||b.size()<=16){
		poly c(deg,0); 
		for(int i=0; i<(int)a.size(); i++)
		for(int j=0; j<(int)b.size(); j++)
		Add(c[i+j],mul(a[i],b[j])); return c;
	}
	init(deg); a.resize(up); b.resize(up); NTT(a); NTT(b);
	for(int i=0; i<up; i++) Mul(a[i],b[i]); NTT(a,-1);
	a.resize(deg); return a;
}
poly inv(poly a, int deg){
	poly b(1,ksm(a[0],Mod-2)),c;
	for(int lim=2; (lim>>1)<deg; lim<<=1){
		c=a; c.resize(lim); init(lim<<1);
		c.resize(up); b.resize(up); NTT(c); NTT(b);
		for(int i=0; i<up; i++) Mul(b[i],dec(2,mul(b[i],c[i])));
		NTT(b,-1); b.resize(lim); 
	} b.resize(deg); return b;
}
poly integ(poly a){
	a.pb(0);
	for(int i=a.size()-1; i>=1; i--) a[i]=mul(a[i-1],iv[i]);
	a[0]=0; return a;
}
poly deriv(poly a){
	for(int i=0; i+1<(int)a.size(); i++) a[i]=mul(a[i+1],i+1);
	a.pop_back(); return a;
}
poly ln(poly a, int deg){
	a=integ(inv(a,deg)*deriv(a)); a.resize(deg);
	return a;
}
void work(int n, int m){
	poly f; for(int i=0; i<=n; i++){
		int w=((ll)i*(i-1)>>1)%(Mod-1);
		f.pb(mul(ksm(2,w),ifc[i]));
	} F[0]=f; f=ln(f,n+1);
	poly df=inv(deriv(f),n+1);
	for(int i=1; i<=m; i++){
		F[i]=f*deriv(F[i-1]); F[i].resize(n+1);
		F[i]=F[i]*df;
	} 
}
int main(){
	freopen("dark.in","r",stdin);
	freopen("dark.out","w",stdout);
	NTT_init(); 
	Fac_init(3e4+5);
	work(3e4+3,15);
	T=read(); while(T--){
		n=read(), m=read();
		assert(n<(int)F[m].size());
		cout<<mul(fc[n],F[m][n])<<'\n';
	}
}