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

题目链接：7-18 二分法求多项式单根

一. 题目

1. 题目

2. 输入输出格式

3. 输入输出样例

4. 限制

二、代码

1. 代码实现

#include <stdio.h>
#include <math.h>

// 返回多项式在x点的值
float func(float *an, float x) {
    float x2 = x * x;
    float x3 = x2 * x;
    return an[3] * x3 + an[2] * x2 + an[1] * x + an[0];
}

float getX(float *an, float a, float b) {
    float x, fx, fa, fb;

    fa = func (an, a);
    if (fa == 0) { // 左端点为0
        return a;
    }

    fb = func (an, b);
    if (fb == 0) { // 右端点为0
        return b;
    }

    x = (a + b) / 2; // 防止此时已经满足精确度，但x未赋值的情况

    while (fabs(a - b) > 0.005) { // 根据精确度要求, 小于0.01/2 即可
        x = (a + b) / 2;
        fx = func (an, x);

        if (fa * fx > 0) { // f((a+b)/2)与f(a)同号
            fa = fx;
            fb = func (an, b);
            a  = x;
        } else if (fb * fx > 0) { // f((a+b)/2)与f(b)同号
            fa = func (an, a);
            fb = fx;
            b  = x;
        } else { // 中间 f((a+b)/2)正好为0
        	// 思考：为什么不会是左端点或者右端点为0
            break;
        }
    }

    return x;
}

int main(void) {
    float an[4];
    float a, b;

    if (scanf("%f %f %f %f", &an[3], &an[2], &an[1], &an[0]) != 4) {
        return -1;
    }

    if (scanf("%f %f", &a, &b) != 2) {
        return -1;
    }

    printf ("%.2f\n", getX(an, a, b));

    return 0;
}

2. 提交结果