HDU 2199 Can you solve this equation?

这道题主要是考察“二分法”，是二分法的经典应用（要求函数连续且递增）。层层逼近最佳值。

1.先定义一个函数f(double x) 用来表示求把x带入函数f的值。f就是 8*x^4 + 7*x^3 + 2*x^2 + 3*x + 6 。

2.由于函数递增，容易得到y < f(0)时没有结果，且y > f(100)时没有结果。此时输出“No solution!”。

3.当有结果可以输出时。用二分法，如果f(mid) > y ,则right = mid, 如果f(mid) < y , 则left = mid。（一直到right - left <= epsilon为止 ，epsilon是一个很小的数）。

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;
const double eps = 1e-8;
double flag;
double f(double x) {
    return 8 * pow(x, 4) + 7 * pow(x, 3) + 2 * pow(x, 2) + 3 * x + 6;
}
double cal(double y) {
    if(f(0) > y || f(100) < y) {
        flag = false;
        return -1;
    }
    double l = 0, r = 100;
    while(r - l > eps) {
        //cout << l << "\t" << r << endl;
        double mid = (l + r) / 2;
        if(f(mid) > y) {
            r = mid;
        } else {
            l = mid;
        }
    }
    return l;
}
int main() {
    int T;
    double y;
    cin >> T;
    while(T--) {
        flag = true;
        cin >> y;
        double solve = cal(y);
        if(!flag) {
            printf("No solution!\n");
        } else {
            printf("%.4f\n", solve);
        }
    }
    return 0;
}

Note：精度的问题，刚开始我的精度是1e-6,发现1.6151，而结果是1.6152，再设置到1e-8时可以通过。我听一个老师说如果小数点精确到3位，eps就设为1e-6,4位就设为1e-8，这是根据经验得来的。.