这段代码来自https://github.com/erlang/otp/blob/master/erts/emulator/sys/common/erl_mseg.c
static const int debruijn[32] = {
0, 1, 28, 2, 29, 14, 24, 3, 30, 22, 20, 15, 25, 17, 4, 8,
31, 27, 13, 23, 21, 19, 16, 7, 26, 12, 18, 6, 11, 5, 10, 9
};
#define LOW1BIT(X) (debruijn[((unsigned int)(((X) & -(X)) * 0x077CB531U)) >> 27])
分析:
首先计算(x) & (-x),根据反码计算规则-x就是x按位取反结果再加1,假设为1的最低位是第k位,则-x中的第31~(k+1)位与x中的31~(k+1)位每一位都是相反的,-x和x中的(k-1)~0位都是0,而两个的第k位都是1,因此x和-x做与运算会得到2^k。
此程序比较巧妙的位置是乘和右移,我猜作者的本意是想找到一个映射关系,使得集合A中的元素x=2^i(其中i=0,1,...31),与集合B中的元素(0,1,...31)形成双射关系,能满足这个条件的映射会有很多种,而此程序中找到的映射关系应该是其中比较简洁的一种。0x077CB531U没有什么特别之处,通过下面的程序可以找到4096个能替代它的值,但是可能需要修改debruijn数组。
#include <stdio.h>
int main()
{
unsigned int k = 1 ,t = 1 ,s;
int i;
unsigned int bitmap[32];
unsigned int pow2[32];
memset(bitmap ,0 ,sizeof(bitmap));
for(i = 0 ;i < 32 ;i ++)
pow2[i] = 1<<i;
while(k != 0xFFFFFFFF) {
for(i = 31 ; i >= 0 ;i --) {
s = (pow2[i] * k) >> 27;
if(bitmap[s] < t) {
bitmap[s] = t;
} else {
break;
}
}
if(i < 0) {
printf("0x%XU\n" ,k);
// break;
}
++ k;
++ t;
}
return 0;
}