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

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

二分法求函数根
输出格式：
在一行中输出该多项式在该区间内的根，精确到小数点后2位。
输入样例：
3 -1 -3 1
-0.5 0.5
输出样例：
0.33

题目链接：https://pta.patest.cn/pta/test/14/exam/4/question/798

题解：

没什么好解的，浮点数的二分法，就是计算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;
}