UVA 1494 Qin Shi Huang National Road System（类次小生成树）

先找出最小生成树，在求最小生成树的过程中，把uv两点路径上的最大权值边存到d[u][v]上，因为当枚举点uv时，uv就连接了，那么就有n条边了，必须去掉一条边，如果去掉uv路径之外的边，则整棵树就被分开了，所以去掉uv路径上的一条边，同时还要保证去掉之后的树的权值之后最小，所以就要去掉uv路径上的权值最大的边得到B，然后uv两城市的人口相加得到A。就可以求出A/B了，不断枚举uv求出最大的即可。为了方便求d[ ][ ]数组，用prim算法比较好，kruskal很难求。我之前先写的kruskal，因为发现求d数组非常麻烦，才写的prim。用prim就非常方便了，只需要在扩展点之前操作即可，具体详见代码：

//
//  main.cpp
//  Richard
//
//  Created by 邵金杰 on 16/8/17.
//  Copyright © 2016年 邵金杰. All rights reserved.
//



#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
#include<cmath>
using namespace std;
struct Point{
    int x,y,p;
    Point(int x,int y,int p): x(x),y(y),p(p) {}
};
struct edge{
    int u,v;
    double dist;
    edge(int u,int v,double dist): u(u),v(v),dist(dist) {}
    bool operator < (const edge &e) const {
        return dist>e.dist;
    }
};
const int maxn=1000+10;
const int INF=100000000;
int vis[maxn];
double d[maxn];
vector<Point> vp;
vector<vector<edge> > G(maxn);
int n;
double mindist;
double max_val[maxn][maxn];
double dist(Point p,Point q)
{
    double x=p.x-q.x;
    double y=p.y-q.y;
    return hypot(x,y);
}
void prim()
{
    memset(vis,0,sizeof(vis));
    memset(max_val,0,sizeof(max_val));
    for(int i=0;i<n;i++) d[i]=INF;
    priority_queue<edge> pq;
    pq.push(edge(0,0,0));
    edge nd(0,0,0);
    int ncount=0;
    while(ncount<n&&!pq.empty())
    {
        do{
            nd=pq.top();
            pq.pop();
        }while(vis[nd.v]==1&&!pq.empty());
        if(vis[nd.v]==0)
        {
            vis[nd.v]=1;
            ncount++;
            mindist+=nd.dist;
            max_val[nd.u][nd.v]=nd.dist;
            max_val[nd.v][nd.u]=nd.dist;
            for(int i=0;i<n;i++)
            {
                if(vis[i]&&i!=nd.v)
                {
                    max_val[i][nd.v]=max(max_val[i][nd.u],max_val[nd.u][nd.v]);
                    max_val[nd.v][i]=max(max_val[nd.u][i],max_val[nd.v][nd.u]);
                }
            }
            for(int i=0;i<G[nd.v].size();i++)
            {
                int v=G[nd.v][i].v;
                if(vis[v]==0)
                {
                    double distance=G[nd.v][i].dist;
                    if(d[v]>distance)
                    {
                        d[v]=distance;
                        pq.push(edge(nd.v,v,distance));
                    }
                }
            }
        }
    }
}
void work()
{
    double ab=0;
    for(int i=0;i<vp.size();i++)
    {
        for(int j=i+1;j<vp.size();j++)
        {
            ab=max(ab,(vp[i].p+vp[j].p)/(mindist-max_val[i][j]));
        }
    }
    printf("%.2f\n",ab);
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int x,y,p;
        scanf("%d",&n);
        vp.clear();
        mindist=0;
        for(int i=0;i<n;i++)
            G[i].clear();
        for(int i=0;i<n;i++)
        {
            scanf("%d%d%d",&x,&y,&p);
            vp.push_back(Point(x,y,p));
        }
        for(int i=0;i<vp.size();i++)
        {
            for(int j=i+1;j<vp.size();j++)
            {
                G[i].push_back(edge(i,j,dist(vp[i],vp[j])));
                G[j].push_back(edge(j,i,dist(vp[i],vp[j])));
            }
        }
        prim();
        work();
    }
    return 0;
}