一元三次方程求解

一元三次方程求解

题目描述

有形如: a x 3 + b x 2 + c x + d = 0 ax^3+bx^2+c^x+d=0 ax3+bx2+cx+d=0一元三次方程。给出该方程中各项的系数 ( a a a, b b b, c c c, d d d 均为实数 ),并约定该方程存在三个不同实根 (根的范围在 − 100 -100 −100至 100 100 100之间 ),且根与根之差的绝对值 ≤ 1 ≤1 ≤1。要求由小到大依次在同一行上输出这三个实根。

提示

记方程 f ( x ) = 0 f(x)=0 f(x)=0,若存在两个数 x 1 x_1 x1和 x 2 x_2 x2,且 x 1 < x 2 x_1 < x_2 x1<x2, f ( x 1 ) × f ( x 2 ) < 0 f(x_1)×f(x_2)<0 f(x1)×f(x2)<0,则在 ( x 1 , x 2 ) (x_1,x_2) (x1,x2)之间一定有一个根。

输入输出格式

输入格式

输入 a , b , c , d a,b,c,d a,b,c,d

输出格式

三个实根(根与根之间留有空格)

输入输出样例

输入样例

1 -5 -4 20

输出样例

-2.00 2.00 5.00


算法分析

这是一道有趣的解方程题。为了便于求解,设方程 f ( x ) = a x 3 + b x 2 + c x + d = 0 f(x)=ax^3+bx^2+cx+d=0 f(x)=ax3+bx2+cx+d=0,设根的值域( − 100 -100 −100至 100 100 100之间)中有 x x x, 其左右两边相距 0.0005 0.0005 0.0005的地方有 x 1 x_1 x1和 x 2 x_2 x2两个数,即 x 1 = x − 0.0005 x_1=x-0.0005 x1=x−0.0005, x 2 = x + 0.0005 x_2=x+0.0005 x2=x+0.0005。 x 1 x_1 x1和 x 2 x_2 x2间的距离( 0.001 0.001 0.001)满足精度要求(精确到小数点后 2 2 2位)。


A C 代码

暴力出奇迹~~

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
double a,b,c,d;
int main()
{
	cin >>a >>b >>c >>d;
	for(double i=-100;i<=100;i+=0.001)
	{
		double j=i+0.001;
		double l=a*i*i*i+b*i*i+c*i+d;
		double r=a*j*j*j+b*j*j+c*j+d;
		if(l*r<=0)
		{
			printf("%.2lf ",(i+j)/2);
		}
	}
	return 0;
}
相关推荐
GGBondlctrl6 分钟前
【leetcode】递归,回溯思想 + 巧妙解法-解决“N皇后”,以及“解数独”题目
算法·leetcode·n皇后·有效的数独·解数独·映射思想·数学思想
武子康18 分钟前
大数据-276 Spark MLib - 基础介绍 机器学习算法 Bagging和Boosting区别 GBDT梯度提升树
大数据·人工智能·算法·机器学习·语言模型·spark-ml·boosting
武子康21 分钟前
大数据-277 Spark MLib - 基础介绍 机器学习算法 Gradient Boosting GBDT算法原理 高效实现
大数据·人工智能·算法·机器学习·ai·spark-ml·boosting
成工小白1 小时前
【C++ 】智能指针:内存管理的 “自动导航仪”
开发语言·c++·智能指针
sc写算法1 小时前
基于nlohmann/json 实现 从C++对象转换成JSON数据格式
开发语言·c++·json
SunkingYang1 小时前
C++中如何遍历map?
c++·stl·map·遍历·方法
Andrew_Xzw1 小时前
数据结构与算法(快速基础C++版)
开发语言·数据结构·c++·python·深度学习·算法
库库的里昂1 小时前
【C++从练气到飞升】03---构造函数和析构函数
开发语言·c++
momo卡1 小时前
MinGW-w64的安装详细步骤(c_c++的编译器gcc、g++的windows版,win10、win11真实可用)
c语言·c++·windows