OpenJudge | 用二分法求方程的根

总时间限制: 1000ms 内存限制: 65536kB

描述

用二分法求下面方程在(-10, 10)之间的一个根。
2x3- 4x2+ 3x- 6 = 0

输入

一个小于1的非负实数e，它的值表示所能允许的误差

输出

一个实数，其值为求得的一个根，要求精确到小数点后8位。
若该区间上没有根，则输出“No Solution”

样例输入

0

样例输出

2.00000000

提示

对于一个连续函数f(x)，若f(a)*f(b) <= 0，则f(x)在区间[a, b]内至少有一个根。
特别的，对于一个单调的连续函数，上述定理得逆定理也成立
若[a, b]上有根，则可进一步考察根是否在 [a, (a+b)/2]内，或者在[(a+b)/2, b]内。

若b-a <= e 则可终止迭代，并认为(a+b)/2是一个近似解，将它输出

请使用double类型！

Code

#include <bits/stdc++.h>
using namespace std;double f(double x) {return 2.0*x*x*x-4.0*x*x+3.0*x-6.0;
}void fun(double *left, double *mid, double *right) {if(f(*left)*f(*mid) <= 0.0) {*right = *mid;*mid = ((*left)+(*right))/2;} else {*left = *mid;*mid = ((*left)+(*right))/2;}
}int main() {double left = -10, right = 10, mid = (left+right)/2;double e;scanf("%lf", &e);while(1) {fun(&left, &mid, &right);if(right-left <= e || right-left <= 0.000000001) {printf("%.8lf", mid);goto out;}}printf("No Solution");out:;
}