69.x的平方根

实现 int sqrt(int x) 函数。

计算并返回 x 的平方根，其中 x 是非负整数。

由于返回类型是整数，结果只保留整数的部分，小数部分将被舍去。

示例 1：输入: 4，输出: 2

示例 2：输入: 8，输出: 2（说明: 8 的平方根是 2.82842..., 由于返回类型是整数，小数部分将被舍去。）

 

1.二分法：就是大了就往小的找，小了就往大了找

int mySqrt(int x) 
{
    if(x == 1)
        return 1;
    int min = 0;
    int max = x;
    while(max-min>1)
    {
        int m = (max+min)/2;
        if(x/m<m)
            max = m;
        else
            min = m;
    }
    return min;
}

2.根据数学公式求：exp(0.5 * log(x))就是x的平方根的意思

int mySqrt(int x)
{
    int answer;
    if(x == 0)
        return 0;
    if(x == 1)
        return 1;
    
    answer = exp(0.5 * log(x));

    if(x >=(long long) (answer + 1) * (answer + 1))
    {
        return (long long)(answer + 1);
    }
    else
    {
        return (long long)answer;
    }
}

3.牛顿迭代法

int mySqrt(int x){
    if(x==0 || x==1) return x;
    long xn = x;
    while(xn > x / xn){
        xn = (xn + x/xn)/2;
    }
    return xn;
}

4.卡马克算法：作者是Nvidia的Gary Tarolli：http://www.matrix67.com/data/InvSqrt.pdf（改天折腾一下）

再加一个数学编程论坛：https://www.gamedev.net/

int InvSqrt(float x)
{
    float xhalf = 0.5f*x;
    int i = *(int*)&x; // get bits for floating VALUE 
    i = 0x5f375a86- (i>>1); // gives initial guess y0
    x = *(float*)&i; // convert bits BACK to float
    x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
    x = x*(1.5f-xhalf*x*x); // Newton step, repeating increases accuracy
    return (int)1/x;
}  

int mySqrt(int x){
    float p = (float)x;
    int ans=InvSqrt(p);
    if (ans*ans > x) return ans-1;
    return ans;
}

 

这题俨然进入数学范畴，头疼，自己摸着边用个二分法搞定。