hdu 4801 最小比率生成树 次小生成树

有n坐城市，输入每坐城市的坐标和人口。

　　现在要在所有城市之间修路，保证每个城市都能相连，并且保证A/B 最大，所有路径的花费和最小

　　A是某条路i两端城市人口的和

　　B表示除路i以外所有路的花费的和（路径i的花费为0）

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string>
#include <cstring>
#include <cmath>
#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <functional>
#include <sstream>
#include <iomanip>
#include <cmath>
#include <cstdlib>
#include <ctime>
typedef long long ll;
//#pragma comment(linker, "/STACK:1024000000,1024000000")  //手动扩栈
#define INF 1e9
#define maxn 1100
#define maxm 100086+10
#define mod 1000000007
#define eps 1e-7
#define PI acos(-1.0)
#define rep(i,n) for(int i=0;i<n;i++)
#define rep1(i,n) for(int i=1;i<=n;i++)
#define scan(n) scanf("%d",&n)
#define scanll(n) scanf("%I64d",&n)
#define scan2(n,m) scanf("%d%d",&n,&m)
#define scans(s) scanf("%s",s);
#define ini(a) memset(a,0,sizeof(a))
#define out(n) printf("%d\n",n)
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
double Dis(double x1,double y1,double x2, double y2) {
    return sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
}
double g[maxn][maxn];
bool used[maxn][maxn];
double maxCost[maxn][maxn]; //maxCost[i][j]表示i与j最小生成树路径
double dis[maxn]; //dis[i]表示未标记节点i加入最小生成树的最小权值
double p[maxn][2];
double val[maxn];
bool v[maxn];
int pa[maxn]; //p[i]记录最小生成树中的节点i的父亲，方便记录路径
double Prim(int n) {
	memset(maxCost,-1,sizeof(maxCost));
	ini(used);
	ini(pa);
	rep1(i,n) dis[i] = INF;
    dis[1] = 0;
	ini(v);
	double ans = 0;
	rep1(i,n)
	{
		int u = -1;
		rep1(j,n)
		{
			if(!v[j])
			{
				if(u == -1 || dis[j] < dis[u]) u = j;
			}
		}
		used[u][pa[u]] = used[pa[u]][u] = 1;
		if(u == -1) break;
		v[u] = 1;
		ans += g[pa[u]][u];
		rep1(j,n)
		{
			if(!v[j])
			{
				if(dis[j] > g[u][j])
				{
					dis[j] = g[u][j];
					pa[j] = u;
				}
			}
			else
			{
				if(j == u) continue;
				maxCost[j][u] = maxCost[u][j] = max(dis[u],maxCost[pa[u]][j]);
			}
		}
	}
	return ans;
}

int main() {
#ifndef ONLINE_JUDGE
	freopen("in.txt","r",stdin);
	//  freopen("out.txt","w",stdout);
#endif
	int T;
	cin>>T;
	int n;
	while(T--)
	{
		cin>>n;
		rep1(i,n)
		{
			scanf("%lf%lf%lf",&p[i][0],&p[i][1],&val[i]);
		}
		rep1(i,n) rep1(j,n)
		{
			g[i][j] = Dis(p[i][0],p[i][1],p[j][0],p[j][1]);
		}
		double B = Prim(n);
		double ans = -1;
		rep1(i,n) //枚举要用魔法的两点
		{
			rep1(j,n)
			{
				if(used[i][j]) ans = max(ans, (val[i]+val[j])/(B-g[i][j]));
				else ans = max(ans, (val[i]+val[j])/(B-maxCost[i][j])); //减去路径中最长的边
			}
		}
		printf("%.2lf\n",ans);
	}
	return 0;
}