一元三次方程求解

原题链接：

https://www.luogu.com.cn/problem/P1024

题目描述

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

提示

记方程f(x) = 0, 若存在2个数x₁和x₂ ，且x₁ < x₂ ,f(x₁) × f(x₂) < 0,则在（x₁ , x₂）之间一定有一个根。

输入格式

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

输出格式

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

输入输出样例：

输入#1

1 -5 -4 20

输出#1

-2.00 2.00 5.00

法一：暴力出奇迹！

思路方法：

根据提示和题目中对精度的要求，我们可以固定x₁和x₂的间隔为0.001，然后遍历区间[-100, 100], 将满足f(x₁) × f(x₂) < 0条件的x₁和x₂取平均值输出即可。

AC代码：

#include<bits/stdc++.h>
using namespace std;
int main(){
	double a, b, c, d;
	cin>>a>>b>>c>>d;
	for(double i = -2; i <= 100; i += 0.001){
		double j = i + 0.001;
		double l = a * i * i * i + b * i * i + c * i + d;
		double r = a * j * j * j + b * j * j + c * j + d;
		if(l * r <= 0){
			printf("%.2lf ", (i + j) / 2);
		}
	}
	return 0;
} 

法二：二分查找

思路方法：

总体思路跟上面的暴力枚举差不多，但二分查找做了一些优化。这里我们可以充分利用题目给定的条件：根于根之差的绝对值>=1;根据这一点，我们可以扩大遍历区间[-100, 100]时的增量，找到满足条件的区间后，再根据二分查找法将查找的精度限定在题目要求的范围之内。

AC代码：

#include<bits/stdc++.h>
using namespace std;

double a, b, c, d;
double func(double x){
	return a * x * x * x + b * x * x + c * x + d;
}
int main(){
	cin>>a>>b>>c>>d;
	int count = 0;
	double l, r, x1, x2;
	for(int i = -100; i <= 100; i++){
		l = i;
		r = i + 1;
		x1 = func(l);
		x2 = func(r);
		if(!x1){//左端点是零点直接输出结果 
			printf("%.2lf ",l);
			count++;
		}
		if(x1 * x2 < 0){
			double mid;
			while(r - l >= 0.001){//二分查找，不断缩小区间 
				mid = (l + r) / 2;
				if(func(mid) * func(l) <= 0){
					r = mid;
				}else{
					l = mid;
				}
			}
			printf("%.2lf ",mid);
			count++;
		}
		if(count == 3){//若找到三个零点提前结束循环 
			break;
		}
	}
	return 0;
}