7-18 二分法求多项式单根

二分法求函数根的原理为：如果连续函数f(x)在区间[a,b]的两个端点取值异号，即f(a)f(b)<0，则它在这个区间内至少存在1个根r，即f®=0。
二分法的步骤为：


输入格式：
输入在第1行中顺序给出多项式的4个系数a​3​​、a​2​​、a​1​​、a​0​​，在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。

输出格式：
在一行中输出该多项式在该区间内的根，精确到小数点后2位。

输入样例：
3 -1 -3 1
-0.5 0.5

输出样例：
0.33

错误代码：

#include <iostream>
using namespace std;
double f(double x,double a3,double a2,double a1,double a0)
{
      return x*x*x*a3+x*x*a2+a1*x+a0;
}
int main()
{
    double a0,a1,a2,a3;
    double start,end;
    cin>>a3>>a2>>a1>>a0;
    cin>>start>>end;
    double mid = start+(end-start)/2;
    while(f(mid,a3,a2,a1,a0)!=0)
    {
        if(f(mid,a3,a2,a1,a0)>0)
        {
            end = mid;
        }
        else{
            start = mid;
        }
        mid = start+(end-start)/2;
    }
    printf("%.2lf",mid);
}

问题分析：

1. 程序员又又又不需要眼睛了，题目并没有指出f(x)为递增函数，所以用＞＜0判断是错误的
2. 按照上述代码会陷入死循环

修改代码
注意到错误后，我修改了代码，用示例运行成功，但提交后提示中点的f值很小但不一定是根说明可能需要求近似根

看了发现题目二分法的第一步已经给了提示，看来我的眼睛真的可以捐了

正确代码

#include <iostream>
using namespace std;
double f(double x,double a3,double a2,double a1,double a0)
{
      return x*x*x*a3+x*x*a2+a1*x+a0;
}
int main()
{
    double a0,a1,a2,a3;
    double start,end;
    cin>>a3>>a2>>a1>>a0;
    cin>>start>>end;
    double mid = start+(end-start)/2;
    while(f(mid,a3,a2,a1,a0)!=0&&end-start>0.0001)
    {
        if(f(mid,a3,a2,a1,a0)*f(end,a3,a2,a1,a0)>0)
        {
            end = mid;
        }
        else{
            start = mid;
        }
        mid = start+(end-start)/2;
    }
    printf("%.2lf",mid);
}

总结
眼睛不需要捐给需要的人