POJ 2728 Desert King

最新推荐文章于 2021-01-21 20:49:31 发布
最新推荐文章于 2021-01-21 20:49:31 发布
阅读量570 收藏
点赞数
分类专栏： POJ - 最优比率生成树 ######## 图论 ####### +++ 最小生成树 +++ 文章标签： 最优比率生成树 Prim 二分 迭代 最大生成树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013569656/article/details/38588687
版权

题目大意：

        沙漠中有N个城市编号1 ~ N（2 ≤ N ≤ 1,000），首都是1号城市，只有1号城市有水源，现在需要修管道使首都的水能引导其它所哟城市，要求其它所有城市到首都的路径唯一，但是各个城市位于不同的海拔，而水管只能是水平铺设，因此连接一个水管的两个城市还需要一个抽水泵以便将水从低处导向高处或从高处引向低处，现保证输入的各个城市在水平平面内不会出现三点共线的情况，并且各个城市之间海拔不相同。

        现有多个测例，每个测例都给出N个城市的坐标信息"x y z"，表示该城市的坐标为(x, y, z)，其中0 ≤ x, y < 10,000，0 ≤ z < 10,000,000，都为整型数，其中管道长度即为水平平面内的欧几里得距离（即xOy平面），抽水泵的费用为两个城市之间的高度差，现在要求出如何铺设管道才能使抽水泵费用总和比上管道总长最小，并输出该最小比率并四舍五入至3为小数，输入以N = 0表示结束。

背景：

        由于各点到1号点路径唯一，则必定是求一棵生成树，可以把各海拔差看做各点的点权值，若以1号点作为参考海拔，则1号点的点权值即为0，其余点的权值即为高度差，而边权值即为两点之间的欧几里得距离。

        即求r = (C1+C2+…+Cn-1)/(L1+L2+…+Ln-1)，其中r为比率，C为点权值，L为边权值，共N - 1个点（1号点权值为0）N - 1条边（生成树），并且要求r最小。

        现在设f(r) = r*L1-C1 + r*L2-C2 +…+ r*Ln-1-Cn-1，则f(r_min) ≤ 0且f(r_max) ≥ 0，并且f(r)是单调递增的，当然r_min可以取无穷小但是这是没有意义的，当r_min为无穷小时f(r)必定<0，因此可以先大致定一个较大的范围[r_min, r_max]，使得f(r_min) < 0且f(r_max) > 0，并用二分法使r逼近f(r)刚好≥0的点。

        那现在的问题是该如何求f(r)呢？以r*Li-Ci为边求最大生成树（Prim+堆优化）就行了，因此二分法求最优比率生成树的算法就描述完了。

        但是还有更快的方法，那就是迭代求f(r)，现定一个初始值r0，然后每次用Prim时都会得到一个对应的原图的生成树，prim返回两个值，一个是f(r)还有一个就是原图对应的生成树的r（r由prim直接返回，f(r)可以通过静态变量返回），因此迭代算法为for ( r = r0; |f(r)| < ε; r = prim(r) )即可，迭代法求出的f(r)无限向0逼近，并且逼近时围绕0上下都抖动。

题目链接

二分：

注释代码：

/*                                   
 * Problem ID : POJ 2728 Desert King 
 * Author     : Lirx.t.Una                                   
 * Language   : C++                       
 * Run Time   : 2469 ms                                   
 * Run Memory : 26568 KB                                   
*/ 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>

//城市最大数量
#define	MAXN	1000
//二分精度
#define	EPS		1E-5

using namespace std;

struct	Point {//各个城市坐标信息

	int		x, y, z;

	friend istream &
	operator>>(istream &is, Point &p) {
	
		is >> p.x >> p.y >> p.z;
		return is;
	}

	double
	POW(int x) { return x * x; }

	double
	operator^(Point &oth) {//平面上的距离
	
		return sqrt(POW( x - oth.x ) + POW( y - oth.y ));
	}

