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

P1024 [NOIP2001 提高组] 一元三次方程求解 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

#include <bits/stdc++.h>
using namespace std;
double a, b, c, d;
double f[10];

double jisuan(double x) {
	return a * x * x * x + b * x * x + c * x + d;
}

int main() {
	//double a,b,c,d;
	int ans = 0;
	cin >> a >> b >> c >> d;
	for (int i = -100; i <= 100; i++) {
		double l = i;
		double r = i + 1;
	//	cout << l << " " << r << endl;
		double x1 = jisuan(l);
		double x2 = jisuan(r);
		if (!x1) {
			f[ans] = l;
			ans++;
		}
		/*if(!x2){//判断一个端点就好，另一个端点会遍历到
			f[ans]=1;
			ans++;
		}*/
		if (x1 * x2 < 0) {
			while (r - l >= 0.001) {//二分边界超过答案精度
				double mid = (r + l) / 2;
				if (jisuan(mid)*jisuan(r) > 0)
					r = mid;//范围往值小（右端点左移）的方向移动
				else
					l = mid;//范围往大（左端点右移）的方向移动
			}
			f[ans] = r;
			ans++;
		}
		if (ans == 3)
			break;
	}
	sort(f, f + 3);
	for (int i = 0; i < 3; i++) {
		if (i)
			cout << " ";
		cout << fixed << setprecision(2) << f[i];
	}
	return 0;
}