心血来潮去报了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还是满惶恐的。希望能够正确实现这一题啦~