位運算之美——用+,-和位運算實現(xiàn)整數(shù)除法和取模(一)

??? 今天看了一位師兄去年的筆經(jīng)總結(jié)，其中有一題是“不許用%和/來實現(xiàn)求任意數(shù)除以3的余數(shù)”，我想考官的目的應(yīng)該是想考察學生對位運算的熟悉程度吧，于是我把題目擴展成“只能用+,-和位運算實現(xiàn)整數(shù)除法(/)和取模(%)”，注意：這里不能使用其它的庫例程來輔助計算，如log,log10等。在思考這道題目的過程中，我又涉及到了許多二進制相關(guān)的題目，如：
??? 判斷給定的整數(shù)是不是2的整數(shù)次冪
??? 判斷給定的整數(shù)是不是4的整數(shù)次冪
??? 求給定整數(shù)的二進制表示中1的個數(shù)
??? 求給定整數(shù)的二進制表示中0的個數(shù)
??? 求給定整數(shù)的二進制表示中最高位1的位置
??? 求大于等于給定整數(shù)的最小的2的整數(shù)次冪
??? 求給定整數(shù)的二進制表示的有效位數(shù)
??? ...
??? 這些題目都是經(jīng)典老題，頻繁出現(xiàn)于各類筆試面試題中，除了能考察位運算外，還能考察應(yīng)聘者能否給出創(chuàng)新的算法來更好地解決問題?？梢哉f這些題目都不難，如果使用32位的int來表示整數(shù)的話，蠻力法都可以比較好地完成任務(wù)，但是如果想盡可能地提高效率，那就需要動一番腦經(jīng)了。下面給出我對這些問題的整理和C++實現(xiàn)，并在此基礎(chǔ)上給出原題(只能用+,-和位運算實現(xiàn)整數(shù)除法(/)和取模(%)，下文都稱為原題)的實現(xiàn)。
??? 當然，從某種意義上講，特別是從充分利用底層硬件的計算能力(利用特殊的cpu指令)來看，這些解法肯定不是最優(yōu)的，希望大俠們多多指點。
??? 還要說明的是，下面各題的順序是按照我在思考原題時的思維過程來安排的，在給出原題的實現(xiàn)時會詳細說明。
??? 判斷給定的整數(shù)是不是2的整數(shù)次冪
??? 這應(yīng)該是最簡單的，利用最高位是1，其后所有位為0的特性，常數(shù)時間解決問題：

????求給定整數(shù)的二進制表示中1的個數(shù)
????考慮到n-1會把n的二進制表示中最低位的1置0并把其后的所有0置1，同時不改變此位置前的所有位，那么n&(n-1)即可消除這個最低位的1。這樣便有了比順序枚舉所有位更快的算法：循環(huán)消除最低位的1，循環(huán)次數(shù)即所求1的個數(shù)。此算法的時間復(fù)雜度為O(n的二進制表示中的1的個數(shù))，最壞情況下的復(fù)雜度O(n的二進制表示的總位數(shù))。
??? 既然有了求給定整數(shù)的二進制表示中1的個數(shù)的辦法，那么想要求給定整數(shù)的二進制表示中0的個數(shù)就很簡單了。事實上，在二進制中，完全可以把0和1看作是對稱的兩個對象，取反操作(~)可以任意的切換這兩個對象，只要先對n進行一次取反，然后再用上述算法即能得到二進制表示中0的個數(shù)。首先看下面的代碼：
??? 不知大家有沒有看出問題來？是的~操作符會把所有高位的都取反，而不是只把有效位取反，所以我們需要一個能保持高位不變的位取反操作，下面是我的實現(xiàn)，時間復(fù)雜度和求二進制表示中1的個數(shù)的算法相同，都與二進制表示中1的個數(shù)有關(guān)：
??? 有了這個特殊的取反操作，求給定整數(shù)的二進制表示中0的個數(shù)的辦法就簡單了：
??? 看到這里，聰明的讀者肯定看出問題來了，其實我干了一件很蠢的事情?？纯瓷鲜鏊惴ǖ臅r間復(fù)雜度，negate_bits花了O(n的二進制表示中1的個數(shù))，while循環(huán)計算取反后的n的二進制表示中1的個數(shù)，事實上就是O(n的二進制表示中0的個數(shù))，兩部分加起來其實就是二進制表示總的有效位數(shù)，換句話說，這個算法是線性的，而事實上，我們完全可以先線性地求出這個總的有效位數(shù)，然后減去1的位數(shù)，即得到0的位數(shù)，根本不用費那么大勁去整個保持高位的取反操作，兩者的時間復(fù)雜度在漸進意義上也是相同的。所以，我犯傻了，但是這里又引出另一個問題：

