【poj2728】Desert King
【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的更多相关文章
- 【POJ2728】Desert King(分数规划)
[POJ2728]Desert King(分数规划) 题面 vjudge 翻译: 有\(n\)个点,每个点有一个坐标和高度 两点之间的费用是高度之差的绝对值 两点之间的距离就是欧几里得距离 求一棵生成 ...
- 【POJ2728】Desert King - 01分数规划
Description David the Great has just become the king of a desert country. To win the respect of his ...
- 【POJ2728】Desert King 最优比率生成树
题目大意:给定一个 N 个点的无向完全图,边有两个不同性质的边权,求该无向图的一棵最优比例生成树,使得性质为 A 的边权和比性质为 B 的边权和最小. 题解:要求的答案可以看成是 0-1 分数规划问题 ...
- 【POJ 2728 Desert King】
Time Limit: 3000MSMemory Limit: 65536K Total Submissions: 27109Accepted: 7527 Description David the ...
- 【HDOJ】【1512】Monkey King
数据结构/可并堆 啊……换换脑子就看了看数据结构……看了一下左偏树和斜堆,鉴于左偏树不像斜堆可能退化就写了个左偏树. 左偏树介绍:http://www.cnblogs.com/crazyac/arti ...
- POJ2728:Desert King——题解
http://poj.org/problem?id=2728 题目大意:求一棵生成树使得路费用和/路长之和最小(路的费用是两端点的高度差) 最小比率生成树. 我们还是01分数规划的思想将边权变为路费用 ...
- poj2728 Desert King【最优比率生成树】【Prim】【0/1分数规划】
含[最小生成树Prim]模板. Prim复杂度为$O(n^2),适用于稠密图,特别是完全图的最小生成树的求解. Desert King Time Limit: 3000MS Memory Li ...
- POJ2728 Desert King 【最优比率生成树】
POJ2728 Desert King Description David the Great has just become the king of a desert country. To win ...
- 【BZOJ】1087: [SCOI2005]互不侵犯King
[算法]状态压缩型DP [题解]http://www.cnblogs.com/xtx1999/p/4620227.html (orz) https://www.cnblogs.com/zbtrs/p/ ...
随机推荐
- (3)redis队列功能
Redis队列功能介绍 List 常用命令: Blpop删除,并获得该列表中的第一元素,或阻塞,直到有一个可用 Brpop删除,并获得该列表中的最后一个元素,或阻塞,直到有一个可用 Brpoplpus ...
- bmp格式图片文件读取
C++读取bmp图片 #include <windows.h> #include <stdio.h> #include <stdlib.h> #include &l ...
- java中为什么byte的取值范围是-128到+127
概念:java中用补码表示二进制数,补码的最高位是符号位,最高位为“0”表示正数,最高位为“1”表示负数.正数补码为其本身:负数补码为其绝对值各位取反加1:例如:+21,其二进制表示形式是000101 ...
- CodeForces 496B Secret Combination
Secret Combination Time Limit:2000MS Memory Limit:262144KB 64bit IO Format:%I64d & %I64u ...
- Create,Insert
创建表 create table people ( id int ,name ) ) create table toys ( id int ,name ) ,people_id int ) CREAT ...
- Java中正则表达式的使用
public class Test{ public static void main(String args[]) { String str="@Shang Hai Hong Qiao Fe ...
- CUBRID学习笔记 30 复制表结构 cubrid教程
语法 CREATE {TABLE | CLASS} <new_table_name> LIKE <old_table_name> 如下 CREATE TABLE a_tbl( ...
- Cocos2d-x优化中纹理优化
转自 http://blog.csdn.net/tonny_guan/article/details/41016241 Cocos2d-x优化中纹理优化 1.纹理像素格式纹理优化工作的另一重要的指标是 ...
- SpringMVC后缀
<!-- 在这里,使用*.html为后缀的URL都能被baobaotao Servlet截获,进而转由SpringMVC框架进行处理.在 Struts框架中,一般将URL后缀配置为*.do:在w ...
- Hibernate各种主键生成策略与配置详解《转》
1.assigned 主键由外部程序负责生成,在 save() 之前必须指定一个.Hibernate不负责维护主键生成.与Hibernate和底层数据库都无关,可以跨数据库.在存储对象前,必须要使用主 ...