二分查找 适用于单调函数中逼近求解某点的值。
如果遇到凸性或凹形函数时,可以用三分查找求那个凸点或凹点。
下面的方法应该是三分查找的一个变形。
如图所示,已知左右端点L、R,要求找到白点的位置。
思路:通过不断缩小 [L,R] 的范围,无限逼近白点。
做法:先取 [L,R] 的中点 mid,再取 [mid,R] 的中点 mmid,通过比较 f(mid) 与 f(mmid) 的大小来缩小范围。
当最后 L=R-1 时,再比较下这两个点的值,我们就找到了答案。
1、当 f(mid) > f(mmid) 的时候,我们可以断定 mmid 一定在白点的右边。
反证法:假设 mmid 在白点的左边,则 mid 也一定在白点的左边,又由 f(mid) > f(mmid) 可推出 mmid < mid,与已知矛盾,故假设不成立。
所以,此时可以将 R = mmid 来缩小范围。
2、当 f(mid) < f(mmid) 的时候,我们可以断定 mid 一定在白点的左边。
反证法:假设 mid 在白点的右边,则 mmid 也一定在白点的右边,又由 f(mid) < f(mmid) 可推出 mid > mmid,与已知矛盾,故假设不成立。
同理,此时可以将 L = mid 来缩小范围。
先增再减的模型
double three_devide(double low,double up)
{
double m1,m2;
while(up-low>=eps)
{
m1=low+(up-low)/3;
m2=up-(up-low)/3;
if(f(m1)<=f(m2))
low=m1;
else
up=m2;
}
return (m1+m2)/2;
}
先减再增
double three_devide(double low,double up)
{
double m1,m2;
while(up-low>=eps)
{
m1=low+(up-low)/3;
m2=up-(up-low)/3;
if(f(m1)<=f(m2))
up=m2;
else
low=m1;
}
return (m1+m2)/2;
}
P3382 【模板】三分法
题目描述
如题,给出一个N次函数,保证在范围[l,r]内存在一点x,使得[l,x]上单调增,[x,r]上单调减。试求出x的值。
输入输出格式
输入格式:
第一行一次包含一个正整数N和两个实数l、r,含义如题目描述所示。
第二行包含N+1个实数,从高到低依次表示该N次函数各项的系数。
输出格式:
输出为一行,包含一个实数,即为x的值。四舍五入保留5位小数。
输入输出样例
输入样例#1: 复制
3 -0.9981 0.5 1 -3 -3 1
输出样例#1: 复制
-0.41421
说明
时空限制:50ms,128M
数据规模:
对于100%的数据:7<=N<=13
样例说明:
如图所示,红色段即为该函数f(x)=x^3-3x^2-3x+1在区间[-0.9981,0.5]上的图像。
当x=-0.41421时图像位于最高点,故此时函数在[l,x]上单调增,[x,r]上单调减,故x=-0.41421,输出-0.41421。
(Tip.l&r的范围并不是非常大ww不会超过一位数)
事实证明,对于单峰函数,二分还是够的
实际上,就是一个爬坡的过程——
发现前面大,就往前面爬;发现后面大,就往后面爬。爬到尖顶了,往刚才的方向一看:哎呀!要掉下去了!连忙回缩——无法回缩了,则当前即为最高点
#include<bits/stdc++.h>
using namespace std;
const double eps=1e-6;
int n;
double a[20];
double f(double x)
{
double sum=0;
for(int i=n;i>=0;i--){
sum=sum*x+a[i];
}
return sum;
}
int main()
{
scanf("%d",&n);
double l,r;
scanf("%lf %lf",&l,&r);
for(int i=n;i>=0;i--) scanf("%lf",&a[i]);
double m1,m2;
while(r-l>=eps){
double mid=(l+r)/2;
if(f(mid-eps)>f(mid)) r=mid;
else l=mid;
}
printf("%.5f\n",l);
return 0;
}
二分效率为log2,三分效率为2*log3
能用二分还是二分吧,简洁好写还效率高
对比三分
#include<bits/stdc++.h>
using namespace std;
const double eps=1e-6;
int n;
double a[20];
double f(double x)
{
double sum=0;
for(int i=n;i>=0;i--){
sum=sum*x+a[i];
}
return sum;
}
int main()
{
scanf("%d",&n);
double l,r;
scanf("%lf %lf",&l,&r);
for(int i=n;i>=0;i--) scanf("%lf",&a[i]);
double m1,m2;
while(r-l>=eps){
m1=l+(r-l)/3;
m2=r-(r-l)/3;
if(f(m1)<=f(m2)) l=m1;
else r=m2;
}
printf("%.5f\n",l);
return 0;
}