这篇博客是个极其简单的东西。
如题,给出一个N次函数,保证在范围[l,r]内存在一点x,使得[l,x]上单调增,[x,r]上单调减。试求出x的值。
输入格式
第一行一次包含一个正整数N和两个实数l、r,含义如题目描述所示。
第二行包含N+1个实数,从高到低依次表示该N次函数各项的系数。
输出格式
输出为一行,包含一个实数,即为x的值。四舍五入保留5位小数。
输入
3 -0.9981 0.5
1 -3 -3 1
输出
-0.41421
三分主要就是用于求函数极值,当然只用于求凸函数或者凹函数的极值。
尽管求这个东西可以通过对函数求导,再二分。。。嘤嘤嘤,但是对于有些麻烦的函数,求导是很不方便的,而且也没有复杂度的优势。
对于上面这个凸函数,我们每次取三分之一点 m1 和 三分之二点 m2 然后我们通过比较两点之间的函数值然后让 l ,r 逼近极值即可。
double f_max(){
while(fabs(l-r)>1e-7){
double mid=(l+r)/2;
if(f(mid-eps)<f(mid+eps)) l=mid;
else r=mid;
}
return l;
}
为什么上面的代码没有取三分之一和三分之二点呢?其实三分之二和三分之一只是理论点,但是在实际情况中,我们一般取中点附件更快哦!!!
AC代码:
#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define int long long
using namespace std;
const double eps=1e-8;
int n;
double l,r,a[20];
double f(double x){
double res=0;
for(int i=n;i>=0;i--) res=res*x+a[i];
return res;
}
double f_max(){
while(fabs(l-r)>1e-7){
double mid=(l+r)/2;
if(f(mid-eps)<f(mid+eps)) l=mid;
else r=mid;
}
return l;
}
signed main(){
cin>>n>>l>>r;
for(int i=n;i>=0;i--) cin>>a[i];
printf("%.5lf\n",f_max()+0.0000005);
return 0;
}