【uva10341】二分法解方程,注意单调性

博客内容讲述了利用二分法解决0到1区间内的方程,通过求导判断单调性。若f(0)<0和f(1)>0,则方程无解。输出解时需使用`fixed`和`setprecision(4)`来控制精度。代码中,精度设置为10的-4次方不够精确,需提升至-9次方才能正确AC。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

p*e-xq*sin(x) + r*cos(x) + s*tan(x) + t*x2 + u = 0
        where 0 <= x <= 1.


好久没碰高数,差点连求导都换给老师了,还好哥数学基础比较扎实,求这个程度的导数还是轻轻松松的。


求个一次导数发现恒小于0,说明其在【0,1】上是单调递减的。


提前判断f(0)是否<0 和 f(1)是否 >0,如有上述两种情况直接输出无解。原因不用说了吧


输出的时候注意

cout<<fixed<<setprecision(4)<<x<<endl;

这样写是固定小数点后4位,不加fixed的话是保留四位有效数字,是不一样的。记得include<iomanip>

下面这个代码是ac不了的,控制循环精度10的-4次方太不精确,看网上有人说-9才可以,试了下果然,明明只保留四位,不明白为什么要精确到-9次方呢???

#include<iostream>
#include<cmath>
#include<iomanip>
#include<stdio.h>
using namespace std;

int p,q,r,s,t,u;


double fx(double x){
	return p*exp(-x)+ q*sin(x) + r*cos(x) + s*tan(x) + t*x*x + u;
}
int main()
{
	
	while(scanf("%d%d%d%d%d%d",&p,&q,&r,&s,&t,&u)!= EOF){
		if ((fx(1)>0)||fx(0)<0) cout<<"No solution"<<endl;
		else{
			double l = 0;
			double h = 1;
			double x = (h+l)/2;
			while(h-l>1e-4){
				if (fx(x)<0) h=x;
				else l=x;
				x=(h+l)/2;
			}
			cout<<fixed<<setprecision(4)<<x<<endl;
		}
		
	}
	return 0;

}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值