一元三次方程求解(二分法)

题目:[NOIP2001]一元三次方程求解 (nowcoder.com)

#include<bits/stdc++.h>
#define int long long
#define PII pair<int,int>
#define endl '\n'
#define ll long long

using namespace std;
//const int N=1e5+10,M=5010,mod=998244353,INF=0x3f3f3f;
double a,b,c,d;
int ans=0;
inline double f(double x)
{
	return a*x*x*x+b*x*x+c*x+d;
}

signed main()
{
	scanf("%lf %lf %lf %lf",&a,&b,&c,&d);
	for(double i=-100;i<=99;i++)
	{
//题目中说了只有3个,所以当我们已经找到三个的时候可以直接退出
		if(ans>=3) break;
		if(f(i)*f(i+1)<=0)
		{
			if(f(i)==0) printf("%.2f ",i),ans++;
			else if(f(i+1)==0) printf("%.2f ",i+1),i++,ans++;
//为了确保不重复计算这个根,并且能够继续寻找其他的根,需要将 i 值递增到下一个位置,即 i++。这样在//下一轮循环中,i 的值就从下一个位置开始继续搜索其他的根。
			else{ 
 //说明在区间 [i, i+1] 内存在一个根,使用二分法进行精确求解。
				double l=i;
				double r=i+1;
				double mid=(l+r)/2;
//精度为2位小数,说明这个区域范围要小于0.01
				while(r-l>0.01)
				{
					if(f(l)*f(mid)<=0)
					{
					r=mid;	
					 } 
					else {
					l=mid;	
					}
					mid=(l+r)/2;
				}
				printf("%.2f ",mid);
				ans++;
			}
		}
	}
	return 0;
}

从这道题学会的:

  •  二分法的实际运用
  • 题意理解!!
  • 注意输出格式
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值