计算机酷学（七）

数字（五）

编译器如何优化乘法操作

      在大多数机器上，整数乘法指令比其他整数运算(如加减、位级和移位)慢，因此编译器通过使用移位和加减运算的组合来代替乘以常熟因子的乘法。
      首先我们从乘以2的幂出发。假设0 ≤ k < ω，x << k可以表示x × 2^k，但是会有溢出的可能性。比方说一个数14，它可以拆分成2³ + 2² + 2¹，也可以拆分成2⁴ - 2¹，因此x × 14可以写成(x << 3) + (x << 2) + (x << 1)，也可以写成(x << 4) - (x << 1)。
      对于一组从位位置n到位位置m的连续的1(n≥m)，可以有两种形式，(x << n) + (x << (n - 1)) + … + (x << m)和(x << (n + 1)) - (x << m)。乘以其他值可以根据这两种形式去推，这里就不多说了。

编译器如何优化除法操作

      在C语言中，整数除法总是舍入到零的。先讲一下除以2的幂的无符号除法，这可以看成x >> k，而且舍入到零，这没有问题。但是除以2的幂的有符号除法会遇到一个问题，就是被除数小于0时，仅通过移位会向下舍入，而不是向零舍入。
      那么要怎么编译器要怎么解决这个问题呢？首先编译器在移位前给这个值进行偏置，偏置的值的大小为2^k-1.为什么是这个值呢？
      首先设x=qy+r，其中0≤r<y。(x + y - 1) / y = q + (r + y - 1) / y。如果r = 0，那么这个结果就是q；如果r != 0，那么这个结果就是q + 1。因此该除法可以写成(x + (1 << k) - 1) >> k。
      C语言代码相当于(x < 0 ? x + (1 << k) - 1 : x) >> k。但是这不能推广到除以任意整数值，太可惜了。

写一个函数div1024

      要求是不使用除法、模运算、乘法、任何条件语句、任何比较运算符或任何循环，对于整数x返回x/1024的值。

int div1024(int x){
	int bias = (x >> 31) & 0x3FF;	// 确定x的正负，如果为正，则为0，如果为负，则为1023
	return (x + bias) >> 10;
}