总时间限制: 

1000ms

内存限制: 

65536kB

描述

有形如：ax³+bx²+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

来源

NOIP2001复赛 提高组 第一题

1 #include
       
         2 using namespace std; 3 const double eps=0.0001; 4 double a,b,c,d; 5 inline double calc(double x){ 6     double h=a*x*x*x+b*x*x+c*x+d; 7     return h; 8 } 9 double ans[5];10 inline void find(double,double);11 int sum;12 int main(){13     scanf("%lf%lf%lf%lf",&a,&b,&c,&d);14     for(int i=-101;i<=101;i++){15         double tmp1=calc((double)i);16         double tmp2=calc(double(i+1));17         if(tmp1==0){18             ans[++sum]=(double)i;19         } 20         else if(tmp2==0){21             ans[++sum]=(double)(i+1);22             i++;23         } 24         else if(tmp1*tmp2<0){25             find((double)i,(double)i+1);26         }27         if(sum==3){28             printf("%.2lf %.2lf %.2lf",ans[1],ans[2],ans[3]);29             return 0;30         }31     }32     33     return 0;34 }35 inline void find(double l,double r){36     if(r-l<=eps){37         ans[++sum]=(l+r)/(double)2;38         return ;39     }40     double mid=(l+r)/(double)2;41     double tmp=calc(mid);42     double tmpl=calc(l); 43     double tmpr=calc(r);44     if(tmp==0){45         ans[++sum]=tmp;46         return ;47     }48     else if(tmp*tmpl<0){49         find(l,mid);50     }51     else if(tmp*tmpr<0){52         find(mid,r);53     }54 }