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

最新推荐文章于 2019-08-07 11:25:17 发布

原创最新推荐文章于 2019-08-07 11:25:17 发布 · 468 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#浮点二分

二分法专栏收录该内容

2 篇文章

订阅专栏

利用二分法寻找多项式在指定区间的根，要求结果精确到小数点后2位。输入包含多项式的系数和区间，输出该根的值。例如，对于多项式3x^3 - x^2 - 3x + 1，在区间[-0.5, 0.5]内，根为0.33。" 122310674,9936310,2021年电工中级试题解析与实操视频,"['电工考试', '职业技能', '电气工程', '理论知识', '实操技巧']

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;
}