Codeforces 512D

题面传送门

题意：
给出一张无向图，每次你可以选择一个度数 (leq 1) 的点并将其删除。
问对于 (k=0,1,2,dots,n)，有多少个删除 (k) 个点的序列，答案模 (10^9+9)。
(1leq nleq 100)

很容易发现，环上的点是永远不可能被删除的。于是我们对原图做一遍拓扑排序，求出哪些点可能被删除。
显然这些能被删除的点一定会形成一棵森林。
而森林中的每个连通块又可以分为两类：一种是该连通块在原图中属于的连通块中有环相连。另一种是该连通块在原图中本身就是一棵树。
对于第一种情况，一定有且仅有一个点 (x) 与环相邻。此时 (x) 一定是连通块中最后一个被删除的。我们就以 (x) 为根跑一遍树上背包。具体来说，设 (dp_{i,j}) 表示 (i) 子树中选了 (j) 个点的方案数。由于我们这是有根树，必须选完 (x) 子树内所有点才能选 (x)。转移就正常求个卷积，再乘上个组合数的系数即可。由于任意两点之间最多被合并一次，复杂度 (n^2)。
对于第二种情况，我们并没有给这棵树定根。所以我们考虑以树上每个点为根都跑一次树上背包。但这样会重复计算。对于选了 (i) 个点的方案，在以未被选择 (s-i) 个点为根情况下都被算了一次。故答案除以 (s-i)。
然后最后再总体跑一遍背包即可。时间复杂度 (mathcal O(n^3))

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define fz(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define ffe(it,v) for(__typeof(v.begin()) it=v.begin();it!=v.end();it++)
#define fill0(a) memset(a,0,sizeof(a))
#define fill1(a) memset(a,-1,sizeof(a))
#define fillbig(a) memset(a,63,sizeof(a))
#define pb push_back
#define ppb pop_back
#define mp make_pair
template<typename T1,typename T2> void chkmin(T1 &x,T2 y){if(x>y) x=y;}
template<typename T1,typename T2> void chkmax(T1 &x,T2 y){if(x<y) x=y;}
typedef pair<int,int> pii;
typedef long long ll;
template<typename T> void read(T &x){
	x=0;char c=getchar();T neg=1;
	while(!isdigit(c)){if(c=='-') neg=-1;c=getchar();}
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	x*=neg;
}
const int MAXN=100;
const int MAXM=4950;
const int MOD=1e9+9;
int qpow(int x,int e){int ret=1;for(;e;x=1ll*x*x%MOD,e>>=1) if(e&1) ret=1ll*ret*x%MOD;return ret;}
int n,m,deg[MAXN+5],to[MAXM*2+5],ec=0,nxt[MAXM*2+5],hd[MAXN+5];
void adde(int u,int v){to[++ec]=v;nxt[ec]=hd[u];hd[u]=ec;}
bool vis[MAXN+5];
int from[MAXN+5],cmp=0,siz[MAXN+5],rt[MAXN+5],con[MAXN+5],c[MAXN+5][MAXN+5];
void dfs0(int x){
	if(from[x]) return;from[x]=cmp;siz[cmp]++;
//	printf("%d %d
",x,cmp);
	for(int e=hd[x];e;e=nxt[e]) if(vis[to[e]]) dfs0(to[e]);
}
int dp[MAXN+5][MAXN+5],s[MAXN+5],tmp[MAXN+5],dpc[MAXN+5][MAXN+5];
void dfs(int x,int f){
//	printf("%d %d
",x,f);
	s[x]=1;dp[x][0]=1;
	for(int e=hd[x];e;e=nxt[e]){
		int y=to[e];if(y==f) continue;dfs(y,x);
		memset(tmp,0,sizeof(tmp));
		for(int i=0;i<=s[x];i++) for(int j=0;j<=s[y];j++)
			tmp[i+j]=(tmp[i+j]+1ll*dp[x][i]*dp[y][j]%MOD*c[i+j][i]%MOD)%MOD;
		s[x]+=s[y];
		for(int i=0;i<=s[x];i++) dp[x][i]=tmp[i];
	} dp[x][s[x]]=dp[x][s[x]-1];
//	for(int i=0;i<=s[x];i++) printf("%d %d %d
",x,i,dp[x][i]);
}
int iv[MAXN+5],ans[MAXN+5]; 
int main(){
	scanf("%d%d",&n,&m);
	for(int i=0;i<=MAXN;i++){
		c[i][0]=1;
		for(int j=1;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%MOD;
	}
	iv[0]=1;for(int i=1;i<=MAXN;i++) iv[i]=qpow(i,MOD-2);
	for(int i=1;i<=m;i++){
		int u,v;scanf("%d%d",&u,&v);
		adde(u,v);adde(v,u);deg[u]++;deg[v]++;
	}
	queue<int> q;
	for(int i=1;i<=n;i++) if(deg[i]<=1) q.push(i);
	while(!q.empty()){
		int x=q.front();q.pop();vis[x]=1;//printf("%d
",x);
		for(int e=hd[x];e;e=nxt[e]){
			int y=to[e];deg[y]--;
			if(deg[y]<=1&&!vis[y]) q.push(y);
		}
	}
	for(int i=1;i<=n;i++) if(vis[i]&&!from[i]) cmp++,dfs0(i);
	for(int i=1;i<=n;i++) for(int e=hd[i];e;e=nxt[e]) if(vis[i]&&!vis[to[e]]) con[i]=to[e],rt[from[i]]=i;
//	for(int i=1;i<=cmp;i++) printf("%d
",rt[i]);
	for(int i=1;i<=n;i++) if(vis[i]){
		if(rt[from[i]]&&i==rt[from[i]]){
			memset(dp,0,sizeof(dp));memset(s,0,sizeof(s));
			dfs(i,con[i]);
			for(int j=0;j<=siz[from[i]];j++) dpc[from[i]][j]=dp[i][j];
		} else if(!rt[from[i]]){
			memset(dp,0,sizeof(dp));memset(s,0,sizeof(s));
			dfs(i,0);//printf("%d
",i);
//			for(int j=1;j<=siz[from[i]];j++) printf("%d ",dp[i][j]);printf("
");
			for(int j=0;j<=siz[from[i]];j++) dpc[from[i]][j]=(dpc[from[i]][j]+dp[i][j])%MOD;
		}
	}
	for(int i=1;i<=cmp;i++) if(!rt[i]){
		for(int j=0;j<=siz[i];j++) dpc[i][j]=1ll*dpc[i][j]*iv[siz[i]-j]%MOD;
	} ans[0]=1;
//	for(int i=1;i<=cmp;i++) for(int j=0;j<=siz[i];j++) printf("%d %d %d
",i,j,dpc[i][j]); 
	for(int i=1;i<=cmp;i++) for(int j=n;j;j--) for(int k=1;k<=min(siz[i],j);k++)
		ans[j]=(ans[j]+1ll*ans[j-k]*dpc[i][k]%MOD*c[j][k])%MOD;
	for(int i=0;i<=n;i++) printf("%d
",ans[i]);
	return 0;
}