hdu4081次小生成树

最新推荐文章于 2021-04-09 08:11:34 发布

转载最新推荐文章于 2021-04-09 08:11:34 发布 · 54 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/acjiumeng/p/7122488.html

先求一遍最小生成树，然后遍历所有边，如果这条边在最小生成树中就直接减去这条边的距离，如果不在最小生成树中，那么就构成了一个环，此时需要减去最小生成树中最大的边，即求次小生成树时的maxx，

有一点要注意当求maxx最大值时j！=k，虽然不知道原理是什么。。。。如果有大佬知道求告知

好像求次小生成树时没有这个条件，据说还能prim+dfs做。。。。

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<vector>
#include<cstdio>
#include<cassert>
#include<iomanip>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define pi acos(-1)
#define ll long long
#define mod 1000000007
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;

const double g=10.0,eps=1e-9;
const int N=1000+10,maxn=111117,inf=0x3f3f3f;

struct edge{
    double x,y,p;
}e[N];
double c[N][N],d[N];
double maxx[N][N];
int pre[N],n;
bool vis[N],used[N][N];
double dis(edge a,edge b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
double prim()
{
    memset(vis,0,sizeof vis);
    memset(used,0,sizeof used);
    memset(maxx,0,sizeof maxx);
    for(int i=2;i<=n;i++)
    {
        pre[i]=1;
        d[i]=c[1][i];
    }
    vis[1]=1;
    pre[1]=0;
    d[1]=0;
    double ans=0;
    for(int i=1;i<n;i++)
    {
        double mind=inf;
        int k;
        for(int j=1;j<=n;j++)
        {
            if(!vis[j]&&mind>d[j])
            {
                mind=d[j];
                k=j;
            }
        }
        vis[k]=1;
        ans+=mind;
        used[k][pre[k]]=used[pre[k]][k]=1;
        for(int j=1;j<=n;j++)
        {
            if(vis[j]&&k!=j)maxx[j][k]=maxx[k][j]=max(maxx[j][pre[k]],d[k]);
            if(!vis[j]&&d[j]>c[k][j])
            {
                d[j]=c[k][j];
                pre[j]=k;
            }
        }
    }
    return ans;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout<<setiosflags(ios::fixed)<<setprecision(2);
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        for(int i=1;i<=n;i++)
            cin>>e[i].x>>e[i].y>>e[i].p;
        for(int i=1;i<=n;i++)
        {
            c[i][i]=0;
            for(int j=i+1;j<=n;j++)
            {
                c[i][j]=c[j][i]=dis(e[i],e[j]);
            }
        }
        double B=prim(),ans=-1;
        for(int i=1;i<=n;i++)
        {
            for(int j=i+1;j<=n;j++)
            {
                if(used[i][j])ans=max(ans,(e[i].p+e[j].p)/(B-c[i][j]));
                else ans=max(ans,(e[i].p+e[j].p)/(B-maxx[i][j]));
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}

View Code

转载于:https://www.cnblogs.com/acjiumeng/p/7122488.html