【图论+二分】【CQBZOJ 1635】【图论专项赛】无线通讯网

问题 G(1635): 【图论专项赛】无线通讯网
时间限制: 1 Sec 内存限制: 64 MB

题目描述
国防部计划用无线网络连接若干个边防哨所。2种不同的通讯技术用来搭建无线网络：每个边防哨所都要配备无线电收发器；有一些哨所还可以增配卫星电话。
任意两个配备了一条卫星电话线路的哨所均可以通话，无论它们相距多远。而只通过无线电收发器通话的哨所之间的距离不能超过D，这是受收发器的功率限制。收发器的功率越高，通话距离D会更远，但同时价格也更贵。
收发器需要统一购买和安装，所以全部哨所只能选择安装一种型号的收发器。换句话说，每一对哨所之间的通话距离都是同一个D。
你的任务是确定收发器必须的最小通话距离D，使得每一对哨所之间至少有一条通话路径(直接的或者间接的)。这段话很重要，附上原文：
Your job is to determine the minimum D required for the transceivers. There must be at least one communication path (direct or indirect) between every pair of outposts.

输入
第1行：2个整数S(1 <= S <= 100)和P(S < P <= 500)，S表示可安装的卫星电话的线路数，P表示边防哨所的数量。
接下来P行，每行描述一个哨所的平面坐标(x,y)，以km为单位，整数，0<=x,y<=10,000

输出
第1行：1个实数D，表示无线电收发器的最小传输距离。精确到小数点后2位。

样例输入
2 4
0 100
0 300
0 600
150 750

样例输出
212.13

翻译一下题意，要求出一个电话线的最短长度，使得每一个哨所都相连，又因为有卫星电话的缘故，所以连完电话线后保证图中至多有S个联通块就行了。

于是有两种做法，
一是用并查集，从小到大加入边，当连通块数量等于卫星电话数时，便找到答案。
二是二分电话线长度，然后并查集判断连通块数量。

（这里是二分）

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;

#define MAXN 500
#define MAXM
#define INF 0x3f3f3f3f
typedef long long int LL;

int getint()
{
    int rn=0;
    char c=getchar();

    while(c<'0'||'9'<c)
        c=getchar();

    while('0'<=c&&c<='9')
    {
        rn=rn*10+c-'0';
        c=getchar();
    }

    return rn;
}

int x[MAXN+10],y[MAXN+10];
double dist[MAXN+10][MAXN+10];
int N,M;

double cal(int i,int j)
{
    return sqrt((double)(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
}

int fa[MAXN+10];
int root(int x)
{
    if(fa[x]==-1)return x;
    else return fa[x]=root(fa[x]);
}

bool join(int a,int b)
{
    int rta=root(a),rtb=root(b);

    if(rta!=rtb)
    {
        fa[rta]=rtb;
        return true;
    }
    else return false;
}

bool check(double D)
{
    memset(fa,-1,sizeof(fa));

    int ans=N;
    for(int i=1;i<=N;++i)
        for(int j=i+1;j<=N;++j)
            if(dist[i][j]<=D)ans-=join(i,j);

    return ans<=M;
}

int main()
{
    scanf("%d%d",&M,&N);
    int i,j;

    for(i=1;i<=N;++i)
        x[i]=getint(),y[i]=getint();

    for(i=1;i<=N;++i)
        for(j=1;j<=N;++j)
            dist[i][j]=cal(i,j);

    double l=0,r=20000,mid;
    while(r-l>=0.0001)
    {
        mid=(l+r)/2;
        if(check(mid))
            r=mid;
        else l=mid+0.0001;
    }

    printf("%0.2lf\n",r);
}