AcWing 348. 沙漠之王 (最优比率生成树)

(1)01分数规划问题
(2)因为要求比率最大所以求最小生成树,由于是完全图,所以使用prim


int n;//村庄数目
int x[1111],y[1111],h[1111];
double a[1111][1111],b[1111][1111];
//长度:距离 / 成本:高度差(保证不为0)最大

int vis[1111];
double dis[1111];//到生成树的距离
int check(double k)//最小生成树
{
    rpp(i,n) dis[i]=a[1][i]-k*b[1][i];
    memset(vis,0,sizeof(vis));
    double mint=0.0;
    vis[1]=1;
    rpp(_,n-1)//每轮挑出一个点加入树
    {
        int id=1;
        double mn=-1;
        rpp(i,n) if(vis[i]==0&&(mn==-1||mn>dis[i])) mn=dis[i],id=i;
        mint+=mn,vis[id]=1;
        rpp(i,n) if(vis[i]==0) dis[i]=min(dis[i],a[id][i]-k*b[id][i]);
    }
    return mint>=0.0;
}
signed main()
{
    while(cin>>n)
    {
        if(n==0) break;
        rpp(i,n)  cin>>x[i]>>y[i]>>h[i];
        rpp(i,n) 
        {
            for(int j=i+1;j<=n;++j)
            {
                a[i][j]=a[j][i]=fabs(h[i]-h[j]);
                b[i][j]=b[j][i]=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
            }
        }
        double l=0.0,r=10000001.0;//
        while(r-l>1e-6)
        {
            double mid=(l+r)/2.0;
            if(check(mid)) l=mid;
            else r=mid;
        }
        printf("%.3f\n",l);
    }
    stop;
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值