一元三次方程求解（暴力二分）

BreaCham

于 2023-01-09 21:39:43 发布

阅读量218

点赞数

分类专栏：刷题感悟文章标签：算法

本文链接：https://blog.csdn.net/BreaCham/article/details/128621563

版权

刷题感悟专栏收录该内容

2 篇文章 0 订阅

订阅专栏

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

题目描述

有形如： $a x^3 + b x^2 + c x + d = 0$ 这样的一个一元三次方程。给出该方程中各项的系数（ $a, b, c, d$ 均为实数），并约定该方程存在三个不同实根（根的范围在 $- 100$ 至 $100$ 之间），且根与根之差的绝对值 $\ge 1$ 。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 $2$ 位。

提示：记方程 $f (x) = 0$ ，若存在 $2$ 个数 $x_1$ 和 $x_2$ ，且 $x_1 < x_2$ ， $f(x_1) \times f(x_2) < 0$ ，则在 $x_1, x_2)$ 之间一定有一个根。

输入格式

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

输出格式

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

样例 #1

样例输入 #1

1 -5 -4 20

样例输出 #1

-2.00 2.00 5.00

提示

【题目来源】

NOIP 2001 提高组第一题

题干中这句话：

根的范围在 $- 100$ 至 $100$ 之间，且根与根之差的绝对值 $\ge 1$

就说明了可以暴力遍历-100~100 之间200个长度为1的小区间（一开始真没想到）。

注意细节

参考代码：

#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() {
	cin >> a >> b >> c >> d;
	int cnt = 0;
	for (double i = -100; i < 100; ++i) {//直接暴力遍历（题目给的区间很小）
		if (f(i) == 0) {
			++cnt;
			printf("%.2lf ", i);
			if (cnt >= 3)
				break;
			continue;
		}
		if (f(i) * f(i + 1) > 0)
			continue;
		double l = i, r = i + 1;
		if (f(i) * f(i + 1) < 0) {
			while (fabs(f(l)) > 1e-4) {//记得绝对值
				double mid = (r + l) / 2;
				if (f(l) * f(mid) < 0)//别忘了f
					r = mid;
				else
					l = mid;
			}
			printf("%.2lf ", l);
			++cnt;
		}
	}
	return 0;
}