HDU 4081（转自YouthDance的CSDN博客）

Qin Shi Huang's National Road System

题目链接：Click Here~

题目重述：

    题目说秦始皇要修路，把N个城市用N-1条边连通。且他希望花费最小，但是这时候有一个多管闲事的道士出来说他有魔法可以帮助秦始皇变成一条路，但是只能变出一条。但是，两个人对修路的法案存在歧义，道士希望修路可以给更多的百姓带来福利，而秦始皇希望修路要尽量使花费小。最后，秦始皇拿出了一个公式A/B，A表示两个城市的人数，B表示出了用魔法变出来的路外，最短的总距离。现在要你求出A/B的最大值。

思路分析：

   从题目可以看出是一个此小生成树问题，但是也有人有DP AC了。具体思路就看网上别人的解释吧。

[html] view plaincopyprint?

1. #include <iostream>  
2. #include <algorithm>  
3. #include <cstdio>  
4. #include <cstring>  
5. #include <cmath>  
6. using namespace std;  
7.   
8. const int N = 1e3 + 5;  
9. int n,m;  
10. double path[N][N],G[N][N];  
11. bool used[N][N],hash[N];  
12. struct Point  
13. {  
14.   double x,y,cost;  
15. }p[N];  
16.   
17. double Dist(Point a,Point b)  
18. {  
19.     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));  
20. }  
21.   
22. double Prim()  
23. {  
24.     int pre[N];  
25.     double sum = 0,mincost[N];  
26.     memset(hash,0,sizeof(hash));  
27.     memset(used,0,sizeof(used));  
28.     memset(path,0,sizeof(path));  
29.     hash[1] = 1;  
30.     for(int i = 1;i <= n;++i){  
31.        mincost[i] = G[1][i];  
32.        pre[i] = 1;  
33.     }  
34.     for(int i = 1;i < n;++i)  
35.     {  
36.         int u = -1;  
37.         for(int j = 1;j <= n;++j)if(!hash[j]){  
38.            if(u == -1||mincost[j] < mincost[u])  
39.              u = j;  
40.         }  
41.         hash[u] = 1;  
42.         used[u][pre[u]] = used[pre[u]][u] = 1;                          //判断两点之间是否在生成树中  
43.         sum += G[pre[u]][u];  
44.         for(int j = 1;j <= n;++j)  
45.         {  
46.             if(hash[j]&&u != j){  
47.                path[j][u] = path[u][j] = max(mincost[u],path[j][pre[u]]);  //生成树中的最大权值  
48. //               path[j][u] = path[u][j] = mincost[u];  
49.             }  
50.             if(!hash[j]){  
51.                if(mincost[j] > G[u][j]){  
52.                   mincost[j] = G[u][j];  
53.                   pre[j] = u;  
54.                }  
55.             }  
56.         }  
57.     }  
58.     return sum;  
59. }  
60. int main()  
61. {  
62.     int T;   
63.     scanf("%d",&T);  
64.     while(T--)  
65.     {  
66.         scanf("%d",&n);  
67.         for(int i = 1;i <= n;++i){  
68.            scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].cost);  
69.         }  
70.         memset(G,0,sizeof(G));  
71.         for(int i = 1;i <= n;++i)  
72.           for(int j = i+1;j <= n;++j){  
73.              G[i][j] = G[j][i] = Dist(p[i],p[j]);  
74.           }  
75.        double mincost = Prim();  
76.        double ratio = -1;  
77.        for(int i = 1;i <= n;++i)  
78.          for(int j = 1;j <= n;++j)  
79.            printf("path[%d][%d] = %lf\n",i,j,path[i][j]);  
80.        for(int i = 1;i <= n;++i)  
81.          for(int j = i+1;j <= n;++j){  
82.              if(used[i][j]){  
83.                  ratio = max(ratio,(p[i].cost+p[j].cost)/(mincost-G[i][j]));  
84.              }else  
85.              {    //删除生成树中的最大的一条边，生成一个类此小生成树(因为，不是严格的此小生成树)  
86.                  ratio = max(ratio,(p[i].cost+p[j].cost)/(mincost-path[i][j]));     
87.              }  
88.          }  
89.        printf("%.2lf\n",ratio);  
90.     }  
91.     return 0;  
92. }