【ZJOI2017】仙人掌

题目大意：给你一张联通图，问你加边使得他是个仙人掌的方案数

解析：

首先如果一开始就不是个仙人掌，那方案数显然是0

否则，我们发现由于每条边都只能被至多覆盖一次（仙人掌的性质），所以多出来的边必定可以把这棵树分割成一些独立的树，那么我们把问题转化为对于一颗树来说，求出他加边使得他是个仙人掌的方案数

那么接下来我们思考一下dp方程

设dp[i][0]表示当前这个点i的子树内部不存在向上连的边的方案数

dp[i][1]表示当前这个点i的子树内部存在向上连的边的方案数

那么子树之间是可以两两配对的所以我们还需要预处理出对于i棵子树来说他们任意配对的方案数，不妨设他为g[i]

那么我们考虑一下

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll mod=998244353;
const int N=5e5;
struct node{
	int u,to;
}edge[N<<2];
int head[N+10],f[N+10],dep[N+10];
bool b[N+10],c[N+10];
ll fac1[N+10],dp[N+10][2];
int T,n,m,x,k,y,sum;
ll ans;
void add(int x,int y)
{
    edge[k].u=y;
    edge[k].to=head[x];
    head[x]=k++; 
}
void init()
{
	for (int i=1;i<=n;i++) head[i]=-1;
	for (int i=1;i<=n;i++) f[i]=0;
	for (int i=1;i<=n;i++) b[i]=0;
}
void la()
{
	/*
	fac[0]=1; fac[1]=0;//full
	for (int i=2;i<=N;i++)
		 if (i&1) fac[i]=0;
		else fac[i]=(fac[i-2]*(ll)(i-1))%mod;
	*/
	fac1[0]=1;	fac1[1]=1;//not full
    for (int i=2;i<=N;i++) 
    {
      fac1[i]=(fac1[i-1]+(ll)(i-1)*fac1[i-2]%mod)%mod;
    }
    // cout << fac1[N] << endl;
}  
/******************************************************************/
void dfs(int now,int fa,int d)
{
	b[now]=1; dep[now]=d;
	for (int i=head[now];i!=-1;i=edge[i].to)
	{
		int u=edge[i].u;
		if (u==fa) continue;
        if (b[u]&&dep[u]<dep[now]) {//
        	f[now]++; f[u]--;
        } else if (!b[u]) dfs(u,now,d+1),f[now]+=f[u];//
	}
}
/******************************************************************/
void solve(int now,int fa)
{
	b[now]=1;
	ll ans=1;
	int sz=0;
	for (int i=head[now];i!=-1;i=edge[i].to)
	{
		int u=edge[i].u;
		if (b[u]||f[u]) continue;
        solve(u,now);
        sz++;
        ans=(ans*((dp[u][0]+dp[u][1])%mod)%mod);
	}
    if (sz) dp[now][1]=(ans*(fac1[sz-1]*(ll)sz%mod)%mod); else dp[now][1]=0;
    dp[now][0]=(ans*fac1[sz])%mod;
}
void dfs1(int now,int fa)
{
	if (!b[now]) {solve(now,0); ans=(ans*dp[now][0])%mod;}
	sum++;
	c[now]=1;
	for (int i=head[now];i!=-1;i=edge[i].to)
	{
		int u=edge[i].u;
		if (u==fa||c[u]) continue;
		dfs1(u,now);
	}
}
int main()
{
	scanf("%d",&T);
	la();
	while (T--){
		scanf("%d%d",&n,&m);
		k=0;
		init();
		for (int i=1;i<=m;i++) {
			scanf("%d%d",&x,&y);
			add(x,y); add(y,x);
		}
		dfs(1,0,1); 
		bool p=1;
		for (int i=1;i<=n;i++) if (f[i]>1) {
			printf("%d\n",0); p=0; break;
		}
		if (!p) continue;
		ans=1;
		for (int i=1;i<=n;i++) b[i]=0,c[i]=0;
		dfs1(1,0);
	    printf("%lld\n",ans);
	}
}