hdu 4081 次小生成树

/*
题意：每个城市有人口，现在有一棵最小生成树，将其中一条边割开，然后选出两个城市的人口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\n",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\n",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\n",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\n",maxx);
    }
return 0;}