BZOJ 1016 最小生成树计数（矩阵树定理）

我们把边从小到大排序，然后依次插入一种权值的边，然后把每一个联通块合并。
然后当一次插入的边不止一条时做矩阵树定理就行了。算出有多少种生成树就行了。
剩下的交给乘法原理。
实现一不小心就会让程序变得很丑

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
#define int long long
const int mod=31011;
const int N=110;
int fa[N],a[N][N][N],n,m,b[1010],cnt[N],id[N],w[N],ans[N],mmp[N];
struct edge{
	int u,v,w;
}e[1010];
bool cmp(edge a,edge b){
	return a.w<b.w;
}
int find(int x){
	if(fa[x]==x)return x;
	else return fa[x]=find(fa[x]);
}
int gauss(int x,int n){
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			a[x][i][j]=(a[x][i][j]+mod)%mod;
	int f=1,ans=1;
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			int A=a[x][i][i],B=a[x][j][i];
			while(B){
				int t=A/B;A%=B;swap(A,B);
				for(int k=i;k<=n;k++)a[x][i][k]=(a[x][i][k]-t*a[x][j][k]%mod+mod)%mod;
				for(int k=i;k<=n;k++)swap(a[x][i][k],a[x][j][k]);
				f=-f;
			}
		}
		ans=ans*a[x][i][i]%mod;
	}
	memset(a[x],0,sizeof(a[x]));
	return (ans*f+mod)%mod;;
}
int read(){
	int sum=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){sum=sum*10+ch-'0';ch=getchar();}
	return sum*f;
}
void init(){
	for(int j=1;j<=n;j++)cnt[j]=0;
	for(int j=1;j<=n;j++)fa[j]=j;
}
signed main(){
	n=read(),m=read();
	for(int i=1;i<=m;i++)e[i].u=read(),e[i].v=read(),e[i].w=read(),b[i]=e[i].w;
	sort(b+1,b+1+m);
	int num=unique(b+1,b+1+m)-b-1;
	for(int i=1;i<=m;i++)e[i].w=lower_bound(b+1,b+1+num,e[i].w)-b;
	sort(e+1,e+1+m,cmp);
	int now=1;
	for(int i=1;i<=n;i++)ans[i]=1;
	for(int i=1;i<=num;i++){
		init();
		int line=now,tmp=0;
		while(line<=m&&e[line].w==i){
			if(e[line].u==e[line].v){line++;continue;}
			int x=find(e[line].u),y=find(e[line].v);
			if(x!=y)fa[x]=y;
			line++;
		}
		for(int j=1;j<=n;j++)fa[j]=find(j);
		for(int j=1;j<=n;j++)id[j]=++cnt[fa[j]];
		for(int j=now;j<=line-1;j++){
			if(e[j].u==e[j].v)continue;
			a[fa[e[j].u]][id[e[j].u]][id[e[j].u]]++;
			a[fa[e[j].v]][id[e[j].v]][id[e[j].v]]++;
			a[fa[e[j].u]][id[e[j].u]][id[e[j].v]]--;
			a[fa[e[j].u]][id[e[j].v]][id[e[j].u]]--;
		}
		for(int j=1;j<=n;j++)w[j]=1;
		for(int j=1;j<=n;j++)w[fa[j]]=w[fa[j]]*ans[j]%mod;
		for(int j=1;j<=n;j++)
			if(cnt[j]){
				mmp[j]=++tmp;
				if(cnt[j]==1)ans[tmp]=w[j];
				else ans[tmp]=w[j]*gauss(j,cnt[j]-1)%mod;
			}
		for(int j=line;j<=m;j++)e[j].u=mmp[fa[e[j].u]],e[j].v=mmp[fa[e[j].v]];
		now=line;n=tmp;
	}
	if(n>1)printf("0");
	else printf("%lld",ans[1]);
	return 0;
}