codevs 1038 一元三次方程求解 NOIP2001提高组

题目链接：http://codevs.cn/problem/1038/

题解：

　　嗯，exm？才知道二分隶属搜索专题……

　　对-100到100枚举，按照题目中的提示，当当f_i*f_i+1<0时，二分深搜，直到精度达到小数点后4位为止（保守起见），当f_i*f_i+1=0时，判定i是否为解，否则跳过（如果i和i+1都为解，下一次循环时还会再判断i+1）

#include<cstdio>
#include<cmath>
double a,b,c,d;
double f(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
void search(double l,double r)
{
    if(r-l<0.001)
    {
        printf("%.2lf ",l);
        return;
    }
    if(f(l)*f(r)<0)
    {
        double mid=(l+r)/2;
        if(f(l)*f(mid)==0)
        {
            printf("%.2lf ",l);
            return;
        }
        else if(f(l)*f(mid)<0)
        {
            search(l,mid);
            }
        else if(f(mid)*f(r)<0)
        {
            search(mid,r);
        }
    }
}
int main()
{
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    for(int i=-100;i<=100;i++)
    {
        if(f(i)*f(i+1)<0)
        {
            search(i,i+1);
        }
        else if(f(i)*f(i+1)==0)
        {
            if(f(i)==0)printf("%.2lf ",(double)i);
        }
    }
    return 0;
}