hdu 4081 Qin Shi Huang's National Road System（次小生成树）

下面让我们来分析一下这个题 
秦始皇统一中国之后要在全国修公路连接各个城市，抠门皇帝只想修成最小生成树（距离最小，不考虑人力），一个道士说自己可以不花人力物力修一条路，经过两方妥协，选择max（两个城市人口/（生成树长度-这条路的长度））的路让他变，求这个比值最大值。 先说说次小生成树的形成过程，在prim最小生成树的算法基础上，记录每两个点之间路径上最长的边长，如果要求次小生成树的数值的话，遍历每一条不在最小生成树中的边，用来替换一条这两个点在最小生成树上那条路径的最长的边，记录当前生成树的总长度，遍历一圈找最小的。 这个题没问次小生成树的值，但是用到了路径最长边这个数组，我叫他Max[][]，题中问的是（两个城市人口/（生成树长度-这条路的长度）），我们首先考虑分母最小，因为对于某个边来说分子是定值啊，分母分成两种情况：1.该边就在最小生成树里，则直接减掉。2.该边在最小生成树外，那么减掉的是构成的环里面除了这条边以外的最大的边。

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>

using namespace std;
#define MAXN 1005
#define INF 0x3f3f3f3f

int n;
double mp[MAXN][MAXN];
double MAX[MAXN][MAXN];
struct Node{
	int x,y;
	int pop;
}poi[1005];


double disten(Node a,Node b){
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int pre[MAXN],used[MAXN][MAXN];
double prim(int s){
	
	double lowdis=0.0;
	double  dis[MAXN];
	int vis[MAXN];
	
	memset(MAX,0,sizeof(MAX));
	memset(used,0,sizeof(used));
	memset(vis,0,sizeof(vis));
	
	for(int i=0;i<n;i++){
		dis[i]=mp[s][i];
		pre[i]=s;
	}
	pre[s]=-1;
	vis[s]=1;
	for(int i=1;i<n;i++){
		double minc=INF*1.0;
		int p=-1;
		for(int j=0;j<n;j++){
			if(!vis[j]&&minc>dis[j]){
				minc=dis[j];
				p=j;
			}
		}
		
		if(minc==INF) return -1;
		lowdis+=minc;
		vis[p]=1;
		used[p][pre[p]]=used[pre[p]][p]=1;
		for(int j=0;j<n;j++){
			
			if(vis[j]&&j!=p){
				MAX[j][p]=MAX[p][j]=max(MAX[j][pre[p]],dis[p]);
			}
			
			if(!vis[j]&&dis[j]>mp[p][j]){
				dis[j]=mp[p][j];
				pre[j]=p;	
			}	
		}
	}
	return lowdis;
}

int main(){
	int t;
	scanf("%d",&t);
	while(t--){
	
	
		scanf("%d",&n);
		
		for(int i=0;i<n;i++)
		scanf("%d%d%d",&poi[i].x,&poi[i].y,&poi[i].pop);
		
		for(int i=0;i<n;i++)
		{
			mp[i][i]=0;
			for(int j=i+1;j<n;j++){
				mp[i][j]=mp[j][i]=disten(poi[i],poi[j]);
			}
		}
		double dis=prim(0);
		double sum=-1;
		
		for(int i=0;i<n;i++){
			for(int j=i+1;j<n;j++){
				if(used[i][j])
				sum=max(sum,1.0*(poi[i].pop+poi[j].pop)/(dis-mp[i][j]));
				else sum=max(sum,1.0*(poi[i].pop+poi[j].pop)/(dis-MAX[i][j]));
			}
		}
		printf("%.2lf\n",sum);
	}  
	return 0;
}