三分法（笔记）

三分法

mid=(Left+Right)/2
midmid=(mid+Right)/2
如果mid接近极值点Right=midmid否则Left=mid

例题

三角形相似 H/(D+y)=x/y
y=x*D/(H-x) r
三角形相似 h/(L-x+y)=x/y
带入y得 L=(h-x)*D/(H-x)+x
易知x>=0&&x<=h

#include<iostream>
#include<cstdio>
using namespace std;
double H,h,D;
double calce(double x)
{
    //cout<<(h-x)*D/(H-x)<<endl;
    return x+(((h-x)*D)/(H-x));
}
int main()
{
    int num;
    scanf("%d",&num);
    while(num--)
    {
        scanf("%lf%lf%lf",&H,&h,&D);
        double l=0,r=h;
        while(l+0.00001<r)//使精确到小数点后三位
        {
            double mid=(l+r)/2.0;
            double midmid=(mid+r)/2.0;
            if(calce(mid)>calce(midmid))
                r=midmid;
            else l=mid;
        }
        printf("%.3f\n",calce(r));
    }
    return 0;
}

二分法实在单调递增或递减数据基础上不断缩小范围（通过判断大于小于查询值）选取合适的数据，三分法是在一段有单个极值的数据上，缩小范围（通过极值的性质取较小或较大值），取到极值。