AT4352-[ARC101C] Ribbons on Tree【dp,容斥】

正题

题目链接:
https://www.luogu.com.cn/problem/AT4352
https://atcoder.jp/contests/arc101/tasks/arc101_c

题目大意

$n$个点之间两两配对，要求配对点之间的路径覆盖整棵树，求方案数

解题思路

考虑容斥，我们钦定有$l$条边没有路径覆盖，就有$dp$状态$f_{i,j,l}$表示$i$的子树中，目前该子树的联通块大小为$j$，已经切断了$l$条边，我们发现该状态已经是$O(n^3)$的，显然无法通过。

考虑优化，每个$dp$状态的容斥系数是$(-1{)}^l$，所以我们可以不用记录$l$这一维度，之间用状态表示乘上了容斥系数的值。

那么我们就有$dp$方程，定义$calc(x)$表示$x$个点两两匹配的值
$$(-1)\ast f_{x,i}\ast f_{y,j}\ast calc(j)\to f_{x,i}$$
$$f_{x,i}\ast f_{y,j}\to f_{x,i+j}$$
该方程的复杂度就是每个子树的乘积，可以理解为每个点堆之间进行一次贡献，时间复杂度$O(n^2)$

$calc(x)={\prod }_{i=1}^{x-1}[i\&1==1]i$

$code$

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=5100,XJQ=1e9+7;
struct node{
	ll to,next;
}a[N*2];
ll n,tot,ls[N],siz[N],f[N][N],ans,g[N],fac[N];
void addl(ll x,ll y){
	a[++tot].to=y;
	a[tot].next=ls[x];
	ls[x]=tot;
}
void dp(ll x,ll fa){
	siz[x]=f[x][1]=1;
	for(ll i=ls[x];i;i=a[i].next){
		ll y=a[i].to;
		if(y==fa)continue;
		dp(y,x);
		for(ll i=1;i<=siz[x];i++)
			g[i]=f[x][i],f[x][i]=0;
		for(ll j=1;j<=siz[x];j++){
			for(ll k=1;k<=siz[y];k++){
				(f[x][j+k]+=f[y][k]*g[j]%XJQ)%=XJQ;
				if((k&1)==0)
					(f[x][j]+=XJQ-f[y][k]*g[j]%XJQ*fac[k-1]%XJQ)%=XJQ;
			}
		}
		siz[x]+=siz[y];
	}
	return;
}
int main()
{
	scanf("%lld",&n);
	for(ll i=1;i<n;i++){
		ll x,y;
		scanf("%lld%lld",&x,&y);
		addl(x,y);addl(y,x);
	}
	fac[0]=fac[1]=1;
	for(ll i=2;i<=n;i++)
		fac[i]=fac[i-2]*i%XJQ;
	dp(1,1);
	for(ll i=2;i<=n;i++)
		ans=(ans+f[1][i]*fac[i-1]%XJQ)%XJQ;
	printf("%lld",ans);
}