洛谷P4047 [JSOI2010]部落划分题解

矫情吗；* 2021-09-29 23:40 333阅读 0赞

## 洛谷P4047 \[JSOI2010\]部落划分题解 ##

## 题目描述 ##

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

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

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

例如，下面的左图表示了一个好的划分，而右图则不是。请你编程帮助聪聪解决这个难题。

![30573.png][]

## 输入输出格式 ##

输入格式：

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

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

输出格式：

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

## 输入输出样例 ##

输入样例\#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
    
    解析：

输入所有的点的坐标,然后就枚举记录每一个点到另外其他点的距离,注意要用double类型.  
然后就排一个序,我这里是用结构体重载运算符,先从最小的开始合并并查集,并把这个边设为一个极大值,注意要把其他的点也要合并,就是说for循环一遍把所有在一个部落里的点与点之间的距离也设为一个极大值.  
最后就在所有的边中取一个最小值,该值即为所值.

![ContractedBlock.gif][] ![ExpandedBlockStart.gif][]

1 #include <cstdio>
     2 #include <iostream>
     3 #include <cstdlib>
     4 #include <queue>
     5 #include <cstring>
     6 #include <stack>
     7 #include <cmath>
     8 #include <algorithm>
     9 #include <iomanip>
    10 #include <ctime>
    11 #define INF 0x7fffffff
    12 #define Max 1000025
    13 inline int read()
    14 {
    15     int s=0,f=1;
    16     char ch=getchar();
    17     while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar();}
    18     while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
    19     return s*f;
    20 }
    21 int n,k,pa[Max],x[Max],y[Max],ID[Max];
    22 int find(int x)
    23 {
    24     if(x!=pa[x])
    25         pa[x]=find(pa[x]);
    26     return pa[x];
    27 }
    28 void add(int x,int y)
    29 {
    30     x=find(x);y=find(y);
    31     pa[y]=x;
    32 }
    33 #define D double
    34 D d(int x_1,int y_1,int x_2,int y_2)
    35 {
    36     return (D)sqrt((x_1-x_2)*(x_1-x_2) + (y_1-y_2)*(y_1-y_2));
    37 }
    38 struct point
    39 {
    40     int x,y;
    41     D dis;
    42     friend bool operator<(point a,point b) {
    43         return a.dis<b.dis;
    44     }
    45 }t[Max];
    46 int main()
    47 {
    48 //    freopen("Group.in","r",stdin);
    49 //    freopen("Group.out","w",stdout);
    50     n=read();k=read();
    51     for(int i = 1 ; i <= n ; ++ i)
    52         pa[i]=i;
    53     for(int i = 1 ; i <= n ; ++ i)
    54         x[i]=read(),y[i]=read(),ID[i]=i;
    55     int num=0;
    56     for(int i = 1 ; i <= n ; ++ i)
    57         for(int j = 1 ; j <= n ; ++ j)
    58             if(i!=j) {
    59                 t[++num].dis=d(x[i],y[i],x[j],y[j]);
    60                 t[num].x=ID[i];t[num].y=ID[j];
    61             }
    62     int m=n;
    63     D Min=INF;
    64     std::sort(t+1,t+1+num);
    65     for(int i = 1 ; i <= num ; ++ i) {
    66         int x=find(t[i].x),y=find(t[i].y);
    67         if(x!=y)
    68             add(x,y),m--,t[i].dis=INF*1.00;
    69         if(m==k)
    70             break;
    71     }
    72     for(int i = 1 ; i <= num ; ++ i)
    73         if(find(t[i].x)==find(t[i].y))
    74             t[i].dis=INF*1.0;
    75     for(int i = 1 ; i <= num ; ++ i)
    76         Min=std::min(Min,t[i].dis);
    77     printf("%.2lf",Min);
    78     return 0;
    79 }

AC代码

转载于:https://www.cnblogs.com/handsomegodzilla/p/11181848.html

[30573.png]: /images/20210725/77f1142fad8945e98ec9d985ed30bcb4.png
[ContractedBlock.gif]: https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif
[ExpandedBlockStart.gif]: /images/20210725/f2d94f036ac64423bbe4daadadc7bec3.png