	int
	operator&(Point &oth) {//海拔差
	
		return abs(z - oth.z);
	}
};

struct	Node {//Prim+堆优化的结点

	int		u;
	double	d;

	Node(void) {}

	Node(int uu, double dd) : u(uu), d(dd) {}

	bool//大顶堆求最大生成树
	operator<(const Node &oth)
	const {
	
		return d < oth.d;
	}
};

double	g[MAXN + 1][MAXN + 1];//r*Li-Ci图中各点之间的边的权值
Point	p[MAXN + 1];//各个点
//原图上个点之间的平面欧几里得距离以及海拔差
double	dist[MAXN + 1][MAXN + 1];
int		high[MAXN + 1][MAXN + 1];

//Prim数据域
bool	vis[MAXN + 1];
double	d[MAXN + 1];

int		n;//城市个数

double
prim(double r) {

	int		i;
	int		u, v;
	int		ne;

	double	ans;//最大生成树返回的f(r)

	Node	node;

	priority_queue<Node>	heap;

	memset(vis, 0, sizeof(vis));

	for ( u = 1; u <= n; u++ )
		for ( v = u + 1; v <= n; v++ )//Prim需要用到全边g[u][v]、g[v][u]
			g[u][v] = g[v][u] = r * dist[u][v] - high[u][v];//r*Li-Ci构造新图

	vis[1] = true;
	for ( i = 2; i <= n; i++ ) {
	
		d[i] = g[1][i];
		heap.push(Node( i, d[i] ));
	}

	ans = 0.0;
	ne	= 1;
	while (true) {
	
		while (true) {
		
			node = heap.top();
			heap.pop();

			if ( !vis[ u = node.u ] ) {

				vis[u] = true;
				ans += node.d;
				ne++;

				break;
			}
		}

		if ( ne == n ) break;

		for ( v = 2; v <= n; v++ )
			if ( !vis[v] && g[u][v] > d[v] ) {
			
				d[v] = g[u][v];
				heap.push(Node( v, d[v] ));
			}
	}

	return ans;
}

int
main() {

	int		i, j;

	double	lft, rht, mid;//二分
	double	r;//二分比率

	while ( scanf("%d", &n), n ) {
	
		for ( i = 1; i <= n; i++ ) cin >> p[i];
		for ( i = 1; i <= n; i++ )//只需要用到半边	
			for ( j = i + 1; j <= n; j++ ) {
			
				dist[i][j] = p[i] ^ p[j];
				high[i][j] = p[i] & p[j];
			}

		//初始二分区间（经验值）
		lft = 0.0;
		rht = 100.0;

		while ( rht - lft > EPS ) {
		
			mid = ( lft + rht ) * 0.5;
			if ( prim(mid) >= 0.0 ) rht = mid;
			else lft = mid;
		}
		printf("%.3lf\n", rht);//一定要求≥0的一端
	}

	return 0;
}
无注释代码： 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>

#define	MAXN	1000
#define	EPS		1E-5

using namespace std;

struct	Point {

	int		x, y, z;

	friend istream &
	operator>>(istream &is, Point &p) {
	
		is >> p.x >> p.y >> p.z;
		return is;
	}

	double
	POW(int x) { return x * x; }

	double
	operator^(Point &oth) {
	
		return sqrt(POW( x - oth.x ) + POW( y - oth.y ));
	}

	int
	operator&(Point &oth) {
	
		return abs(z - oth.z);
	}
};

struct	Node {

	int		u;
	double	d;

	Node(void) {}

	Node(int uu, double dd) : u(uu), d(dd) {}

	bool
	operator<(const Node &oth)
	const {
	
		return d < oth.d;
	}
};

double	g[MAXN + 1][MAXN + 1];
Point	p[MAXN + 1];
double	dist[MAXN + 1][MAXN + 1];
int		high[MAXN + 1][MAXN + 1];

bool	vis[MAXN + 1];
double	d[MAXN + 1];

