整数快速开平方算法

//平方根square root
unsigned int __stdcall i32_sqrt(unsigned long x)
{   unsigned int r, i;  // 结果、循环计数
    unsigned long t, k;  // 临时变量
    if (x == 0 || x== 1) // 被开方数，0,1开方结果为自身
        return x;
 
    r = 0;
 
    t = (x >> 30);          // 获取最高位
    x <<= 2;
    if (t > 1)              // 最高位为1
    {
        r ++;    // 结果当前位为1，否则为默认的0
        t -= r;
    }
 
    for (i=15; i>0; i--) // 求剩余的15位
    {
        r <<= 1;   // 左移一位
  
        t <<= 2;
        t += (x >> 30);  // 假设
  
        k = r;
        k = (k<<1)+1;
  
        x <<= 2;
        if (t >= k)   // 假设成立
        {
            t -= k;
            r ++;
        }
  
    }
 
    return r;
};
//平方根square root
unsigned long __stdcall i64_sqrt(unsigned __int64 x)
{ unsigned int i;   //循环计数
    unsigned long r;  // 结果
    unsigned __int64 t, k; //临时变量
    if (x == 0 || x== 1) // 被开方数，0,1开方结果为自身
        return (unsigned long)x;
 
    r = 0;
 
    t = (x >> 62);   // 获取最高位
    x <<= 2;
    if (t > 1)    // 最高位为1
    {
        r ++;    // 结果当前位为1，否则为默认的0
        t -= r;
    }
 
    for (i=31; i>0; i--) // 求剩余的31位
    {
        r <<= 1;   // 左移一位
  
        t <<= 2;
        t += (x >> 62);  // 假设
  
        k = r;
        k = (k<<1)+1;
  
        x <<= 2;
        if (t >= k)   // 假设成立
        {
            t -= k;
            r ++;
        }
  
    }
 
    return r;
};

 

//
// 计算32位整数的平方根
//
int __stdcall u32_sqrt(int val)
{
  int r = 0;
  for(int shift=0;shift<32;shift+=2)
  { int x=0x40000000l >> shift;
   
   if(x + r <= val)
   { val -= x + r;
    r = (r >> 1) | x;
   } else
   { r = r >> 1;
   }
  }
        //if(r < val) ++r;

        return r;
};

 

__int64 __stdcall u64_sqrt(__int64 val)
{
  __int64 r = 0;
  for(int shift=0;shift<64;shift+=2)
  { __int64 x=0x4000000000000000l >> shift;
   
   if(x + r <= val)
   { val -= x + r;
    r = (r >> 1) | x;
   } else
   { r >>= 1;
   }
  }

        //if(r < val) ++r;

        return r;
};