hdu2675二分解方程

#include<iostream>  
#include<cmath>  
#include<cstdio>
#define EXP 1e-8 
#define  f(x) log(x)/x
using namespace std;
const double e = 2.718281828459;
int main()
{
    //连边取对数 化简
    double Y;
    while (cin >> Y)
    {
        if (!(fabs(Y - 1)<EXP))//有第二个解 这个值在1到e之间(增区间)
        {
            double left = 1, right = e;
            double ey = e*Y;
            while ((right - left)>EXP)
            {
                double mid = (left + right) /2;

                if (f(mid) - f(ey)>0.000000000)
                {
                    right = mid;
                }
                else
                {
                    left = mid;
                }
            }
            printf("%0.5lf ", left);
        }
        printf("%0.5lf\n", e*Y);
    }
    return 0;
}