【poj2728】Desert King

题意

最优比率生成树。

http://blog.csdn.net/ophunter_lcm/article/details/10113817

分析

Dinkelbach算法,通过迭代的思想进行参数搜索。

每一次的参数搜索,求出来的最小生成树又会生成一个更优的参数,然后用新的参数继续搜...以此类推。

为什么生成的参数会更优?

这是因为我们要最小化\(\sum C\over \sum D\),我们设答案为\({\sum C\over \sum D}<L\),这等价于\(\sum(C-DL)<0\)

我们满足以上等式,就是说明存在\(C,D\)满足\(\sum(C-DL)<0\),所以逆推回来就有\(L>{\sum C\over \sum D}\),且\(\sum C\over \sum D\)是可取的,所以就进一步的逼近了。

可以证明逼近最优解的次数为\(O(\log NC)\),即是一个\(\log\)级别的次数。

所以可以大胆地去用。

它跑得一般比二分还快。

但是适用的范围没有二分广。

总之具体问题具体分析。

为什么我的姿势这么挫,跑得这么慢...

#include <cstdio>
#include <cstring>
#include <cctype>
#include <climits>
#include <cmath>
#include <algorithm>
using namespace std;

#define rep(i,a,b) for (int i=(a);i<=(b);i++)

const int N=1024;

const int INF=INT_MAX;

const double MAX=1e10;
const double EPS=1e-4;

int n;
double x[N],y[N],z[N];

double d[N][N];
double h[N][N];

double ans;
double mp[N][N];
int vis[N]; double dis[N]; int src[N];
double sumD,sumH;

inline int cmp(double a,double b)
{
    if (fabs(a-b)<EPS) return 0;
    return a<b?-1:1;
}

inline double GetDis(double xi,double yi,double xj,double yj)
{
    double dx=xi-xj,dy=yi-yj;
    return sqrt(pow(dx,2)+pow(dy,2));
}

double Prim(void)
{
    double sum=0;

    memset(mp,0,sizeof mp);
    memset(vis,0,sizeof vis); memset(dis,0,sizeof dis); memset(src,0,sizeof src);
    sumD=sumH=0;

    rep(i,1,n) rep(j,1,n)
        mp[i][j]=h[i][j]-ans*d[i][j];
    vis[1]=1; dis[1]=0; rep(i,2,n) dis[i]=mp[1][i],src[i]=1;

    rep(tms,1,n-1)
    {
        double t=MAX; int frm=0,to=0;
        rep(i,1,n) if (!vis[i])
            if (dis[i]<t) t=dis[i],frm=src[i],to=i;

        vis[to]=1;
        sumD+=d[frm][to],sumH+=h[frm][to];
        sum+=t;

        rep(i,1,n) if (!vis[i])
            if (mp[to][i]<dis[i]) dis[i]=mp[to][i],src[i]=to;
    }

    return sum;
}

int main(void)
{
    #ifndef ONLINE_JUDGE
    freopen("poj2728.in","r",stdin);
    freopen("poj2728.out","w",stdout);
    #endif

    rep(tms,1,INF)
    {
        scanf("%d",&n); if (!n) break;
        memset(x,0,sizeof x); memset(y,0,sizeof y); memset(z,0,sizeof z);
        rep(i,1,n) scanf("%lf%lf%lf",x+i,y+i,z+i);

        memset(d,0,sizeof d); memset(h,0,sizeof h);
        rep(i,1,n) rep(j,1,n)
        {
            d[i][j]=GetDis(x[i],y[i],x[j],y[j]);
            h[i][j]=fabs(z[i]-z[j]);
        }

        ans=MAX;
        while (1)
        {
            double t=Prim();
            if (cmp(t,0)!=0)
                ans=sumH/sumD;
            else break;
        }
        printf("%0.3lf\n",ans);
    }

    return 0;
}