int		n;

double
prim(double r) {

	int		i;
	int		u, v;
	int		ne;

	double	ans;

	Node	node;

	priority_queue<Node>	heap;

	memset(vis, 0, sizeof(vis));

	for ( u = 1; u <= n; u++ )
		for ( v = u + 1; v <= n; v++ )
			g[u][v] = g[v][u] = r * dist[u][v] - high[u][v];

	vis[1] = true;
	for ( i = 2; i <= n; i++ ) {
	
		d[i] = g[1][i];
		heap.push(Node( i, d[i] ));
	}

	ans = 0.0;
	ne	= 1;
	while (true) {
	
		while (true) {
		
			node = heap.top();
			heap.pop();

			if ( !vis[ u = node.u ] ) {

				vis[u] = true;
				ans += node.d;
				ne++;

				break;
			}
		}

		if ( ne == n ) break;

		for ( v = 2; v <= n; v++ )
			if ( !vis[v] && g[u][v] > d[v] ) {
			
				d[v] = g[u][v];
				heap.push(Node( v, d[v] ));
			}
	}

	return ans;
}

int
main() {

	int		i, j;

	double	lft, rht, mid;
	double	r;

	while ( scanf("%d", &n), n ) {
	
		for ( i = 1; i <= n; i++ ) cin >> p[i];
		for ( i = 1; i <= n; i++ )	
			for ( j = i + 1; j <= n; j++ ) {
			
				dist[i][j] = p[i] ^ p[j];
				high[i][j] = p[i] & p[j];
			}

		lft = 0.0;
		rht = 100.0;

		while ( rht - lft > EPS ) {
		
			mid = ( lft + rht ) * 0.5;
			if ( prim(mid) >= 0.0 ) rht = mid;
			else lft = mid;
		}
		printf("%.3lf\n", rht);
	}

	return 0;
}

迭代：

注释代码：

/*                                    
 * Problem ID : POJ 2728 Desert King  
 * Author     : Lirx.t.Una                                    
 * Language   : C++                        
 * Run Time   : 391 ms                                    
 * Run Memory : 26424 KB                                    
*/ 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>

#define	MAXN	1000
//迭代的精度可以更加低一点
#define	EPS		1E-4

using namespace std;

struct	Point {

	int		x, y, z;

	friend istream &
	operator>>(istream &is, Point &p) {
	
		is >> p.x >> p.y >> p.z;
		return is;
	}

	double
	POW(int x) { return x * x; }

	double
	operator^(Point &oth) {
	
		return sqrt(POW( x - oth.x ) + POW( y - oth.y ));
	}

	int
	operator&(Point &oth) {
	
		return abs(z - oth.z);
	}
};

struct	Node {

	int		u;
	double	d;

	Node(void) {}

	Node(int uu, double dd) : u(uu), d(dd) {}

	bool
	operator<(const Node &oth)
	const {
	
		return d < oth.d;
	}
};

double	g[MAXN + 1][MAXN + 1];
Point	p[MAXN + 1];
double	dist[MAXN + 1][MAXN + 1];
int		high[MAXN + 1][MAXN + 1];
int		pre[MAXN + 1];

bool	vis[MAXN + 1];
double	d[MAXN + 1];

int		n;

double 	ans;//r*Li-Ci图的最大生成树的边权之和

