【Kruskal】 洛谷1991 无线通讯网

本文主要介绍了如何运用Kruskal算法解决一道关于无线通讯网的问题。题目要求确保所有哨所之间至少存在一条通话路径。通过分析,为使最大通话距离最小,应选择距离最大的边连接卫星电话覆盖不到的哨所。在Kruskal算法实现中,优先连接P-S个哨所,然后用卫星电话覆盖剩下的S-1个哨所,以此达到最优解。
摘要由CSDN通过智能技术生成

题目描述

国防部计划用无线网络连接若干个边防哨所。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。


分析

这道题几乎可以算作是裸的Kruskal

而因为有P个哨所,其中有S个哨所可以装卫星电话,要使剩余通话距离的最大值最小,因而要使装上卫星电话的哨所之间的距离尽可能的大。而我们又可以发现,如果装上电话的几个哨所连成的边是所有边中边权最大的边,就有可能浪费这几个减小距离的机会。因此只要在Kruskal的过程中,连通P-S条边以后,剩余的S-1条边便可以通过电话来连接。附上代码。


#include<algorithm>
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#define union UNION //防止union和库函数重名
#define MAXN 5005
#define MAXM 200000
using namespace std;
struct Edge
{
    int from;
    int to;
    double weight;
}edges[MAXM];
int cnt;//cnt表示边的条数
int f[MAXN];
int r[MAXN],c[MAXN];//r数组、c数组分别为横纵坐标
int S,P;//S个点可装卫星电话 P个哨所
double ans=0;
bool operator <(const Edge &a,const Edge &b)//重载小于号
{
    return a.weight<b.weight;
}
double cal(int x1,int y1,int x2,int y2)//计算两点间的距离
{
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int find(int x)//find函数与union函数均为并查集
{
    if(f[x]==x) return x;
    else return f[x]=find(f[x]);
}
void union(int x,int y)
{
    x=find(x),y=find(y);
    if(x!=y) f[x]=y;
}
void Kruskal()//Kruskal最小生成树
{
    int t=0;//t表示当前已经连的边数目
    for(int i=1;i<=P;i++) f[i]=i;
    sort(edges,edges+cnt);
    for(int i=0;i<cnt;i++)
    {
        int x=edges[i].from;
        int y=edges[i].to;
        if(find(x)!=find(y))
        {
            union(x,y);
            t++;
            ans=edges[i].weight;
            //printf("%.2lf",edges[i].weight);
        }
        if(t==P-S) break;//如果当前连上的边达到了P-S条 就不用再往下找了
    }
}
int main()
{
    scanf("%d%d",&S,&P);
    for(int i=1;i<=P;i++)
        scanf("%d%d",&r[i],&c[i]);
    for(int i=1;i<=P;i++)
        for(int j=1;j<=P&&j!=i;j++)
        {
            double dist=cal(r[i],c[i],r[j],c[j]);
            edges[cnt].from=i;
            edges[cnt].to=j;
            edges[cnt].weight=dist;
            cnt++;//统计边的条数
        }
    Kruskal();
    printf("%.2lf",ans);
    return 0;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值