Sqrt(x)

Sqrt(x)

Implement int sqrt(int x).

Compute and return the square root of x.

解题技巧：

方法1：采用二分搜索法

方法2：牛顿切线法

计算x² = n的解，令f(x)=x²-n，也就是求f(x)=0的解

首先取x₀，如果x₀不是解，做经过点(x₀,f(x₀))的切线，与x轴的交点为x₁；

同样地，如果x₁不是解，做经过点(x₁,f(x₁))的切线，与x轴的交点为x₂；

以此类推，以这样的方式得到的x_i会无限趋近于f(x)=0的解。

判断x_i是否是f(x)=0的解有两种方法：1. 直接计算f(x_i)的值判断是否为0；2. 判断前后两个解x_i和x_i-1是否无限接近

经过点(x_i,f(x_i))的切线方程为f(x) = f(x_i) + f'(x_i)(x - x_i)，其中f'(x)为f(x)的导数，本题中为2x。令切线方程等于0，可得x_i+1=x_i - f(x_i) / f'(x_i)。

简化：x_i+1=x_i - (x_i² - n) / (2x_i) = x_i - x_i / 2 + n / (2x_i) = x_i / 2 + n / 2x_i = (x_i + n/x_i) / 2。

最后得到迭代公式：x_i+1= (x_i + n/x_i) / 2

代码：

int mySqrt(int x)
{
    unsigned long long bg = 0, ed = x/2 + 1, mid, tmp;
    while(bg <= ed)
    {
        mid = (bg + ed) / 2;
        tmp = mid * mid;
        if(tmp == x) return mid;
        else if(tmp < x) bg = mid+1;
        else ed = mid-1;
    }
    return ed;
}

int mySqrt(int x)
{
    if (x ==0)  return 0;
    double pre;
    double cur = 1;
    do
    {
        pre = cur;
        cur = x / (2.0 * pre) + pre / 2.0;
    } while (abs(cur - pre) > 0.00001);
    return int(cur);
}