深入理解计算机系统的数值类型

前言

在日常编程中,数值类型numeric types)是我们打交道最多的类型,可能没有之一。除了最熟悉的 int,还有 longfloatdouble 等。正因太熟悉,我们往往不会深究它们的底层原理。因为平时的工作中,知道个大概,也够用了。

但,在某些业务场景下,比如金融业务,数值运算不准确会带来灾难性的后果。这时,你就必须清楚数值类型的二进制表示、截断、转型等原理,否则很难保证运算结果的正确性。

另外,数值类型也是一个容易被黑客攻击的点,考虑如下一段代码:

// C++
/* Declaration of library function memcpy */
void *memcpy(void *dest, void *src, size_t n);
/* Kernel memory region holding user-accessible data */
#define KSIZE 1024
char kbuf[KSIZE];
/* Copy at most maxlen bytes from kernel region to user buffer */
int copy_from_kernel(void *user_dest, int maxlen) {
    /* Byte count len is minimum of buffer size and maxlen */
    int len = KSIZE < maxlen ? KSIZE : maxlen;
    memcpy(user_dest, kbuf, len);
    return len;
}

如果你熟悉数值类型的原理,一定会敏锐察觉出第 10 行存在 intsize_t 的类型转换。在 64 位系统中,size_t 通常被定义为 unsigned long 类型,如果攻击者在调用 copy_from_kernel 时,特意传入一个负数的 maxlen,转型到 memcpy 中的 n 将会是一个很大的正数,从而导致了内存拷贝的越界!

数值类型是计算机编程的基础,用的很多,也很重要,理解它的底层原理,有助于写出正确的代码,避免一些意料之外的错误

每个计算机系统都有固定的 word size,也即常说的 xx 位,它也是指针的大小,跟 虚拟内存 相关,比如一个 w 位系统上的应用程序,最多能够访问 2 w 2^{w} 2w byte 大小的虚拟内存。

最常用的是 32 位 和 64 位 系统,某些数值类型在它们之上会有些差异,比如 long 类型 在 32 位系统上是 32 bit 大小,在 64 位系统上是 64 bit 大小。考虑如今 64 位系统逐渐成为主流,本文会以它作为基础,进行数值类型的介绍

整数

在计算机系统中,整数可以分成 无符号unsigned)整数 和 有符号signed)整数 两大类,这之下,按照类型表示的 bit 位大小,又可细分成 8 位的 char/byte/int8 、16 位的 short/innt16、32 位的 int/int32 和 64 位的 long/int64,它们的取值范围如下:

类型最小值最大值
[signed] char-128127
unsigned char0255
short-32,76832,767
unsigned short065,535
int−2,147,483,6482,147,483,647
unsigned int04,294,967,295
long−9,223,372,036,854,775,8089,223,372,036,854,775,807
unsigned long018,446,744,073,709,551,615

死记这个表不容易,下面我们将试图从二进制编码层面去理解它。

二进制编码

整数在计算机系统上都是以二进制存储的,对于一个 w 位的整数 x x x,它的二进制表示写成这样:
[ x w − 1 , x w − 2 , . . . , x 0 ] b [x_{w-1}, x_{w-2}, ..., x_0]_b [xw1,xw2,...,x0]b
其中, x i x_i xi 取值 0 0 0 1 1 1

无符号编码(Unsigned Encodings)

在二进制表示的基础上,无符号编码 是这样:
x = ∑ i = 0 w − 1 x i ⋅ 2 i x = \sum_{i=0}^{w-1}x_i \cdot 2^i x=i=0w1xi2i
比如,w = 4 场景下的一些例子:
[ 0001 ] b = 0 ⋅ 2 3 + 0 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 = 0 + 0 + 0 + 1 = 1 [0001]_b=0\cdot2^3+0\cdot2^2+0\cdot2^1+1\cdot2^0=0+0+0+1=1 [0001]b=023+022+021+120=0+0+0+1=1

[ 0101 ] b = 0 ⋅ 2 3 + 1 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 = 0 + 4 + 0 + 1 = 5 [0101]_b=0\cdot2^3+1\cdot2^2+0\cdot2^1+1\cdot2^0=0+4+0+1=5 [0101]b=023+122+021+120=0+4+0+1=5

[ 1011 ] b = 1 ⋅ 2 3 + 0 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 = 8 + 0 + 2 + 1 = 11 [1011]_b=1\cdot2^3+0\cdot2^2+1\cdot2^1+1\cdot2^0=8+0+2+1=11 [1011]b=123+022+121+120=8+0+2+1=11

[ 1111 ] b = 1 ⋅ 2 3 + 1 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 = 8 + 4 + 2 + 1 = 15 [1111]_b=1\cdot2^3+1\cdot2^2+1\cdot2^1+1\cdot2^0=8+4+2+1=15 [1111]b=123+122+121+120=8+4+2+1=15

由上述可知,无符号编码无法表示负数,因此只能表示无符号整数。为了表示有符号整数,还要探寻另一种编码方式。

原码编码(True Form Encodings)

为了区分正数和负数,很容易想到使用一个 bit 位作为符号位 0 0 0 表示正数, 1 1 1 表示负数。在无符号编码的基础上,使用最高位作为符号位,其他位含义不变,得出 原码编码 形式:
x = ( − 1 ) x w − 1 ⋅ ∑ i = 0 w − 2 x i 2 i x = (-1)^{x_{w-1}} \cdot \sum_{i=0}^{w-2}x_i2^i x=(1)xw1i=0w2xi2i
比如,w = 4 场景下的一些例子:
[ 0001 ] b = ( − 1 ) 0 ⋅ ( 0 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 ) = 1 ⋅ ( 0 + 0 + 1 ) = 1 [0001]_b=(-1)^0\cdot(0\cdot2^2+0\cdot2^1+1\cdot2^0)=1\cdot(0+0+1)=1 [0001]b=(1)0(022+021+120)=1(0+0+1)=1

[ 0101 ] b = ( − 1 ) 0 ⋅ ( 1 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 ) = 1 ⋅ ( 4 + 0 + 1 ) = 5 [0101]_b=(-1)^0\cdot(1\cdot2^2+0\cdot2^1+1\cdot2^0)=1\cdot(4+0+1)=5 [0101]b=(1)0(122+021+120)=1(4+0+1)=5

[ 1011 ] b = ( − 1 ) 1 ⋅ ( 0 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 ) = − 1 ⋅ ( 0 + 2 + 1 ) = − 3 [1011]_b=(-1)^1\cdot(0\cdot2^2+1\cdot2^1+1\cdot2^0)=-1\cdot(0+2+1)=-3 [1011]b=(1)1(022+121+120)=1(0+2+1)=3

[ 1111 ] b = ( − 1 ) 1 ⋅ ( 1 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 ) = − 1 ⋅ ( 4 + 2 + 1 ) = − 7 [1111]_b=(-1)^1\cdot(1\cdot2^2+1\cdot2^1+1\cdot2^0)=-1\cdot(4+2+1)=-7 [1111]b=(1)1(122+121+120)=1(4+2+1)=7

虽然原码编码方式简单直观,但它还存在两个问题:

(1) 0 0 0 存在两种编码形式

原码编码方式下, 0 0 0 存在两种编码形式, [ 0 , 0 , . . . , 0 ] b [0,0,...,0]_b [0,0,...,0]b [ 1 , 0 , . . . , 0 ] b [1,0,...,0]_b [1,0,...,0]b。同一个整数值,却有两种编码,这对计算机系统来说没什么意义,反而是一种浪费。

(2)带负数的加法运算不正确

原码编码方式下,两个正数的加法没问题,一旦带上负数,结果就出错了:

所以,原码编码方式,注定不会被使用。

补码编码(Two’s-complement Encodings)

