一元三次方程求解

一元三次方程求解

题目描述

有形如: 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;
}
相关推荐
虾球xz10 分钟前
游戏引擎学习第282天:Z轴移动与摄像机运动
c++·学习·游戏引擎
.小墨迹33 分钟前
Apollo学习——planning模块(3)之planning_base
linux·开发语言·c++·学习·自动驾驶
feifeigo12334 分钟前
基于粒子群算法的配电网重构
算法·重构
龙湾开发1 小时前
计算机图形学编程(使用OpenGL和C++)(第2版)学习笔记 10.增强表面细节(一)过程式凹凸贴图
c++·笔记·学习·3d·图形渲染
德亦周1 小时前
如何在Mac电脑上的VScode去配置C/C++环境
c++·vscode·macos
XiaoyaoCarter1 小时前
每日一道leetcode(新学数据结构版)
数据结构·c++·算法·leetcode·职场和发展·哈希算法·前缀树
晴空闲雲1 小时前
数据结构与算法-线性表-单链表(Linked List)
数据结构·算法·链表
zm1 小时前
服务器连接多客户端
java·javascript·算法
八月的雨季 最後的冰吻2 小时前
SIP协议栈--osip源码梳理
linux·服务器·网络·c++·网络协议
fancy1661662 小时前
搜索二维矩阵 II
c++·算法·矩阵