PAT 5-18 二分法求多项式单根 浮点数二分法

5-18 二分法求多项式单根 (20分)

题解：

没什么好解的，浮点数的二分法，就是计算l的函数值lf，中点m的函数值mf，右边r的
函数值rf，mf*lf<0那根就在左边，令r=m;  mf*rf<0那根就在右边，令l=m;
但是只写了上述过程只能得15分，因为漏了对端点的判断。abs(lf-0)<eps和
abs(rf-0)<eps(如果直接采用==0的形式判定只能得18分)

好心塞啊，因为卡在了之前那个“一帮一”的简单题目上，导致整个人变蠢，一直以为自己
这题已经做了，最后三分钟才发现漏做这道二分法 = = 不过，这也是一次不错的经验吧！
下次加油咯！

#include <iostream>
#include <string.h>
#include <algorithm>
#include <stdio.h>
#include <math.h>
#include <stack>
#include <queue>
#include <vector>
#define INF 0xffffffff
using namespace std;

double a[4];
double l, r;

double f(double x){
return a[0]*x*x*x+a[1]*x*x+a[2]*x+a[3];
}

int main(){
while(scanf("%lf%lf%lf%lf%lf%lf", &a[0], &a[1], &a[2], &a[3], &l, &r) == 6){
double lf = 0;
double rf = 0;
double m = 0;
while(l<r){
lf = f(l);
rf = f(r);
if(abs(lf-0)<0.00000001){
m = l;
break;
}
if(abs(rf-0)<0.00000001){
m = r;
break;
}
m = (l+r)/2;
double mf = f(m);
if(lf*mf<0){
r = m;
}
else{
l = m;
}
}
printf("%.2lf\n", m);
}
return 0;
}