*循环-08. 二分法求多项式单根*

/*
 * Main.c
 * C8-循环-08. 二分法求多项式单根
 *  Created on: 2014年7月26日
 *      Author: Boomkeeper
 *****部分通过********
 */
#include <stdio.h>
#include <math.h>

float a3 = 0, a2 = 0, a1 = 0, a0 = 0;

double func(double x) {
    return (a3 * pow(x, 3) + a2 * pow(x, 2) + a1 * x + a0);
}


int main(void) {
    float a = 0, b = 0; //区间
    float midvalue = 0; //midvalue=(a+b)/2
    float threshold = 0;

    scanf("%f%f%f%f%f%f", &a3, &a2, &a1, &a0, &a, &b);

    do {
        midvalue = (a + b) / 2; //计算中间值
        if (func(a) * func(b) < 0) {
            if (func(midvalue) == 0) {
                printf("%.2f
", midvalue);
                return 0;
            } else {
                if (func(midvalue) * func(a) > 0) {
                    //同号，说明偏左了，应当往右偏
                    a = midvalue;
                    threshold = fabs(a - b); //计算区间阀值
                    continue;
                }
                if (func(midvalue) * func(b) > 0) {
                    //同号，说明偏右了，应当往左偏
                    b = midvalue;
                    threshold = fabs(a - b); //计算区间阀值
                    continue;
                }
            }
        }
    } while (threshold > 0.0001);
    /**
     * 因为threshold小到一定程度后，说明a、b可以看做成一点，也就是题目中f(r)=0的r点。
     */
    printf("%.2f
",(a+b)/2);

    return 0;
}

有更好方法的大虾们，希望指点，我试过将 while (threshold > 0.0001); 提高精度，测试结果仍然一样，我想还是算法上可能还有问题。

题目链接：

http://pat.zju.edu.cn/contests/basic-programming/%E5%BE%AA%E7%8E%AF-08