double
prim(double r) {

	int		i;
	int		u, v;
	int		ne;

	double	we, wv;//原图中生成树的边权和点权

	Node	node;

	priority_queue<Node>	heap;

	memset(vis, 0, sizeof(vis));

	for ( u = 1; u <= n; u++ )
		for ( v = u + 1; v <= n; v++ )
			g[u][v] = g[v][u] = r * dist[u][v] - high[u][v];

	vis[1] = true;
	for ( i = 2; i <= n; i++ ) {
	
		d[i] = g[1][i];
		pre[i] = 1;
		heap.push(Node( i, d[i] ));
	}

	ans = 0.0;
	ne	= 1;
	we = 0.0;
	wv = 0.0;
	while (true) {
	
		while (true) {
		
			node = heap.top();
			heap.pop();

			if ( !vis[ u = node.u ] ) {

				vis[u] = true;
				ans += node.d;
				we += dist[ pre[u] ][u];
				wv += high[ pre[u] ][u];
				ne++;

				break;
			}
		}

		if ( ne == n ) break;

		for ( v = 2; v <= n; v++ )
			if ( !vis[v] && g[u][v] > d[v] ) {
			
				d[v] = g[u][v];
				pre[v] = u;
				heap.push(Node( v, d[v] ));
			}
	}

	return wv / we;//返回原图对应生成树的比率
}

int
main() {

	int		i, j;
	double	r;

	while ( scanf("%d", &n), n ) {
	
		for ( i = 1; i <= n; i++ ) cin >> p[i];
		for ( i = 1; i <= n; i++ )	
			for ( j = i + 1; j <= n; j++ ) {//由于prim中需要计算原图的r，因此要计算全边
			
				dist[i][j] = dist[j][i] = p[i] ^ p[j];
				high[i][j] = high[j][i] = p[i] & p[j];
			}

		//迭代初值取经验值9，ans只要>EPS就行了
		for ( r = 9, ans = 1; fabs(ans) > EPS; r = prim(r) );
		printf("%.3lf\n", r);
	}

	return 0;
}
无注释代码： 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>

#define	MAXN	1000
#define	EPS		1E-4

using namespace std;

struct	Point {

	int		x, y, z;

	friend istream &
	operator>>(istream &is, Point &p) {
	
		is >> p.x >> p.y >> p.z;
		return is;
	}

	double
	POW(int x) { return x * x; }

	double
	operator^(Point &oth) {
	
		return sqrt(POW( x - oth.x ) + POW( y - oth.y ));
	}

	int
	operator&(Point &oth) {
	
		return abs(z - oth.z);
	}
};

struct	Node {

	int		u;
	double	d;

	Node(void) {}

	Node(int uu, double dd) : u(uu), d(dd) {}

	bool
	operator<(const Node &oth)
	const {
	
		return d < oth.d;
	}
};

double	g[MAXN + 1][MAXN + 1];
Point	p[MAXN + 1];
double	dist[MAXN + 1][MAXN + 1];
int		high[MAXN + 1][MAXN + 1];
int		pre[MAXN + 1];

bool	vis[MAXN + 1];
double	d[MAXN + 1];

int		n;

double 	ans;

double
prim(double r) {

	int		i;
	int		u, v;
	int		ne;

	double	we, wv;

	Node	node;

	priority_queue<Node>	heap;

	memset(vis, 0, sizeof(vis));

	for ( u = 1; u <= n; u++ )
		for ( v = u + 1; v <= n; v++ )
			g[u][v] = g[v][u] = r * dist[u][v] - high[u][v];

	vis[1] = true;
	for ( i = 2; i <= n; i++ ) {
	
		d[i] = g[1][i];
		pre[i] = 1;
		heap.push(Node( i, d[i] ));
	}

	ans = 0.0;
	ne	= 1;
	we = 0.0;
	wv = 0.0;
	while (true) {
	
		while (true) {
		
			node = heap.top();
			heap.pop();

			if ( !vis[ u = node.u ] ) {

				vis[u] = true;
				ans += node.d;
				we += dist[ pre[u] ][u];
				wv += high[ pre[u] ][u];
				ne++;

				break;
			}
		}

		if ( ne == n ) break;

		for ( v = 2; v <= n; v++ )
			if ( !vis[v] && g[u][v] > d[v] ) {
			
				d[v] = g[u][v];
				pre[v] = u;
				heap.push(Node( v, d[v] ));
			}
	}

	return wv / we;
}

