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

参考大佬

https://www.luogu.com.cn/blog/user55918/solution-h1024

题目描述

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

提示：记方程 f(x) = 0f(x)=0，若存在 22 个数 x_1x1​ 和 x_2x2​，且 x_1 < x_2x1​<x2​，f(x_1) \times f(x_2) < 0f(x1​)×f(x2​)<0，则在 (x_1, x_2)(x1​,x2​) 之间一定有一个根。

输入格式

一行，44 个实数 a, b, c, da,b,c,d。

输出格式

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

输入输出样例

输入 #1复制

1 -5 -4 20

输出 #1复制

-2.00 2.00 5.00

说明/提示

【题目来源】

NOIP 2001 提高组第一题

 

#include <cstdio>
double a, b, c, d;
int ans = 0;

double check(double n) {
	return a*n*n*n + b*n*n + c*n + d;;
}


int main() {
	double l, r, m, x1, x2;
	int ans=0;
	scanf("%lf%lf%lf%lf", &a, &b, &c, &d);
	//因为两根差大于等于1，所以区间长为一的区间最多只有一个根，我们可以枚举不重叠的区间，然后各个区间二分加精度求根 
	for (int i=-100; i<100; i++) {//枚举各个区间 
		l = i;
		r = i+1;
		x1 = check(l);
		x2 = check(r);
		if(!x1) { //若为0，则直接输出根 不能判断x2 
			printf("%.2lf ",l);
			ans++;
		}
		if (x1*x2<0) {
			while (r-l >= 0.001) {//保留两位，算到第三位，这样子前两位才是确定值 
				m = l + (r-l) / 2;
				if (check(m)*check(r) <= 0) {
					l = m;
				}
				else {
					r = m;
				}	
			}
			printf("%.2lf ", l);//输出哪个都一样，在保留两位的情况下值相同 
			ans++;
		}
		if (ans==3) 
        break;       
	}
	return 0;
}