????求給定整數(shù)的二進制表示的有效位數(shù)
??? 上面提到了線性地求這個位數(shù)(下文記為m)，即“循環(huán)右移1位，記錄右移次數(shù)”，時間復(fù)雜度O(m)。但是我想，一看到這個題目，所有人的第一反應(yīng)應(yīng)該是floor(log₂(n))+1吧，但是請注意，本文在一開始就規(guī)定了“不能使用庫例程”，那么在這個限制下該怎么做呢？有沒有比線性時間更好的算法呢？其實到目前為止我也沒有什么特別好的算法，希望誰有什么精妙的算法能指點一下，不要打我。。。

?1//求給定整數(shù)的二進制表示的位數(shù)，線性算法
?2int?count_bit_1(unsigned int?n)
?3{
?4????int?r?=?0;
?5????while(n)
?6????{
?7????????n>>=1;
?8????????r++;
?9????}
10????return?r;
11}

:???求大于等于給定整數(shù)的最小的2的整數(shù)次冪
????首先是最簡單的思路：求出n的二進制表示的總位數(shù)m，于是1<<m即為所求值，當然這里要排除n自身就是2的整數(shù)次冪的情況，復(fù)雜度O(m)，實現(xiàn)如下：
??? 事實上這就涉及到上面求二進制表示位數(shù)的問題，所以目前為止在此基礎(chǔ)上的算法都是線性時間的。???
??? 那有沒有不用計算位數(shù)m，從而效率更好的算法呢，能不能像在計算二進制表示中1的個數(shù)時那樣根據(jù)1的個數(shù)來設(shè)計算法呢？回到那一題中，“n-1會把n的二進制表示中最低位的1置0并把其后的所有0置1”，那么n|=n-1就把n的二進制表示中最低位1后的所有0置1，再加上1，那么就把最低位1左移了一位。于是，便有了更好的算法：循環(huán)左移最低位的1，直到n是2的整數(shù)次冪。該算法跟二進制表示中的1個數(shù)和位置有關(guān)，最壞時間復(fù)雜度還是O(二進制表示位數(shù))，但是比起上一個實現(xiàn)，這個算法在多數(shù)情況下都比上一個算法快。實現(xiàn)如下：
????
??? 最后來一個簡單的擴充題目：
??? 判斷給定的整數(shù)是不是4的整數(shù)次冪
????觀察4的整數(shù)次冪的特征，容易發(fā)現(xiàn)除了滿足n&(n-1)==0外，唯一的1位后的0的個數(shù)是偶數(shù)，這從4^x=2^2k也能簡單地得到。這就很直觀地衍生出一個簡單的算法：
????算法很直觀，但是比起is_2exp的常數(shù)時間is_4exp的線性時間總讓我覺得不能接受，不過無奈還是沒有想出好辦法來，哎。。。求大牛指點啊

??? 說明：寫這篇文章，已經(jīng)三次丟失全文了，把我快搞瘋了，firefox下好像有點問題，先把文章發(fā)上來，過會兒換到IE下繼續(xù)。。。
??? 再說明：換了IE后就沒再出問題了，不過寫著寫著發(fā)現(xiàn)寫了好久，先歇會兒，得看書補習功課了
??? 最后的說明：下次會基于上面的內(nèi)容，給本文最初提出的問題(只能用+,-和位運算實現(xiàn)整數(shù)除法(/)和取模(%))的實現(xiàn)