[模板]三分法

模板第二弹

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

double l,r,g[14],_1,_2;int n;#define eps 1e-7double 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=2l+r3
mid2=l+2r3
#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