B - Can you solve this equation?

Description
现在，给出等式8* X^4+ 7* X^3+ 2* X^2+ 3 * X +6= Y，请找出他在0和100之间的解（包含0和100）。
现在，请你试试运气。。。。

Input

输入的第一行包含一个整数T（1 <= T <=100），表示测试用例的数目。接下来T个数字，每一行都有一个实数Y（abs（Y）<=10^10）;

Output
对于每个测试用例，如果有解，你应该输出一个实数（精确到小数点后4位，四舍五入），如果在0到100之间无解，就输出“No solution!”。

Sample Input

2

100

-4

Sample Output

1.6152

No solution!

题目很简单，直接用二分搜索法查找满足条件的x，排除y值可以取的范围，同时要注意结果的精度，四舍五入。
代码如下：

#include <iostream>
#include <iomanip>
#include <math.h>
#include <algorithm>
using namespace std;

const double EPS = 1e-7;

int f(double x, double y)
{
    if(8 * pow(x,4) + 7 * pow(x,3) + 2 * pow(x,2) + 3 * x + 6 < y)
        return 1;
    return 0;
}

double bs(double y)
{
    double lo = 0, hi = 100;
    double mi, ans=hi+1;
    while(lo + EPS  < hi)
    {
        mi = (lo + hi) / 2.0;
        if(f(mi,y))
        {
            lo = mi;
            if(ans == hi + 1)
                ans = mi;
            ans = max(ans,mi);
        }
        else
        {
            hi = mi;
            ans = min(ans,mi);
        }
    }
    if(ans == hi + 1)
        ans = -1;
    return ans;
}

int main()
{
    int t;
    double y;
    const long long max_num = 8 * pow(100,4) + 7 * pow(100,3) + 2 * pow(100,2) + 3 * 100 + 6;
    cin >> t;
    while(t--)
    {
        cin >> y;
        if(y > max_num || y < 6)
        {
            cout << "No solution!" << endl;
        }
        else
        {
            if(y == 6)
            {
                cout << fixed << setprecision(4) << 0.0000 << endl;
                continue;
            }
            if(y == max_num)
            {
                cout << fixed << setprecision(4) << 100.0000 << endl;
                continue;
            }
            double ans = bs(y);
            if(ans != -1)
                cout << fixed << setprecision(4) << (int)(ans * 10000 + 0.5) / 10000.0 << endl;
            else
                cout << "No solution!" << endl;
        }
    }
    return 0;
}