秦始皇的国家道路系统

题目描述
有n个城市，秦始皇要修用n-1条路把它们连起来，要求从任一点出发，都可以到达其它的任意点。秦始皇希望这所有n-1条路长度之和最短。然后徐福突然有冒出来，说是他有魔法，可以不用人力、财力就变出其中任意一条路出来。

秦始皇希望徐福能把要修的n-1条路中最长的那条变出来，但是徐福希望能把要求的人力数量最多的那条变出来。对于每条路所需要的人力，是指这条路连接的两个城市的人数之和。

最终，秦始皇给出了一个公式，A/B，A是指要徐福用魔法变出的那条路所需人力， B是指除了徐福变出来的那条之外的所有n-2条路径长度之和，选使得A/B值最大的那条。

希望A/B的值最大，那么最大是多少？

输入
第一行输入T,表示测试组数

每组测试首先输入n，表示城市

接下来n行，每行输入x,y,p，分别表示第i个城市的坐标和人数

输出
输出最大的A/B，保留两位小数

样例输入
2
4
1 1 20
1 2 30
200 2 80
200 1 100
3
1 1 20
1 2 30
2 2 40
样例输出
65.00
70.00
提示
T<=10

2 < n <=1000, 0 <=x, y <=1000, 0 <=p <=100000

任意两个城市不在同一坐标

题解
先用prim求最小生成树，然后枚举边，形成环中找一条最大边删除求最大值。

代码

#include<bits/stdc++.h>
#define inf 2147483647
#define N 400005
#define M 1000005
#define ll long long
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int tot,n,Head[1005],ret[2005],Next[2005],p[1005],x[1005],y[1005],rt;
double ans,sum,len[2005],mn[1005],a[1005][1005],dis[1005],c[1005][1005];
bool vis[1005],flag[1005][1005];
struct node{double dis;int now,pre;};
struct cmp
{
    bool operator()(node a,node b)
    {
        return a.dis>b.dis;
    }
};
inline void ins(int u,int v,double l)
{
    ret[++tot]=v;len[tot]=l;
    Next[tot]=Head[u];Head[u]=tot;
}
void prim()
{
    priority_queue<node,vector<node>,cmp>q;
    for (int i=1;i<=n;i++) mn[i]=inf;
    mn[1]=0;q.push((node){0,1,0});
    while (!q.empty())
    {
        node now=q.top();q.pop();
        if (vis[now.now]) continue;vis[now.now]=1;sum+=mn[now.now];
        if (now.pre)ins(now.pre,now.now,now.dis),ins(now.now,now.pre,now.dis);
        for (int i=1;i<=n;i++)
        {
            if (now.now==i) continue;
            if (a[now.now][i]<mn[i])
            {
                mn[i]=a[now.now][i];
                q.push((node){mn[i],i,now.now});
            }
        }
    }
}
void dfs(int u,int f)
{
    for (int i=Head[u];i;i=Next[i])
    {
        flag[u][ret[i]]=1;
        if (ret[i]==f) continue;
        ans=max(ans,(double)(p[u]+p[ret[i]])/(sum-a[u][ret[i]]));
        c[rt][ret[i]]=max(c[rt][u],len[i]);
        dfs(ret[i],u);
    }
}
int main()
{
    int Case=read();
    while (Case--)
    {
        n=read();tot=0;sum=ans=0;
        for (int i=1;i<=n;i++)
        {
            vis[i]=Head[i]=0;
            for (int j=1;j<=n;j++)
                flag[i][j]=0;
        }

        for (int i=1;i<=n;i++)
        {
            x[i]=read();y[i]=read();p[i]=read();
        }
        for (int i=1;i<n;i++)
            for (int j=i+1;j<=n;j++)
                a[i][j]=a[j][i]=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])); 
        prim();
        for (int i=1;i<=n;i++)
        {
            rt=i;dfs(i,0);
        }
        for (int i=1;i<n;i++)
            for (int j=i+1;j<=n;j++)
            {
                if (flag[i][j]) continue;
                ans=max(ans,(double)(p[i]+p[j])/(sum-c[i][j]));
            }
        printf("%.2lf\n",ans);
    }
    return 0;
}