【题解】洛谷 P4047

---

题目

原题链接

描述

聪聪研究发现，荒岛野人总是过着群居的生活，但是，并不是整个荒岛上的所有野人都属于同一个部落，野人们总是拉帮结派形成属于自己的部落，不同的部落之间则经常发生争斗。只是，这一切都成为谜团了——聪聪根本就不知道部落究竟是如何分布的。

不过好消息是，聪聪得到了一份荒岛的地图。地图上标注了 n 个野人居住的地点（可以看作是平面上的坐标）。我们知道，同一个部落的野人总是生活在附近。我们把两个部落的距离，定义为部落中距离最近的那两个居住点的距离。聪聪还获得了一个有意义的信息——这些野人总共被分为了 k 个部落！这真是个好消息。聪聪希望从这些信息里挖掘出所有部落的详细信息。他正在尝试这样一种算法：

对于任意一种部落划分的方法，都能够求出两个部落之间的距离，聪聪希望求出一种部落划分的方法，使靠得最近的两个部落尽可能远离。

输入格式

• 输入文件第一行包含两个整数 n 和 k，分别代表了野人居住点的数量和部落的数量。

• 接下来 n 行，每行包含两个整数 x，y，描述了一个居住点的坐标。

输出格式

• 输出一行一个实数，为最优划分时，最近的两个部落的距离，精确到小数点后两位。

输入输出样例

输入1

4 2
0 0
0 1
1 1
1 0

输出 1

1.00

输入 2

9 3
2 2
2 3
3 2
3 3
3 5
3 6
4 6
6 2
6 3

输出 2

2.00

说明/提示

数据规模与约定
对于 100% 的数据，保证 $2\le k\le n\le 10^3，0\le x,y\le 10^4$。

---

解题思路

分析

我们先算出边权，然后从小到大排序，然后把小的边划分到一个部落里，这样任意两个部落间的最小距离就一定大于等于这条边。
如图：

当我们把红色分成一堆时，ans=红色线段<绿色线段。
这里需要自己理解。

代码实现

1. 首先先算出两点间的距离，不会的学一下勾股定理。

   edge[l].w=sqrt(pow(x[j]-x[i],2)+pow(y[i]-y[j],2));
   

2. 然后，我们开始最小生成树。
3. 当连接了n-k+1条边后，我们就可以退出了。
4. 注意：这里是 $ans=edge[i].w$

code

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e6+1;
inline char gc()
{
  	static char now[1<<20],*S,*T;
  	if(T==S)
  	{
		T=(S=now)+fread(now,1,1<<20,stdin);
    	if(T==S) 
			return EOF;
  	}
  return *S++;
}
template <typename T>
inline void Read(T&x) 
{
  	x=0;
  	char c=gc();
  	while(c<'0'||c>'9') 
  		c=gc();
  	x=c-'0';
  	while((c=gc())>='0'&&c<='9') 
  		x=x*10+c-'0';
}
template <typename T, typename... Args>
inline void Read(T&x,Args&...args)
{
  	Read(x);
  	Read(args...);
}
struct fy
{
	int u,v;
	double w;
}edge[N<<1];
int fa[N],n,m,cnt,x[N],y[N],l;
double ans;
bool cmp(fy a,fy b)
{
    return a.w<b.w;
}
int find(int x)
{
    while(x!=fa[x]) 
		x=fa[x]=fa[fa[x]];
    return x;
}
void kr()
{
    for(int i=1;i<=l;i++)
    {
        int fx=find(edge[i].u),fy=find(edge[i].v);
        if(fx==fy)
            continue;
        ans=edge[i].w;
        fa[fx]=fy;
        cnt++;
        if(cnt==n-m+1)
        	return ;
    }
}
signed main()
{
    Read(n,m);
    for(int i=1;i<=n;i++)
        fa[i]=i;
    for(int i=1;i<=n;i++)
    {
        Read(x[i],y[i]);
        for(int j=1;j<i;j++)
        	edge[++l].u=i,edge[l].v=j,edge[l].w=sqrt(pow(x[j]-x[i],2)+pow(y[i]-y[j],2));
    }
    sort(edge+1,edge+l+1,cmp);
    kr();
    printf("%.2f",ans);
    return 0;
}

更多方法

更多方法