三分查找

一 问题

在一个平面上有N个点，每个点的坐标已经给出，现在要求在X轴上找一个点，使得这个点到所有点中最大的距离最小。

二 分析

不知道是否有基于策略选择的算法，也许是DP也许是Greedy。本文通过数值计算来解决此问题。

问题抽象：

设f(x) = max(node1_distance, node2_distance, ……nodeN_distance), a<=x<=b;

求x*使min(f(x)) = f(x*),a<=x*<=b;

如果存在x*使目标函数极小，那么对任意的a<=x1<x2<=b,当x2<=x*时，f(x1)>f(x2);当x1>=x*时，f(x1)<f(x2). 也就是f(x)是一个单峰函数。

问题就转化为求单峰函数的极值。

三 算法

网上主要有两种搜索算法，虽然都称之为三分搜索(或称查找)。核心思想都是通过“去坏留好”的原则不断压缩解空间来求极值。

流传算法1：

void Solve(void){    double Left, Right;    double mid, midmid;    double mid_value, midmid_value;    Left = MIN; Right = MAX;    while (Left + EPS < Right)    {        mid = (Left + Right) / 2;        midmid = (mid + Right) / 2;        mid_area = Calc(mid);        midmid_area = Calc(midmid);        // 假设求解最小极值.        if (mid_area >= midmid_area) Left = mid;        else Right = midmid;    }}

流传算法2：

double ternary_search(double l,double r){ //三分查找求最小值while(l + EPS < r){double left = (2*l + r)/3;double right = (2*r + l)/3;double tmp1 = cal(left);double tmp2 = cal(right);if(tmp1 > tmp2) l = left;else     r = right;}return l;}