zoukankan      html  css  js  c++  java
  • 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
    ",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;}


  • 相关阅读:
    【转载】以太网帧类型速查表
    oracle 无法解析指定的连接标识符
    最近观影有感~
    myeclipse 8.5 安装svn插件
    Symbian S60 设备配置 IMAP
    VS 2008 自动缩进修改
    RGB YUV [转]
    NoteExpress 无法访问国外数据库
    小s坐月子
    2011124 code
  • 原文地址:https://www.cnblogs.com/thefirstfeeling/p/4410559.html
Copyright © 2011-2022 走看看