【题解】洛谷P1991

---

题目

原题链接

描述

国防部计划用无线网络连接若干个边防哨所。2 种不同的通讯技术用来搭建无线网络；
每个边防哨所都要配备无线电收发器；有一些哨所还可以增配卫星电话。

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

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

输入格式

• 第一行，2 个整数 S 和 P，S 表示可安装的卫星电话的哨所数，P 表示边防哨所的数量。
• 接下来
  P 行，每行两个整数 x，y 描述一个哨所的平面坐标 (x,y)，以 km 为单位。

输出格式

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

输入/输出例子

输入

2 4
0 100
0 300
0 600
150 750

输出

212.13

提示

数据范围 $1\le S\le 100,S<P\le 500，0\le x,y\le 10000$。

---

解题思路

分析

1. 首先，我们先不考虑装卫星电话，则此题就是一道最小生成树。
2. 当我们考虑安装卫星电话时，考虑最优方案，我们把距离大的两点间安装卫星电话。
   注意：在考虑在哪安装卫星电话时，安装时得保证最小生成树合法。
3. 经过模拟，我们可以得出连接的边数=P-S
4. 熟悉以上之后，我们就可以开始敲代码了

代码实现

首先最小生成树必须学会，不会的戳一戳
基于刚刚的分析，我们按照正常最小生成树的方法去做，但是：当边连接到P-S条时，就可以得出答案了。
同时，这里的ans应该是 $ans=edge[i].w$

for(int i=1;i<=l;i++)
    {
        int fx=find(edge[i].u),fy=find(edge[i].v);
        if(fx==fy)
            continue;
        ans=edge[i].w;
        fa[fx]=fy;
        cnt++;
        if(cnt>=m-n)
        	return ;
    }

code

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e6+1;
inline char gc()
{
  	static char now[1<<20],*S,*T;
  	if(T==S)
  	{
		T=(S=now)+fread(now,1,1<<20,stdin);
    	if(T==S) 
			return EOF;
  	}
  return *S++;
}
template <typename T>
inline void Read(T&x) 
{
  	x=0;
  	char c=gc();
  	while(c<'0'||c>'9') 
  		c=gc();
  	x=c-'0';
  	while((c=gc())>='0'&&c<='9') 
  		x=x*10+c-'0';
}
template <typename T, typename... Args>
inline void Read(T&x,Args&...args)
{
  	Read(x);
  	Read(args...);
}
struct fy
{
	int u,v;
	double w;
}edge[N<<1];
int fa[N],n,m,cnt,x[N],y[N],l;
double ans;
bool cmp(fy a,fy b)
{
    return a.w<b.w;
}
int find(int x)
{
    while(x!=fa[x]) 
		x=fa[x]=fa[fa[x]];
    return x;
}
void kr()
{
    for(int i=1;i<=l;i++)
    {
        int fx=find(edge[i].u),fy=find(edge[i].v);
        if(fx==fy)
            continue;
        ans=edge[i].w;
        fa[fx]=fy;
        cnt++;
        if(cnt>=m-n)
        	return ;
    }
}
signed main()
{
    Read(n,m);
    for(int i=1;i<=m;i++)
        fa[i]=i;
    for(int i=1;i<=m;i++)
    {
        Read(x[i],y[i]);
        for(int j=1;j<i;j++)
        	edge[++l].u=i,edge[l].v=j,edge[l].w=sqrt(pow(x[j]-x[i],2)*1.0+pow(y[i]-y[j],2)*1.0);
    }
    sort(edge+1,edge+l+1,cmp);
    kr();
    printf("%.2f",ans);
    return 0;
}

更多方法

更多方法