【51nod1806】wangyurzee的树（prufer序列）（容斥原理）

传送门

---

题解：

关于度数的限制我们可以考虑放到prufer序列上。

一个点的度数就是它在prufer序列中出现次数+1。

所以问题就变成了一个长度为$n-2$的序列，每个位置的数可以是$[1,n]$中间的任意整数，限制某些数出现次数不能为某个值，求方案数。

很显然的容斥，考虑用总方案数减去所有不合法的方案，就是求满足$0,1,2\cdots m$个限制的所有方案数。

假设我们现在满足$|S|$中的限制，设$x=|S|，sum=\sum _{i=1}^x{S}_i$，考虑从$n-2$个里面选出$sum$个来$\left(\genfrac{}{}{0ex}{}{n-2}{sum}\right)$，然后这$sum$个位置要对$S$做可重集合的排列$\frac{sum!}{\prod S_i!}$，剩下$n-2-sum$可以随便填剩下$n-x$个数$(n-x{)}^{n-2-sum}$。

所以方案数为：$$\frac{(n-2)!}{(n-2-sum)!\prod S_i!}(n-x{)}^{n-2-sum}$$

然后容斥dfs就行了，注意可能有对于同一个点的限制，要在dfs的时候注意处理一下。

复杂度：$O(2^m)$

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;
using pii=std::pair<int,int>;
#define fi first
#define se second

cs int mod=1e9+7;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){ll r=(ll)a*b;;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline int Inc(int &a,int b){return (a+=b)>=mod&&(a-=mod),a;}
inline int Dec(int &a,int b){return (a-=b)<0&&(a+=mod),a;}

cs int P=1e6+6;
int fac[P],ifac[P];

int n,m,cnt;
int ans=0;

inline void init(){
	fac[0]=ifac[0]=1;
	for(int re i=1;i<=n;++i)fac[i]=mul(fac[i-1],i);
	ifac[n]=power(fac[n],mod-2);
	for(int re i=n-1;i;--i)ifac[i]=mul(ifac[i+1],i+1);
}

std::vector<int> d[20];
pii limit[20];

void dfs(int cur,int r1,int r2,int down,int coef){
	if(r1<0)return ;
	if(cur>cnt){
		Inc(ans,mul(coef,mul(mul(fac[n-2],down),power(r2,r1,ifac[r1]))));
		return ;
	}
	dfs(cur+1,r1,r2,down,coef);
	for(int re v:d[cur])dfs(cur+1,r1-v,r2-1,mul(down,ifac[v]),mod-coef);
}

signed main(){
//	freopen("tree.in","r",stdin);
	scanf("%d%d",&n,&m);init();if(n==1)cout<<1,exit(0);
	for(int re i=1;i<=m;++i)scanf("%d%d",&limit[i].fi,&limit[i].se);
	std::sort(limit+1,limit+m+1);
	for(int re i=1;i<=m;++i){
		if(limit[i].fi!=limit[i-1].fi)++cnt;
		d[cnt].push_back(limit[i].se-1);
	}
	dfs(1,n-2,n,1,1);
	cout<<ans<<"\n";
	return 0;
}