UVA - 10341 Solve It (数学--二分法求解)

大体题意:

给你一个函数,给你定义域,求出是否有解,有解输出解,无解输出 No solution

思路:

分析式子和输入  知道这个式子是一个单调递减的,那么就可以用二分了!

当两个端点 的函数值乘积大于0时,说明同号,则无解!

否则有唯一解!  直接二分找就可以了!

教训:

这种求解方程的二分,涉及到小数的二分 ,最好用循环多少来解,不要用l < r  容易出错!!

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
double p,q,r,s,t,u;
const double e =  2.718281828459045;
const double eps = 1e-10;
double f(double x){
	return p*pow(e,-x) + q*sin(x) + r * cos(x) + s*tan(x) + t * x * x + u;
}
int main(){
	while(scanf("%lf %lf %lf %lf %lf %lf",&p, &q, &r, &s, &t, &u) == 6){
		if (fabs(f(1)) < eps)printf("1.0000\n");
		else if (fabs(f(0)) < eps)printf("0.0000\n");
		else if (f(1) * f(0) > eps)printf("No solution\n");
		else {
			double l = 0;
			double r = 1;
			for (int i = 0; i < 100; ++i){
				double mid = (l + r) / 2.0;
				if (f(mid) > 0)l = mid;
				else r = mid;
			}
			printf("%.4lf\n",(l+r)/2.0);		
		}
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值