【poj2728】Desert King的更多相关文章

  1. 【POJ2728】Desert King(分数规划)

    [POJ2728]Desert King(分数规划) 题面 vjudge 翻译: 有\(n\)个点,每个点有一个坐标和高度 两点之间的费用是高度之差的绝对值 两点之间的距离就是欧几里得距离 求一棵生成 ...

  2. 【POJ2728】Desert King - 01分数规划

    Description David the Great has just become the king of a desert country. To win the respect of his ...

  3. 【POJ2728】Desert King 最优比率生成树

    题目大意:给定一个 N 个点的无向完全图,边有两个不同性质的边权,求该无向图的一棵最优比例生成树,使得性质为 A 的边权和比性质为 B 的边权和最小. 题解:要求的答案可以看成是 0-1 分数规划问题 ...

  4. 【POJ 2728 Desert King】

    Time Limit: 3000MSMemory Limit: 65536K Total Submissions: 27109Accepted: 7527 Description David the ...

  5. 【HDOJ】【1512】Monkey King

    数据结构/可并堆 啊……换换脑子就看了看数据结构……看了一下左偏树和斜堆,鉴于左偏树不像斜堆可能退化就写了个左偏树. 左偏树介绍:http://www.cnblogs.com/crazyac/arti ...

  6. POJ2728:Desert King——题解

    http://poj.org/problem?id=2728 题目大意:求一棵生成树使得路费用和/路长之和最小(路的费用是两端点的高度差) 最小比率生成树. 我们还是01分数规划的思想将边权变为路费用 ...

  7. poj2728 Desert King【最优比率生成树】【Prim】【0/1分数规划】

    含[最小生成树Prim]模板. Prim复杂度为$O(n^2),适用于稠密图,特别是完全图的最小生成树的求解.   Desert King Time Limit: 3000MS   Memory Li ...

  8. POJ2728 Desert King 【最优比率生成树】

    POJ2728 Desert King Description David the Great has just become the king of a desert country. To win ...

  9. 【BZOJ】1087: [SCOI2005]互不侵犯King

    [算法]状态压缩型DP [题解]http://www.cnblogs.com/xtx1999/p/4620227.html (orz) https://www.cnblogs.com/zbtrs/p/ ...

随机推荐

  1. Bitmap的分析与使用.md

    https://github.com/GeniusVJR/LearningNotes/blob/master/Part1/Android/Bitmap的分析与使用.md

  2. 关于C指针

    地址概念:内存基本单元是一个字节,一个字节8个位.在定义变量的时候,如int a=10:系统为变量a分配2个字节空间:1000~1001,并在1001~1002中存有数据10.内存中没有变量a,只有1 ...

  3. java properties读取与设值

    import java.io.FileInputStream;import java.io.FileNotFoundException;import java.io.FileOutputStream; ...

  4. VBA中操作XML

    OFFICE2007之后使用了OpenXml标准(伟大的改变),定制文本级的Ribbon可以通过修改压缩包内的xml文件来实现. 先学习一下VBA中操作XML的方法 先引用Microsoft XML ...

  5. ABAP Enhancement:第一部分

    声明:原创作品,转载时请注明文章来自SAP师太技术博客( 博/客/园www.cnblogs.com):www.cnblogs.com/jiangzhengjun,并以超链接形式标明文章原始出处,否则将 ...

  6. JVM 1.类的加载、连接、初始化

    Java类的加载是由类加载器来完成的,过程如下: 首先,加载是把硬盘.网络.数据库等的class文件中的二进制数据加载到内存的过程,然后会在Java虚拟机的运行时数据区的堆区创建一个Class对象,用 ...

  7. NOJ 1063 生活的烦恼

    描述 生活的暑假刚集训开始,他要决心学好字典树,二叉树,线段树和各种树,但生活在OJ上刷题的时候就遇到了一个特别烦恼的问题.那当然就是他最喜欢的二二叉树咯!题目是这样的:给你一颗非空的二叉树,然后再给 ...

  8. 使用Select命令创建菜单

    创建文本菜单的一半功夫都花在了创建菜单布局和获取输入的字符上.bash shell提供了一个很容易上手的小工具来自动完成这些工作select命令允许从单个命令行创建菜单,然后在提取输入的答案并自动处理 ...

  9. 如何删除github里面的文件夹?

    按照以下步骤即可(本地删除) 1. git pull you git url2. git checkout 3. rm -r dirName4. git add --all5. git commit  ...

  10. SAP FI/CO凭证不一致的解决办法

    First, use program RKACOR20 to delete the incorrect CO documents. OKBA - Transfer FI Documents to CO ...