1238:一元三次方程求解
时间限制: 1000 ms 内存限制: 65536 KB
【题目描述】
形如:ax3+bx2+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
1. #include <stdio.h> 2. #include <iostream> 3. #define N 100002 4. using namespace std; 5. double a,b,c,d,x1,x2,midx; 6. double jg(double x) { 7. return a*x*x*x+b*x*x+c*x+d; 8. } 9. int main(int argc, char *argv[]) 10. { 11. scanf("%lf %lf %lf %lf",&a,&b,&c,&d); 12. for(int i=-100;i<=100;i++){ 13. x1=i; 14. x2=x1+1; 15. if(jg(x1)==0) printf("%.2lf ",x1); 16. else if(jg(x1)*jg(x2)<0){ 17. while(x2-x1>=0.001){ 18. midx=(x1+x2)/2; 19. if(jg(x1)*jg(midx)<=0)x2=midx; 20. else x1=midx; 21. } 22. printf("%.2lf ",x1); 23. } 24. } 25. return 0; 26. }