- 链接 : 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(){
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即可。