[P4208][JSOI2008]最小生成树计数（矩阵树定理）

按照kruskal尽量选权值小的边的方法，相同权值的边，不管怎么选，对连通块的改变情况一定是一样的。所以从小到大，对于每种权值的边做一次矩阵树定理，乘积就是答案。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int mod=31011,N=105,M=1005;
struct edge{
	int x,y,w;
}d1[M];
int n,m,f[N],a[N][N],d[N][N],c[N][N],num1,ans,cnt,l1[N],l2[N],flag;
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;
}
int getfa(int x){
	if(f[x]!=x)f[x]=getfa(f[x]);
	return f[x];
}
bool cmp(edge i,edge j){
	return i.w<j.w;
}
inline int gauss(int n1){
	int ans1=1;
	for(int i=1;i<=n1;++i){
		for(int j=i+1;j<=n1;++j)
			while(c[j][i]){
				int x=c[i][i]/c[j][i];
				for(int k=i;k<=n1;++k)c[i][k]=(c[i][k]-x*c[j][k]%mod)%mod;
				swap(c[i],c[j]),ans1=-ans1;
			}
		ans1=(ans1*c[i][i]%mod+mod)%mod;
	}
	return ans1;
}
int main(){
	n=read(),m=read();
	for(int i=1;i<=n;++i)f[i]=i;
	for(int i=1;i<=m;++i)d1[i].x=read(),d1[i].y=read(),d1[i].w=read();
	sort(d1+1,d1+m+1,cmp);
	d1[0].w=d1[m+1].w=-1,ans=1,cnt=0,num1=0,flag=0,memset(l1,0,sizeof l1);
	for(int l=0,i=1;i<=m+1;++i){
		if(d1[i].w!=d1[i-1].w){
			for(int j=l;j<i;++j){
				int xx=getfa(d1[j].x),yy=getfa(d1[j].y);
				if(xx!=yy)f[xx]=yy,++cnt;
			}
			for(int j=1;j<=num1;++j)for(int k=1;k<=num1;++k)c[j][k]=d[j][k]-a[j][k];
			for(int xx,j=1;j<=num1;++j){xx=getfa(l2[j]);if(xx==l2[j])++c[l1[xx]][l1[xx]];}
			ans=ans*gauss(num1)%mod;
			if(cnt==n-1){flag=1;break;}
			l=i,num1=0,memset(a,0,sizeof a),memset(d,0,sizeof d),memset(l1,0,sizeof l1);
		}
		int xx=getfa(d1[i].x),yy=getfa(d1[i].y);
		if(xx==yy)continue;
		if(!l1[xx])l1[xx]=++num1,l2[num1]=xx;
		if(!l1[yy])l1[yy]=++num1,l2[num1]=yy;
		++d[l1[xx]][l1[xx]],++d[l1[yy]][l1[yy]],++a[l1[xx]][l1[yy]],++a[l1[yy]][l1[xx]];
	}
	flag?printf("%d",ans):printf("0");
	return 0;
}