新博文地址:[leetcode]Sqrt(x)
http://oj.leetcode.com/problems/sqrtx/
Compute and return the square root of x.
初看这道题目的时候,想的太简单了,以为参数是int类型,返回值是int类型,就直接用了时间复杂度为O(n)的做法,从 x/2 开始往后遍历,知道找出满足 i * i == x的数,并返回。
后来勉强想到了二分查找,但是因为int类型发生的数组越界问题,好一顿纠结。最后看了别人的算法,看到了最合口味的牛顿迭代法:
算法思想也是很简单的,用到了极限的概念(大学数学木有白学啊,看到之后还是感觉自然而然的T_T)
不罗嗦:
为了方便理解,就先以本题为例:
计算x2 = n的解,令f(x)=x2-n,相当于求解f(x)=0的解,如左图所示。
首先取x0,如果x0不是解,做一个经过(x0,f(x0))这个点的切线,与x轴的交点为x1。
同样的道理,如果x1不是解,做一个经过(x1,f(x1))这个点的切线,与x轴的交点为x2。
以此类推。
以这样的方式得到的xi会无限趋近于f(x)=0的解。
判断xi是否是f(x)=0的解有两种方法:
一是直接计算f(xi)的值判断是否为0,二是判断前后两个解xi和xi-1是否无限接近。
经过(xi, f(xi))这个点的切线方程为f(x) = f(xi) + f’(xi)(x - xi),其中f'(x)为f(x)的导数,本题中为2x。令切线方程等于0,即可求出xi+1=xi - f(xi) / f'(xi)。
继续化简,xi+1=xi - (xi2 - n) / (2xi) = xi - xi / 2 + n / (2xi) = xi / 2 + n / 2xi = (xi + n/xi) / 2。
有了迭代公式xi+1= (xi + n/xi) / 2,程序就好写了。
public int sqrt(int x) {
if(x <= 1){
return x;
}
return newTon(x);
}
private int newTon(int x){
int result = 0;
double pre = x / 2 + 1;
double post = (pre + x / pre)/2;
while(Math.abs(pre - post) > 0.0000001){
pre = post ;
post = (pre + x / pre)/2;
}
result = (int) (pre + post) / 2;
return result;
}
至于二分法,资源很多,就不再啰嗦了。