一个笔试题

并不是我去笔试，在其他地方看来的

求一个整数n的开方，不能使用库函数sqrt()，要求精确到0.001

这个说实话，可能就是那种，知道就知道，不知道就搞不来的题。

牛顿迭代法

#include "stdafx.h"
#include <iostream>

#define ABS(val) val>0?val:-val

double mysqrt(double x){
	double d1,d2;
	if(x==0)
		return 0;
	d1=x/2;
	d2=(d1+2)/2;
	while(ABS(d2-d1)>0.001){
		d1=d2;
		d2=(d1+x/d1)/2;
	}
	return d2;
}

int _tmain(int argc, _TCHAR* argv[])
{
	std::cout<<mysqrt(4)<<std::endl;

	return 0;
}

思路很巧妙，主要是个思路，当然答案不是很理想，比如mysqrt(1)，其答案就搞笑了