于是,补码编码 被发明,它也是建立在无符号编码的基础上,仍然取最高位为符号位,编码方式是这样:
x = − x w − 1 ⋅ 2 w − 1 + ∑ i = 0 w − 2 x i 2 i x = -x_{w-1} \cdot 2^{w-1} + \sum_{i=0}^{w-2}x_i2^i x=xw12w1+i=0w2xi2i
它与无符号编码的唯一区别是,最高位的取值从 x w − 1 ⋅ 2 w − 1 x_{w-1} \cdot 2^{w-1} xw12w1 变成了 − x w − 1 ⋅ 2 w − 1 -x_{w-1} \cdot 2^{w-1} xw12w1

比如,w = 4 场景下的一些例子:
[ 0001 ] b = − 0 ⋅ 2 3 + 0 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 ) = 0 + 0 + 0 + 1 = 1 [0001]_b=-0\cdot2^3 +0\cdot2^2+0\cdot2^1+1\cdot2^0)=0+0+0+1=1 [0001]b=023+022+021+120)=0+0+0+1=1

[ 0101 ] b = − 0 ⋅ 2 3 + 1 ⋅ 2 2 + 0 ⋅ 2 1 + 1 ⋅ 2 0 = 0 + 4 + 0 + 1 ) = 5 [0101]_b=-0\cdot2^3+1\cdot2^2+0\cdot2^1+1\cdot2^0=0+4+0+1)=5 [0101]b=023+122+021+120=0+4+0+1)=5

[ 1011 ] b = − 1 ⋅ 2 3 + 0 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 = − 8 + 0 + 2 + 1 = − 5 [1011]_b=-1\cdot2^3+0\cdot2^2+1\cdot2^1+1\cdot2^0=-8+0+2+1=-5 [1011]b=123+022+121+120=8+0+2+1=5

[ 1111 ] b = − 1 ⋅ 2 3 + 1 ⋅ 2 2 + 1 ⋅ 2 1 + 1 ⋅ 2 0 = − 8 + 4 + 2 + 1 = − 1 [1111]_b=-1\cdot2^3+1\cdot2^2+1\cdot2^1+1\cdot2^0=-8+4+2+1=-1 [1111]b=123+122+121+120=8+4+2+1=1

补码编码很巧妙地解决了原码编码的两个问题:

首先,0 在补码编码下只有一种编码形式, [ 0 , 0 , . . . , 0 ] b [0,0,...,0]_b [0,0,...,0]b

此外,带负数的加法运算,也正确了。

因为补码编码的简单和正确性,目前,几乎所有的计算机系统,都采用补码编码来表示有符号整数

位运算

位运算主要包含 取反异或移位 等几种,我们在业务开发时用得比较少,但如果你有阅读开源代码的习惯,就会经常发现它们的踪迹。如果碰巧对位运算不熟悉,那么阅读这些代码,就同读天书一般。

取反(~)、与(&)、或(|)、异或(^)的规则比较简单:

移位运算,可以分成 左移右移 两种,其中,右移又可分为 逻辑右移算术右移

左移(<<)运算,是对二进制整数,向左移若干位,高位丢弃,低位补零。也即,对 [ x w − 1 , x w − 2 , . . . , x 0 ] b [x_{w-1}, x_{w-2}, ..., x_0]_b [xw1,xw2,...,x0]b 左移 k k k 位,得到 [ x w − k − 1 , x w − k − 2 , . . . , x 0 , 0 , . . . , 0 ] b [x_{w-k-1}, x_{w-k-2}, ..., x_0, 0, ..., 0]_b [xwk1,xwk2,...,x0,0,...,0]b

比如,对 int i = -1 左移 10 位,会得到 i = -1024 的结果:

