一元三次方程求解

一元三次方程求解点击打开链接

总时间限制:  

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

提示

记方程f(x)=0，若存在2个数x1和x2，且x1 < x2，f(x1)*f(x2) < 0，则在(x1，x2)之间一定有一个根。 源码

#include<cstdio>  
#include<iostream>   
#include<algorithm>  
#include<cstring>
using namespace std;
double a,b,c,d;
double f(double x)
{
	return x*x*x*a+b*x*x+x*c+d;
}//算出当前值
int main()
{  
	scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
	for(double x=-100;x<=100;x++)
	{
		double x1=x,x2=x+1;
		if(f(x1)==0)//处理x1就为最小值的情况
			printf("%.2lf ",x1);
		else
			if(f(x1)*f(x2)<0)
			{
				while(x2-x1>=1e-3)
				{
					double xx=(x2+x1)/2;
					if((f(x1)*f(xx))<=0)
						x2=xx;
					else
						x1=xx;
				}//分治处理
				printf("%.2lf ",x1);
			}
	}
	printf("\n");
	return 0;  
}  

谢谢观赏！！！