NOI2.4 一元三次方程求解

7891:一元三次方程求解
查看提交统计提问
总时间限制: 1000ms 内存限制: 65536kB
描述
有形如：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
提示
记方程f(x)=0，若存在2个数x1和x2，且x1 < x2，f(x1)*f(x2) < 0，则在(x1，x2)之间一定有一个根。
来源
NOIP2001复赛 提高组 第一题

#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
int main()
{
    double a,b,c,d;
    cin>>a>>b>>c>>d;
    if(a<0)
    {
        a=-a;
        b=-b;
        c=-c;
        d=-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);
    double mid=(-100+x1)/2;
    double low=-100,high=x1;
    while(abs(a*mid*mid*mid+b*mid*mid+c*mid+d)>=0.001)
    {
        if((a*mid*mid*mid+b*mid*mid+c*mid+d)<0)
        {
            low=mid;
            mid=(low+high)/2;
        }
        else
        {
            high=mid;
            mid=(low+high)/2;
        }

    }
    cout<<fixed<<setprecision(2)<<mid<<" ";

    mid=(x1+x2)/2;
    low=x1,high=x2;
    while(abs(a*mid*mid*mid+b*mid*mid+c*mid+d)>=0.001)
    {
        if((a*mid*mid*mid+b*mid*mid+c*mid+d)>0)
        {
            low=mid;
            mid=(low+high)/2;
        }
        else
        {
            high=mid;
            mid=(low+high)/2;
        }

    }
    cout<<fixed<<setprecision(2)<<mid<<" ";

    mid=(x2+100)/2;
    low=x2,high=100;
    while(abs(a*mid*mid*mid+b*mid*mid+c*mid+d)>=0.001)
    {
        if((a*mid*mid*mid+b*mid*mid+c*mid+d)<0)
        {
            low=mid;
            mid=(low+high)/2;
        }
        else
        {
            high=mid;
            mid=(low+high)/2;
        }

    }
    cout<<fixed<<setprecision(2)<<mid;
}