BZOJ 1821 [JSOI2010]Group 部落划分 Group 题解与分析

1821: [JSOI2010]Group 部落划分 Group

Time Limit: 10 Sec   Memory Limit:64 MB
Submit: 825   Solved: 386

Description

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

 

Input

第一行包含两个整数N和K(1<=N<=1000,1<K<=N)，分别代表了野人居住点的数量和部落的数量。 接下来N行，每行包含两个正整数x,y，描述了一个居住点的坐标(0<="x," y<="10000)。" < div>

Output

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

Sample Input

4 2
0 0
0 1
1 1
1 0

Sample Output

1.00

HINT

Source

JSOI2010第二轮Contest1

 

【分析】：

          将各点间连一条边，对这些边从小到大排序，然后将前N-K条边归为一个集合<加边的条件为边的起始点不在同一集合>，保证它们不参与答案贡献，剩下K个直接单独放同一集合，这样的贪心就保证了答案尽可能的大

【代码】：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<vector>
#include<stack>
#include<queue>
#include<cmath>
using namespace std;
#define MAX 1001
#define MAXM 1000001
struct POINT{int x,y;};
struct EDGE{double v;int f,t;};
POINT a[MAX];
EDGE b[MAXM];
int N,K,tot=0,now=0,fa[MAX];
double dist(int x,int y)
{
      double x1=(double)a[x].x*1.0,y1=(double)a[x].y*1.0,x2=(double)a[y].x*1.0,y2=(double)a[y].y*1.0;
      return (sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)));
}
bool cmp(EDGE x,EDGE y){return x.v<y.v;}
int get(int x){return (fa[x]==x ? x : (fa[x]=get(fa[x])));}
int main()
{
      //freopen("input.in","r",stdin);
	  //freopen("output.out","w",stdout); 
	  scanf("%d%d",&N,&K);
	  for(int i=1;i<=N;i++)
	        scanf("%d%d",&a[i].x,&a[i].y);
      for(int i=1;i<=N;i++)
            fa[i]=i;
      //now=N;
      for(int i=1;i<=N;i++)
            for(int j=i+1;j<=N;j++)
            {
                  b[++tot].f=i;
                  b[tot].t=j;
                  b[tot].v=dist(i,j);
            }
      sort(b+1,b+1+tot,cmp);
      for(int i=1;i<=tot;i++)
      {
            if(get(b[i].f)!=get(b[i].t))
            {
                  fa[get(b[i].f)]=get(b[i].t);
                  now++;
            }
            if(now==N-K+1)
            {
                  printf("%.2lf\n",b[i].v);
                  return 0;
            }
      }
	  //system("pause");
      return 0;
}

 

转载注明出处：http://blog.csdn.net/u011400953