#10015. 「一本通 1.2 练习 2」扩散（Floyd&&Kruskal加并查集）

题目描述

一个点每过一个单位时间就会向四个方向扩散一个距离，如图。

两个点a、b连通，记作e(a,b),当且仅当a、b的扩散区域有公共部分。连通块的定义是块内的任意两个点u、v都必定存在路径e(u,a0),e(a0,a1),…,e(ak,v)。给定平面上的n给点，问最早什么时刻它们形成一个连通块。

输入格式

第一行一个数n，以下n行，每行一个点坐标。

【数据规模】

对于20%的数据，满足1≤N≤5; 1≤X[i],Y[i]≤50;

对于100%的数据，满足1≤N≤50; 1≤X[i],Y[i]≤10^9。

输出格式

一个数，表示最早的时刻所有点形成连通块。

样例 #1

样例输入 #1

2
0 0
5 5
1
2
3

样例输出 #1

5

$$--\ast \ast Floyd\ast \ast --$$

大致思路

• 两个按坐标轴扩散的点，很容易就能想到曼哈顿距离
  $$d(i,j)=|X1-X2|+|Y1-Y2|$$.
  反映在代码上就是

inline int dist(node x,node y){
	int a=abs(x.x-y.x),b=abs(x.y-y.y);
	return  (a+b+1)/2;
}

注意此处是两个点两头扩散，因此曼哈顿距离为一半

• 建图，保存这个点到其他点的距离

for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			tmp=dist(a[i],a[j]);
			dis[i][j]=tmp;
		}
	}

3.主要部分Floyd

for(int k=1;k<=n;k++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				dis[i][j]=min(max(dis[i][k],dis[k][j]),dis[i][j]);
			}
		}
	}

有关其算法我不做太多解释，大致为每个点保存最小路径

全部代码

#include<bits/stdc++.h>
using namespace std;
int n,mid,k1,k2,ans=-1e9,tmp;
int dis[9999][9999];
struct node{
	int x,y;
}a[19999];
inline int dist(node x,node y){
	int a=abs(x.x-y.x),b=abs(x.y-y.y);
	return  (a+b+1)/2;
}
int main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i].x>>a[i].y;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			tmp=dist(a[i],a[j]);
			dis[i][j]=tmp;
		}
	}
	for(int k=1;k<=n;k++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				dis[i][j]=min(max(dis[i][k],dis[k][j]),dis[i][j]);
			}
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			ans=max(dis[i][j],ans);
		}
	}
	cout<<ans;
	return 0;
}

$$--\ast \ast Kruskal加并查集做法\ast \ast --$$

首先介绍一下Kruskal

他属于最小生成树的算法，大致过程为：

• 将所有边从小到大排序；
• 将所有点分别放入各自的并查集中
• 选择所有未使用边中边权最小的
• 若该边所连接的两个点已经联通，则舍去，否则合并这两个并查集
• 重复上述两步，直到所有点都已经被标记或所有未标记的点都不存在连边与已标记的点直接相连

Kruskal需要使用到并查集

int find(int x){
	return aa[x]==x?x:aa[x]=find(aa[x]);
	}

且解决这个问题同样需要曼哈顿距离

int dist(node x,node y){
	int a=abs(x.x-y.x),b=abs(x.y-y.y);
	return  (a+b+1)/2;
}

AC node

#include<bits/stdc++.h>
using namespace std;
int n,mid,k1,k2,ans;
struct node{
	int x,y;
}a[19999];
int u[100000],v[100000],w[100000],r[100000],aa[100000];
int cnt=1,tmp=0;
bool cmp(int i,int j){
	return w[i]<w[j];
	}
int find(int x){
	return aa[x]==x?x:aa[x]=find(aa[x]);
	}
int dist(node x,node y){
	int a=abs(x.x-y.x),b=abs(x.y-y.y);
	return  (a+b+1)/2;
}
int main(){
	cin>>n;
	ans=0;
	for(int i=1;i<=n;i++){
		cin>>a[i].x>>a[i].y;
	}
	for(int i=1;i<=n;i++){
		aa[i]=i;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(i==j)continue;
			tmp=dist(a[i],a[j]);
			u[cnt]=i;v[cnt]=j;
			w[cnt]=tmp;r[cnt]=cnt;
			cnt++;
		}
	}
	sort(r+1,r+1+cnt,cmp);
	for(int i=1;i<=cnt;i++){
		int k=r[i];
		if(find(u[k])!=find(v[k])){
			ans=w[k];
			aa[find(u[k])]=find(v[k]);
		}
	}
	cout<<ans;
	return 0;
}