UVA - 10341 Solve It (高精度 二分法解方程)

题目链接:点击打开链接


题意 找一个0 到 1之间的数使得方程满足,若没有输出 No solution;

题解:首先方程根据给出的 数 的范围可知:方程为单调递减,这样就可以运用二分解方程了

代码:

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#include<math.h>
#define EE   1e-8    // 不是l,是数字1啊; 
double p,q,r,s,t,u;
double 	j_fang(double x)
{
	double k;
	k = p*exp(-x)+q*sin(x)+r*cos(x)+s*tan(x)+t*x*x+u;
	return k;
}

double er_feng()    // 我这个方程竟然1e-8 才过; 
{
	double l = 0,r = 1,mid;
	while(r - l>EE)
	{
		mid = (l+r)/2;
		double k = j_fang(mid);
		if(fabs(k)<=EE) return mid;
		else if(k<0)
			r = mid - EE;
		else
			l = mid + EE;
				 
	}
	return l;
}
int main()
{
	int i,j; 
	while(~scanf("%lf%lf%lf%lf%lf%lf",&p,&q,&r,&s,&t,&u))
	{
		if(fabs(j_fang(0.0))<=EE)
			printf("0.0000\n");
		else if(fabs(j_fang(1.0))<=EE)
			printf("1.0000\n");
		else if(j_fang(0)*j_fang(1)>0)
			printf("No solution\n");
		else
			printf("%.4f\n",er_feng());
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值