一元三次方程求解（二分法）

链接：https://ac.nowcoder.com/acm/problem/16694
来源：牛客网
 

题目描述

有形如：ax3+bx2+cx+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 -5 -4 20

1 -5 -4 20

输出

复制-2.00 2.00 5.00

-2.00 2.00 5.00

具体可以参考二分法的基本思想，根据题目可知根的范围是-100~100，我们可以以一个整数区间查找根的近似值，如果根刚好是整数，可以直接输出，如果是小数则需要步步逼近。还要注意，防止重复计算。

#include "bits/stdc++.h"
typedef long long int ll;
using namespace std;
double a,b,c,d;
double def(double x)
{
    double y=((a*x+b)*x+c)*x+d; //秦九韶算法
    return y;
}
int main()
{

    cin>>a>>b>>c>>d;
    for(double i=-100;i<=100;i++)
    {
        double f1=def(i);
        double f2=def(i+1);
        if(f1*f2<=0){
            if(f1==0)
                printf("%.2f ",i);
            else if(f2==0)
            {
                 printf("%.2f ",i+1);
                 i++;   //防止重复判断
            }
            else
            {
                double l=i,r=i+1;
                double mid;
                //在整数区间内求近似值，精度在0.001即可
                while(r-l>1e-3){
                    mid=l+(r-l)/2.0;

                    //二分法的技巧
                    if(def(l)*def(mid)<=0)
                        r=mid;
                    else
                        l=mid;
                }
                 printf("%.2f ",mid);
            }
        }
    }
   return 0;
}