二分求方程的解

浮点二分

浮点二分的一般写法；

//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，其实这里浮点二分是没什么的，但是这里你看不同的求方程的解，它在处理上是不同的；