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

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;

}



你可能感兴趣的:(二分法,uva,有效数字,解方程,单调)