发现一个很牛的开方算法

本文介绍了一种计算整数平方根的方法,通过位移和迭代实现高效运算,无论输入数值大小,均可在固定步骤内得出结果。此外,还探讨了一种快速求取浮点数倒数平方根的算法,该算法使用位操作和一次牛顿迭代,极大地提高了计算效率。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

static int SQRT (int nRoot)
 {
     int nSqrt = 0;

     for(int i = 0x10000000; i != 0; i >>= 2)
     {
         int nTemp = nSqrt + i;
         nSqrt >>= 1;
         if(nTemp <= nRoot)
         {
                 nRoot -= nTemp;
                 nSqrt += i;
          }
  }
    return nSqrt;
 }


用到了 (a+b)(a+b) = a*a + b*b + 2ab

测试了下, 由于 for 循环执行完成需16次循环。因此, 无论给定什么数, 都在16步给出结果。



Another way:


/*
** float q_rsqrt( float number )
*/  
float q_rsqrt( float number )  
{  
   long i;  
   float x2, y;  
   const float threehalfs = 1.5F;  
 
   x2 = number * 0.5F;  
   y  = number;  
   i  = * ( long * ) &y;                       // evil floating point bit level hacking  
   i  = 0x5f3759df - ( i >> 1 );               // what the fuck?  
   y  = * ( float * ) &i;  
   y  = y * ( threehalfs - ( x2 * y * y ) );   // 1st iteration  
//  y  = y * ( threehalfs - ( x2 * y * y ) );   // 2nd iteration, this can be removed  
 
#ifdef __linux__  
   assert( !isnan(y) ); // bk010122 - FPE?  
#endif  
   return y;  
}  



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值