P1024一元三次方程求解题解

题目分析

这是一道有趣的解方程题。为了便于求解,设方程 f(x)=ax3+bx2+cx+d=0f(x)=ax3+bx2+cx+d=0f(x)=ax3+bx2+cx+d0,设根的值域(-100至100之间)中有 xxx, 其左右两边相距 0.00050.00050.0005 的地方有 x1x1x1x2x2x2
两个数,即 x1=x−0.0005x1=x-0.0005x1=x0.0005x2=x+0.0005x2=x+0.0005x2=x+0.0005。只要x1x1x1x2x2x2 间的距离满足精度要求(精确到小数点后 222 位)。因此我们就可以使用暴力枚举

AC代码

#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
int main() {
	cin>>a>>b>>c>>d;
	for(double i=-100;i<=100;i+=0.001){
		double j=i+0.001;
		double y1=a*i*i*i+b*i*i+c*i+d;
		double y2=a*j*j*j+b*j*j+c*j+d;
		if(y1*y2<=0){
			cout<<fixed<<setprecision(2)<<i<<" ";
		}
	}
}
### NOIP2001 提高组 一元三次方程求解 题解 #### 实数域二分查找的应用 针对该问题,实数域上的二分查找用于解决函数零点问题。具体实现方式如下: 当给定区间 \([l, r]\),通过不断缩小范围来逼近根的位置。每次迭代计算中间值 \(mid\) 并评估其对应的多项式的符号,以此决定更新左边界还是右边界直到满足精度要求。 ```cpp #include <iostream> #include <cmath> using namespace std; double a, b, c, d; const int N = 100; // 假设最多有N个解 pair<double, double> ans[N]; int cnt = 0; // 计算f(x)=ax^3+bx^2+cx+d的值 double f(double x) { return (((a * x + b) * x + c) * x + d); } void solve(double l, double r) { if (fabs(f(l)) < 1e-8 || fabs(f(r)) < 1e-8) { // 如果端点本身就是解,则记录并返回 printf("%.2lf\n", min(fabs(f(l)) < 1e-8 ? l : r)); return; } while (r - l >= 1e-6) { // 控制精度至小数点后六位 double mid = (l + r) / 2; if (f(mid) == 0) break; else if (f(mid) * f(l) > 0) l = mid; else r = mid; } printf("%.2lf ", (l + r) / 2); // 输出平均作为近似解 } ``` 此段代码实现了对指定区间的扫描,并利用二分法逐步锁定可能存在的实根位置[^1]。 #### 处理多个解的情况 由于可能存在多于一个实根,在遍历整个定义域时需特别注意区分不同间隔内的变化趋势。每当检测到相邻两点间存在异号情况时启动一次新的二分解过程,从而确保不会遗漏任何潜在的真实交点。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值