[hdu-4081] Qin Shi Huang's National Road System题解

题目传送门
题意解析:题目的要求就是给你一张完全图，然后让你找出一棵生成树，并且在生成树上找出一条边，使这条边上的两个点的点权之和除以这棵生成树去掉这条边之后的边权之和最大。

---

My opinion:我一开始看懂题意后是一脸懵逼的。什么鬼啊，这题目。于是我开始简化题目的意思，先想到了如果这棵生成树给你的后该怎么办？很显然，如果我们已经知道生成树之后，只要枚举每条边，然后找最大的答案就好了。然而我们并不知道这棵树。从答案入手，答案是max(A/B),那么显然我们需要让A尽可能的大，而B尽可能的小。如果我们固定了其中一个值，剩下的就可以枚举了。显然Ｂ的值是最好固定的，只要求一遍最小生成树，Ｂ的值算是固定了，但是还有Ａ，如果Ａ所在的这条边不在生成树上怎么办？当然是强行让它连在树上好了。显然我们可以开始枚举每一条边，这时候就会分成两种情况，这条边是不是在树上，如果在树上，那一切好说，但是如果不在树上呢？那么我们通过随便一画图，因为这是一棵树，所以如果再次添加一条边的话，就会形成一个环，我们只要把环上的最大的那条边去掉，留下这条边，那么B就会成为存在这条边的情况下的最小值，符合题意。
总结:结题步骤如下：
1、先求出最小生成树
2、枚举每条边
3、判断这条边是否属于这棵树，分两种情况处理
4、找出最大值并输出

---

我那总感觉很丑的代码:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#define rep(i,a,n) for (int i=a;i<=n;i++)
#define per(i,a,n) for (int i=a;i>=n;i--)
#define Clear(a,x) memset(a,x,sizeof(a))
#define sqr(x) (x)*(x)
#define ll long long
#define db double
#define INF 2000000000
#define eps 1e-8
using namespace std;
ll read(){
    ll x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9') f=ch=='-'?-1:f,ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int maxn=1005;
int n;
db A,B;
double map[maxn][maxn],dis[maxn],a[maxn][maxn];
int x[maxn],y[maxn];
int point[maxn];
bool flag[maxn][maxn],vis[maxn];
int fa[maxn];
double calc(int x1,int y1,int x2,int y2){
    return sqrt((db)sqr(x1-x2)+(db)sqr(y1-y2));
}
void prim(){
    A=0,B=0;
    Clear(flag,0);
    Clear(vis,0);
    Clear(a,0);
    vis[1]=1;
    rep(i,1,n){
        dis[i]=map[1][i];
        fa[i]=1;
    }
    rep(i,1,n-1){
        int u=-1;
        rep(j,1,n)
            if (!vis[j]&&(u==-1||dis[j]<dis[u])) u=j;
        flag[u][fa[u]]=flag[fa[u]][u]=1;
        B+=map[fa[u]][u];
        vis[u]=1;
        rep(j,1,n){
            if (vis[j]&&j!=u) a[u][j]=a[j][u]=max(a[j][fa[u]],dis[u]);
            if (!vis[j]&&dis[j]>map[u][j]){
                dis[j]=map[u][j];
                fa[j]=u;
            }
        }
    }
}
int main(){
    int T=read();
    while (T--){
        n=read();
        rep(i,1,n){
            x[i]=read(),y[i]=read();
            point[i]=read();
        }
        Clear(map,0);
        rep(i,1,n)
            rep(j,1,n)
                if (i!=j) map[i][j]=map[j][i]=calc(x[i],y[i],x[j],y[j]);
        prim();
        db ans=-1;
        rep(i,1,n)
            rep(j,1,n)
                if (i!=j)
                    if (flag[i][j]) ans=max(ans,(point[j]+point[i])/(B-map[i][j]));
                        else ans=max(ans,(point[i]+point[j])/(B-a[i][j]));
        printf("%.2f\n",ans);
    }
    return 0;
} 

附上AC记录