int
main() {

	int		i, j;
	double	r;

	while ( scanf("%d", &n), n ) {
	
		for ( i = 1; i <= n; i++ ) cin >> p[i];
		for ( i = 1; i <= n; i++ )	
			for ( j = i + 1; j <= n; j++ ) {
			
				dist[i][j] = dist[j][i] = p[i] ^ p[j];
				high[i][j] = high[j][i] = p[i] & p[j];
			}

		for ( r = 9, ans = 1; fabs(ans) > EPS; r = prim(r) );
		printf("%.3lf\n", r);
	}

	return 0;
}
单词解释：

prime：adj, 首席的

lifter：n, 升降机，起重机

vertical：adj, 垂直的

altitude：n, 海拔，高度

elegant：adj, 高雅的

ruler：n, 统治者

dominate：vt, 控制，支配

tyrant：n, 暴君

Terra Cotta Warriors：n, 兵马俑

terracotta：n, 陶瓦，吃土色

life-sized：adj, 原物大小的

mausoleum：n, 陵墓

gigantic：adj, 巨大的，庞大的

undertake：vt, 承担，从事

imperial：adj, 帝国的，皇帝的

unified：adj, 统一的

Warring States Period：n, 战国时代

warring：adj, 交战的，敌对的

__Lirx_t_Una__
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【01规划】POJ 2728 Desert King
AsUs.的博客
05-02 325
POJ 2728 Desert King 题意：给出 n 个点的坐标和它的高度，求一颗生成树使得树上所连边的两点高度差之和除以距离之和最小。 思路：同样构造f(l)方程， 令∑hight / ∑dis = l，那么∑hight = l*∑dis, 令f(l) = ∑(hight - l*dis) ,题目所求是最小值 那么∑hight / ∑dis ≤ l 说明存在更优解，即f(l) ≤ 0 ，二分 l 得出答案 关于f(l)和二分边界判断 ： 我们设 cost = hight - l*dis ,把它
poj2728 Desert King
dfssbqyrq96616258的博客
01-29 122
大概题意： 每两个点中的边权有两个：一个是两点坐标的欧几里得距离（horizontal distance），暂且成为ai,第二个是两点的海拔之差，称为bi.然后需要一个生成树使sum(ai)\sum(bi)最小。 这里可以引入分数规划：我们设ai\bi=k,那么ai-bi*k=0 我们只需要二分一个值mid，当ai-bi*mid=0时，这时的mid便是最优值。 对于每...
POJ2728 Desert King
forezxl的博客
05-23 182
01分数规划 题目传送门 题目大意：有nnn个点，两点之间都有道路，道路有两个权值：len=len=len=两点之间距离，cost=cost=cost=两点高度之差。求使∑cost∑len∑cost∑len\frac {\sum cost}{\sum len}最小的生成树的该比值。 这题就是求最优比例生成树。二分比值midmidmid，把每条边的权值改为cost[i]−mid∗len[i]c...
poj2728 desert king
Elijahqi
05-04 230
http://www.elijahqi.win/archives/1136 Description David the Great has just become the king of a desert country. To win the respect of his people, he decided to build channels all over his country to...
poj2728 Desert King（最优比率生成树【Prim）
Coco_T的博客
03-04 258
题目链接 分析： 最小化：∑cost∑len∑cost∑len\sum cost \over \sum len 01分数规划：最优比率生成树 渠道的长度是两个村庄之间的水平距离，通道的成本是升降机的高度 题目约定： He just needs to build the necessary channels to bring water to all the villages, whic...
POJ2728 Desert King - （0/1）分数规划
Zolrk的博客
08-07 213
题目求一种方案，使得图全连通并且所有边费用与距离之商最小 ∑i∈ecosti∑i∈ecosti\sum_{i∈e}cost_i除以∑i∈edisi∑i∈edisi\sum_{i∈e}dis_i最小 可以考虑二分求解 可以假设这个值小于等于L时存在一个解，然后检查是否存在这个解，如果不存在说明L取小了 问题是为什么要假设“存在”，事实上如果假设“任意”，那么就要检查每种可能都要小于，就很麻烦...
POJ2728Desert King
AcerMoOi之路
07-04 685
最优比率生成树问题，就是01分数规划在树上的应用是一类求解在树上选一些边，是这些边的花费比长度最大的问题，多次prime就好代码//By AcerMo #include&lt;cmath&gt; #include&lt;cstdio&gt; #include&lt;cstring&gt; #include&lt;iostream&gt; #include&lt;algorithm&gt; usin...
poj 2728 Desert King
weixin_30691871的博客
08-27 87
Description David the Great has just become the king of a desert country. To win the respect of his people, he decided to build channels all over his country to bring water to every village. Vil...
Poj 2728 Desert King
aichangjie8078的博客
07-28 97
http://poj.org/problem?id=2728 题意：一个无向图，每条边有两个边权，cost和len,求出一个最小生成树，使所有边的cost/len和最小。 题解：最优比率生成树：http://www.cnblogs.com/lotus3x/archive/2009/03/21/1418480.html 　　　使用迭代法，代码抄的http://www.cnblogs....
POJ2728 Desert King 最优比率生成树
Luowaterbi的博客
01-21 218
一开始没理解最优比率生成树为什么要用最小生成树。 书上说的是最大生成树，我一开始觉得是最大生成树，但这样求出来是合法值的最大值，也不是题目上的最小值啊。后来我发现，是在0/1规划二分的意义上理解出错了。 在一般的二分答案中，比如函数是单调递减的，>0的都不是解，<0的都是解，=0的是我们可能要求的解里面的最大值。但是对于0/1规划，二分的是唯一解，>0、<0的都不是解，只有=0的是解。 我们设二分的比率是kkk，树用xxx表示，最小生成树用xminx_{min}xmin​表示，最大生
poj 2728 Desert King (最优比率生成树)
越看越喜欢啊
02-25 356
poj 2728 Desert King (最优比率生成树)
poj 2728 Desert King（最小比率生成树 / 0-1分数规划 / 二分）
繁凡さん的博客
09-11 224
二分答案，我们要找最小的答案，如果有更小的答案说明 ∑W−Z∗∑L<=0∑W−Z∗∑L <= 0∑W−Z∗∑L<=0。 #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> using namespace std; const int N = 500007, INF = 0x3f3f3f3f; const double eps = 1e-6; do..
[poj2728]Desert King
StaroForgin的博客
09-23 161
Desert King 题解 最优比率生成树的一道板子题。 很容易发现，我们要求的是最小的生产树上边集，使得最大。 我们可以考虑0/1分数规划。 记表示第条边是否属于集合，答案变成了，其中表示这个图的边集。 设表示这个比值，构造函数。 其表示将边权修改为后，最小生成树的值。 而当使得时，我们就得到了答案。 由于本题是一个稠密图，，所以需采用prim求最小生成树。 先二分，通过生成最小生成树来求出当前的值，二分左右边界不能直接设为最大值，会T，需设为与。 最后输出答案需要用%f，%lf会
[POJ 2728] Desert King
C202044zxy的博客
09-23 154
一、题目 点此看题
POJ 1797 Heavy Transportation
源代码
03-02 857
题目大意： 注释代码： /* * Problem ID : POJ 1797 Heavy Transportation * Author : Lirx.t.Una * Language : GCC * Run Time : 279 ms * Run Memory
POJ 2485 Highways
源代码
03-07 675
题目大意： 注释代码： 无注释代码： 单词解释：
POJ 1860 Currency Exchange
源代码
08-20 644
题目大意： 题目链接 注释代码：
POJ 2502 Subway
源代码
08-22 632
题目大意： 题目链接 注释代码： 无注释代码： #include #include #include #include #define INF 1E15 #define MAXN 202 #define V_W 10 #define V_S 40 using namespace std; struct Point { int x, y; friend i
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值