luogu1024 一元三次方程求解（NOIP2001提高组第1题）

luogu1024  一元三次方程求解（NOIP2001提高组第1题）

时空限制    1000ms/128MB

题目描述

有形如：ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a，b，c，d 均为实数)，并约定该方程存在三个不同实根(根的范围在-100至100之间)，且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。

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

输入输出格式

输入格式：

一行，4个实数A，B，C，D。

输出格式：

一行，三个实根，并精确到小数点后2位。

输入输出样例

输入样例#1：

1 -5 -4 20

输出样例#1：

-2.00 2.00 5.00

代码

#include<cstdio>
using namespace std;
double a,b,c,d;

double f(double x){
	return a*x*x*x+b*x*x+c*x+d;
}

int main(){
	scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
	int cnt=0;
	for (double x=-100,x1,x2,xx; x<=100; x+=0.9999){	//枚举可能有答案的区间
		x1=x; x2=x+0.9999;
		if (f(x1)==0)
			++cnt==3?printf("%.2lf\n",x1):printf("%.2lf ",x1);
		else if (f(x1)*f(x2)<=0){
			while (x2-x1>=0.001){
				xx=(x1+x2)/2;
				if (f(x1)*f(xx)<=0) x2=xx;
				else x1=xx;
			}
			++cnt==3?printf("%.2lf\n",x1):printf("%.2lf ",x1);
		}
	}
	return 0;
}