hdu4081 Qin Shi Huang's National Road System MST+树形DP

        平面内有n座城市，每座城市有一个人口数，现在可以预先不用任何代价在两个城市之间连一条边，并且把被连的两个城市的人口数之和记为A，在添加了这一条边的基础上，求一颗最小生成树，权值和记为B，求最大的A/B。 n不超过1000，可以先用普利姆求一个最小生成树，然后以每个城市为起点，做一个DFS，遍历的过程中，记录下起点到当前节点的路径上最长的一条边ds，那么遍历到当前城市k时， 起点和k的人口数之和/（最小生成树的权值和 - ds）的最大值就是要的答案了。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <memory.h>
using namespace std;
typedef long long ll;
const double inf=99999999999.9999;
double gs[1010][1010];
double g[1010][1010];
double dis[1010];
bool vis[1010];
int id[1010];
double ans,res;

struct nd
{
    double w;
    int d,next;
}edge[2200];
int gt[2200];
int ne;
int st;
struct node
{
    double x,y;
    int p;
}ct[1010];
inline double len(node a,node b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
void add(int x,int y,double len)
{
    edge[ne].d=y;
    edge[ne].w=len;
    edge[ne].next=gt[x];
    gt[x]=ne;
    ne++;
}
int n,tt;
int m,p,q;
double b[1010];
void dfs(int i,double ds)
{
    int j=gt[i];
    while(j!=-1)
    {
        if (vis[edge[j].d])
        {
            j=edge[j].next;
            continue;
        }
        else
        {
            res=max(res,(ct[st].p+ct[edge[j].d].p)/(ans-max(ds,edge[j].w)));
            vis[edge[j].d]=true;
            dfs(edge[j].d,max(edge[j].w,ds));
        }

    }
}
int main()
{
//    freopen("in.txt","r",stdin);
    scanf("%d",&tt);
    while(tt--)
    {
        scanf("%d",&n);
        for (int i=1; i<=n; i++)
         for (int j=1; j<=n; j++)
         gs[i][j]=inf;
        for (int i=1; i<=n; i++)
        {
            scanf("%lf%lf%d",&ct[i].x,&ct[i].y,&ct[i].p);
        }

        for (int i=1; i<=n; i++)
         for (int j=1; j<=n; j++)
         g[i][j]=g[j][i]=len(ct[i],ct[j]);

        memset(vis,false,sizeof vis);
        vis[1]=true;
        dis[1]=0;
        id[1]=0;
        for (int i=2; i<=n; i++)
        dis[i]=g[1][i],id[i]=1;

        ans=0;
        for (int i=2; i<=n; i++)
        {
            double minn=inf;
            int k;
            for (int j=1; j<=n; j++)
            {
                if (!vis[j] && dis[j]<minn)
                {
                    minn=dis[j];
                    k=j;
                }
            }
            vis[k]=true;
            ans+=minn;
            gs[id[k]][k]=gs[k][id[k]]=g[id[k]][k];

            for (int i=0; i<=n; i++)
            gt[i]=-1;

            for (int j=1; j<=n; j++)
            {
                if (!vis[j] && dis[j]>g[k][j])
                {
                    dis[j]=g[k][j];
                    id[j]=k;
                }
            }
        }
//        for (int i=1; i<=n; i++)
//        {
//         for (int j=1; j<=n; j++)
//         if (gs[i][j]<inf) printf("%.2lf ",gs[i][j]);
//         else printf("x ");
//         puts("");
//        }

        ne=0;
        for (int i=1; i<=n; i++)
         for (int j=1; j<=n; j++)
         {
             if (gs[i][j]<inf-1)
             {
                 add(i,j,gs[i][j]);
             }
         }
        res=0;
        for (int i=1; i<=n; i++)
        {
            st=i;
            memset(vis,false,sizeof vis);
            vis[i]=true;
            dfs(i,0);
        }
        printf("%.2lf\n",res);
    }
    return 0;
}