码迷,mamicode.com
首页 > 其他好文 > 详细

[日常摸鱼]三分法

时间:2018-01-24 15:30:41      阅读:139      评论:0      收藏:0      [点我收藏+]

标签:bsp   一个   res   clu   三分法   --   log   can   int   

翻到一个三分法的模板发现没有写掉…今天干脆写掉算了…(luogu3382)

跟二分基本差不多…

#include<cstdio>
typedef double dl;
const double eps=1e-7;
const int N=20;
dl l,r,a[N];
int n;
inline dl calc(dl x)
{
    dl res=0,tx=1;
    for(register int i=0;i<=n;i++)
    {
        res+=tx*a[i];tx*=x;
    }
    return res;
}
int main()
{
    scanf("%d%lf%lf",&n,&l,&r);
    for(register int i=n;i>=0;i--)scanf("%lf",&a[i]);
    while(r-l>eps)
    {
        dl x1=l+(r-l)/3.0,x2=l+2.0*(r-l)/3.0;
        dl t1=calc(x1),t2=calc(x2);
        if(t1>t2)r=x2;
        else l=x1;
    }
    printf("%.5lf",l);
    return 0;
}

 

[日常摸鱼]三分法

标签:bsp   一个   res   clu   三分法   --   log   can   int   

原文地址:https://www.cnblogs.com/yoooshinow/p/8341365.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!