一元三次方程求解

题目描述
形如：ax³ + bx² + cx + d = 0 这样的一个一元三次方程。

• 给出该方程中各项的系数 (a，b，c，d 均为实数)，
• 并约定该方程存在三个不同实根 (根的范围在 −100 至 100 之间)，
• 且根与根之差的绝对值 ≥ 1。

要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 2 位。

输入格式
一行，包含四个实数 a，b，c，d，相邻两个数之间用单个空格隔开。

输出格式
一行，包含三个实数，为该方程的三个实根，按从小到大顺序排列，相邻两个数之间用单个空格隔开，精确到小数点后 2 位。

输入样例
1.0 -5.0 -4.0 20.0

输出样例
-2.00 2.00 5.00

---

题解
实数二分：

解题步骤：

1. 先对原方程求导，即 3ax² + 2bx + c = 0;
2. 利用求根公式求出导函数的两个根 x1 x2，即原方程的两个极值；
3. 最后三个实根分别在 (-100, x1) (x1, x2) (x2, 100) 之间；

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

double a, b, c, d;
	
double check(double x)
{
	return a*x*x*x + b*x*x + c*x + d;
}

double cut(double l, double r)
{
	bool flag = check(l) > 0? true : false;
	
	if(flag) 
	{
		while(r - l > 1e-3)
		{			
			double mid = (l + r) / 2;
			if(check(mid) > 0) l = mid;
			else r = mid;
		}
	}
	else
	{
		while(r - l > 1e-3)
		{
			double mid = (l + r) / 2;
			if(check(mid) > 0) r = mid;
			else l = mid;
		}
	}
	
	return l;
}

int main()
{
	cin >> a >> b >> c >> d;
	
	double x1 = (-2*b + sqrt(4*b*b - 12*a*c)) / 6*a;
	double x2 = (-2*b - sqrt(4*b*b - 12*a*c)) / 6*a;
	
	if(x1 > x2) swap(x1, x2);
	 
	printf("%.2f %.2f %.2f", cut(-100, x1), cut(x1, x2), cut(x2, 100));
	return 0;
}

ps：今天我才知道 double 输出时不能用 %lf 😨