一元三次方程求解

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

输入格式
一行，4 个实数 A,B,C,D。

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

Sample Input
1 -5 -4 20
Sample Output
-2.00 2.00 5.00

---

对方称求导后算出导数等于零的点即为驻点，驻点两边单调性不同，对每个单调区间进行求解即可。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>

using namespace std ;

const int N = 1e5 + 100 ;
double a , b , c , d ;
double f(double x)
{
    return a*x*x*x+b*x*x+c*x+d ;
}
double g(double m , double n)
{
    while(1)
    {
        double mid = (m+n)/2 ;
        if(fabs(m - n) < 0.00001)
            return mid ;
        else
        {
            if(f(mid) == 0)
                return mid ;
            else if(f(m)*f(mid) < 0)
                n = mid ;
            else
                m = mid ;
        }
    }
}
int main()
{
    double x1 , x2 , x3 ;
    double  x , y;
    scanf("%lf%lf%lf%lf" ,&a ,&b ,&c ,&d) ;
    //(-2b+-sqrt(4b*b-4*3a*c))/6a ;
    double m = (-2*b+sqrt(4*b*b-12*a*c))/(6*a) ;
    double n = (-2*b-sqrt(4*b*b-12*a*c))/(6*a) ;
    if(m < n)
        swap( m , n ) ;//m为大
    //m , n 为驻点
    //cout << m << ' ' << n << endl;
    x1 = g(-100 , n) ;
    x2 = g(n , m) ;
    x3 = g(m , 100) ; //每个单调区间分别求解
    printf("%.2f %.2f %.2f",x1 ,x2 ,x3) ;
    return 0 ;
}