【C++算法讲解】实数二分

上期我们讲了如何用整数二分解决问题，那么今天我们就来学习实数二分。

其实实数二分和整数二分没有什么不同，只是要注意一下浮点数的问题。

那么废话不多说，我们来看例题。

题目：P1024 [NOIP2001 提高组] 一元三次方程求解

题目描述

有形如：ax^3+bx^2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数（a,b,c,d 均为实数），并约定该方程存在三个不同实根（根的范围在 −100 至 100 之间），且根与根之差的绝对值 ≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 2 位。

提示：记方程 f(x)=0，若存在 22 个数 x1​ 和 x2​，且 x1​<x2 ， f(x1​)×f(x2​)<0，则在 (x1​,x2​) 之间一定有一个根。

输入格式

一行，4 个实数 a,b,c,d。

输出格式

一行，3 个实根，从小到大输出，并精确到小数点后 2 位。

输入输出样例

输入 #1

1 -5 -4 20

输出 #1

-2.00 2.00 5.00

这道题如果用枚举的话，则需要把小数点后两位也枚举出来。这时候，我们就可以用二分来。

不过这个时候，我们就要建一个常量，来控制精度。把它赋值为1e-4（也就是0.0001）。

题目上的提示信息也写了：记方程 f(x)=0，若存在 22 个数 x1​ 和 x2​，且 x1​<x2 ， f(x1​)×f(x2​)<0，则在 (x1​,x2​) 之间一定有一个根。

#include<bits/stdc++.h>
using namespace std;
#define ok 1e-4
double l,r,mid;
double A,B,C,D;
double f(double x){                        //一元三次方程
	return (A*x*x*x+B*x*x+C*x+D);
}
int main(){
	cin>>A>>B>>C>>D;
	for(int i=-100;i<=100;i++){            //二分
		l=i,r=i+1,mid;
		if(fabs(f(l))<ok){
			printf("%.2f ",l);
		}else if(fabs(f(r))<ok){           //控制精度
			continue;
		}else if(f(l)*f(r)<0){
			while(r-l>ok){
				mid=(l+r)/2;
				if(f(mid)*f(r)<0) l=mid;
				else r=mid;
			}
			printf("%.2f ",l);
		}
	}
	return 0;
}

总结来说，实数二分只需要控制浮点数的精度就好了，所以也非常的简单。