P1024 [NOIP2001 提高组] 一元三次方程求解

题目描述

有形如:a x^3 + b x^2 + c x + 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​) 之间一定有一个根。

输入格式

一行,4 个实数 a,b,c,d。

输出格式

一行,3 个实根,从小到大输出,并精确到小数点后 2 位。

输入输出样例

输入 #1

1 -5 -4 20

输出 #1

-2.00 2.00 5.00

说明/提示

【题目来源】

NOIP 2001 提高组第一题

【分析】

用二分查找找答案,再按公式判断是否正确。

【AC代码】

#include<iostream>
#include<cmath>
#include<cstdio>

using namespace std;

double a,b,c,d;

double f(double x) {
	return a * pow(x, 3) + b * pow(x, 2)+ c * x + d; 
}

int main() {
	cin >> a >> b >> c >> d;
	for (int i = -100; i <= 100; i++) {
		double x1 = i, x2 = i + 1;
		if (f(x1) == 0) { 
			printf("%.2lf ", x1);
		} else if (f(x1) * f(x2) < 0) {
			while (x2 - x1 > 0.0000000000001) {
				double mid = (x1 + x2) / 2;
				if  (f(mid) == 0) {
					x1 = mid;
					break;
				}
				if (f(x1) * f(mid) < 0) { 
					x2=mid;
				} else {
					x1=mid;
				}
			}
			printf("%.2lf ", x1);
		}
	}
	return 0;
}

上一篇:C# winform 生成两位随机数


下一篇:冬季实战营学习报告