该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
考虑我们那 32-bit unsigned int 的性质,它做加法、乘法只会保留低 32-bit 的结果。
实际上,它是模 2^32 的剩余系。我们可以找到这样一个数 p 使得 3*p = 1 (mod 2^32)
p 可以取 0xAAAAAAAB,3*p = 0x200000001 = 2^33 + 1
对于给出的数 a,只要乘上 p(0xAAAAAAAB),如果 a 能被3整除,那么 a*p*3 = a (mod 2^32)
如果 a 能被3整除,那么 a*p mod 2^32 就是所求的商
如果 a 不能被3整除呢?a*p*3 = a (mod 2^32) 仍然成立。如何分辨呢?
设 a = 3q + r(保证 0 <= r
ap = (2^33+1)q + pr,其中 pr 最大为 2*0xAAAAAAAB
如果计算 ap 时,保留 64-bit 结果,那么右移33位就得到了 q,q 就是商,余数 r 可以是0~2
计算 q*3 就可以验证 a 是否是3的倍数
gcc/g++ 32位无符号数除以3的优化实现用到了以上方法
其实,gcc/g++ 对很多常数除法做了优化……
另外,(a+r)*p(规定 0 <= r
-----
考虑数 4*a+b 和 a+b 模3同余,所以:c 和 (c >> 2) + (c & 3) 模3同余
对于 c >= 4,(c >> 2) + (c & 3) 至少比 c 小1
反复令 c = (c >> 2) + (c & 3),当 c <= 3 时,判断 c 是否为0或3,是则原始的 c 被3整除
-----
设 c = 2*a+b = 3*a+b-a,那么如果 b-a 和 c 同余
设 c = 4*a+2*b+d = 3*a+3*b + a-b+d,a-b+d 和 c 同余
设 c = 8a+4b+2d+e = 9a+3b+3d -a +b -c +d,-a +b -c +d 和 c 同余
……
二进制表示中,偶数位和 减去 奇数位和,与原数同余