[P3317][SDOI2014]重建（矩阵树定理）

矩阵树定理应用：

给定一个边带权的图，求其所有生成树的边权积之和。 
A为边权矩阵， D为发出边权和矩阵。 
即对于一条边(u,v)，边权为w， A[u,v]=A[v,u]=w。 
如果有重边， A[u,v]为所有边(u,v)的权值和。 
D[u,u]等于由u发出的所有边的权值之和。 
对于任意 u≠v，D[u,v]=0。 
C=D−A。 
所有生成树的边权积之和等于C的任意余子式M[i,i]的行列式。 

矩阵树定理可以求出所有生成树的边权积之和，但这道题，每个生成树对答案的贡献是树边的概率与(1-非树边的概率)的乘积。考虑所有生成树树边概率积之和可求，要是把这个值乘上原图中(1-每条边概率)的乘积，就可以转化成类似所求的形式了。但这样得到的是，所有生成树树边的概率*(1-概率)与非树边的(1-概率)的乘积之和。所以我们在使用矩阵树定理时，把边权设为该边的概率/(1-概率)，就可以得到所求了。

#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;
}