hdu 4081 最小生成树变形

题目意思是给你一个完全图，每个点有一个权值。

求一棵生成树。你可以使其中一条边权值变成0.使得该边两点的权值和除以生成树的边劝和最小。

思路：

    具有最优解的树一定有这么一个性质。除了权值变成0的边外，其余的边都必须是最小生成树的边。

    基于这点，我们就可以先求出最小生成树，然后枚举的删边。删了一条边之后就变成了两个独立集，

    而我们要在这两个集合中找一条边，使得它具有最优解。注意到，分母是固定的，所以只需在两边找到具有最大权值的点即可。

    最小生成树的复杂度n^2 枚举删边加上找点的复杂度也是n^2.最后复杂度是n^2

#include <iostream>

#include <stdio.h>

#include <string.h>

#include <math.h>

using namespace std;

const int maxn=1111,inf=99999999;

double d[maxn][maxn],dist;

int n;

double ans;

struct

{

    int x,y,p;

}a[maxn];

double f(int t,int s)

{

    double x=a[t].x-a[s].x;

    x*=x;

    double y=a[t].y-a[s].y;

    y*=y;

    return(sqrt(x+y));

}

int text[maxn],fa[maxn];

double key[maxn];

int min()

{

    int txt=inf,ted;

    for(int i=1;i<=n;i++)

    if(key[i]<txt&&!text[i])

    {

        txt=key[i];

        ted=i;

    }

    return(ted);

}

int head[maxn],lon;

struct

{

    int to,next;

    double w;

}e[111111];

int edgemake(int from,int to,double w)

{

    e[++lon].w=w;

    e[lon].to=to;

    e[lon].next=head[from];

    head[from]=lon;

}

int prim()

{

    memset(key,111,sizeof(key));

    memset(text,0,sizeof(text));

    text[1]=1;

    key[1]=0;

    for(int i=1;i<=n;i++)

    {

        key[i]=d[1][i];

        fa[i]=1;

    }

    for(int k=2;k<=n;k++)

    {

        int u=min();

        edgemake(u,fa[u],key[u]);

        edgemake(fa[u],u,key[u]);

        text[u]=1;

        dist+=key[u];

        for(int i=1;i<=n;i++)

        if(d[u][i]<key[i])

        {

            key[i]=d[u][i];

            fa[i]=u;

        }

    }

}

int dfs(int t,int s,int sum)

{

    sum=sum>a[t].p?sum:a[t].p;

    for(int k=head[t];k!=-1;k=e[k].next)

    {

        int u=e[k].to;

        if(u==s) continue;

        if(text[u]) continue;

        text[u]=1;

        int txt=dfs(u,s,sum);

        sum=sum>txt?sum:txt;

        text[u]=0;

    }

    return(sum);

}

int main()

{

    int tcase;

    scanf("%d",&tcase);

    while(tcase--)

    {

        memset(head,-1,sizeof(head));

        lon=0;

        scanf("%d",&n);

        for(int i=1;i<=n;i++)

        {

            int x,y,p;

            scanf("%d %d %d",&x,&y,&p);

            a[i].x=x;

            a[i].y=y;

            a[i].p=p;

        }

        for(int i=1;i<=n;i++)

        for(int j=1;j<=n;j++)

        d[i][j]=f(i,j);

        dist=0;

        prim();

        ans=0;

        for(int i=1;i<=n;i++)

        {

            memset(text,0,sizeof(text));

            for(int k=head[i];k!=-1;k=e[k].next)

            {

                int txt1=dfs(i,e[k].to,0);

                int txt2=dfs(e[k].to,i,0);

                double tmp=(txt1+txt2)/(dist-e[k].w);

                ans=ans>tmp?ans:tmp;

            }

        }

        printf("%.2f\n",ans);

    }

    return 0;

}