这几天好多次都遇到了二分法求方程的解。今天就花了点时间学了一下。其实,说起来也没有什么。高中的时候思想都是有的,只是当时,由于计算繁杂,不作为学习的重点。可是现在不一样了,我们现在,不需要计算,只需要告诉计算机怎么算,来让计算机计算就可以了。所以,现在这种方法也是很受用的吧。
二分法求方程的解的前提是,该函数是单调函数。
思想估计都没有什么问题。关键是代码的问题。
NYOJ 503 解方程为例点击打开链接
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const double eps=1e-6;//eps用来保证数据近似解的精度.当然,得到的解是越精确越好.但是,也就越花费时间.
double ans;
double f(double x){
return (((8*x-7)*x+2)*x+3)*x+6;
}
double solve(double left,double right){
double mid;
while(fabs(right-left)>eps){
mid=(right+left)/2;
if(f(mid)==ans) return mid;
if(f(mid)>ans) right=mid;
else left=mid;
}
return (left+right)/2;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
double max_val=f(100),min_val=f(0);
scanf("%lf",&ans);
if(ans>=min_val&&ans<=max_val)
printf("%.4lf\n",solve(0,100));
else printf("No solution!\n");
}
return 0;
}
NYOJ 965 Crossed Ladders
这题,也是可以用二分法求方程的近似解。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const double eps=1e-6;
double x,y,c;
double f(double ans){
double ans1=(y*y-ans*ans);
double ans2=(x*x-ans*ans);
return sqrt(ans1*ans2)/(sqrt(ans1)+sqrt(ans2));
}
double solve(double low,double high){
double mid;
while(high-low>eps){
mid=(high+low)/2;
double temp=f(mid);
// printf("%.3lf %.3lf %.3lf %.3lf\n",low,high,temp,c);
if(temp==c) return temp;
if(temp>c) low=mid;
else high=mid;
}
return (high+low)/2;
}
int main()
{
// freopen("output.txt","r",stdin);
// freopen("output1.txt","w",stdout);
while(scanf("%lf%lf%lf",&x,&y,&c)!=EOF)
{
double low=min(x,y);
printf("%.3lf\n",solve(0,low));
}
return 0;
}