图论复习【Kruskal】

K r u s k a l Kruskal Kruskal

【算法思想】
一种贪心算法,将边按权值排序,每次从剩下的边集里选择权值最小且两个端点不在统一集合的边加入生成树中,反复操作,直到加入了 n − 1 n-1 n1条边。
【算法步骤】

  1. 将原图中边按权值从小到大快排。
  2. 按照权值从小到大选边。若当前选取的边加入后使生成树 T T T形成环,则舍弃当前边;否则标记当前边并计数。
  3. 重复 2. 2. 2.操作,直到生成树 T T T中包含 n − 1 n-1 n1条边。否则当遍历完所有的边之后仍然选取不到 n − 1 n-1 n1条边表示 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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值