HDU-2675 Equation Again

#include <iostream>
#include <cmath>
#include <cstdio>
using namespace std;

const double e = 2.718281828459;
const double eps = 1e-7;
int main()
{
    double y;
    while(cin >> y)
    {
        double temp = log(e*y)/(e*y);
        double l = 0, r = e;
        while(r - l > eps)
        {
            double mid = (l + r) / 2;
            double ans = log(mid) / mid;
            if(ans >= temp) r = mid - eps;
            else l = mid;
        }
        if(fabs(y-1) > eps)
        printf("%.5f %.5f\n", r, e*y);
        else
        printf("%.5f\n", e*y);
    }
    return 0;
}

题意：给这个公式X (eY) == (eY) x，Y值确定，求X的值。

题解：是水题  还是水题。然而我还是纠结了2个小时。（感谢阿信叠的友情出演..╮(╯▽╰)╭）一道数学题，然而脑残了..两边取对数，得：lnx / x = ln(ey) / (ey) 再对lnx / x求导，发现（0，e）单调递增，（e，+∞）单调递减。然而题目定y>=1,则ey>=e 。然后就这么陷进去了...........（当时想法 既然>=e,不就单调递减么 单调函数么，然后就x = ey 不就一一对应了么？还算什么，不就x = ey么？还什么二分查值）....现在回想...有点醉。跟阿信叠说歪理也一本正紧 他也陷进去了.... 然而ln(ey) /(ey)，应该看做一整个常数，这样ey只是大于e的一个解 还有一个(0,e)的解。接下来只要取l = 0, r = e进行二分查值就好。这里是浮点数的二分查值，必须注意精度(eps)。还是得注意while(>) ，r = mid - ，l = mid。举个实例吧。一开始while(r - 1 > eps) r = mid - eps l =mid  输出 r +eps 是不过的。（注意只是>并不是>=）所以 应该只要输出r。（不需要+eps）（l也是不需要+eps的）。

总结浮点数形式应该是while(r - l > eps)  if(>=) r = mid - eps;else l = mid ; pirntf(r); （尝试了 while(r - l >= eps) if(>=) r = mid - eps;else l = mid + eps; pirntf(r + eps)是错的）。