[模板]三分法

模板第二弹

题目选自：Luogu P3382
显然三分。
标题就是三分。
先给个代码。
每次都给那么长的代码，好像有点不太好，所以这次只给主代码（想要前面的很有意义莫名其妙的代码的可以翻我以前的博客）。

double l,r,g[14],_1,_2;
int n;
#define eps 1e-7
double f(double x)
{
    double y=g[0];
    fr(i,1,n)
        y+=g[i]*pow(x,i);
    rt y;
}
int main(){
#ifndef ONLINE_JUDGE
    freopen("","r",stdin);
    freopen("","w",stdout);
#endif
    n=read();
    scanf("%lf%lf",&l,&r);
    fd(i,n,0)
        scanf("%lf",&g[i]);
    while(l+eps<r)
    {
        _1=(l+l+r)/3;
        _2=(r+_1)/2;
        if(f(_1)<=f(_2))
            l=_1;
        else
            r=_2;
    }
    printf("%.5lf",l);
    rt 0;
}

可以参考一下:

注：
#1:
$mid1= \dfrac{2l+r}{3}$
$mid2= \dfrac{l+2r}{3}$
#2
$f(mid1)<=f(mid2)$
#3
$l+eps<r$
#4
$l=mid1$
#5
$r=mid2$
eps
$eps=10^{-7}$
解释：

#1将$mid1,mid2$设置为$l,r$的两个三分点

#2判断两点函数值大小

#3判断$l,r$值是否不足够接近（区间是否足够大）

#4舍去$[l,mid1)$

#5舍去$(mid2,r]$

eps可以看做处理精度的标准

思路大概清晰了一些吧。
来梳理一下：
1. 该搞什么初始化就搞什么初始化
2. 重复执行3~5，直到精度足够
3. 算出等分点
4. 分别求出亮等分点的函数值
5. 比较，并做相应处理（看是舍小还是舍大）
6. 得出结果
以上就是本博客的精华部分，下期博客，再见。

最后，祝大家：

rp=(unsigned long long)-1