题目描述
形如:ax3 + bx2 + 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
题解
实数二分:
解题步骤
:
- 先对原方程求导,即 3ax2 + 2bx + c = 0;
- 利用求根公式求出导函数的两个根
x1 x2
,即原方程的两个极值; - 最后三个实根分别在
(-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
😨