HDU 2899 Strange fuction——三分搜索

• 链接 ： http://acm.hdu.edu.cn/showproblem.php?pid=2899
• 题意 ： 给定一个函数，判断其最小值。
• 分析：
  首先回顾一下，单调函数查找特定的值，是通过二分搜索得到的。
  即 l = min, r = max. mid = (l + r) /2 此处是下界。 mid = l + (l+r)/2 是上界。
  三分搜索，则针对的是凹凸函数。
  如图：
  
  其中 mid = l + (r - l)/3
  midmid = r - (r - l)/3。
  如果mid更靠近最大值，则舍弃区间(midmid, r)；反之舍弃(l, mid)
  算法的正确性 ： （1）若是mid 和midmid在最大值一边，因为凹凸函数在最大值一边是单调的，所以舍弃另一个远离最大值的区间则不会有任何影响。 （2）若是mid和midmid在两边，由于最大值在中间的区间，所以舍弃两边的区间不会有任何影响。
• 代码

#include "bits/stdc++.h"
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define mk make_pair
typedef long long ll;
const int maxn = 1e5 + 7;
int t,x,y;
double fct(double x){         //计算函数值
    return 6*pow(x,7) + 8*pow(x,6) + 7*pow(x,3) + 5*pow(x,2) - x*y;
}
int main(){
    // freopen("1.txt","r",stdin);
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin>>t;
    while (t--) {
        cin>>y;
        double l = 0, r = 100;
        while (r - l > 1e-6) {
            double midl = l + (r - l)/3;
            double midr = r - (r - l)/3;
            if (fct(midl) > fct(midr)) {
                l = midl;
            }else{
                r = midr;
            }
        }
        printf("%.4lf\n",fct(r));
    }
    return 0;
}

• 遇到的问题：
  跳出循环的条件，因为精确到1e-4，所以一开始我设置的是 r - l > 1e-5，不过wa了，设置成 1e-6即可。