三分搜索

二分搜索方程的根大家应该都明白，即在给定的函数是单调函数的情况下，我们可以不断的二分来搜索方程的根

单是对于凸函数来说二分就不行了

什么是凸函数了，大家看图就明白了

当然开口向上也是 

算法比较简单，我们直接讲解模板

<span style="font-size:18px;">#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const double eps = 0.0000001;
double thirdsearch(double y)
{
    double left=0.0;
    double right = 100.0;
    double mid,mmid;
    while(right-left>=eps)
    {
       mid = (left+right)/2;
       mmid = (mid+right)/2;
        if(cal(mid,y)<cal(mmid,y))//cal函数根据需要定义
         right=mmid;
         else
          left = mid;
    }
    return (left+right)/2;
}</span>

只看一下什么时候更新left什么时候更新right

M1=(L+R)/2;

M2=(M+R)/2;

其中f(M1)比f(M2)更接近target

故R = M1;

反之L = M2;

显然三分每次都把区间缩小到原来的3/4

下面推荐一个三分搜索的入门题

http://acm.hdu.edu.cn/showproblem.php?pid=2899