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

传送门
思路：
把关于度数的限制放到prufer序列上。
发现可以容斥一波。
我们强制一些点一定不满足限制，然后记得特判一个点被加了多个限制的情况。
然后就是简单的组合计数。
我才不会告诉你我因为阶乘逆元预处理错了wa了十几次呢
代码：

#include<bits/stdc++.h>
#define ri register int
#define fi first
#define se second
using namespace std;
const int rlen=1<<18|1;
inline char gc(){
	static char buf[rlen],*ib,*ob;
	(ib==ob)&&(ob=(ib=buf)+fread(buf,1,rlen,stdin));
	return ib==ob?-1:*ib++;
}
inline int read(){
	int ans=0;
	char ch=gc();
	while(!isdigit(ch))ch=gc();
	while(isdigit(ch))ans=((ans<<2)+ans<<1)+(ch^48),ch=gc();
	return ans;
}
typedef long long ll;
const int mod=1e9+7;
inline int add(const int&a,const int&b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(const int&a,const int&b){return a>=b?a-b:a-b+mod;}
inline int mul(const int&a,const int&b){return (ll)a*b%mod;}
inline void Add(int&a,const int&b){a=a+b>=mod?a+b-mod:a+b;}
inline void Dec(int&a,const int&b){a=a>=b?a-b:a-b+mod;}
inline void Mul(int&a,const int&b){a=(ll)a*b%mod;}
inline int ksm(int a,int p){int ret=1;for(;p;p>>=1,a=mul(a,a))if(p&1)Mul(ret,a);return ret;}
const int N=20,M=1e6+5;
int n,m,fac[M],ifac[M];
typedef pair<int,int> pii;
pii upd[N];
inline void init(int n){
	fac[0]=fac[1]=ifac[0]=ifac[1]=1;
	for(ri i=2;i<=n;++i)fac[i]=mul(fac[i-1],i),ifac[i]=mul(ifac[mod-mod/i*i],mod-mod/i);
	for(ri i=2;i<=n;++i)Mul(ifac[i],ifac[i-1]);
}
bool vis[N];
inline int C(int n,int m){return n>=m?mul(fac[n],mul(ifac[m],ifac[n-m])):0;}
inline int solve(int sta){
	ll sum=n-2;
	int cnt=0,ret=1;
	memset(vis,0,sizeof(vis));
	for(ri i=1;i<=m;++i){
		if(!((sta>>(i-1))&1))continue;
		++cnt;
		if(vis[upd[i].fi])return 0;
		vis[upd[i].fi]=1;
		Mul(ret,C(sum,upd[i].se-1));
		sum-=upd[i].se-1;
		if(sum<0)return 0;
	}
	Mul(ret,ksm(n-cnt,sum));
	return cnt&1?mod-ret:ret;
}
int main(){
	n=read(),m=read();
	for(ri i=1;i<=m;++i)upd[i].fi=read(),upd[i].se=read();
	if(n==1)return cout<<(m==1?0:1),0;
	init(n);
	int ans=0;
	for(ri i=0,up=1<<m;i<up;++i)Add(ans,solve(i));
	cout<<ans;
	return 0;
}