一元三次方程求解——二分

题目描述

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

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

输入格式

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

输出格式

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

样例 1

样例输入 1

1 -5 -4 20

样例输出 1

-2.00 2.00 5.00

高中数学的零点存在性定理：如果函数y= f(x)在区间[a,b]上的图象是连续不断的一条曲线，并且有f(a)·f(b)<0,那么，函数y= f(x)在区间(a,b)内有零点，即至少存在一个c∈(a,b),使得f(c)=0,这个c也就是方程f(x)= 0的根。

那么此题我们就可以枚举每个长度为1的区间，来判断给区间是否存在零点

高中的二分法：对于在区间{a,b}上连续不断，且满足f(a)f(b)<0的函数y=f(x)，通过不断地把函数f(x)的零点所在的区间二等分，使区间的两个端点逐步逼近零点，进而得到零点近似值的方法叫做二分法。

我们就可以用二分来寻找满足f(a)*f(b)<0条件的区间的零点了。

#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;
double solve(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    cin>>a>>b>>c>>d;
    for (double i=-100.0;i<=100.0;i++)
    {
        double x1=i,x2=i+1.0;
        if (solve(x1)==0.0) printf("%.2lf ",i);
        if (solve(x1)*solve(x2)<0)
        {
            while (x2-x1>=0.001)     //此题要求结果保留2位小数，取0.001就行
            {
                double x=(x1+x2) / 2.0;
                if (solve(x1)*solve(x)<=0.0) x2=x;
                else x1=x;
            }
            printf("%.2lf ",x1);
        }
    }
    return 0;
}