[NOIP2001]一元三次方程求解

输入描述:

一行，4个实数A,B,C,D。

输出描述:

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

1 -5 -4 20

-2.00 2.00 5.00

#include<bits/stdc++.h>
using namespace std;
double a, b, c, d;

double f(double x)
{
    return a*x*x*x + b*x*x + c*x + d;
}
void half(double l, double r)
{
    if (r - l <= 0.001)
    {
        printf("%.2lf ", l);
        return;
    }
    double mid = (l + r) / 2;
    double ans_l, ans_r;
    ans_l = f(l)*f(mid);
    ans_r = f(r)*f(mid);
    if (f(mid) == 0)
        printf("%,2lf ", mid);
    if (f(r) == 0)
        printf("%.2lf ", r);
    if (ans_l < 0)
        half(l, mid);
    else if (ans_r < 0)
        half(mid, r);
}

int main(){
    cin >> a >> b >> c >> d;
    for (double i = -100; i <= 99; ++i)
    {
        if (f(i)*f(i + 1) <= 0)
        {
            half(i, i + 1.0);
        }
    }
    cout << endl;
}