// Java语言
public static void main(String[] args) {
    int i = -1;
    System.out.println("Before << , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
    i <<= 10;
    System.out.println("After << , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
}
// 输出结果:
Before << , i's value is -1
i's binary string is 11111111111111111111111111111111
After << , i's value is -1024
i's binary string is 11111111111111111111110000000000

在 C/C++ 中,两种右移操作符都是 >>,对无符号整数用的是逻辑右移,对有符号整数用的是算术右移;在 Java 中,逻辑右移的操作符是 >>>,算术右移的操作符是 >>。为了方便区分,下文统一用 Java 的表示方法。

逻辑右移(>>>)运算,是对二进制整数,向右移若干位,高位补零,低位丢弃。也即,对 [ x w − 1 , x w − 2 , . . . , x 0 ] b [x_{w-1}, x_{w-2}, ..., x_0]_b [xw1,xw2,...,x0]b 逻辑左移 k 位,得到 [ 0 , . . . , 0 , x w − 1 , x w − 2 , . . . , x k ] b [0, ..., 0,x_{w-1}, x_{w-2}, ..., x_k]_b [0,...,0,xw1,xw2,...,xk]b

比如,对 int i = -1 逻辑右移 10 位,会得到 i = 4194303 的结果:

// Java语言
public static void main(String[] args) {
    int i = -1;
    System.out.println("Before >>> , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
    i >>>= 10;
    System.out.println("After >>> , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
}
// 输出结果:
Before >>> , i's value is -1
i's binary string is 11111111111111111111111111111111
After >>> , i's value is 4194303
i's binary string is 1111111111111111111111

算术右移(>>)运算,是对二进制整数,向右移若干位,高位补符号位,低位丢弃。也即,对 [ x w − 1 , x w − 2 , . . . , x 0 ] b [x_{w-1}, x_{w-2}, ..., x_0]_b [xw1,xw2,...,x0]b 逻辑左移 k 位,得到 [ x w − 1 , . . . , x w − 1 , x w − 1 , x w − 2 , . . . , x k ] b [x_{w-1}, ..., x_{w-1},x_{w-1}, x_{w-2}, ..., x_k]_b [xw1,...,xw1,xw1,xw2,...,xk]b

比如,对 int i = -1 算术右移 10 位,仍会得到 i = -1 的结果:

// Java语言
public static void main(String[] args) {
    int i = -1;
    System.out.println("Before >> , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
    i >>= 10;
    System.out.println("After >> , i's value is " + i);
    System.out.println("i's binary string is " + Integer.toBinaryString(i));
}
// 输出结果:
Before >> , i's value is -1
i's binary string is 11111111111111111111111111111111
After >> , i's value is -1
i's binary string is 11111111111111111111111111111111

目前为止,介绍移位运算的原理时,我们都默认 k < w,如果 k >= w 会怎样

比如, [ x w − 1 , x w − 2 , . . . , x 0 ] b [x_{w-1}, x_{w-2}, ..., x_0]_b [xw1,xw2,...,x0]b 左移 w 位,结果会是 [ 0 , 0 , . . . , 0 ] b [0, 0, ..., 0]_b [0,0,...,0]b 吗:

// Java语言
public static void main(String[] args) {
    int i1 = -1;
    System.out.println("Before << 31, i1's value is " + i1);
    System.out.println("i1's binary string is " + Integer.toBinaryString(i1));
    i1 <<= 31;
    System.out.println("After << 31, i1's value is " + i1);
    System.out.println("i1's binary string is " + Integer.toBinaryString(i1));

    int i2 = -1;
    System.out.println("Before << 32, i2's value is " + i2);
    System.out.println("i2's binary string is " + Integer.toBinaryString(i2));
    i2 <<= 32;
    System.out.println("After << 32, i2's value is " + i2);
    System.out.println("i2's binary string is " + Integer.toBinaryString(i2));
}
// 输出结果:
Before << 31, i1's value is -1
i1's binary string is 11111111111111111111111111111111
After << 31, i1's value is -2147483648
i1's binary string is 10000000000000000000000000000000
Before << 32, i2's value is -1
i2's binary string is 11111111111111111111111111111111
After << 32, i2's value is -1
i2's binary string is 11111111111111111111111111111111

上述例子中, w = 32,我们发现 k = 31 时,结果还符合预期;当 k = 32 时,结果不是 0,而是 -1,也即相当于 k = 0 时的结果。

原因是这样,w 位整数 x,当执行 x << k 时,实际执行的是 x << (k % w)。所以,当 i2 << 32 时,实际是 i2 << 32 % 32 = i2 << 0

右移操作也遵循同样的规则,也即 x >> k = x >> (k % w)x >>> k = x >>> (k % w)

算术运算

整数的算术运算也不总符合我们的常识,比如下面这个例子:

// Java语言
public static void main(String[] args) {
    int i1 = 2147483647;
    int i2 = 1;
    System.out.printf("2147483647+1=%d\n", i1+i2);
}
// 输出结果:
2147483647+1=-2147483648

两个正数 21474836471 相加,实际运行结果却是一个负数 -2147483648。原因是,运算结果超出了 int 类型的表示范围,我们把这种现象称作 溢出

接下来,我们将介绍整数类型常见的几种算术运算,包括当运算出现溢出时,计算机系统的处理方式。

加法

(1)无符号整数加法运算

对两个 w 位的无符号整数 x , y ∈ [ 0 , 2 w − 1 ] x, y \in[0, 2^w-1] x,y[0,2w1],有 x + y ∈ [ 0 , 2 w + 1 − 2 ] x+y\in[0, 2^{w+1}-2] x+y[0,2w+12];如果仍用一个 w 位的整数表示相加结果,那么当结果在 [ 2 w , 2 w + 1 − 2 ] [2^w, 2^{w+1}-2] [2w,2w+12] 范围时,也即结果需要 w + 1 位表示时,运算就产生溢出了。

如果出现溢出,系统会对结果进行截断,只保留低 w 位

比如,w = 4 场景下的一些例子:

因此,如果用 + w u +_w^{u} +wu 表示 w 位无符号整数的加法运算,那么它有如下规则:
x + w u y = { x + y , x + y < 2 w   正常 x + y − 2 w , 2 w ≤ x + y < 2 w + 1   溢出 x+_w^{u}y= \begin{cases} x+y, &x+y<2^w\ &正常\\ x+y-2^w, &2^w\le x+y<2^{w+1}\ &溢出 \end{cases} x+wuy={x+y,x+y2w,x+y<2w 2wx+y<2w+1 正常溢出
其中,溢出场景因为对结果的进行截断,舍去了 w + 1 w+1 w+1 位,所以结果需要减去 2 w 2^w 2w

(2)有符号整数加法运算

对两个 w 位的有符号整数 x , y ∈ [ − 2 w − 1 , 2 w − 1 − 1 ] x, y \in[-2^{w-1}, 2^{w-1}-1] x,y[2w1,2w11],有 x + y ∈ [ − 2 w , 2 w − 2 ] x+y \in[-2^{w}, 2^{w}-2] x+y[2w,2w2];如果仍用一个 w 位的整数表示相加结果,那么当结果在 [ − 2 w , − 2 w − 1 − 1 ] [-2^{w}, -2^{w-1}-1] [2w,2w11] [ 2 w − 1 , 2 w − 2 ] [2^{w-1}, 2^{w}-2] [2w1,2w2] 范围时,运算就产生溢出了。

比如,w = 4 场景下的一些例子:

由上述例子可知,两个 w 位的负数相加,即使结果需要 w + 1 位表示,也可能是正常场景。这取决于截断位 w 位后的最高位,为 1 则还属于正常场景,为 0 则属于溢出场景

因此,如果用 + w t +_w^{t} +wt 表示 w 位有符号整数的加法运算,那么它有如下规则:
x + w t y = { x + y − 2 w , 2 w − 1 ≤ x + y   正溢出 x + y , − 2 w − 1 ≤ x + y < 2 w − 1   正常 x + y + 2 w , x + y < − 2 w − 1   负溢出 x+_w^{t}y= \begin{cases} x+y-2^w, &2^{w-1}\le x+y\ &正溢出\\ x+y, &-2^{w-1}\le x+y <2^{w-1}\ &正常\\ x+y+2^w, &x+y<-2^{w-1}\ &负溢出 \end{cases} x+wty= x+y2w,x+y,x+y+2w,2w1x+y 2w1x+y<2w1 x+y<2w1 正溢出正常负溢出

  • 正溢出场景下,原本 x + y = 2 w − 1 + ∑ i = 0 w − 2 x i ⋅ 2 i x+y=2^{w-1}+\sum_{i=0}^{w-2}x_i\cdot2^i x+y=2w1+i=0w2xi2i,但按照补码编码方式,实际变成了 x + w t y = − 2 w − 1 + ∑ i = 0 w − 2 x i ⋅ 2 i x+_w^{t}y=-2^{w-1}+\sum_{i=0}^{w-2}x_i\cdot2^i x+wty=2w1+i=0w2xi2i,所以就有 x + w t y = x + y − 2 w x+_w^{t}y=x+y-2^w x+wty=x+y2w
  • 负溢出场景下,w 位一定是 0,原本 x + y = − 2 w + ∑ i = 0 w − 2 x i ⋅ 2 i x+y=-2^{w}+\sum_{i=0}^{w-2}x_i\cdot2^i x+y=2w+i=0w2xi2i,结果截断之后,实际变成了 x + w t y = ∑ i = 0 w − 2 x i ⋅ 2 i x+_w^{t}y=\sum_{i=0}^{w-2}x_i\cdot2^i x+wty=i=0w2xi2i,所以就有 x + w t y = x + y + 2 w x+_w^{t}y=x+y+2^w x+wty=x+y+2w

为什么负溢出场景下,w 位一定是 0

负溢出场景出现的条件是 x + y < − 2 w − 1 x+y<-2^{w-1} x+y<2w1,也即 :
x + y = − 2 w + ∑ i = 0 w − 1 x i ⋅ 2 i < − 2 w − 1 = > ∑ i = 0 w − 1 x i ⋅ 2 i < 2 w − 2 w − 1 = 2 w − 1 = > x w − 1 ⋅ 2 w − 1 + ∑ i = 0 w − 2 x i ⋅ 2 i < 2 w − 1 x+y=-2^w+\sum_{i=0}^{w-1}x_i\cdot 2^i < -2^{w-1}\\ =>\sum_{i=0}^{w-1}x_i\cdot 2^i < 2^w-2^{w-1}=2^{w-1}\\ =>x_{w-1}\cdot 2^{w-1} + \sum_{i=0}^{w-2}x_i\cdot 2^i < 2^{w-1}\\ x+y=2w+i=0w1xi2i<2w1=>i=0w1xi2i<2w2w1=2w1=>xw12w1+i=0w2xi2i<2w1
如果要使上述不等式成立, x w − 1 x_{w-1} xw1 必然为 0.

前面很多公式,但记住这个就行,无符号整数和有符号整数的加法运算规则,本质都是一样的,可以拆解成 4 步

  1. 先计算出真实加法运算的结果值。
  2. 对结果值用 w + 1 位二进制表示。
  3. 再将 w + 1 位的二进制结果截断,保留低 w 位。
  4. 将截断后的 w 位二进制值转换回十进制整数,得到最终结果。无符号整数用无符号编码,有符号整数用补码编码。
取反

取反,也即求相反数,给定一个整数 x x x,那么它的相反数 y = − x y = -x y=x,也即满足 x + y = 0 x + y = 0 x+y=0

(1)无符号整数的取反

对一个 w 位的无符号整数 x ∈ [ 0 , 2 w − 1 ] x \in [0, 2^w-1] x[0,2w1],对它取反,也即找到一个整数 y ∈ [ 0 , 2 w − 1 ] y \in [0, 2^w-1] y[0,2w1],使得 x + y = 0 x + y = 0 x+y=0

  • x = 0 x = 0 x=0,那么很容易得出 y = 0 y = 0 y=0
  • 当 x > 0,只有在溢出场景才会存在 x + y = 0 x+y=0 x+y=0 的可能。由前文可知,无符号整数加法溢出场景下, x + w u y = x + y − 2 w = 0 x+_w^uy=x+y-2^w=0 x+wuy=x+y2w=0 时,有, y = 2 w − x y=2^w-x y=2wx

因此,如果用 − w u -_w^{u} wu 表示 w 位无符号整数的取反运算,那么它有如下规则:
− w u x = { x , x = 0 2 w − x , x > 0 -_w^{u}x= \begin{cases} x, &x=0\\ 2^w-x, &x>0 \end{cases} wux={x,2wx,x=0x>0
比如,w = 8 场景下的例子:

// C++
int main() {
    uint8_t i1 = 0;
    uint8_t i2 = -i1;
    printf("i1=%u\n", i1);
    printf("-i1=%u\n", i2);

    uint8_t i3 = 100;
    uint8_t i4 = -i3;
    printf("i3=%u\n", i3);
    printf("-i3=%u\n", i4);
    printf("2^8-i3=256-%u=%u\n", i3, 256-i3);
    return 0;
}
// 输出结果
i1=0
-i1=0
i3=100
-i3=156
2^8-i3=256-100=156

(2)有符号整数的取反

对一个 w 位的有符号整数 x ∈ [ − 2 w − 1 , 2 w − 1 − 1 ] x \in [-2^{w-1}, 2^{w-1}-1] x[2w1,2w11],对它取反,也即找到一个整数 y ∈ [ − 2 w − 1 , 2 w − 1 − 1 ] y \in [-2^{w-1}, 2^{w-1}-1] y[2w1,2w11],使得 x + y = 0 x + y = 0 x+y=0

  • x > − 2 w − 1 x > -2^{w-1} x>2w1,很容易得出 y = − x y = -x y=x,因为此时 − x ∈ [ − 2 w − 1 + 1 , 2 w − 1 − 1 ] -x\in [-2^{w-1}+1, 2^{w-1}-1] x[2w1+1,2w11] ,仍是有效的范围。
  • x = − 2 w − 1 x = -2^{w-1} x=2w1,因为 2 w − 1 2^{w-1} 2w1 已经不再有效范围内,所以只能是溢出场景。而 − 2 w − 1 − 2 w − 1 = − 2 w -2^{w-1}-2^{w-1}=-2^w 2w12w1=2w,截断为 w 之后,刚好为 0。所以,此时 − x = − 2 w − 1 -x = -2^{w-1} x=2w1

因此,如果用 − w t -_w^{t} wt 表示 w 位有符号整数的取反运算,那么它有如下规则:
− w t x = { − 2 w − 1 , x = − 2 w − 1 − x , x > − 2 w − 1 -_w^{t}x= \begin{cases} -2^{w-1}, &x=-2^{w-1}\\ -x, &x>-2^{w-1} \end{cases} wtx={2w1,x,x=2w1x>2w1
比如,w = 8 场景下的例子:

// C++
int main() {
    int8_t i1 = -128;
    int8_t i2 = -i1;
    printf("i1=%d\n", i1);
    printf("-i1=%d\n", i2);

    int8_t i3 = 100;
    int8_t i4 = -i3;
    printf("i3=%d\n", i3);
    printf("-i3=%d\n", i4);
    return 0;
}
// 输出结果
i1=-128
-i1=-128
i3=100
-i3=-100
乘法

前面介绍加法时说过,无符号整数和有符号整数的加法运算都可以拆成 4 步,这对乘法运算也适用。

对于无符号整数 x , y ∈ [ 0 , 2 w − 1 ] x,y \in [0, 2^w-1] x,y[0,2w1],那么 x ⋅ y ∈ [ 0 , ( 2 w − 1 ) 2 ] = [ 0 , 2 2 w − 2 w + 1 ] x\cdot y \in [0, (2^w-1)^2]=[0,2^{2w}-2^w+1] xy[0,(2w1)2]=[0,22w2w+1],即无符号整数的乘法运算结果最多需要 2w 位来表示。

比如,w = 8 时,无符号整数的例子:

// C++
int main() {
    uint8_t i1 = 100;
    uint8_t i2 = 2;
    uint8_t i3 = i1 * i2;
    printf("normal: i1 * i2 = %d * %d = %d\n", i1, i2, i3);

    uint8_t i4 = 100;
    uint8_t i5 = 3;
    uint8_t i6 = i4 * i5;
    printf("overflow: i4 * i5 = %d * %d = %d\n", i4, i5, i6);
    return 0;
}
// 输出结果
normal: i1 * i2 = 100 * 2 = 200
overflow: i4 * i5 = 100 * 3 = 44

对于有符号整数 x , y ∈ [ − 2 w − 1 , 2 w − 1 − 1 ] x,y \in [-2^{w-1}, 2^{w-1}-1] x,y[2w1,2w11],那么 x ⋅ y ∈ [ − 2 w − 1 ⋅ ( 2 w − 1 − 1 ) , ( − 2 w − 1 ) 2 ] = [ − 2 2 w − 2 + 2 w − 1 , 2 2 w − 2 ] x\cdot y \in [-2^{w-1}\cdot(2^{w-1}-1), (-2^{w-1})^2]=[-2^{2w-2}+2^{w-1},2^{2w-2}] xy[2w1(2w11),(2w1)2]=[22w2+2w1,22w2],即有符号整数的乘法运算结果最多也需要 2w 位来表示。

比如,w = 8 时,有符号整数的例子:

// C++
int main() {
    int8_t i1 = -50;
    int8_t i2 = 2;
    int8_t i3 = i1 * i2;
    printf("normal: i1 * i2 = %d * %d = %d\n", i1, i2, i3);

    int8_t i4 = -128;
    int8_t i5 = 127;
    int8_t i6 = i4 * i5;
    printf("overflow: i4 * i5 = %d * %d = %d\n", i4, i5, i6);
    return 0;
}
// 输出结果
normal: i1 * i2 = -50 * 2 = -100
overflow: i4 * i5 = -128 * 127 = -128

大部分机器上,乘法运算消耗 3 ~ 10 个 CPU 时钟周期,而加法运算和位运算只消耗 1 个时钟周期,所以,追求极致性能的程序都会想办法通过加法运算和位运算来替代乘法运算。

如果不考虑截断,对 x = [ x w − 1 , x w − 2 , . . . , x 0 ] b x=[x_{w-1}, x_{w-2}, ..., x_0]_b x=[xw1,xw2,...,x0]b 左移 k 位,会得到:
[ x w − 1 , x w − 2 , . . . , x 0 , 0 , . . . , 0 ] b = ∑ i = 0 w − 1 x i ⋅ 2 i + k = [ ∑ i = 0 w − 1 x i ⋅ 2 i ] ⋅ 2 k = x ⋅ 2 k [x_{w-1}, x_{w-2}, ..., x_0,0,...,0]_b=\sum_{i=0}^{w-1}x_i\cdot 2^{i+k}=[\sum_{i=0}^{w-1}x_i\cdot 2^i]\cdot 2^k=x\cdot 2^k [xw1,xw2,...,x0,0,...,0]b=i=0w1xi2i+k=[i=0w1xi2i]2k=x2k
也即,对 x x x 左移 k 位 相当于乘以 2 k 2^k 2k

如果考虑截断,就会存在溢出场景,对于 w 位的整数,左移 k 位效果也等同于 x ∗ w u 2 k x*_w^u2^k xwu2k x ∗ w t 2 k x*_w^t2^k xwt2k

比如,w = 8 时,无符号整数的例子:

// C++
int main() {
    uint8_t i1 = 100;
    uint8_t i2 = 4;
    uint8_t i3 = i1 * i2;
    printf("i1 * i2 = %d * %d = %d\n", i1, i2, i3);

    uint8_t k = 2;
    uint8_t i4 = i1 << k;
    printf("i1 << k = %d << %d = %d\n", i1, k, i4);
    return 0;
}
// 输出结果
i1 * i2 = 100 * 4 = 144
i1 << k = 100 << 2 = 144

有符号整数的例子:

// C++
int main() {
    int8_t i1 = -50;
    int8_t i2 = 4;
    int8_t i3 = i1 * i2;
    printf("i1 * i2 = %d * %d = %d\n", i1, i2, i3);

    int8_t k = 2;
    int8_t i4 = i1 << k;
    printf("i1 << k = %d << %d = %d\n", i1, k, i4);
    return 0;
}
// 输出结果
i1 * i2 = -50 * 4 = 56
i1 << k = -50 << 2 = 56

那么,对于与任意常数 K 的相乘,有没可能转换为移位操作

任意常数 K,可以表示成 [ ( 0...0 ) ( 1...1 ) ( 0...0 ) ⋅ ⋅ ⋅ ( 1...1 ) ] b [(0...0)(1...1)(0...0)\cdot\cdot\cdot(1...1)]_b [(0...0)(1...1)(0...0)(1...1)]b 的形式,也即,由一系列连续的 0 和 连续的 1 组成,比如 14 可以表示成 [ ( 0...0 ) ( 111 ) ( 0 ) ] b [(0...0)(111)(0)]_b [(0...0)(111)(0)]b

假设只存在一个连续的 1 序列,从高到低, 位于 n 到 m 位,比如 14 中,n = 3,m = 1,那么:
x ⋅ K = ( x < < n ) + ( x < < n − 1 ) + . . . + ( x < < m ) x\cdot K=(x<< n)+(x<< n-1)+...+(x<< m) xK=(x<<n)+(x<<n1)+...+(x<<m)

x ⋅ K = ( x < < ( n + 1 ) ) − ( x < < m ) x\cdot K=(x<< (n+1))-(x<< m) xK=(x<<(n+1))(x<<m)

比如, x ⋅ 14 x\cdot14 x14 就可以表示成 ( x < < 3 ) + ( x < < 2 ) + ( x < < 1 ) (x<< 3)+(x<< 2)+(x<< 1) (x<<3)+(x<<2)+(x<<1) 或者 ( x < < 4 ) − ( x < < 1 ) (x<< 4)-(x<< 1) (x<<4)(x<<1)

// C++
int main() {
    int8_t x = 5;
    int8_t K = 14;
    printf("x * 14 = %d\n", x*K);
    printf("(x<<3) + (x<<2) + (x<<1) = %d\n", (x<<3)+(x<<2)+(x<<1));
    printf("(x<<4) - (x<<1) = %d\n", (x<<4)-(x<<1));
    return 0;
}
// 输出结果
x * 14 = 70
(x<<3) + (x<<2) + (x<<1) = 70
(x<<4) - (x<<1) = 70

同理,当 K 的二进制表示,存在多个连续的 1 序列时,也成立。

除法

除法运算比乘法运算更慢,通常需要 30 个 CPU 时钟以上。同理, x / 2 k x/2^k x/2k 也可以转换成右移运算,注意结果的取整:
⌊ x / 2 k ⌋ = x > > k   o r   x > > > k \lfloor x/2^k\rfloor = x >> k\ or\ x >>> k x/2k=x>>k or x>>>k

⌈ x / 2 k ⌉ = ( x + ( 1 < < k ) − 1 ) > > k   o r   ( x + ( 1 < < k ) − 1 ) > > > k \lceil x/2^k \rceil = (x+(1<< k)-1)>> k\ or\ (x+(1<< k)-1)>>> k x/2k=(x+(1<<k)1)>>k or (x+(1<<k)1)>>>k

// C++
int main() {
    int8_t x = -50;
    int8_t K = 4;
    printf("x / 4 = %d\n", x/K);
    printf("x >> 2 = %d\n", x>>2);
    printf("(x+(1<<2))>>2 = %d\n", (x+(1<<2))>>2);

    uint8_t y = 50;
    uint8_t Z = 4;
    printf("y / 4 = %d\n", y/Z);
    printf("y >>> 2 = %d\n", y>>2);
    printf("(y+(1<<2))>>>2 = %d\n", (y+(1<<2))>>2);
    return 0;
}
// 输出结果
x / 4 = -12
x >> 2 = -13
(x+(1<<2))>>2 = -12

y / 4 = 12
y >>> 2 = 12
(y+(1<<2))>>>2 = 13

浮点数

考虑如下一段代码:

// Java
public class Example {
    public static void main(String[] args) {
        int i = 12345;
        float f = 12345.0F;
        System.out.printf("binary str of i: %s\n", int2BinaryStr(i));
        System.out.printf("binary str of f: %s\n", float2BinaryStr(f));
    }
    // 将int转为32位的二进制表示
    private static String int2BinaryStr(int i) {
        String str = Integer.toBinaryString(i);
        for (int len = str.length(); len < 32; len++) {
            str = "0" + str;
        }
        return str;
    }
    // 将float转为32位的二进制表示
    private static String float2BinaryStr(float f) {
        return int2BinaryStr(Float.floatToRawIntBits(f));
    }
}

我们将整数 i = 12345 和 浮点数 f = 12345.0 转成二进制表示,结果如下:

binary str of i: 00000000000000000011000000111001
binary str of f: 01000110010000001110010000000000

虽然从十进制上看,值都是 12345,但是两种类型的二进制表示却差别很大,说明,计算机系统对整数和浮点数的处理是两套不同的机制

大部分编程语言支持的整数最大是 64 位,但很多场景下,我们需要更大的取值范围,或者是小数运算,这些都是 64 位整数无法满足的场景。为了解决这些,计算机系统引入了 浮点数Floating Point)类型。

浮点数类型可以分为 32 位单精度 float/float32 类型和 64 位双精度 double/float64 类型,取值范围如下:

类型最小值最大值
float-3.40282347E+383.40282347E+38
double-1.79769313486231570E+3081.79769313486231570E+308

虽然浮点数的取值范围很广,但它只能精确表示其中一小部分,其他都是近似表示

下面,我们将深入介绍浮点数的二进制表示和近似规则。

简单浮点数编码

对于十进制数 d d d,可以表示成 [ d m d m − 1 ⋅ ⋅ ⋅ d 1 d 0 . d − 1 d − 2 ⋅ ⋅ ⋅ d − n + 1 d − n ] d [d_md_{m-1}\cdot\cdot\cdot d_1d_0.d_{-1}d_{-2}\cdot\cdot\cdot d_{-n+1}d_{-n}]_d [dmdm1d1d0.d1d2dn+1dn]d,其中,以小数点 . . . 为界,左边为整数部分,右边为小数部分,那么:
d = ∑ i = − n m 1 0 i ⋅ d i d=\sum_{i=-n}^m10^i\cdot d_i d=i=nm10idi
比如, 123.4 5 d = 1 0 2 ⋅ 1 + 1 0 1 ⋅ 2 + 1 0 0 ⋅ 3 + 1 0 − 1 ⋅ 4 + 1 0 − 2 ⋅ 5 123.45_d=10^2\cdot1+10^1\cdot2+10^0\cdot3+10^{-1}\cdot4+10^{-2}\cdot5 123.45d=1021+1012+1003+1014+1025

同理,对于二进制数 b b b,也可以表示成 [ b m b m − 1 ⋅ ⋅ ⋅ b 1 b 0 . b − 1 b − 2 ⋅ ⋅ ⋅ b − n + 1 b − n ] b [b_mb_{m-1}\cdot\cdot\cdot b_1b_0.b_{-1}b_{-2}\cdot\cdot\cdot b_{-n+1}b_{-n}]_b [bmbm1b1b0.b1b2bn+1bn]b,同样以小数点 . . . 分割整数和小数部分,那么:
b = ∑ i = − n m 2 i ⋅ b i b=\sum_{i=-n}^m2^i\cdot b_i b=i=nm2ibi
比如, 101.1 1 b = 2 2 ⋅ 1 + 2 1 ⋅ 0 + 2 0 ⋅ 1 + 2 − 1 ⋅ 1 + 2 − 2 ⋅ 1 = 5.7 5 d 101.11_b=2^2\cdot1+2^1\cdot0+2^0\cdot1+2^{-1}\cdot1+2^{-2}\cdot1=5.75_d 101.11b=221+210+201+211+221=5.75d

计算机系统的这种编码方式,注定只能表示 M ∗ 2 E M*2^E M2E 形式的数,其他的,只能近似表示。

比如,无法精确表示 0.2:

二进制表示分数十进制表示
0.00/20.0
0.011/40.25
0.0102/80.25
0.00113/160.1875
0.001106/320.1875
0.00110113/640.203125
0.001101026/1280.203125
0.0011001151/2560.19921875

浮点数除了表示小数之外,还须表示大数,按照这种二进制编码思路,如果要表示 5 ⋅ 2 100 5\cdot2^{100} 52100 ,需要 102 位,受限于计算机系统的编码长度,这显然不能接受。

其实,对 M ⋅ 2 E M\cdot2^E M2E 形式的数,我们只须存储 M 和 E 即可,再来一个符号位 s,即可表示这种形式的数: V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E

对 w 位浮点数,可以用 1 位表示 s,用 k 位表示 E,用 n 位表示 M,那么就有 w = 1 + k + n:

(1)s 的编码

用 0 表示正数,1 表示负数。

(2)E 的编码

二进制表示为 [ e k − 1 e k − 2 . . . e 1 e 0 ] [e_{k-1}e_{k-2}...e_1e_0] [ek1ek2...e1e0],因为 E 可以是正,也可以是负,因此,可以直接用 补码 进行编码。

比如,k = 4 时,E 的取值是这样的:

eE
00000
00011
01117
1000-8
1001-7
1111-1

但是,这样从 [ 0000 ] b [0000]_b [0000]b [ 1111 ] b [1111]_b [1111]b 并不是递增的趋势。另一种方法是,对 e 用无符号编码,令 E = e − B i a s E = e - Bias E=eBias,其中, B i a s = 2 k − 1 − 1 Bias=2^{k-1}-1 Bias=2k11。比如,k = 4 时, B i a s = 7 Bias=7 Bias=7,那么 E 的取值是这样的:

eBiasE
00007-7
00017-6
011170
100071
100172
111178

(3)M 的编码

二进制表示为 [ m n − 1 m m − 2 . . . m 1 m 0 ] [m_{n-1}m_{m-2}...m_1m_0] [mn1mm2...m1m0],M 不涉及负数,因此可以只用最高位表示整数,其他表示小数,即 [ m n − 1 . m m − 2 . . . m 1 m 0 ] [m_{n-1}.m_{m-2}...m_1m_0] [mn1.mm2...m1m0]。比如,当 n = 3 时,M 的取值如下:

mM十进制
0000.000.00 (0/4)
0010.010.25 (1/4)
0100.100.5 (2/4)
0110.110.75 (3/4)
1001.001.00 (4/4)
1011.011.25 (5/4)
1101.101.5 (6/4)
1111.111.75 (7/4)

但是,这种编码方式会导致 0 有很多种表示。比如 w = 8,k = 4,n = 3 时,因为 V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E,所以 [ 00000000 ] b [00000000]_b [00000000]b [ 00001000 ] b [00001000]_b [00001000]b [ 00010000 ] b [00010000]_b [00010000]b 等的浮点数值都是 0。

可以这么改动,令 m = [ 0. m n − 1 m m − 2 . . . m 1 m 0 ] m=[0.m_{n-1}m_{m-2}...m_1m_0] m=[0.mn1mm2...m1m0],当 e = 0 e = 0 e=0 时, M = m M = m M=m;当 e ≠ 0 e\neq0 e=0 时, M = 1 + m M=1+m M=1+m

emM十进制
00000000.0000.000 (0/8)
00000010.0010.125 (1/8)
00001110.1110.875 (7/8)
00010001.0001.000 (8/8)
00011.0011.125 (9/8)

这样,既解决了 0 的表示问题,又让 M 的精度更大了

IEEE 浮点数编码

上述这些编码方法,就是 IEEE 754 Floating-Point Representation 标准中定义的浮点数编码方式的基本思路,标准会在此基础上做了一些调整,比如新增了 Infinity 和 NaN 的表示。

IEEE 标准将浮点数分成 单精度双精度 两种,分别用 32 位和 64 位表示,它们的 k 值和 n 值都不同:

在此基础上,根据 e 的取值不同,又分为 3 种场景,Denormalized ValuesNormalized ValuesSpecial Values

(1)Denormalized Values

此场景下,e 取值为全 0,用来表示 0 值,以及接近 0 的小数

V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E 的表示下, M = m M = m M=m E = 1 − B i a s = 1 − ( 2 k − 1 − 1 ) E = 1 - Bias = 1 - (2^{k-1}-1) E=1Bias=1(2k11)

为什么 ,而不是 E = -Bias

根据前文分析,应该有 E = e - Bias,但 Denormalized Values 中,确是 E = 1 - Bias,而不是 E = 0 - Bias。

这主要考虑到,从 Denormalized Values 到 Normalized Values 的平滑过度,让 Largest Denormalized Value 和 Smallest Normalized Value 更接近。后面的举例可以看到。

这里,0 分为了 -0.0 和 +0.0,在一些科学计算场景下,它们会表示不同的含义。比如 1 / ( + 0.0 ) = + ∞ 1/(+0.0) =+\infty 1/(+0.0)=+ 1 / ( − 0.0 ) = − ∞ 1/(-0.0) =-\infty 1/(0.0)=

(2) Normalized Values

此场景下,e 取值不是全 0,也不是全 1,是最常见的场景。

V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E 的表示下, M = 1 + m = 1. m n − 1 m n − 2 . . . m 0 M = 1 + m = 1.m_{n-1}m_{n-2}...m_0 M=1+m=1.mn1mn2...m0 E = e − B i a s = e − ( 2 k − 1 − 1 ) E = e - Bias = e - (2^{k-1}-1) E=eBias=e(2k11)

(3)Special Values

此场景下,e 取值为全 1,用来表示 Infinity 和 NaN。

  • m 取值全 0,表示 Infinity:1)s = 0 时,为 + ∞ +\infty +;2)s = 1 时,为 − ∞ -\infty
  • m 取值不是全 0,表示 NaN(Not a Number),比如 − 1 \sqrt{-1} 1 或者 ∞ − ∞ \infty-\infty

比如,w = 8,k = 4(此时, B i a s = 2 k − 1 − 1 = 7 Bias=2^{k-1}-1=7 Bias=2k11=7),n = 3 时,IEEE 标准浮点数的编码例子如下:

场景二进制eE 2 E 2^E 2EmM 2 E ⋅ M 2^E\cdot M 2EM十进制
Denormalized0 0000 0000-61/640/80/80/5120.0
0 0000 0010-61/641/81/81/5120.001953
0 0000 1110-61/647/87/87/5120.013672
Normalized0 0001 0001-61/640/88/88/5120.015625
0 0001 0011-61/641/89/89/5120.017578
0 1110 1111471287/815/81920/8240.0
Infinity0 1111 000------ + ∞ +\infty +
NaN0 1111 001------NaN
NaN

从上表可以看出,因为 Denormalized 场景下令 E = 1 - Bias,从 Denormalized 到 Normalized 的过度,也即 7/512 到 8/512,更平滑了。如果令 E = - Bias,则是从 7/1024 到 8/512,跨度太大。

到这里,我们已经深入介绍了 IEEE 浮点数编码格式,回到本节最开始的例子,如何从 int i = 12345 的二进制表示,推断出 float f = 12345.0 的二进制表示?

  1. 首先,12345 整数的二进制表示为 1100000011100 1 b 11000000111001_b 11000000111001b ,也可以表示成 1.100000011100 1 b ⋅ 2 13 1.1000000111001_b \cdot 2^{13} 1.1000000111001b213 ,对应到 V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E 的形式,可以确认 s = 0,E = 13,M = 1.1000000111001。
  2. 另外,12345 属于 Normalized 场景,而 float 类型中 k = 8,有, E = e − ( 2 8 − 1 − 1 ) = 13 E = e - (2^{8-1}-1) = 13 E=e(2811)=13,得出 e = 140,按 k 位无符号编码表示为 [ 10001100 ] b [10001100]_b [10001100]b
  3. 同理,由 M = 1 + 0. m n − 1 m n − 2 . . . m 0 = 1.1000000111001 M = 1 + 0.m_{n-1}m_{n-2}...m_0=1.1000000111001 M=1+0.mn1mn2...m0=1.1000000111001,float 类型中,n = 23,所以,得出 m 的二进制表示为 [ 10000001110010000000000 ] b [10000001110010000000000]_b [10000001110010000000000]b,注意,低位补零
  4. 最后将 s、e、m 按照 float 单精度的编码格式组合起来,就是 [ 0   10001100   10000001110010000000000 ] b [0\ 10001100\ 10000001110010000000000]_b [0 10001100 10000001110010000000000]b,也即 12345.0 的二进制编码。

近似规则(Rounding)

前面说过,浮点数只能表示 M ∗ 2 E M*2^E M2E 形式的数,其他的,只能近似表示。

近似规则,我们最熟悉的是 “四舍五入”,比如,要保留 2 位小数,那么 1.234、1.235、1.236 近似之后分别是 1.23、1.24、1.24。

IEEE 浮点数标准中,并没采用 “四舍五入” 法,定义了如下 4 种近似规则:

  • Round-to-even,往更近的方向靠,如果向上和向下的距离一样,则往偶数的方向靠。
  • Round-toward-zero,往 0 的方向靠。
  • Round-down,往更小的方向靠。
  • Round-up,往更大的方向靠。

比如,下面的例子,需要近似为整数:

规则1.401.601.502.50-1.50
Round-to-even1222-2
Round-toward-zero1112-1
Round-down11-2
Round-up223-1

这些规则对二进制也生效,比如,Round-to-even 规则,在二进制中,0 代表偶数,1 表示奇数。下面例子中,需要按照 Round-to-even 近似保留 2 位小数:

二进制(十进制)向下中点向上Round-to-even
10.00011 (67/32)10.00 (64/32)68/3210.01 (72/32)10.00 (64/32)
10.00110 (70/32)10.00 (64/32)68/3210.01 (72/32)10.01 (72/32)
10.11100 (23/8)10.11 (22/8)23/811.00 (24/8)11.00 (24/8)
10.10100 (21/8)10.10 (20/8)21/810.11 (22/8)10.10 (20/8)

上述例子中,前 2 个按照就近原则近似;后 2 个处于中点位置,往偶数方向靠。

注意,IEEE 标准规定 Round-to-even 为默认的近似规则

浮点数运算

浮点数的运算规则比较简单,令 ⊕ \oplus 指代后一类的运算符,x 和 y 为浮点数类型,那么 x ⊕ y x \oplus y xy 的运算结果为 R o u n d ( x ⊕ y ) Round(x \oplus y) Round(xy),也即对真实计算结果进行近似。

注意,算术运算的结合律分配率在浮点数运算下是不生效的,比如:

// Java
public static void main(String[] args) {
    // 结合律不满足示例
    float f1 = 3.14F;
    float f2 = 1e10F;
    float f3 = (f1 + f2) - f2;
    float f4 = f1 + (f2 - f2);
    System.out.printf("(3.14 + 1e10) - 1e10 = %f\n", f3);
    System.out.printf("3.14 + (1e10 - 1e10) = %f\n", f4);

    // 分配律不满足示例
    float f5 = 1e20F;
    float f6 = f5 * (f5 - f5);
    float f7 = f5 * f5 - f5 * f5;
    System.out.printf("1e20 * (1e20 - 1e20) = %f\n", f6);
    System.out.printf("1e20 * 1e20 - 1e20 * 1e20 = %f\n", f7);
}
// 输出结果
(3.14 + 1e10) - 1e10 = 0.000000
3.14 + (1e10 - 1e10) = 3.140000
1e20 * (1e20 - 1e20) = 0.000000
1e20 * 1e20 - 1e20 * 1e20 = NaN

结合律例子中,(3.14+1e10)-1e10 结果是 0,因为 3.14 在近似时,精度丢失了,也即 3.14+1e10=1e10;类似,分配律例子中,1e20*1e20 的结果超出了 float 类型的表示范围,得到 Infinity,而 Infinity - Infinity 的结果是 NaN

注意,浮点数运算结果,如果超出了浮点数表示范围,会得到 Infinity/-Infinity。这与整数运算的溢出机制有所区别。

数值类型转换

数值类型间的转换,可以分成 2 类:宽转换(Widening Conversion)和 窄转换(Narrowing Conversion)。

宽转换指往表示范围更广的类型转换,比如从 int 到 long、从 long 到 float;窄转换则相反。

整型间转换

(1)宽转换

整型间的宽转换不会产生溢出,无符号整数场景,高位补零;有符号整数场景,高位补符号位。

// C++
int main() {
    int8_t i1 = 100;
    cout << "int8_t  i1: " << bitset<8>(i1) << endl;
    cout << "int16_t i1: " << bitset<16>((int16_t) i1) << endl;

    int8_t i2 = -100;
    cout << "int8_t  i2: " << bitset<8>(i2) << endl;
    cout << "int16_t i2: " << bitset<16>((int16_t) i2) << endl;

    uint8_t i3 = 200;
    cout << "uint8_t  i3: " << bitset<8>(i3) << endl;
    cout << "uint16_t i3: " << bitset<16>((uint16_t) i3) << endl;
    return 0;
}
// 输出结果
int8_t  i1: 01100100
int16_t i1: 0000000001100100
int8_t  i2: 10011100
int16_t i2: 1111111110011100
uint8_t  i3: 11001000
uint16_t i3: 0000000011001000

(2)窄转换

整型间的窄转换直接进行高位截断,只保留低 n 位。比如 16 位的 int16 转换为 8 位的 int8,直接保留 int16 类型值的低 8 位作为转换结果。

// C++
int main() {
    int16_t i1 = 200;
    cout << "int16_t i1: " << bitset<16>(i1) << endl;
    cout << "int8_t  i1: " << bitset<8>((int8_t) i1) << endl;

    int16_t i2 = -200;
    cout << "int16_t i2: " << bitset<16>(i2) << endl;
    cout << "int8_t  i2: " << bitset<8>((int8_t) i2) << endl;

    uint16_t i3 = 300;
    cout << "uint16_t i3: " << bitset<16>(i3) << endl;
    cout << "uint8_t  i3: " << bitset<8>((uint8_t) i3) << endl;
    return 0;
}
// 输出结果
int16_t i1: 0000000011001000
int8_t  i1: 11001000
int16_t i2: 1111111100111000
int8_t  i2: 00111000
uint16_t i3: 0000000100101100
uint8_t  i3: 00101100

(3)无符号整数与有符号整数间的转换

无符号整数与有符号整数间的转换规则是:

  • 如果两者二进制位数一致,比如 int8uint8 的转换,则二进制数值不变,只是改变编码方式;
  • 如果位数不一致,比如 int16uint8 的转换,则二进制数值,先按照宽转换或窄转换规则转换,再改变编码方式。
// C++
int main() {
    uint8_t i1 = 200;
    cout << "uint8_t i1, decimal: " << +i1 << ", binary: " << bitset<8>(i1) << endl;
    cout << "int8_t  i1, decimal: " << +(int8_t) i1 << ", binary: " << bitset<8>((int8_t) i1) << endl;

    int16_t i2 = -300;
    cout << "int16_t i2, decimal: " << +i2 << ", binary: " << bitset<16>(i2) << endl;
    cout << "uint8_t i2, decimal: " << +(uint8_t) i2 << ",  binary: " << bitset<8>((uint8_t) i2) << endl;
    return 0;
}
// 输出结果
uint8_t i1, decimal: 200, binary: 11001000
int8_t  i1, decimal: -56, binary: 11001000
int16_t i2, decimal: -300, binary: 1111111011010100
uint8_t i2, decimal: 212,  binary: 11010100

整数与浮点数间转型

(1)宽转换

整型到浮点数类型的转换这一方向,为宽转换:

  • 如果浮点数的精度,能够表示整数,则正常转换。
  • 如果浮点数精度,无法表示整数,则需要近似,会导致精度丢失。
// Java
public static void main(String[] args) {
    int i1 = 1234567;
    System.out.printf("int i1: %d, float i1: ", i1);
    System.out.println((float) i1);

    int i2 = 123456789;
    System.out.printf("int i2: %d, float i2: ", i2);
    System.out.println((float) i2);
}
// 输出结果
int i1: 1234567, float i1: 1234567.0
int i2: 123456789, float i2: 1.23456792E8

上述例子中,i2=123456789 超过 float 类型能够表示的精度,所以为近似后的结果 1.23456792E8

那么,为什么 123456789 会近似为 1.23456792E8

要解释该问题,首先要把它们转换成二进制表示:

public static void main(String[] args) {
    ...
    System.out.println("int i2:   " + int2BinaryStr(i2));
    System.out.println("float i2: " + float2BinaryStr((float) i2));
}
// 输出结果
int i2:   00000111010110111100110100010101
float i2: 01001100111010110111100110100011

接下来,我们根据 IEEE 浮点数的编码规则,尝试将 int i2 转换成 float i2

  1. int i2 的二进制 0000011101011011110011010001010 1 b 00000111010110111100110100010101_b 00000111010110111100110100010101b,可以写成 1.11010110111100110100010101 ⋅ 2 26 1.11010110111100110100010101\cdot2^{26} 1.11010110111100110100010101226,对应到 V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E 的形式,可以确认 s = 0,E = 26,M = 1.11010110111100110100010101。
  2. float 类型中 k = 8,有, E = e − ( 2 8 − 1 − 1 ) = 26 E = e - (2^{8-1}-1) = 26 E=e(2811)=26,得出 e = 153,按 k 位无符号编码表示为 [ 10011001 ] b [10011001]_b [10011001]b
  3. 同理,由 M = 1 + 0. m n − 1 m n − 2 . . . m 0 = 1.11010110111100110100010101 M = 1 + 0.m_{n-1}m_{n-2}...m_0=1.11010110111100110100010101 M=1+0.mn1mn2...m0=1.11010110111100110100010101,但由于 float 类型的 n = 23,而 m 一共有 26 位,因此需要按照 round-to-even 规则,对 0.11010110111100110100010101进行近似,保留 23 位小数,得到 0.11010110111100110100011,所以 m 为 [ 11010110111100110100011 ] b [11010110111100110100011]_b [11010110111100110100011]b
  4. 最后,将 s、e、m 按照 float 单精度的编码格式组合起来,就是 [ 0   10011001   11010110111100110100011 ] b [0\ 10011001\ 11010110111100110100011]_b [0 10011001 11010110111100110100011]b,转换成十进制,就是 1.23456792E8

(2)窄转换

浮点数类型到整型的转换这一方向,为窄转换:

  • 如果浮点数的整数部分,能够用整型表示,则直接舍去小数,保留整数部分。
  • 如果超出了整型范围,则结果为该整型的最大/最小值。
// Java
public static void main(String[] args) {
    float f1 = 12345.123F;
    System.out.print("float f1: ");
    System.out.print(f1);
    System.out.printf(",  int f1: %d\n", (int) f1);

    float f2 = 1.2345E20F;
    System.out.print("float f2: ");
    System.out.print(f2);
    System.out.printf(",  int f2: %d\n", (int) f2);

    float f3 = -1.2345E20F;
    System.out.print("float f3: ");
    System.out.print(f3);
    System.out.printf(", int f3: %d\n", (int) f3);
 }
// 输出结果
float f1: 12345.123,  int f1: 12345
float f2: 1.2345E20,  int f2: 2147483647
float f3: -1.2345E20, int f3: -2147483648

浮点数间转型

(1)宽转换

单精度 float 到 双精度 double 为宽转换,不会出现精度丢失的问题。

对于 V = ( − 1 ) s ⋅ M ⋅ 2 E V=(-1)^s\cdot M\cdot2^E V=(1)sM2E ,规则如下:

  • s 保持不变。
  • 在 E 保持不变的前提下,因为 float 的 k = 8,而 double 的 k = 11,所以两者的 e 会有所不同。
  • 在 M 保持不变的前提下,float 的 n = 23,而 double 的 n =52,所以 m 需要低位补 52 - 23 = 29 个 0。
// Java
public static void main(String[] args) {
    float f1 = 1.2345E20F;
    System.out.print("float  f1: ");
    System.out.print(f1);
    System.out.print(", double f1: ");
    System.out.println((double) f1);

    System.out.println("float  f1: " + float2BinaryStr(f1));
    System.out.println("double f1: " + double2BinaryStr((double) f1));
}
// 输出结果
float  f1: 1.2345E20, double f1: 1.2344999897320129E20
float  f1: 01100000110101100010011011010000
double f1: 0100010000011010110001001101101000000000000000000000000000000000

(2)窄转换

doublefloat 为窄转换,会存在精度丢失问题。

如果 double 值超出了 float 的表示范围,则转换结果为 Infinity

// Java
public static void main(String[] args) {
    double d1 = 1E200;
    System.out.print("double d1: ");
    System.out.println(d1);
    System.out.print("float d1: ");
    System.out.println((float) d1);

    double d2 = -1E200;
    System.out.print("double d2: ");
    System.out.println(d2);
    System.out.print("float d2: ");
    System.out.println((float) d2);
 }
// 输出结果
double d1: 1.0E200
float d1: Infinity
double d2: -1.0E200
float d2: -Infinity

如果 double 值还在 float 的表示范围内,则按照如下转换规则:

  • s 保持不变。
  • 在 E 保持不变的前提下,因为 float 的 k = 8,而 double 的 k = 11,所以两者的 e 会有所不同。
  • 对于 M,因为 float 的 n = 23,而 double 的 n = 52,所以转换到 float 之后,需要进行截断,只保留高 23 位。
// Java
public static void main(String[] args) {
    double d1 = 3.267393471324506;
    System.out.print("double d1: ");
    System.out.println(d1);
    System.out.print("float  d1: ");
    System.out.println((float) d1);
    System.out.println("double d1: " + double2BinaryStr(d1));
    System.out.println("float  d1: " + float2BinaryStr((float) d1));
}
// 输出结果
double d1: 3.267393471324506
float  d1: 3.2673936
double d1: 0100000000001010001000111001111100110000001101000000010101110110
float  d1: 01000000010100010001110011111010

最后

本文花了很长的篇幅,深入介绍了计算机系统对数值类型的编码、运算、转换的底层原理。

数值类型间的转换是最容易出现隐藏 bug 的地方,特别是无符号整数与有符号整数之间的转换。所以,很多现代的编程语言,如 Java、Go 等都不再支持无符号整数,根除了该隐患。

另外,浮点数的编码方式,注定它只能精确表示一小部分的数值范围,大部分都是近似,所以才有了不能用等号来比较两个浮点数的说法。

数值类型虽然很基础,但使用时一定要多加小心。希望本文能够加深你对数值类型的理解,让你写出更健壮的程序。

文章配图

可以在 用Keynote画出手绘风格的配图 中找到文章的绘图方法。

参考

[1] Computer Systems: A Programmer’s Perspective (3rd edition), Randal E. Bryant .etc

[2] The Java® Language Specification (Java SE 18 Edition), James Gosling .etc

[3] IEEE Standard 754 Floating-Point Representation, IEEE

[4] 漫话:为什么计算机用补码存储数据?, 漫话编程

更多文章请关注微信公众号:元闰子的邀请

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值