洛谷P1991 无线通讯网——MST

题目：https://www.luogu.org/problemnew/show/P1991

本题难在对题意的理解，也就是难在卫星电话的放置上。以样例说明，最小生成树MST如下图所示：

可以放置2部卫星电话，所以还要放置2部接发器。

将边BC去掉，得到两个连通块AB与CD，不妨卫星电话放在点A，点C，则连通块AB内接发器功率要求为100，连通块CD内接发器功率为212.13，求出最大功率为212.13。

从而转化为求最小生成树中的第2条边的权值。而2=(4-1)-(2-1)=4-2。

分析：求出最小生成树，这棵树共p-1条边，去掉树中最大的s-1条边后，将得到s个连通子图，每个连通子图放置一部卫星电话。所以具体做法是用Kruskal算法输出第p-s条边的权值。

AC代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
struct edge{
	int from,to,next;
	double value;
};
edge a[500005];
int n,m,x[501],y[501];
int num,head[501],father[501],vis[500005];
int edge_num,ans;
double dis(int i,int j){
	return sqrt( (x[i]-x[j])*(x[i]-x[j]) + (y[i]-y[j])*(y[i]-y[j]) );
}
void join(int from,int to,double value){
	a[++num].from=from;
	a[num].to=to;
	a[num].value=value;
	a[num].next=head[from];
	head[from]=num;
	/*上面两行代码不用写*/ 
}
int find(int x){
	if(father[x]==x) return x;
	return father[x]=find(father[x]);
}
void unionn(int x,int y){
	int r1=find(x),r2=find(y);
	father[r1]=r2;
}
void qsort(int l,int r){
	edge mid=a[(l+r)/2];
	int i=l,j=r;
	while(i<=j){
		while( a[i].value<mid.value )i++;
		while( a[j].value>mid.value )j--;
		if(i<=j){
			edge tmp=a[i];
			a[i]=a[j];
			a[j]=tmp;
			i++;j--;
		}
	}
	if(l<j)qsort(l,j);
	if(i<r)qsort(i,r);	
}
void Kruskal(){
	for(int i=1;i<=num;i++){
		if(!vis[i]){
			vis[i]=1;
			if( find(a[i].from)!=find(a[i].to) ){
				unionn(a[i].from,a[i].to);
				edge_num++;
				if(edge_num==m-n){
					printf("%.2f\n",a[i].value);
					break;
				}
			}
		}
	}
}
int main(){
	memset(head,-1,sizeof(head));
	cin>>n>>m;
	for(int i=1;i<=m;i++)cin>>x[i]>>y[i];
	for(int i=1;i<=m;i++)
		for(int j=1;j<=m;j++)
			if(i!=j){
				double value=dis(i,j);
				join(i,j,value);
				join(j,i,value);
			}
	qsort(1,num);
	for(int i=1;i<=m;i++)father[i]=i;
	Kruskal();
	return 0;
}