/* 题意:每个城市有人口,现在有一棵最小生成树,将其中一条边割开,然后选出两个城市的人口A/剩余的最小生成树中的 边和B,使其比值尽可能大。 次小生成树求环上的最大值 */ #include<stdio.h> #include<string.h> #include<stdlib.h> #include<math.h> #include<queue> using namespace std; #define inf 0x3fffffff #define N 1100 struct info{ int x,y,num; }f[N]; struct node{ int u,v,next; double w; }ff[N*N],bian[N*2]; int head[N],yong; double distan(int i,int j) { return sqrt((f[i].x-f[j].x)*(f[i].x-f[j].x)+(f[i].y-f[j].y)*(f[i].y-f[j].y)); } int cmp(const void *a,const void *b) { return (*(struct node *)a).w>(*(struct node *)b).w?1:-1; } void init() { yong=0; memset(head,-1,sizeof(head)); } void addedge(int u,int v,double w) { bian[yong].v=v; bian[yong].next=head[u]; head[u]=yong++; } int pre[N]; int findd(int u) { if(u!=pre[u]) pre[u]=findd(pre[u]); return pre[u]; } double sum; int fp[N]; void min_tree(int n,int len) { int i,k=0; init(); sum=0; for(i=1;i<=n;i++) pre[i]=i; for(i=0;i<len&&k<n-1;i++) { int aa=findd(ff[i].u); int bb=findd(ff[i].v); if(aa!=bb) { pre[aa]=bb; k++; fp[k]=i; addedge(ff[i].u,ff[i].v,ff[i].w); addedge(ff[i].v,ff[i].u,ff[i].w); sum+=ff[i].w; } } return ; } double ma[N][N],maxd[N][N]; void prime(int n) { int i,j,vis[N],pre[N],stac[N],top=0; double dis[N]; memset(maxd,-1,sizeof(maxd)); for(i=1;i<=n;i++) {//初始化 dis[i]=ma[1][i]; pre[i]=1; } stac[top++]=1;//储存最小生成树的点 memset(vis,0,sizeof(vis)); vis[1]=1; for(i=1;i<n;i++) { double minn=inf; int next=1; for(j=1;j<=n;j++) if(!vis[j]&&minn>dis[j]) {//找到最近的 minn=dis[j]; next=j; } vis[next]=1; //printf("%d ",minn); for(j=0;j<top;j++) maxd[next][stac[j]]=maxd[stac[j]][next]=max(minn,maxd[pre[next]][stac[j]]);//当前点与栈里面的点的最大值为minn和他的父节点与其他点的权值最大值 stac[top++]=next;//加入栈 for(j=1;j<=n;j++) if(!vis[j]&&dis[j]>ma[next][j]) {//更新距离 pre[j]=next;//与j连接的最近的点即为他的父节点 dis[j]=ma[next][j]; } } return ; } void inittwo(int n) { int i,j; for(i=1;i<=n;i++) { ma[i][i]=0; for(j=i+1;j<=n;j++) ma[i][j]=ma[j][i]=inf; } return ; } int main() { int t,n,i,j,len; double maxx; scanf("%d",&t); while(t--) { scanf("%d",&n); for(i=1;i<=n;i++) scanf("%d%d%d",&f[i].x,&f[i].y,&f[i].num); len=0; inittwo(n); for(i=1;i<n;i++) for(j=i+1;j<=n;j++) { ff[len].u=i; ff[len].v=j; ff[len++].w=ma[i][j]=ma[j][i]=distan(i,j); } qsort(ff,len,sizeof(ff[0]),cmp); min_tree(n,len); prime(n); maxx=-1; for(i=1;i<n;i++) for(j=i+1;j<=n;j++) { if(1.0*(f[i].num+f[j].num)/(sum-maxd[i][j])>maxx) maxx=1.0*(f[i].num+f[j].num)/(sum-maxd[i][j]); } printf("%.2f ",maxx); } return 0;}
/* 上一种选定城市人口从边出发,这种选定边从点出发 */ #include<stdio.h> #include<string.h> #include<stdlib.h> #include<math.h> #include<queue> using namespace std; #define N 1100 struct info{ int x,y,num; }f[N]; struct node{ int u,v,next; double w; }ff[N*N],bian[N*2]; int head[N],yong; double distan(int i,int j) { return sqrt((f[i].x-f[j].x)*(f[i].x-f[j].x)+(f[i].y-f[j].y)*(f[i].y-f[j].y)); } int cmp(const void *a,const void *b) { return (*(struct node *)a).w>(*(struct node *)b).w?1:-1; } void init() { yong=0; memset(head,-1,sizeof(head)); } void addedge(int u,int v,double w) { bian[yong].v=v; bian[yong].next=head[u]; head[u]=yong++; } int pre[N]; int findd(int u) { if(u!=pre[u]) pre[u]=findd(pre[u]); return pre[u]; } double sum; int fp[N]; void min_tree(int n,int len) { int i,k=0; init(); sum=0; for(i=1;i<=n;i++) pre[i]=i; for(i=0;i<len&&k<n-1;i++) { int aa=findd(ff[i].u); int bb=findd(ff[i].v); if(aa!=bb) { pre[aa]=bb; k++; fp[k]=i; addedge(ff[i].u,ff[i].v,ff[i].w); addedge(ff[i].v,ff[i].u,ff[i].w); sum+=ff[i].w; } } return ; } int vis[N]; int dfs(int u) { int i,k=f[u].num; for(i=head[u];i!=-1;i=bian[i].next) { int v=bian[i].v; if(!vis[v]) { vis[v]=1; int w=dfs(v); if(k<w) k=w; } } return k; } int main() { int t,n,i,j,len; double maxx; scanf("%d",&t); while(t--) { scanf("%d",&n); for(i=1;i<=n;i++) scanf("%d%d%d",&f[i].x,&f[i].y,&f[i].num); len=0; for(i=1;i<n;i++) for(j=i+1;j<=n;j++) { ff[len].u=i; ff[len].v=j; ff[len++].w=distan(i,j); } qsort(ff,len,sizeof(ff[0]),cmp); min_tree(n,len); // printf("%.0f ",sum); maxx=-1; for(i=1;i<=n-1;i++) { int e=fp[i]; memset(vis,0,sizeof(vis)); vis[ff[e].v]=1; int maxx1=dfs(ff[e].u); memset(vis,0,sizeof(vis)); vis[ff[e].u]=1; int maxx2=dfs(ff[e].v); if(maxx<1.0*(maxx1+maxx2)/(sum-ff[e].w)) maxx=1.0*(maxx1+maxx2)/(sum-ff[e].w); } printf("%.2f ",maxx); } return 0;}