NKOI 1408 【分类练习4.分治法】求方程的根

//蒟蒻的第四篇题解

输入m，n，p，a，b，求方程f(x)＝mx+nx-px＝0在[a,b]内的根。m，n，p，a，b均为整数，且a<b


m，n，p 都大于等于 1。

输入格式

m，n，p，a，b五个正整数。m,n,p大小都不超过1000，a,b不超过100。

输出格式

如果有根，则输出，精确到 1E-10；如果无方程根， 则输出“NO” 。

样例输入

2 3 4 1 2

样例输出

1.5071265916

题目大意

输入正整数m,n, p, a,b，求方程m~＋n一p=0的根。

解析

-二分法求根:

-设当前范围是Lr]，且fOf(r)<0

-设m=(l+r)/2.I如果f(Df(m)<0，则根在L m]中，赋值r = m;

·否则根在[m,r]中，赋值l = m;

·当r -l<A的时候停止，其中△用来控制二分的精度，这道题设为10-11或10-12为宜。

·循环结束后，返回(1+r)/2值。

#include<bits/stdc++.h>
using namespace std;
long long m,n,p,a,b;
double f(double x) {
	double ff;
	ff=pow(m,x)+pow(n,x)-pow(p,x);
	return ff;
}
void ans(double l,double r) {
	double mid;
	while(1e-11<r-l) {
		mid=(l+r)/2;
		if(f(l)*f(mid)<0)	r=mid;
		else	l=mid;
	}
	cout<<fixed<<setprecision(10)<<mid;
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>m>>n>>p>>a>>b;
	if(f(a)*f(b)>0)	cout<<"NO";
	else	ans(a,b);
	return 0;
}//坚持用cin
//czy6666666创作