B1821 [JSOI2010]Group 部落划分 Group 二分答案&&并查集

这个题正解是最小生成树,但是。。。最大值最小?一看就是二分答案啊!不用多想,直接二分答案加暴力验证就行了。

题干:

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

代码:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(int i = a;i <= n;i++)
#define lv(i,a,n) for(int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF = 1 << 30;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = 0;
    while(c = getchar(), c < '0' || c > '9')
        if(c == '-') op = 1;
    x = c - '0';
    while(c = getchar(), c >= '0' && c <= '9')
        x = x * 10 + c - '0';
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < 0) putchar('-'), x = -x;
    if(x >= 10) write(x / 10);
    putchar('0' + x % 10);
}
const db eps = 1e-4;
db x[1010],y[1010];
int fa[1010],n,m;
int get_fa(int x)
{
    if(fa[x] != x)
    return fa[x] = get_fa(fa[x]);
    else
    return x;
}
void un(int x,int y)
{
    x = get_fa(x);
    y = get_fa(y);
    if(x != y)
    fa[x] = y;
}
bool check(db mid)
{
    duke(i,1,n)
    {
        fa[i] = i;
    }
    duke(i,1,n)
    {
        duke(j,1,n)
        {
            if((x[i] - x[j]) * (x[i] - x[j]) + (y[i] - y[j]) * (y[i] - y[j]) < mid)
            {
                un(i,j);
            }
        }
    }
    int tot = 0;
    duke(i,1,n)
    {
        if(fa[i] == i)
        tot++;
    }
    if(tot < m)
    return false;
    else
    return true;
}
int main()
{
    db mx = 0,my = 0;
    read(n);read(m);
    duke(i,1,n)
    {
        read(x[i]);
        read(y[i]);
        mx = max(x[i],mx);
        my = max(y[i],my);
    }
    db l = 0,r = mx * mx + my * my;
    while(r - l > eps)
    {
        db mid = (l + r) / 2;
        if(check(mid))
            l = mid;
        else
            r = mid;
    }
    printf("%.2lf",sqrt(l));
}
原文地址:https://www.cnblogs.com/DukeLv/p/9671138.html