一元三次方程求解

一元三次方程求解  

时间限制：1.0s   内存限制：256.0MB

      

问题描述

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

输入格式

　　四个实数：a，b，c，d

输出格式

　　由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位

样例输入

1 -5 -4 20

样例输出

-2.00 2.00 5.00

数据规模和约定

　　|a|，|b|，|c|，|d|<=10

问题分析：二分加搜索

#include<bits/stdc++.h>
using namespace std;

double a,b,c,d;

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

double dfs(double l,double r)
{
	if (r - l <= 0.001)
		return l;
	
	double mid = (l+r)/2;
	if (mul(l)*mul(mid) <= 0)
		return dfs(l,mid);
	else
		return dfs(mid,r);
}
int main()
{
	cin >> a >> b >> c >> d;
	double ans[4];
	int cnt = 0;
	
	for(int i = -100; i <= 100; i++){
		double l = i,r = i+1;
		if (mul(l) == 0){
			ans[cnt++] = l;
		} 
		else if (mul(l)*mul(r) < 0){
			ans[cnt++] = dfs(l,r);
		}
		if (cnt == 3)
			break;
	}	
	
	printf("%.2lf %.2lf %.2lf",ans[0],ans[1],ans[2]);
	return 0;
	
}