[二分查找]洛谷P1024 一元三次方程求解

Echo♛

已于 2023-09-18 19:26:36 修改

阅读量183

点赞数 2

分类专栏：洛谷练习文章标签：算法

于 2023-09-18 19:23:31 首次发布

本文链接：https://blog.csdn.net/qq_73058399/article/details/132993144

版权

洛谷练习专栏收录该内容

2 篇文章 0 订阅

订阅专栏

文章讲述了如何在给定一个三次多项式函数的情况下，利用分段扫描和二分法求解零点，特别强调了当题目中存在多根且距离大于1时的处理方法。代码展示了如何针对区间判断并执行浮点数二分搜索。

摘要由CSDN通过智能技术生成

二分法最关键的弱点就是对于多根的问题束手无策，在这个题中为了使得二分可用，题目设计了两根的距离大于1，这使得分段扫描成为可能，这也就是分治的思维

首先先从-100到100，以间隔为1进行扫描，找到可能的区间
然后对于区间进行处理，首先要考虑到区间头和区间尾可能为0，为了使得前后不重复处理，我们可以将区间分成[i,i+1)的形式，先判断i时候是否为0，然后再判断i和i+1的函数值是否异号
如果异号且i+1的时候不为0则进行浮点数二分处理

代码如下：

#include<iostream>
#include<cmath>
using namespace std;
float a, b, c, d;
float ans[3];
int k = 0;
float fx(float x)
{
	return a * pow(x, 3) + b * pow(x, 2) + c * x + d;
}
int main()
{
	scanf("%f %f %f %f", &a, &b, &c, &d);
	for (int i = -100; i < 100 && k < 3; i++)
	{
		float f1 = fx(i), f2 = fx(i + 1);
		if (fabs(f1) < 1e-8)
		{
			ans[k++] = i;
			continue;
		}
		else if (f1 * f2 < 0 && fabs(f2)>1e-8)
		{
			float left = i, right = i + 1;
			while (fabs(left - right) >= 0.001)
			{
				float mid = (left + right) / 2;
				if (fabs(fx(mid)) < 1e-8)
				{
					ans[k++] = mid;
					break;
				}
				else if (fx(mid) * fx(left) < 0)
					right = mid;
				else left = mid;
			}
			ans[k++] = left;
		}
	}
	printf("%.2f %.2f %.2f", ans[0], ans[1], ans[2]);
	return 0;
}