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

//新生训练

#include <bits/stdc++.h>
using namespace std;
const double N = 1e-8;
double a, b, c, d;
double ans[3];
int k = 0;
double fx(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 && k < 3; i++)
	{
		double f1 = fx(i), f2 = fx(i + 1);
		if (fabs(f1) < N)
		{
			ans[k++] = i;
			continue;
		}
		else if (f1 * f2 < 0 && fabs(f2) > N)
		{
			double l = i, r = i + 1;
			while (fabs(l - r) >= 0.001)
			{
				double mid = (l + r) / 2;
				if (fabs(fx(mid)) < N)
				{
					ans[k++] = mid;
					break;
				}
				else if (fx(mid) * fx(l) < 0)
					r = mid;
				else
					l = mid;
			}
			ans[k++] = l;
		}
	}
	printf("%.2f %.2f %.2f", ans[0], ans[1], ans[2]);
	return 0;
}

 //一道对于理科生很实用的题 hwh；

//一开始笔者对于浮点型数据的定义是 float ，之后改成了 double ；

//解释一下 pow（） ；是幂函数的调用语句，pow（x，y）等价于 x^y ，也就是x的y次幂；

~~~//仅当笔者个人备忘录使用。