BZOJ 1821 JSOI2010 部落划分 Group prim

Description

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

Input

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

Output

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

Sample Input

Sample Output

1.00

题解

最小生成树求第K大边

代码

 1 #include <cstdio>
 2 #include <cmath>
 3 #include <cstring>
 4 #include <ctime>
 5 #include <iostream>
 6 #include <algorithm>
 7 #include <set>
 8 #include <vector>
 9 #include <queue>
10 #include <typeinfo>
11 #include <map>
12 #include <stack>
13 typedef long long ll;
14 #define inf 0x7fffffff
15 using namespace std;
16 inline ll read()
17 {
18     ll x=0,f=1;
19     char ch=getchar();
20     while(ch<'0'||ch>'9')
21     {
22         if(ch=='-')f=-1;
23         ch=getchar();
24     }
25     while(ch>='0'&&ch<='9')
26     {
27         x=x*10+ch-'0';
28         ch=getchar();
29     }
30     return x*f;
31 }
32 
33 //**************************************************************************************
34 double mp[1111][1111];
35 double lowcost[1111];
36 int used[1111];
37 int n;
38  int k;
39 int kk;
40 double  ans[1111];
41 void prim()
42 {
43    int x=1;
44    double mn;
45    memset(used,0,sizeof(used));
46    memset(lowcost,0,sizeof(lowcost));
47    for(int i=1;i<=n;i++)
48     lowcost[i]=mp[x][i];
49    used[x]=1;
50    for(int i=2;i<=n;i++)
51    {
52        mn=inf*1.0;
53        for(int j=1;j<=n;j++)
54        {
55            if(!used[j]&&lowcost[j]<mn)
56            {
57                x=j;
58                mn=lowcost[j];
59            }
60        }
61        ans[++kk]=mn;
62        used[x]=1;
63        for(int j=1;j<=n;j++)
64        {
65            if(!used[j]&&mp[x][j]<lowcost[j])
66             lowcost[j]=mp[x][j];
67        }
68    }
69 }
70 int main()
71 {
72  double x[1111];
73  double y[1111];
74     scanf("%d%d",&n,&k);
75     for(int i=1;i<=n;i++)
76     {
77         scanf("%lf%lf",&x[i],&y[i]);
78 
79     }
80     for(int i=1;i<=n;i++)
81     {
82         for(int j=1;j<=n;j++)
83         {
84             if(i==j)mp[i][j]=0;
85             else
86             {
87                 mp[i][j]=sqrt((y[i]-y[j])*(y[i]-y[j])+(x[i]-x[j])*(x[i]-x[j]));
88             }
89         }
90     }
91     kk=0;
92     prim();
93     sort(ans+1,ans+1+kk);
94     printf("%.2f
",ans[kk-k+2]);
95     return 0;
96 }