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即可。

猜你喜欢

转载自blog.csdn.net/qq_39763472/article/details/105061044