[P4455][CQOI2018]社交网络（矩阵树定理）

关于矩阵树定理用于有向图生成树计数：

①内向树生成树计数。 
A为邻接矩阵， D为出度矩阵。 
C=D−A。 
以root为根的内向生成树个数为C的余子式M[root,root]的行列式。 
②外向树生成树计数。 
A为邻接矩阵， D为入度矩阵。 
C=D−A。 
以root为根的外向生成树个数为C的余子式M[root,root]的行列式。

然后这就是一道板子题了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=260,mod=10007;
int n,m,a[N][N],c[N][N],d[N][N],ans;
inline int read(){
	int x=0,f=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return f?-x:x;
}
inline void gauss(){
	for(int i=1;i<n;++i){
		for(int j=i+1;j<n;++j)
			while(c[j][i]){
				int x=c[i][i]/c[j][i];
				for(int k=i;k<n;++k)c[i][k]=(c[i][k]-x*c[j][k]%mod)%mod;
				swap(c[i],c[j]);ans=-ans;
			}
		ans=(ans*c[i][i]%mod+mod)%mod;
	}
}
int main(){
	n=read(),m=read();
	memset(a,0,sizeof a),memset(d,0,sizeof d);
	for(int u,v,i=1;i<=m;++i)v=read(),u=read(),++d[v][v],++a[u][v];
	for(int i=2;i<=n;++i)for(int j=2;j<=n;++j)c[i-1][j-1]=d[i][j]-a[i][j];
	ans=1;gauss();
	printf("%d",ans);
	return 0;
}