[P1034][NOIP2001]一元三次方程求解 (二分)

二分

#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
double fc(double x)
{
    return a*pow(x,)+b*pow(x,)+c*x+d;
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    double l=-,r=,ans,x1,x2,mid;
    int s=;
    for(int i=-;i<=;i++){
        l=i;
        r=i+;
        x1=fc(l);
        x2=fc(r);
        if(!x1){
            printf("%.2lf ",l);
            s++;
        }
        if(x1*x2<){
            while(r-l>=0.001)
            {
                mid=(l+r)/;
                if(fc(mid)*fc(r)<=) l=mid;
                else r=mid;
            }printf("%.2f ",r);
        s++;
        }
        if(s==) break;
    }
    return ;
}