洛谷P1024 [NOIP2001 提高组] 一元三次方程求解

这是一道二分题

废话不多说上代码

#include<iostream>
using namespace std;
double a,b,c,d;
int cnt;
double fc(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    
    for(int i=-100;i<100;i++)
    {
        double l=i,r=l+1;
            double x1=fc(l);
            double x2=fc(r);
        if(!fc(l))
        {         
            printf("%.2lf ",l);
            cnt++;
        }
        if(x1*x2<0)
        {
            while(r-l>=0.001)
            {
                double mid=(l+r)/2;
                if(fc(mid)*fc(r)<=0)
                {
                    l=mid;
                }
                else
                r=mid;
            }
            printf("%.2lf ",r);
            cnt++; 
        }
        if(cnt==3)
        return 0;
    }
    return 0;
}