用+,-,移位实现除法运算

本文介绍两种不使用除法运算符实现整数除法的方法:一种是通过循环减法累加商;另一种是利用特殊数值(魔数)进行位操作计算商。详细解释了算法原理,并提供了C++代码示例。

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



如何利用+,-,移位来实现除法.假设求 dividend / divisor

方法一:

首先想到的是用减法来实现.

算法思想:

对dividend 循环减 divisor, 减一次res++, 直到刚好减为0或余数小于divisor.

  1. int integer_div_1(unsigned int dividend, unsigned int divisor)  
  2. {  
  3.     if (divisor == 0)  
  4.     {  
  5.         cout << "除数不能为0" << endl;  
  6.         exit(1);  
  7.     }  
  8.   
  9.     int res = 0;  
  10.     while ((dividend -= divisor) >= 0)  
  11.         ++res;  
  12.   
  13.     return res;  
  14. }  
int integer_div_1(unsigned int dividend, unsigned int divisor)
{
    if (divisor == 0)
    {
        cout << "除数不能为0" << endl;
        exit(1);
    }

    int res = 0;
    while ((dividend -= divisor) >= 0)
        ++res;

    return res;
}


这个算法每次都以一倍的divisor进行叠加,算法效率并不高,可以改进.

以100/3为例.

算法分别比较97, 94, 91, ..., 4,1, -2,最后dividend = -2退出while循环.算法比较了34次.

如果采用每次采用将比较数翻倍的比较方法. 算法会得到优化. 举例如下: k初始化为0,  res = 0;

首先用3与100比,小于. 然后翻倍6, 小于. 12, 24, 48, 96, 192, 因为192 > 100.  退回到 96. 这里共比较了 5次. 每比较一次 k++, res += 1<<k.

100 - 96 = 4 > 除数3. 再用4重做上一步. 先跟3比较, 然后6, 6 > 4. 这次比较2次.

回到3. 4 - 3 = 1 < 除数3. 算法停止. 总共比较了5 + 2 + 1 = 8次, 比原来的34次快了很多.

代码实现如下:

  1. //递归代码  
  2. int integer_div_2(unsigned int dividend, unsigned int divisor)  
  3. {  
  4.     if (divisor == 0)  
  5.     {  
  6.         cout << "除数不能为0" << endl;  
  7.         exit(1);  
  8.     }  
  9.   
  10.     if (dividend < divisor) return 0;  
  11.     unsigned int k = 0, c = divisor, res = 0;  
  12.       
  13.     for ( ; dividend >= c; c <<= 1, k++)  
  14.         if (dividend - c < divisor)  
  15.             return 1<<k;  
  16.   
  17.     return integer_div_1(dividend - (c>>1), divisor) + (1<<(k - 1));  
  18. }  
//递归代码
int integer_div_2(unsigned int dividend, unsigned int divisor)
{
    if (divisor == 0)
    {
        cout << "除数不能为0" << endl;
        exit(1);
    }

    if (dividend < divisor) return 0;
    unsigned int k = 0, c = divisor, res = 0;
    
    for ( ; dividend >= c; c <<= 1, k++)
        if (dividend - c < divisor)
            return 1<<k;

    return integer_div_1(dividend - (c>>1), divisor) + (1<<(k - 1));
}
  1. //非递归算法   
  2. int integer_div_3(unsigned int dividend, unsigned int divisor)  
  3. {  
  4.     if(divisor == 0)  
  5.     {  
  6.         cout << "除数不能为0" << endl;  
  7.         exit(1);  
  8.     }  
  9.   
  10.     if (dividend < divisor)   
  11.         return 0;  
  12.       
  13.     unsigned int k, c, res=0;  
  14.   
  15.     while (dividend > divisor)  
  16.     {  
  17.         for (k = 0,c = divisor; dividend >= c; c <<= 1, k++)  
  18.         {  
  19.             if (dividend - c < divisor)  
  20.             {  
  21.                 res += 1<<k;  
  22.                 break;  
  23.             }          
  24.         }  
  25.         if (dividend - c < divisor)  
  26.             break;  
  27.   
  28.         res += 1<<(k - 1);  
  29.         dividend -= c>>1;  
  30.     }  
  31.   
  32.     return res;  
  33. }  
//非递归算法 
int integer_div_3(unsigned int dividend, unsigned int divisor)
{
    if(divisor == 0)
    {
        cout << "除数不能为0" << endl;
        exit(1);
    }

    if (dividend < divisor) 
        return 0;
    
    unsigned int k, c, res=0;

    while (dividend > divisor)
    {
        for (k = 0,c = divisor; dividend >= c; c <<= 1, k++)
        {
            if (dividend - c < divisor)
            {
                res += 1<<k;
                break;
            }        
        }
        if (dividend - c < divisor)
            break;

        res += 1<<(k - 1);
        dividend -= c>>1;
    }

    return res;
}

方法二:

另一种方法是使用魔数进行运算.

如果你看过编译器生成代码优化或一些大型系统的源码,或者嵌入式代码优化等等性能要求非常高的代码,你就会发现一些称之为magic number的数.

如0xAAAAAAAB、0x66666667、0x24924925、0x51EB851F、0x10624DD3......这些分别是快速除法(3,5,7。。。)的魔幻数字
编译器一般用这个来优化除法,至于怎么得来的自己google吧.

例如求一个数除以3的值可以这样求:

  1. int divby3(int x)  
  2. {  
  3.     return ((__int64)x*0xAAAAAAABUL) >> 33 ;  
  4. }  
int divby3(int x)
{
    return ((__int64)x*0xAAAAAAABUL) >> 33 ;
}

这里推荐一本书,里面有讲关于魔数的算法. <<Hacker's Delight>> By Henry S. Warren, Publisher : Addison Wesley

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值