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;
}