hdu 2899 Strange fuction

Strange fuction 函数求导看凸性用三分法
添加链接描述

//三分法
#include <bits/stdc++.h>
using namespace std;
//F(x) = 6*x^7+8*x^6+7*x^3+5*x^2-y*x (0 <= x <=100)
#define F(x, y) (x*(x*(x*(x*(x*(x*(6*x+8)))+7)+5)-y))
double solve(double l, double r, double y){
    double lm, rm;
    while(r - l > 1e-7){
        lm = (2*l+r)/3;
        rm = (l+2*r)/3;
        if(F(lm, y)>F(rm, y))l = lm;
        else if(F(lm, y)<F(rm, y)) r = rm;
        else l = lm, r = rm;
    }
    return F((l+r)/2, y);
}
int main(){
    int T;
    double y;
    scanf("%d", &T);
    while(T--){
        scanf("%lf", &y);        
        printf("%.4f\n", solve(0, 100, y));
    }
    return 0;
}

给一函数，该函数在任意Y>0的情况下求x在[0,100]内有极小值。
该函数求导，导函数（1.x=0时，F(x)<0；1,当x值很大，F(x)>0）则该函数为凸性。三分法
这种情况先减后增有极小，若lm比rm低（即lm对应的函数值 < rm函数值）则极小点（最低点）肯定在[ left, rm ] ，反之在[ lm, right ]，剩下就跟二分一样根据大小关系调整区间就行了。
那lm和rm取值多少？一个不错的取值是lm为整个区间的1/3点，rm为2/3点，即
**lmid = l + (r - l)/3; rmid = r - (r - l)/3; **
lm低于rm，则极大在[ lm,right ].否则极大在 [ left, rm ]。写代码上就是极小的处理语句反过来就行了。