题目链接: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;
}