网上有一些很数学的证明方法,表示看的挺晕,自己理解了一下后,发表下自己的看法,如果有错误,再进行修改
其实原题就是求 MIN( ∑CiXi / ∑DiXi )Xi∈{0,1} ,对每个生成树,设其比率r=∑CiXi / ∑DiXi ,可得∑CiXi - ∑DiXi * r=0(条件1)
那么对于所有的生成树,显然∑CiXi
- ∑DiXi * min(r) >= 0,当
∑CiXi / ∑DiXi = min(r)时,等号成立。而我们现在不知道min(r)是多少,只好进行枚举,对每个枚举的r ,构建新的权值(Ci-Di*r),然后求最小生成树, 为什么求最小呢? 我的理解就是这是为了寻找使得生成树的总权值为0的可能性,因为只有当其等于0
的时候,才满足了条件1
这个条件, 说明这个r是可行的,并且如果r枚举到值为min(r)时,其最小生成树的的总权值必然恰好等于0,但是如果不能等于0, 比如大于0, 显然是对该r值,所有的生成树上无论如何也满足不了条件1,说明r值就是偏小了。同理如果小于0,r值是偏大的,说明可能存在某些生成树使得满足条件1,而我们的目标是在满足条件1的情况下使得r最小。
根据这个我们可以发现,实际上r的值是可以进行二分查找的。
而也有人给出了更为高效的迭代方法。
代码贴上:
首先是二分法
#include <iostream>
#include <string>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#define MAXN 1005
#define INF 1000000000
#define eps 1e-7
using namespace std;
int n;
double Edge[MAXN][MAXN], lowcost[MAXN];
int nearvex[MAXN];
struct Point
{
int x, y, z;
}p[MAXN];
double cal(int a, int b)
{
return sqrt(1.0 * (p[a].x - p[b].x) * (p[a].x - p[b].x) + 1.0 * (p[a].y - p[b].y) * (p[a].y - p[b].y));
}
double prim(int src, double l)
{
double cost = 0, len = 0;
double sum = 0;
for(int i = 1; i <= n; i++)
{
nearvex[i] = src;
lowcost[i] = abs(p[src].z - p[i].z) - Edge[src][i] * l;
}
nearvex[src] = -1;
for(int i = 1; i < n; i++)
{
double mi = INF;
int v = -1;
for(int j = 1; j <= n; j++)
if(nearvex[j] != -1 && lowcost[j] < mi)
{
v = j;
mi = lowcost[j];
}
if(v != -1)
{
cost += abs(p[nearvex[v]].z - p[v].z);
len += Edge[nearvex[v]][v];
nearvex[v] = -1;
sum += lowcost[v];
for(int j = 1; j <= n; j++)
{
double tmp = abs(p[v].z - p[j].z) - Edge[v][j] * l;
if(nearvex[j] != -1 && tmp < lowcost[j])
{
lowcost[j] = tmp;
nearvex[j] = v;
}
}
}
}
return sum;
}
int main()
{
while(scanf("%d", &n) != EOF && n)
{
for(int i = 1; i <= n; i++)
scanf("%d%d%d", &p[i].x, &p[i].y, &p[i].z);
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
Edge[i][j] = cal(i, j);
double low = 0, high = 10.0; //其实二分20多次已经很足够了
double l = 0.0, r = 100.0, mid;
while(r - l > eps)
{
mid = (l + r) / 2;
if(prim(1, mid) >= 0) l = mid;
else r = mid;
}
printf("%.3f\n", r);
}
return 0;
}
然后是迭代,注意跟二分法不同的就是prim函数的返回值不同
#include <iostream>
#include <string>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#define MAXN 1005
#define INF 1000000000
#define eps 1e-7
using namespace std;
int n;
double Edge[MAXN][MAXN], lowcost[MAXN];
int nearvex[MAXN];
struct Point
{
int x, y, z;
}p[MAXN];
double cal(int a, int b)
{
return sqrt(1.0 * (p[a].x - p[b].x) * (p[a].x - p[b].x) + 1.0 * (p[a].y - p[b].y) * (p[a].y - p[b].y));
}
double prim(int src, double l)
{
double cost = 0, len = 0;
for(int i = 1; i <= n; i++)
{
nearvex[i] = src;
lowcost[i] = abs(p[src].z - p[i].z) - Edge[src][i] * l;
}
nearvex[src] = -1;
for(int i = 1; i < n; i++)
{
double mi = INF;
int v = -1;
for(int j = 1; j <= n; j++)
if(nearvex[j] != -1 && lowcost[j] < mi)
{
v = j;
mi = lowcost[j];
}
if(v != -1)
{
cost += abs(p[nearvex[v]].z - p[v].z);
len += Edge[nearvex[v]][v];
nearvex[v] = -1;
for(int j = 1; j <= n; j++)
{
double tmp = abs(p[v].z - p[j].z) - Edge[v][j] * l;
if(nearvex[j] != -1 && tmp < lowcost[j])
{
lowcost[j] = tmp;
nearvex[j] = v;
}
}
}
}
return cost / len;
}
int main()
{
while(scanf("%d", &n) != EOF && n)
{
for(int i = 1; i <= n; i++)
scanf("%d%d%d", &p[i].x, &p[i].y, &p[i].z);
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
Edge[i][j] = cal(i, j);
double a = 0, b;
while(1)
{
b = prim(1, a);
if(fabs(a - b) < eps) break;
a = b;
}
printf("%.3f\n", b);
}
return 0;
}
分享到:
相关推荐
NULL 博文链接:https://128kj.iteye.com/blog/1705139
NULL 博文链接:https://200830740306.iteye.com/blog/603493
先利用prim算法求出最小生成树,然后通过往MST里加边来判断新生成的最小生成树是否具有最小的权值,POJ上The Unique MST(1679)题是要求判断最小生成树是否唯一,此题其实根本不用这样做,但是为了练习球次小生成树...
POJ上的一道题目,自己写的代码,因为想下载别人的, 所以就放上了。
最小度限制生成树的求法也不是很难,先把度被限制的那个点(s点)去掉,求一次MST,然后把s到各个连通分量的最小权值的边加上,然后继续加边看看能否使生成树权减小(详见解题报告)。
NULL 博文链接:https://128kj.iteye.com/blog/1704752
POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类
北大POJ初级-动态规划 解题报告+AC代码
放炮问题,北大网站 POJ 1185 算法
北大POJ1015-Jury Compromise【动态规划DP】 解题报告+AC代码
poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题...
poj经典动态规划题目解题报告,包括经典的动态规划题目20多道,可以作为学习动态规划系统的资料,包括题目: Pku acm 1179 Polygon Pku acm 1125 Stockbroker Grapevine Pku acm 1160 post office Pku ...
PKU Online Judge上面很全面的动态规划试题总结。动态规划是ACM考点中最重要的一大类算法之一,对于工作人员来说,动态规划也是实际开发中经常会遇到的算法。这是POJ上面很多DP题目的总结与深刻分析。利于算法学习,...
POJ1321棋盘问题 很好两种解法很值得去参考一下 完整的实验报告还有代码希望kan
POJ题解 个人写法 线段树每个人都不一样
北大poj1012-Joseph【经典约瑟夫问题】 poj1012-Joseph【经典约瑟夫问题】
关于poj dp分类,我一直寻找dp的分类,终于找到了,也上传一下
POJ1048,加强版的约瑟夫问题 难度中等
POJ第1861题源码 POJ第1861题源码 POJ第1861题源码
poj 2763 程序 线段树 LCA 2000MS AC