PKU Helper 招新面试一

心血来潮去报了PKU Helper的招新，因为感觉他们都很厉害而且蛮好玩的。没想到居然有面试题，弱渣被吓得瑟瑟发抖23333本着做一做玩一玩的心理决定真的玩一玩啦~

T1
请自己实现一个 sqrt 函数，接受一个浮点型变量，返回其算数平方根，返回结果在 1e-6 的误差范围内相等即可。不得调用语言原生的 sqrt 库函数。

我的第一反应就是写个二分，估计这个精度要求切一切就出来了。然后再考虑下工程上呀特判呀的比如输入负数、NaN之类的怎么办就可以啦~然后check博客的时候惊喜地发现了神奇的“牛顿法”，大概就是小时候听说的方法居然又一次出现感觉非常神奇。于是决定学习下这种神奇的解法啦：）

//From DK
//Ref: http://blog.csdn.net/xusiwei1236/article/details/25657611
double MySqrt(double x)
{
    const double EPS=1e-7;
    const double NaN=0xFFF8000000000000;
    const double Infinity=0x7FF0000000000000;
    const double NegInfinity=0xFFF0000000000000;

    if (x<0) return NaN;
    if (x==NaN || x==Infinity || x==NegInfinity) return NaN;

    double a=x+1,b,point=a;
    for (;;)
    {
        b=(a*a+x)/(2*a);
        if (fabs(b-a)<=EPS) break;
        if (point==b) break;
        point=a;
        a=b;
    }
    return b;
}

开始的时候是把EPS啊NaN啊定义成宏的，后来考虑到交一个函数就都弄成常量塞函数里面啦~邮件里面说暂定最低要求为正确实现第一题，so还是满惶恐的。希望能够正确实现这一题啦~