BZOJ 1821 Group 部落划分 Group

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec   Memory Limit: 64 MB
Submit: 2405   Solved: 1135
[ Submit][ Status][ Discuss]

Description

聪聪研究发现，荒岛野人总是过着群居的生活，但是，并不是整个荒岛上的所有野人都属于同一个部落，野人们总是拉帮结派形成属于自己的部落，不同的部落之间则经常发生争斗。只是，这一切都成为谜团了——聪聪根本就不知道部落究竟是如何分布的。 不过好消息是，聪聪得到了一份荒岛的地图。地图上标注了N个野人居住的地点（可以看作是平面上的坐标）。我们知道，同一个部落的野人总是生活在附近。我们把两个部落的距离，定义为部落中距离最近的那两个居住点的距离。聪聪还获得了一个有意义的信息——这些野人总共被分为了K个部落！这真是个好消息。聪聪希望从这些信息里挖掘出所有部落的详细信息。他正在尝试这样一种算法： 对于任意一种部落划分的方法，都能够求出两个部落之间的距离，聪聪希望求出一种部落划分的方法，使靠得最近的两个部落尽可能远离。 例如，下面的左图表示了一个好的划分，而右图则不是。请你编程帮助聪聪解决这个难题。

Input

第一行包含两个整数N和K(1< = N < = 1000,1< K < = N)，分别代表了野人居住点的数量和部落的数量。

接下来N行，每行包含两个正整数x,y，描述了一个居住点的坐标(0 < =x, y < =10000)

Output

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

Sample Input

4 2
0 0
0 1
1 1
1 0

Sample Output

1.00

题意一眼看上去很诡，但是仔细想想其实可以二分答案，考虑精度，保险起见把精度放在了1e-5，二分的值mid的意思是不能有任何一个长度小于mid的边没有连接两个部落，否则这个边一定是最短距离且小于mid则不合题意。

#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<iostream>
#include<cmath>

using namespace std;

const int MAXN=1010;
struct edge
{
	int u,v;
	double w;
};
edge m[MAXN*MAXN/2];
int fa[MAXN],cnt;
bool go[MAXN];
int cmp(const edge a,const edge b)
{
	return a.w<b.w;
}

struct P
{
	int x,y;
};
P a[MAXN];
int find(int x)
{
	if (fa[x]!=x)
	{
		fa[x]=find(fa[x]);
	}
	return fa[x];
}

int main()
{
	int n,k;
	scanf("%d%d",&n,&k);
	for (int i=1;i<=n;i++)
	{
		scanf("%d%d",&a[i].x,&a[i].y);
		for (int j=1;j<i;j++)
		{
			m[++cnt].u=i;
			m[cnt].v=j;
			m[cnt].w=sqrt((double)pow(a[i].x-a[j].x,2)+(double)pow(a[i].y-a[j].y,2));
		}
	}
	double l=0,r=100000000,mid;
	int tot;
	while (r-l>0.00001)
	{
		memset(go,0,sizeof(go));
		tot=0;
		mid=(l+r)/2.0;
		for (int i=1;i<=n;i++)
		{
			fa[i]=i;
		}
		for (int i=1;i<=cnt;i++)
		{
			if (m[i].w>=mid)
			{
				continue;
			}
			int f1=find(m[i].u),f2=find(m[i].v);
			if (f1!=f2)
			{
				fa[f1]=f2;
			}
		}
		for (int i=1;i<=n;i++)
		{
			if (!go[find(i)]) 
			{
				go[find(i)]=1;
				tot++;
			}
		}
		if (tot<k)
		{
			r=mid;
		}
		else
		{
			l=mid;
		}
	}
	
	printf("%.2lf",l);
    return 0;
}

一开始sqrt只套了一个pow导致wa了好几次。。。