Codeforces 512D

考虑一个简单环上的所有节点都不可以删除。我们可以从所有度$\le 1$的节点开始bfs出所有可以删除的节点。
那么可以删除的节点显然会形成若干个森林，有的树是无根的（最后被删除的节点没有限制），有的树会存在根（连到一个不可被删除的节点）。
考虑一棵树保留的节点显然是可能为空的连通块，那么就可以DP了。对于有根的树，我们设$F[x]$表示删空点$x$为根的子树的方案数，$G[x][i]$表示点$x$保留，在子树中选$i$个点删掉的方案数，转移做个树形背包就好了。对于无根的树，我们保留的连通块可能跟dfs的根不相连，不妨让它在连通块的根处被计数，再算一个删掉上面的方案数即可，还要考虑删空的情况。
时间复杂度是$\mathcal{O}(n^2)$，因为少打一个$\ast f1[x]$和少打一个$/f1[x]$调了1h。。。

#include <bits/stdc++.h>
#define MOD 1000000009
 
using namespace std;
 
typedef long long ll;
 
ll pow_mod(ll x,int k) {
  ll ans=1;
  while (k) {
  	if (k&1) ans=ans*x%MOD;
  	x=x*x%MOD;
  	k>>=1;
  }
  return ans;
}
 
ll C[105][105];
 
void pre(int n) {
  for(int i=0;i<=n;i++) C[i][0]=1;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=i;j++) C[i][j]=(C[i-1][j-1]+C[i-1][j])%MOD;
}
 
vector <int> e[105];
int d[105],fa[105];
bool vis[105];
 
queue <int> q;
 
void bfs(int n) {
  for(int i=1;i<=n;i++)
    if (d[i]<=1) q.push(i);
  while (!q.empty()) {
  	int x=q.front();q.pop();
  	vis[x]=1;
  	if (!d[x]) continue;
  	for(int i=0;i<e[x].size();i++)
  	  if (!vis[e[x][i]]) {
  	  	  int u=e[x][i];
  	  	  fa[x]=u;
  	  	  d[u]--;
  	  	  if (d[u]==1) q.push(u); 
		}
  }
}
 
int size[105];
 
int f1[105],f2[105],g[105][105];
int sum1[105],sum2[105];
 
int cur[105];
 
void dfs1(int x) {
  f1[x]=g[x][0]=1;
  for(int i=0;i<e[x].size();i++)
    if (e[x][i]!=fa[x]&&d[e[x][i]]<=1) {
    	int u=e[x][i];
    	dfs1(u);
    	f1[x]=(ll)f1[x]*f1[u]%MOD*C[size[x]+size[u]][size[u]]%MOD;
    	memset(cur,0,sizeof(cur));
    	for(int j=0;j<=size[x];j++)
    	for(int k=0;k<=size[u];k++)
    	  cur[j+k]=(cur[j+k]+(ll)g[x][j]*g[u][k]%MOD*C[j+k][j])%MOD;
    	memcpy(g[x],cur,sizeof(g[x]));
    	size[x]+=size[u];
	}
  size[x]++;
  g[x][size[x]]=(g[x][size[x]]+f1[x])%MOD;
}
 
int sz;
 
void dfs2(int x) {
  ll s=f2[x];
  int usz=sz-size[x];
  for(int i=0;i<size[x];i++) sum2[usz+i]=(sum2[usz+i]+(ll)f2[x]*g[x][i]%MOD*C[usz+i][i])%MOD;
  for(int i=0;i<e[x].size();i++)
    if (e[x][i]!=fa[x]) {
    	int u=e[x][i];
    	s=s*f1[u]%MOD*C[usz+size[u]][size[u]]%MOD;
    	usz+=size[u];
	}
  sum2[sz]=(sum2[sz]+s)%MOD;
  for(int i=0;i<e[x].size();i++)
    if (e[x][i]!=fa[x]) {
    	int u=e[x][i];
    	f2[u]=s*pow_mod(C[usz][size[u]],MOD-2)%MOD*pow_mod(f1[u],MOD-2)%MOD;
    	dfs2(u);
	}
}
 
int main() {
  int n,m;
  scanf("%d%d",&n,&m);
  pre(n);
  for(int i=1;i<=m;i++) {
  	int x,y;
  	scanf("%d%d",&x,&y);
  	e[x].push_back(y);
  	e[y].push_back(x);
  	d[x]++;d[y]++;
  }
  bfs(n);
  sum1[0]=1;
  int s=0;
  for(int i=1;i<=n;i++)
    if (vis[i]&&!vis[fa[i]]) {
    	memset(sum2,0,sizeof(sum2));
    	if (d[i]==1) {
    		dfs1(i);
    		sz=size[i];
    		for(int j=0;j<=sz;j++) sum2[j]=g[i][j];
		}
		else {
			dfs1(i);
			sz=size[i];
			f2[i]=1;
			dfs2(i);
		}
		memset(cur,0,sizeof(cur));
		for(int j=0;j<=s;j++)
		for(int k=0;k<=sz;k++)
		  cur[j+k]=(cur[j+k]+(ll)sum1[j]*sum2[k]%MOD*C[j+k][j])%MOD;
		memcpy(sum1,cur,sizeof(sum1));
		s+=sz;
	}
  for(int i=0;i<=n;i++) printf("%d\n",sum1[i]);
  return 0;
}