/*首先求出最小生成树大小sum。并记录最小生成树的边。
然后遍历图中所有边。
如果边在最小生成树中.这用边两点的人口和/sum-边长.
如果边不在最小生成树中,这条边一上来就形成了一个环。
例如原来是a->b->c 加入一条边 ac 那么 为了保证sum-边长最小 就得sum-max(ab,bc)
则就是 ac两点的人口和/sum-max(ab,bc)
*/
#include <stdio.h>
#include <string.h>
#include<math.h>
#define MAX 1005
#define MAXCOST 0x7fffffff
#define max(a,b) a>b?a:b
int n;
double graph[MAX][MAX],sum,dis,ans,min;
double city[MAX][3];
double dp[MAX][MAX];
bool adjvex[MAX];//记录顶点是否被访问
double lowcost[MAX];//只管记录的是已是最小生成树集合到各个点的最小值
int pre[MAX],p;//最小生成树中的点的父节点
void Prim()
{
int i,j,k;
for(i=1;i<=n;i++)
{
dp[i][i]=0;
lowcost[i]=graph[1][i]; //从第一个的顶点开始
pre[i]=1;
}
memset(adjvex,false,sizeof(adjvex));
adjvex[1]=true;
for(i=2;i<=n;i++)//这里注意因为第一个点在循环外算了。所以还需进行n-1次
{
min=MAXCOST;
for(j=2;j<=n;j++)//选出当前最小生成树集合到非集合点的最小的边
if(adjvex[j]==false && lowcost[j]<min)
{
min=lowcost[j];
k=j; //记下最小的点
}
p=pre[k];
dp[k][p]=dp[p][k]=graph[p][k];
for(j=1;j<=n;j++)//更新k j上最大的边。
{
if(adjvex[j])
{
dp[k][j]=dp[j][k]=max(dp[j][p],dp[p][k]);//k是刚刚加入生成树的点.他与其他点的联系都父节点建立.
}
}
adjvex[k]=true; //为true表示该权值已经是最小,为flase是还不确定,应继续更新
for(j=2;j<=n;j++)
{
if(adjvex[j]==false && lowcost[j]>graph[k][j]) //更新lowcost
{
pre[j]=k;
lowcost[j]=graph[k][j];
}
}
}
for(i=2;i<=n;i++)//第一个点早选了。从2开始
{
sum+=lowcost[i]; //求最小生成树路径长度
}
}
void getAns()//遍历每一对点
{
ans=0;
for(int i=1;i<n;i++)
{
for(int j=i+1;j<=n;j++)
{
double t1=(city[i][2]+city[j][2])/(sum-dp[i][j]);
ans=max(ans,t1);
}
}
}
int main()
{
freopen("test.txt","r",stdin);
int x,y,z,T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
sum=0;
for(int i=1;i<=n;i++)
{
scanf("%lf%lf%lf",&city[i][0],&city[i][1],&city[i][2]);
}
for(int i=1;i<n;i++)
{
for(int j=i+1;j<=n;j++)
{
dis=sqrt(((city[i][0]-city[j][0])*(city[i][0]-city[j][0])+(city[i][1]-city[j][1])*(city[i][1]-city[j][1])));
graph[i][j]=dis;
graph[j][i]=dis;
dp[i][j]=0;
dp[j][i]=0;
}
}
Prim();
getAns();
printf("%.2lf\n",ans);
}
}
hdu 4081 Qin Shi Huang's National Road System(次小生成树)
最新推荐文章于 2021-08-23 20:47:31 发布