洛谷P3382 【模板】三分法 三分法求 函数最大值

洛谷P3382 【模板】三分法
三分法求 函数最大值
三分 适用于求凸峰函数 或者凹峰函数
但听dalao说不大适合 /
                        ---    
                      /          
这种有平的情况

还有我觉得 f(mid1) == f(mid2 ) 也可特判一下 此时 l = mid1 r = mid2

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
#include <iomanip>
using namespace std ; 
const int maxn = 14 ;
const double eps = 1e-7 ; 
int n ; 
double k[maxn] ; 
double l,r,mid1,mid2 ; 

inline double f(double x ) 
{
    double ans = 0.0 ; 
    double y = 1 ; 
    for(int i=n+1;i>=1;i--) 
    {
        ans = ans + k[ i ] * y ; 
        y*=x ; 
    }
    return ans ; 
} 

int main() 
{
    scanf("%d%lf%lf",&n,&l,&r) ; 
    for(int i=1;i<=n+1;i++) scanf("%lf",&k[ i ]) ; 
    while(abs(r-l)>eps) 
    {
        mid1 = l+(r-l) / 3 ;
        mid2 = r-(r-l) / 3 ; 
        if( f(mid1) < f(mid2) ) 
            l = mid1 ;    
        else 
            r = mid2 ; 
    }
    printf("%.5lf
",l) ; 
    
    return 0 ; 
}