hdu 4081 次小生成树

http://acm.hdu.edu.cn/showproblem.php?pid=4081

我就不说题意了，为了使A/B最大，就应该是B越小，故可以先求出n个点的最小生成树。因此，可以枚举每一条边，假设最小生成树的值是B, 而枚举的那条边长度是edge[i][j],  如果这一条边已经是属于最小生成树上的，那么最终式子的值是A/(B-edge[i][j])。如果这一条不属于最小生成树上的， 那么添加上这条边，就会有n条边，那么就会使得有了一个环，为了使得它还是一个生成树，就要删掉环上的一棵树。 为了让生成树尽量少，那么就要删掉除了加入的那条边以外，权值最大的那条路径。 假设删除的那个边的权值是Max[i][j], 那么就是A/(B-Max[i][j]).

这题的关键也在于怎样求出次小生成树;

先用prim求出最小生成树T.在prim的同时，用一个矩阵max[u][v] 记录 在T中连结任意两点u,v的唯一的路中权值最大的那条边的权值.这是很容易做到的，因为prim是每次增加一个结点s, 而设已经标号了的结点集合为W, 则W中所有的结点到s的路中的最大权值的边就是当前加入的这条边.

View Code

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cmath>
 4 const int N=1010;
 5 const double inf=1e14;
 6 using namespace std;
 7 
 8 struct Point{
 9     int x,y,z;
10 }point[N];
11 
12 int n;
13 double edge[N][N];
14 int nearvex[N];//保存前驱
15 double  lowcost[N]; 
16 double sum;
17 
18 int used[N][N];
19 int visited[N];
20 double  Max[N][N];//用来保存最小生成树中两点之间的权值最大的边
21 
22 
23 void prim(int v0){
24     sum=0;
25     memset(used,0,sizeof(used));
26     memset(visited,0,sizeof(visited));
27     memset(Max,0,sizeof(Max));
28     for(int i=1;i<=n;i++){
29         lowcost[i]=edge[v0][i];
30         nearvex[i]=v0;
31     }
32     visited[v0]=1;
33     for(int i=1;i<n;i++){
34         double min=inf;
35         int v=-1;
36         for(int j=1;j<=n;j++){
37             if(!visited[j]&&lowcost[j]<min){
38                 v=j,min=lowcost[j];
39             }
40         }
41         if(v!=-1){
42             sum+=lowcost[v];
43             used[v][nearvex[v]]=used[nearvex[v]][v]=1;//标记这条边已经是最小使用过//
44             visited[v]=1;
45             for(int k=1;k<=n;k++){
46                 if(visited[k]&&k!=v){
47                     //对于那些已经加入最小生成树的边，只要每次更新所有点到新加入的点之间的边权值最大值即可
48                     Max[v][k]=Max[k][v]=(Max[k][nearvex[v]]>lowcost[v]?Max[k][nearvex[v]]:lowcost[v]);
49                 }
50                 if(!visited[k]&&edge[v][k]<lowcost[k]){
51                     lowcost[k]=edge[v][k];
52                     nearvex[k]=v;
53                 }
54             }
55         }
56     }
57 }
58 
59 
60 int main(){
61     int t;
62     scanf("%d",&t);
63     while(t--){
64         scanf("%d",&n);
65         for(int i=1;i<=n;i++){
66             scanf("%d%d%d",&point[i].x,&point[i].y,&point[i].z);
67         }
68         for(int i=1;i<=n;i++){
69             edge[i][i]=0;
70             for(int j=i+1;j<=n;j++){
71                 double dis=sqrt(pow((point[i].x-point[j].x)*1.0,2)+pow((point[i].y-point[j].y)*1.0,2));
72                 edge[i][j]=edge[j][i]=dis;
73             }
74         }
75         prim(1);
76         double r=-1;
77         for(int i=1;i<=n;i++){
78             for(int j=1;j<=n;j++)if(i!=j){
79                 if(used[i][j]){
80                     r=(r>(point[i].z+point[j].z)*1.0/(sum-edge[i][j])?r:(point[i].z+point[j].z)*1.0/(sum-edge[i][j]));
81                 }else if(!used[i][j]){
82                     r=(r>(point[i].z+point[j].z)*1.0/(sum-Max[i][j])?r:(point[i].z+point[j].z)*1.0/(sum-Max[i][j]));
83                 }
84             }
85         }
86         printf("%.2lf\n",r);
87     }
88     return 0;
89 }