洛谷 P1991 无线通信网 题解

P1991 无线通信网

题目描述

国防部计划用无线网络链接若干个边防哨所。2 种不一样的通信技术用来搭建无线网络；

每一个边防哨所都要配备无线电收发器；有一些哨所还能够增配卫星电话。

任意两个配备了一条卫星电话线路的哨所（两边都ᤕ有卫星电话）都可以通话，不管他们相距多远。而只经过无线电收发器通话的哨所之间的距离不能超过 D，这是受收发器的功率限制。收发器的功率越高，通话距离 D 会更远，但同时价格也会更贵。

收发器须要统一购买和安装，因此所有哨所只能选择安装一种型号的收发器。换句话说，每一对哨所之间的通话距离都是同一个 D。你的任务是肯定收发器必须的最小通话距离 D，使得每一对哨所之间至少有一条通话路径（直接的或者间接的）。

输入格式

从 wireless.in 中输入数据第 1 行，2 个整数 S 和 P，S 表示可安装的卫星电话的哨所数，P 表示边防哨所的数量。接下里 P 行，每行两个整数 x，y 描述一个哨所的平面坐标(x, y)，以 km 为单位。

输出格式

输出 wireless.out 中

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

输入输出样例

输入 #1

2 4
0 100
0 300
0 600
150 750

输出 #1

212.13

说明/提示

对于 20% 的数据：P = 2，S = 1

对于另外 20% 的数据：P = 4，S = 2

对于 100% 的数据保证：1 ≤ S ≤ 100，S < P ≤ 500，0 ≤ x，y ≤ 10000。

【思路】

最小生成树 + 克鲁斯卡尔

【题目大意】

有无线电收发器和卫星电话两种工具
卫星电话不须要考虑距离
无线电收发器有传播距离的限制
求这个限制最小是多少

【题目分析】

卫星电话能够当作免费的
p个哨所须要p-1条边链接起来
而s个卫星电话能够免去s-1条边
因此就只剩下了p-s条边须要找
找最小的
因此最小生成树就很显然了

【核心思路】

将两两之间的匹配方式用结构图储存一下
而后sort排序从最短的边开始试
若是这条边链接的两个点没有被接起来
那就连起来就行了
这样知道用完p-s条边
由于从小到大排的序
因此最后用的那一条边的权值就是最大的
输出就行了

【完整代码】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm> 
#include<cmath>

using namespace std;

int read()
{
    int sum = 0,fg = 1;
    char c = getchar();
    while(c < '0' || c > '9')
    {
        if(c == '-')fg = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9')
    {
        sum = sum * 10 + c - '0';
        c = getchar();
    }
    return sum * fg;
}
const int Max = 503;
int x[Max],y[Max];
int father[Max];
struct node
{
    int xx,yy;
    double l;
}a[Max * Max];

bool cmp(const node aa,const node bb)
{
    return aa.l < bb.l;
}

int find(int xz)
{
    if(father[xz] != xz)father[xz] = find(father[xz]);
    return father[xz];
}

void hebing(int xz,int yz)
{
    xz = find(xz);
    yz = find(yz);
    if(yz != xz)
        father[xz] = yz;
}

int main()
{
    int s = read(),p = read();
    for(register int i = 1;i <= p;++ i)
        father[i] = i;
    for(register int i = 1;i <= p;++ i)
        x[i] = read(),y[i] = read();
    int jj = 0;
    for(register int i = 1;i <= p;++ i)
    {
        for(register int j = i + 1;j <= p;++ j)
        {
            a[++ jj].xx = i;
            a[jj].yy = j;
            a[jj].l = double(sqrt((x[i] - x[j]) * (x[i] - x[j]) * 1.0 + (y[i] - y[j]) * (y[i] - y[j]) * 1.0));
        }
    }
    sort(a + 1,a + jj + 1,cmp);
    int tot = p - s;
    int js = 0;
    double M = 0;
    for(register int i = 1;i <= jj;++ i)
    {
        if(find(a[i].xx) != find(a[i].yy))
        {
            js ++;
            M = a[i].l;
            hebing(a[i].xx,a[i].yy);
        }
        if(js == tot)
            break;
    }
    printf("%.2lf\n",M);
    return 0;
}