矩阵树定理应用:
给定一个边带权的图,求其所有生成树的边权积之和。
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;
}