[蓝桥杯][算法训练VIP]一元三次方程求解

时间限制: 1Sec 内存限制: 128MB 提交: 21 解决: 7

题目描述

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

数据规模和约定 
|a|，|b|，|c|，|d|< =20
 

输入

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

输出

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

样例输入

1  -5  -4  20 

样例输出

-2.00 2.00 5.00

方法一：枚举法

#include<bits/stdc++.h>

using namespace std;

int main()
{
	//freopen("in.txt","r",stdin);
	double a,b,c,d;
	scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
	for (double i=-100;i<=100;i+=0.01)
		if (abs(a*i*i*i+b*i*i+c*i+d)<=0.00001)
			printf("%.2lf ",i);
	printf("\n");
	return 0;
}

方法二：二分法求方程的根。利用根两侧的函数值相乘为负。

#include<bits/stdc++.h>

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()
{
	freopen("in.txt","r",stdin);
	scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
	for (double i=-100;i<=100;i++)
	{
		double x1 = i;
		double x2 = i + 1;
		if (f(x1)==0)
			printf("%.2lf ",x1);
		else if (f(x1)*f(x2)<0)
		{
			while (x2-x1>=0.001)
			{
				double x3 = (x1 + x2) / 2;
				if (f(x1)*f(x3)<=0)
					x2 = x3;
				else
					x1 = x3;
			}
			printf("%.2lf ",x1);
		}
	}
	printf("\n");
	return 0;
}