当前位置: 代码迷 >> 综合 >> bzoj 1016--最小生成树计数
  详细解决方案

bzoj 1016--最小生成树计数

热度:31   发布时间:2023-12-01 22:03:20.0

对着模板一顿乱敲(很无奈)

#pragma comment(linker,"/STACK:1024000000,1024000000")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;typedef long long ll;
const int maxn=105;
const int maxm=1005;struct node{int set[maxn];void init(int n){for(int i=0;i<=n;i++) set[i]=i;}int Find(int x){return x==set[x]?set[x]:set[x]=Find(set[x]);}int Union(int x,int y){int xx=Find(x);int yy=Find(y);if(xx==yy) return -1;set[xx]=yy;return 1;}
}a,b,c;struct Edge{int u,v,len;
}edge[maxm];
int edgenum;void Add(int u,int v,int len){Edge E={u,v,len};edge[++edgenum]=E;
}bool vis[maxn];
vector<int>g[maxn];
ll p[maxn][maxn],deg[maxn][maxn];bool cmp(Edge e1,Edge e2){return e1.len<e2.len;
}ll EDT(ll a[][maxn],int n,ll Mod){ll temp=1;ll t;for(int i=0;i<n;i++){for(int j=0;j<n;j++) a[i][j]%=Mod;}for(int i=1;i<n;i++){for(int j=i+1;j<n;j++){while(a[j][i]){t=a[i][i]/a[j][i];for(int k=i;k<n;k++){a[i][k]-=a[j][k]*t;a[i][k]%=Mod;}for(int k=i;k<n;k++) swap(a[i][k],a[j][k]);temp=-temp;}}temp=temp*a[i][i]%Mod;}return (temp+Mod)%Mod;
}ll cal_MST_count(int n,ll Mod){sort(edge+1,edge+edgenum+1,cmp);int pre=edge[1].len;ll ans=1;a.init(n);b.init(n);memset(vis,false,sizeof(vis));memset(deg,0,sizeof(deg));for(int i=0;i<=n;i++) g[i].clear();for(int t=1;t<=edgenum+1;t++){if(edge[t].len!=pre || t==edgenum+1){for(int i=1;i<=n;i++){if(vis[i]){int k=b.Find(i);g[k].push_back(i);vis[i]=0;}}for(int i=1;i<=n;i++){if(g[i].size()){memset(p,0,sizeof(p));for(int j=0;j<g[i].size();j++){for(int k=j+1;k<g[i].size();k++){int x=g[i][j];int y=g[i][k];p[j][k]=p[k][j]=-deg[x][y];p[j][j]+=deg[x][y];p[k][k]+=deg[x][y];}}ans=ans*EDT(p,g[i].size(),Mod)%Mod;for(int j=0;j<g[i].size();j++) a.set[g[i][j]]=i;}}memset(deg,0,sizeof(deg));for(int i=1;i<=n;i++){b.set[i]=a.Find(i);g[i].clear();}if(t==edgenum+1) break;pre=edge[t].len;}int x=a.Find(edge[t].u);int y=a.Find(edge[t].v);if(x==y) continue;vis[x]=vis[y]=1;b.Union(x,y);deg[x][y]++;deg[y][x]++; }if(!edgenum) return 0;for(int i=2;i<=n;i++){if(b.Find(i)!=b.Find(1)) return 0;}return ans;
}void init(){edgenum=0;
}int n,m,u,v,len;
ll Mod=31011;int main(){cin>>n>>m;init();while(m--){scanf("%d%d%d",&u,&v,&len);Add(u,v,len);}printf("%lld\n",cal_MST_count(n,Mod));
}