浮点二分
浮点二分的一般写法;
//check()函数
bool check()
{
}
double tree(double l, double r)
{
double eps = 1e-6; // eps 表示精度,取决于题目对精度的要求
while (r - l > eps)
{
double mid = (l + r) / 2;
if (check(mid)) r = mid;//check()是一个函数,用于检查是否满足某个性质的;
else l = mid;
}
return l;
}
其实这里没有必要去记忆什么,我们就只要知道浮点二分是没有什么端点+1,-1
的讨论的,因为这里无论怎样,l和r都是更新成mid,不会出现死循环的,然后剩下跟着题意写就行了;
题目一
分析
这题它讲了根在[1.5,2.4]之间并且只有一根,然后由f(1.5)>0
,f(2.4)<0
,可以判断在1.5~2.4
之间是单调递减的;
然后就根据浮点二分,来写求解即可;
#include<bits/stdc++.h>
using namespace std;
//pow(a,b),是求a的b次方
double BinarySearch(double l,double r)
{
double mid;
double f;
while( r-l>1e-11 )
{
mid = (l+r)/2.0;
f = pow(mid,5)-pow(mid,4)*15.0+pow(mid,3)*85.0-pow(mid,2)*225.0+mid*274.0-121.0;
if(f > 0)
{
l = mid;//因为是单调递减的,而此时>0,说明根应在右边找
}
else
{
r = mid;
}
}
return mid;
}
int main()
{
printf("%.6f\n",BinarySearch(1.5,2.4));
return 0;
}
题目2
链接 方程求解
分析
1,这题也是要确定函数根的区间,判断单调性;
2,但是这没给范围,而且也没讲有几个根,那就感觉有点恶心了,那就只能自己猜了;
3,我们尝试求个导看看是不是单调的,通过求导发现,这是一个单调递增的函数,然后,我们再找一下根的区间,这也是靠猜了,应该也不难猜(毕竟题目出的已经够可恶的了);
4,然后可以发现f(0)<0,f(10)>0
,说明在[0,10]必有一根;
#include<bits/stdc++.h>
using namespace std;
double tree(double l,double r)
{
double mid;
while(r-l>1e-11)
{
mid=(l+r)/2;
double f=pow(mid,3)-5*pow(mid,2)+10*mid-80;
if(f<0) l=mid;
else r=mid;
}
return mid;
}
int main()
{
printf("%.9lf\n",tree(0,10));
return 0;
}
题目3
链接 一元三次方程求解
分析
1,这题,与上面相比告诉你有3个根了,并且a,b,c,d还是不确定,那么你求单调性,就不好搞了,因为每次都不一样;
2,由于答案在[-100,100]之间,并且根与根之差的绝对值>=1,那么我们唯一能做的是在[-100,-99],[-99,-98]…[99,100],里面逐个去找是否有根;
3,如何判断该区间有根,我们需要用到高中的零点存在性定理来判断,即f(a)*f(b)<0
那么a,b区间是一定有根的;
4,要注意特判一下端点,因为端点也是有可能的;
#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
double f(double mid)
{
return a*pow(mid,3) + b*pow(mid,2) + c*mid +d;
}
int main()
{
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
for(int i=-100;i<=99;i++)
{
double x1=i*1.0,x2=(i+1)*1.0;
if(f(x1)==0)
printf("%.2lf ",x1);//检测端点
//x1如果不是的话,那么零点可能在(x1,x2],这里不必担心x2,因为下一次循环还是要判断端点的,而此时这个端点就变成x2了呀;
else if(f(x1)*f(x2)<0)
{
while(x2-x1>1e-6)
{
double mid = (x1+x2)/2;
if(f(x1)*f(mid)<0) x2=mid;
else x1=mid;
}
printf("%.2lf ",x1);
}
}
return 0;
}
总结:
1,浮点二分是相对简单的,不会牵涉边界的+1,-1
的问题,所以一般想到怎么写了,那么基本就可以了;
2,其实这里浮点二分是没什么的,但是这里你看不同的求方程的解,它在处理上是不同的;