位运算之美——用+,-和位运算实现正整数除法和取模(二)

作者：翼帆@cppblog 
原文地址：http://www.cppblog.com/xiaoyisnail/archive/2009/09/21/96883.html
本文版权归作者和cppblog共有，欢迎转载，但未经作者同意必须保留此段声明，且在文章页面明显位置给出原文连接，否则保留追究法律责任的权利。

终于有时间写本文的第二部分了。在上一篇文章(下文中称为“文(一)”)中，我提出了一个具体的问题“只能用+,-和位运算实现正整数除法(/)和取模(%)”，并整理了一些和位运算相关的题目和算法，本文将给出上述问题的一个完整的解答思路和实现。

首先思考最简单的除法实现，即循环减除数，减到不能再减为止，所减次数即所求的商，事实上这就是我们初学四则运算时对除法的定义，实现：

 1 int  integer_div_0( int  dividend,  int  divisor)
 2 {
 3    if(divisor == 0)
 4    {
 5        cout<<"非法参数，除零错"<<endl;
 6        exit(1);
 7    }
 8
 9    int res = 0;
10    while((dividend-=divisor)>=0)
11        ++res;
12
13    return res;
14}

很显然，这个简单的实现是非常低效的，求integer_div_0(a,b)的时间复杂度为O(a/b)，当a很大b很小时，计算开销很大，有什么办法能提高效率吗？考虑到循环减除数同时比较差是否大于等于0有点类似于在一个一个数组或是linked list中顺序搜索某个目标值，在较差的情况下近似于搜索整个问题空间，这样的开销必然很大，那么只要找到某个方法来减少需要搜索的问题空间就可以提高性能了。于是，自然可以想到如下算法：以除数为初始测试值，以2的指数为步长来搜索问题空间，当被除数与测试值的差小于除数时便结束搜索，若在这之前测试值大于被除数，则将被除数减去前一个测试值，并重复上述过程直到搜索结束。举个例子，求1200/3：
顺序搜索时，我们要与1200比较的数有：3,6,9,12,15,...,1998,2001，比较次数667次
以2的指数为步长搜索时，与1200比较3,6,12,24,48,96,192,384,768,1536，然后与1200-768=432再进行比较3,6,12,24,48,96,192,384,768，再取432-384=48比较3,6,12,24,48，搜索结束，比较次数共24次，比顺序搜索有很大的提高。你可能会问，为什么要以2的指数为步长来搜索呢？答案是，这样我们就可以使用位操作来进一步提高计算效率了。下面是这个算法的实现：

 1 int  integer_div_1(unsigned  int  dividend, unsigned  int  divisor)
 2 {
 3    if(divisor == 0)
 4    {
 5        cout<<"非法参数，除零错"<<endl;
 6        exit(1);
 7    }
 8
 9    if(dividend < divisor) return 0;
10    unsigned int k=0,c=divisor, res=0;
11    
12    for(;dividend>=c;c<<=1,k++)
13        if(dividend-c < divisor)
14            return 1<<k;
15
16    return integer_div_1(dividend-(c>>1), divisor)+(1<<(k-1));
17}

注意到最后一行的尾递归，再把代码优化为非递归如下：

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

最后，有了整数除法，取模运算就很简单了，从进行整数除法搜索商的最后一步立刻就能得到模除的余数，实现如下，为了方便起见，代码里使用C++ STL中的pair模板以同时返回商和余数：

 1 // 整数除法and取模，返回商和余数
 2 pair < int , int >  integer_div_3(unsigned  int  dividend, unsigned  int  divisor)
 3 {
 4    if(divisor == 0)
 5    {
 6        cout<<"非法参数，除零错"<<endl;
 7        exit(1);
 8    }
 9
10    if(dividend < divisor) 
11        return make_pair(0, dividend);
12    unsigned int k, c, quotient=0, remainder;
13
14    while(dividend > divisor)
15    {
16        for(k=0,c=divisor;dividend>=c;c<<=1,k++)
17        {
18            if(dividend-c < divisor)
19            {
20                quotient += 1<<k;
21                remainder = dividend-c;
22                break;
23            }        
24        }
25        if(dividend-c < divisor)
26            break;
27
28        quotient += 1<<(k-1);
29        dividend -= c>>1;
30    }
31
32    return make_pair(quotient, remainder);
33}

 

 

附上代码：http://www.cppblog.com/Files/xiaoyisnail/bits.rar