K r u s k a l Kruskal Kruskal
【算法思想】
一种贪心算法,将边按权值排序,每次从剩下的边集里选择权值最小且两个端点不在统一集合的边加入生成树中,反复操作,直到加入了
n
−
1
n-1
n−1条边。
【算法步骤】
- 将原图中边按权值从小到大快排。
- 按照权值从小到大选边。若当前选取的边加入后使生成树 T T T形成环,则舍弃当前边;否则标记当前边并计数。
- 重复 2. 2. 2.操作,直到生成树 T T T中包含 n − 1 n-1 n−1条边。否则当遍历完所有的边之后仍然选取不到 n − 1 n-1 n−1条边表示 M S T MST MST不存在。
[SCOI2005]繁忙的都市
注意:
k
r
u
s
k
a
l
kruskal
kruskal合并点是合并当前点指向的祖先点,不是这两个点本身。
#include <bits/stdc++.h>
using namespace std;
const int N=100000+50;
int n,m,fa[N],maxx=-1,tot=0;
struct Edge{
int x,y,w;
}a[N];
int getfa(int x){
return (fa[x]==x)?fa[x]:fa[x]=getfa(fa[x]);
}
inline bool cmp(Edge a,Edge b){
if(a.w!=b.w) return a.w<b.w;
else if(a.x!=b.x) return a.x<b.x;
else return a.y<b.y;
}
inline void kruskal(){
for(int i=1;i<=m;++i){
if(getfa(a[i].x)!=getfa(a[i].y)) {fa[getfa(a[i].x)]=getfa(a[i].y),maxx=max(maxx,a[i].w),++tot;if(tot==n-1) break;}
}
return;
}
inline int read(){
int cnt=0,f=1;char c=getchar();
while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}
while(isdigit(c)){cnt=(cnt<<1)+(cnt<<3)+(c^48);c=getchar();}
return cnt*f;
}
signed main(){
n=read(),m=read();
for(int i=1;i<=m;++i){
a[i].x=read(),a[i].y=read(),a[i].w=read();
}
sort(a+1,a+1+m,cmp);
for(int i=1;i<=n;++i) fa[i]=i;
kruskal();
printf("%d %d",tot,maxx);
return 0;
}
[USACO07DEC]Building Roads S
(这是一句废话但我掉进去了)并查集记得初始化。
关键思想 当最小生成树已经钦定必须选一条边,就把这条边边权赋为
0
0
0,保证这条边一定被选。
根据对权值进行处理来符合对边的要求是常用处理办法
连边的时候不用直接将坐标上的点相连,用编号代替实际坐标上点可以减小编码复杂度。
注意卡精度要不然只有
60
p
t
s
60pts
60pts。
#include <bits/stdc++.h>
#define db double
using namespace std;
const int N=5000100;
int n,m,fa[N],tot=0,cnt=0;
double ans;
struct node{int x,y;}b[N];
struct Edge{int x,y;double w;}a[N];
int getfa(int x){return (fa[x]==x)?fa[x]:fa[x]=getfa(fa[x]);}
inline bool cmp(Edge a,Edge b){
if(a.w!=b.w) return a.w<b.w;
else if(a.x!=b.x) return a.x<b.x;
else return a.y<b.y;
}
inline double calc(int x,int y){return (db)(sqrt((db)(b[x].x-b[y].x)*(b[x].x-b[y].x)+(db)(b[x].y-b[y].y)*(b[x].y-b[y].y)));}
inline void kruskal(){
for(int i=1;i<=tot;++i){
if(getfa(a[i].x)!=getfa(a[i].y)) {
fa[getfa(a[i].x)]=getfa(a[i].y),ans+=a[i].w;
++cnt;
}
if(cnt==n-1) break;
}
return;
}
inline int read(){
int cnt=0,f=1;char c=getchar();
while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}
while(isdigit(c)){cnt=(cnt<<1)+(cnt<<3)+(c^48);c=getchar();}
return cnt*f;
}
signed main(){
// freopen("1.in","r",stdin);
n=read(),m=read();
for(int i=1;i<=n;++i){b[i].x=read(),b[i].y=read();}
for(int i=1;i<=n;++i) fa[i]=i;
for(int i=1;i<=n;++i){
for(int j=i+1;j<=n;++j){
a[++tot].x=i;a[tot].y=j;a[tot].w=calc(i,j);
}
}
for(int i=1;i<=m;++i){a[++tot].x=read();a[tot].y=read();a[tot].w=0.0;}
sort(a+1,a+1+tot,cmp);
kruskal();
printf("%.2lf",ans);
return 0;
}