分治算法问题 C: 一元三次方程求解

题目描述

有形如：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）之间一定有一个根。

输入

a，b，c，d

输出

三个实根（根与根之间留有空格）

样例输入 Copy

1 -5 -4 20

样例输出 Copy

-2.00 2.00 5.00

#include<iostream>
using namespace std;
double a,b,c,d;
double f(double x);
void solve(double x1,double x2);
int main()
{
double i;
cin >> a >> b >> c >> d;
for(i =-100;i <= 99;i++)
{
if(f(i)*f(i+1) <= 0)
{
solve(i,i+1);
}
}
}
double f(double x)
{
return a*x*x*x+b*x*x+c*x+d;
}
void solve(double x1,double x2)
{
if(x2-x1 < 0.001)
{
printf("%.2f ",x1);
return;
}
double x3 = (x1+x2)/2;
if(f(x2)==0)
{
printf("%.2f ",x2);
}
if(f(x3)==0){
printf("%.2f ",x3);
}
if(f(x1)*f(x3)<0)
solve(x1,x3);
if(f(x3)*f(x2)<0)
solve(x3,x2);
}