*循环-08. 二分法求多项式单根*
1 /* 2 * Main.c 3 * C8-循环-08. 二分法求多项式单根 4 * Created on: 2014年7月26日 5 * Author: Boomkeeper 6 *****部分通过******** 7 */ 8 #include <stdio.h> 9 #include <math.h> 10 11 float a3 = 0, a2 = 0, a1 = 0, a0 = 0; 12 13 double func(double x) { 14 return (a3 * pow(x, 3) + a2 * pow(x, 2) + a1 * x + a0); 15 } 16 17 18 int main(void) { 19 float a = 0, b = 0; //区间 20 float midvalue = 0; //midvalue=(a+b)/2 21 float threshold = 0; 22 23 scanf("%f%f%f%f%f%f", &a3, &a2, &a1, &a0, &a, &b); 24 25 do { 26 midvalue = (a + b) / 2; //计算中间值 27 if (func(a) * func(b) < 0) { 28 if (func(midvalue) == 0) { 29 printf("%.2f\n", midvalue); 30 return 0; 31 } else { 32 if (func(midvalue) * func(a) > 0) { 33 //同号，说明偏左了，应当往右偏 34 a = midvalue; 35 threshold = fabs(a - b); //计算区间阀值 36 continue; 37 } 38 if (func(midvalue) * func(b) > 0) { 39 //同号，说明偏右了，应当往左偏 40 b = midvalue; 41 threshold = fabs(a - b); //计算区间阀值 42 continue; 43 } 44 } 45 } 46 } while (threshold > 0.0001); 47 /** 48 * 因为threshold小到一定程度后，说明a、b可以看做成一点，也就是题目中f(r)=0的r点。 49 */ 50 printf("%.2f\n",(a+b)/2); 51 52 return 0; 53 }
有更好方法的大虾们，希望指点，我试过将 1 while (threshold > 0.0001); 提高精度，测试结果仍然一样，我想还是算法上可能还有问题。
题目链接：
http://pat.zju.edu.cn/contests/basic-programming/%E5%BE%AA%E7%8E%AF-08
