DES对称加密算法详解和c++代码实现（带样例和详细的中间数据）

DES对称加密算法详解和c++代码实现（带样例和详细的中间数据）

特点：

1.DES是对称性加密算法，即加密和解密是对称的，用的是同一个密钥

2.DES只处理二进制数据，所以需要将明文转换成为2进制数据

3.DES每次处理64位的数据，所以应该将明文切割成64位的分组，当最后一组数据不足64位的时候，高位补0

4.DES使用64位的密钥，但因为密钥中的每8位会被忽略，所以有效的密钥长度是56位，从而产生16个48位的子密钥（变换过程后面会说明）

5.每64位数据一个块，是DES的永恒组织方式

 

具体样例分析：

（仅以一组64位数据为例分析加密过程）

明文M是：8787878787878787

密钥K是：0E329232EA6D0D73

上面的信息都是16进制的，转换为2进制

明文M是：0000 0001 0010 0011 0100 0101 0110 0111 1000 1001 1010 1011 1100 1101 1110 1111

密钥K是：00010011 00110100 01010111 01111001 10011011 10111100 11011111 11110001

 

第一步：根据密钥生成16个子密钥

1.根据密钥初始置换表将64位的密钥转化为58位的密钥

57 49 41 33 25 17 9

1 58 50 42 34 26 18

10 2 59 51 43 35 27

19 11 3 60 52 44 36

63 55 47 39 31 23 15

7 62 54 46 38 30 22

14 6 61 53 45 37 29

21 13 5 28 20 12 4

由于上表中第一个元素为57，这将使原秘钥的第57位变换为新秘钥K+的第1位。同理，原秘钥的第49位变换为新秘钥的第2位……原秘钥的第4位变换为新秘钥的最后一位。注意原秘钥中只有56位会进入新秘钥，上表也只有56个元素。

原密钥K：00010011 00110100 01010111 01111001 10011011 10111100 11011111 11110001

新密钥K：1111000 0110011 0010101 0101111 0101010 1011001 1001111 0001111

2.将新密钥拆分成C0和D0，每组都有28位

比如新密钥

C0：1111000 0110011 0010101 0101111

D0：0101010 1011001 1001111 0001111

3.根据密钥轮次左移表，左移特定的位数

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16

1  1  2  2  2  2  2  2  1   2    2    2    2    2    2   1

比如第一轮是左移1位，第二轮也是左移1位，第三轮是左移两位

所以

C1：1110000110011001010101011111

D1：1010101011001100111100011110

下面给出C1,D1到C16，D16的数据：

C1 = 1110000110011001010101011111

D1 = 1010101011001100111100011110

C2 = 1100001100110010101010111111

D2 = 0101010110011001111000111101

C3 = 0000110011001010101011111111

D3 = 0101011001100111100011110101

C4 = 0011001100101010101111111100

D4 = 0101100110011110001111010101

C5 = 1100110010101010111111110000

D5 = 0110011001111000111101010101

C6 = 0011001010101011111111000011

D6 = 1001100111100011110101010101

C7 = 1100101010101111111100001100

D7 = 0110011110001111010101010110

C8 = 0010101010111111110000110011

D8 = 1001111000111101010101011001

C9 = 0101010101111111100001100110

D9 = 0011110001111010101010110011

C10 = 0101010111111110000110011001

D10 = 1111000111101010101011001100

C11 = 0101011111111000011001100101

D11 = 1100011110101010101100110011

C12 = 0101111111100001100110010101

D12 = 0001111010101010110011001111

C13 = 0111111110000110011001010101

D13 = 0111101010101011001100111100

C14 = 1111111000011001100101010101

D14 = 1110101010101100110011110001

C15 = 1111100001100110010101010111

D15 = 1010101010110011001111000111

C16 = 1111000011001100101010101111

D16 = 0101010101100110011110001111

需要记住的是：每一对Cn 和 Dn都是由前一对Cn-1 和 Dn-1移位而来！

4.得到Cn，Dn后合并CnDn，然后根据密钥压缩置换表将56位密钥压缩成48位的子密钥

密钥压缩置换表：

14 17 11 24 1 5

3 28 15 6 21 10

23 19 12 4 26 8

16 7 27 20 13 2

41 52 31 37 47 55

30 40 51 45 33 48

44 49 39 56 34 53

46 42 50 36 29 32

每对子秘钥有56位，但PC-2仅仅使用其中的48位。

于是，第n轮的新秘钥Kn 的第1位来自组合子秘钥CnDn的第14位，第2位来自第17位，依次类推，知道新秘钥的第48位来自组合秘钥的第32位。

下面给出子密钥K1到K16的数据：

K1 = 000110 110000 001011 101111 111111 000111 000001 110010

K2 = 011110 011010 111011 011001 110110 111100 100111 100101

K3 = 010101 011111 110010 001010 010000 101100 111110 011001

K4 = 011100 101010 110111 010110 110110 110011 010100 011101

K5 = 011111 001110 110000 000111 111010 110101 001110 101000

K6 = 011000 111010 010100 111110 010100 000111 101100 101111

K7 = 111011 001000 010010 110111 111101 100001 100010 111100

K8 = 111101 111000 101000 111010 110000 010011 101111 111011

K9 = 111000 001101 101111 101011 111011 011110 011110 000001

K10 = 101100 011111 001101 000111 101110 100100 011001 001111

K11 = 001000 010101 111111 010011 110111 101101 001110 000110

K12 = 011101 010111 000111 110101 100101 000110 011111 101001

K13 = 100101 111100 010111 010001 111110 101011 101001 000001

K14 = 010111 110100 001110 110111 111100 101110 011100 111010

K15 = 101111 111001 000110 001101 001111 010011 111100 001010

K16 = 110010 110011 110110 001011 000011 100001 011111 110101

到这里，16个48位的子密钥就生成了

 

第二步：加密数据块

 

1.明文初始变换，仍然是64位变为64位

58 50 42 34 26 18 10 2

60 52 44 36 28 20 12 4

62 54 46 38 30 22 14 6

64 56 48 40 32 24 16 8

57 49 41 33 25 17 9 1

59 51 43 35 27 19 11 3

61 53 45 37 29 21 13 5

63 55 47 39 31 23 15 7

参照上表，M的第58位成为IP的第1位，M的第50位成为IP的第2位，M的第7位成为IP的最后一位。

对明文M：0000 0001 0010 0011 0100 0101 0110 0111 1000 1001 1010 1011 1100 1101 1110 1111

IP：1100 1100 0000 0000 1100 1100 1111 1111 1111 0000 1010 1010 1111 0000 1010 1010

这里M的第58位是1，变成了IP的第1位。M的第50位是1，变成了IP的第2位。M的第7位是0，变成了IP的最后一位。

 

2.数据分组：IP分为左半边的L0和右半边的R0

对上面的IP

L0：1100 1100 0000 0000 1100 1100 1111 1111

R0：1111 0000 1010 1010 1111 0000 1010 1010

 

3.16轮迭代

然后我们定义一个函数f，输入一个32位的数据和一个48位的子密钥，产生一个32位的输出

从n=1到n=16

Ln=Rn-1

Rn=Ln-1^f(Rn-1,Kn)

通俗的说就是当前迭代的L是上一次迭代的L，当前迭代的R是上一次迭代的L异或f（上次迭代的R，当前子密钥）

下面就是要了解f函数是怎么工作的了

f函数第一步：32位的数据扩展成为48位

明文扩展表：

32 1 2 3 4 5

4 5 6 7 8 9

8 9 10 11 12 13

12 13 14 15 16 17

16 17 18 19 20 21

20 21 22 23 24 25

24 25 26 27 28 29

28 29 30 31 32 1

对R0：1111 0000 1010 1010 1111 0000 1010 1010

扩展一下：E(R0) = 011110 100001 010101 010101 011110 100001 010101 010101

f函数第二步：将扩展的结果和48位的子密钥进行异或运算

比如，对K1 , E(R0)，我们有：

K1 = 000110 110000 001011 101111 111111 000111 000001 110010

E(R0) = 011110 100001 010101 010101 011110 100001 010101 010101

K1^E(R0) = 011000 010001 011110 111010 100001 100110 010100 100111.

f函数第三步：S盒替代，48位变成32位

一共有8个S盒，分别位S1到S8，所以48位的数据可以分为8组，每组6个数据，用a1a2a3a4a5a6表示

每组数据对应一个S盒

a1和a6组合成2进制，然后转换为10进制，代表对应S盒的行号

a2,a3,a4,a5组合成2进制，然后转换为10进制，代表对应S盒的列号

从行号和列号可以从对应的S盒中确定一个数X，X是十进制的，转化为4位的二进制输出

所以每个S盒都是输入6位，输出4位，所以最后输出的数据是4*8=32位

                            S1

14 4 13 1 2 15 11 8 3 10 6 12 5 9 0 7

0 15 7 4 14 2 13 1 10 6 12 11 9 5 3 8

4 1 14 8 13 6 2 11 15 12 9 7 3 10 5 0

15 12 8 2 4 9 1 7 5 11 3 14 10 0 6 13

                             S2

15 1 8 14 6 11 3 4 9 7 2 13 12 0 5 10

3 13 4 7 15 2 8 14 12 0 1 10 6 9 11 5

0 14 7 11 10 4 13 1 5 8 12 6 9 3 2 15

13 8 10 1 3 15 4 2 11 6 7 12 0 5 14 9

                            S3

10 0 9 14 6 3 15 5 1 13 12 7 11 4 2 8

13 7 0 9 3 4 6 10 2 8 5 14 12 11 15 1

13 6 4 9 8 15 3 0 11 1 2 12 5 10 14 7

1 10 13 0 6 9 8 7 4 15 14 3 11 5 2 12

                            S4

7 13 14 3 0 6 9 10 1 2 8 5 11 12 4 15

13 8 11 5 6 15 0 3 4 7 2 12 1 10 14 9

10 6 9 0 12 11 7 13 15 1 3 14 5 2 8 4

3 15 0 6 10 1 13 8 9 4 5 11 12 7 2 14

                            S5

2 12 4 1 7 10 11 6 8 5 3 15 13 0 14 9

14 11 2 12 4 7 13 1 5 0 15 10 3 9 8 6

4 2 1 11 10 13 7 8 15 9 12 5 6 3 0 14

11 8 12 7 1 14 2 13 6 15 0 9 10 4 5 3

                            S6

12 1 10 15 9 2 6 8 0 13 3 4 14 7 5 11

10 15 4 2 7 12 9 5 6 1 13 14 0 11 3 8

9 14 15 5 2 8 12 3 7 0 4 10 1 13 11 6

4 3 2 12 9 5 15 10 11 14 1 7 6 0 8 13

                           S7

4 11 2 14 15 0 8 13 3 12 9 7 5 10 6 1

13 0 11 7 4 9 1 10 14 3 5 12 2 15 8 6

1 4 11 13 12 3 7 14 10 15 6 8 0 5 9 2

6 11 13 8 1 4 10 7 9 5 0 15 14 2 3 12

                          S8

13 2 8 4 6 15 11 1 10 9 3 14 5 0 12 7

1 15 13 8 10 3 7 4 12 5 6 11 0 14 9 2

7 11 4 1 9 12 14 2 0 6 10 13 15 3 5 8

2 1 14 7 4 10 8 13 15 12 9 0 3 5 6 11

 

例子：对于第一轮，我们得到这8个S盒的输出：

K1 + E(R0) = 011000 010001 011110 111010 100001 100110 010100 100111.

S1(B1)S2(B2)S3(B3)S4(B4)S5(B5)S6(B6)S7(B7)S8(B8)= 0101 1100 1000 0010 1011 0101 1001 0111

 

f函数第四步：进行P盒置换，32位输入，32位输出

P盒置换表:

16 7 20 21

29 12 28 17

1 15 23 26

5 18 31 10

2 8 24 14

32 27 3 9

19 13 30 6

22 11 4 25

比如，对于8个S盒的输出：

S1(B1)S2(B2)S3(B3)S4(B4)S5(B5)S6(B6)S7(B7)S8(B8)= 0101 1100 1000 0010 1011 0101 1001 0111

我们得到

f = 0010 0011 0100 1010 1010 1001 1011 1011

f函数到此结束

那么，得到当前迭代的R

R1 = L0 ^ f(R0 , K1 )

= 1100 1100 0000 0000 1100 1100 1111 1111
+ 0010 0011 0100 1010 1010 1001 1011 1011
= 1110 1111 0100 1010 0110 0101 0100 0100

在下一轮迭代中，我们的L2 = R1，这就是我们刚刚计算的结果。之后我们必须计算R2 =L1 + f(R1, K2)，一直完成16个迭代。在第16个迭代之后，我们有了区块L16 and R16。

下面给出L1,R1到L16，R16的数据：

L[1]:11110000101010101111000010101010
R[1]:11101111010010100110010101000100

L[2]:11101111010010100110010101000100
R[2]:11001100000000010111011100001001

L[3]:11001100000000010111011100001001
R[3]:10100010010111000000101111110100

L[4]:10100010010111000000101111110100
R[4]:01110111001000100000000001000101

L[5]:01110111001000100000000001000101
R[5]:10001010010011111010011000110111

L[6]:10001010010011111010011000110111
R[6]:11101001011001111100110101101001

L[7]:11101001011001111100110101101001
R[7]:00000110010010101011101000010000

L[8]:00000110010010101011101000010000
R[8]:11010101011010010100101110010000

L[9]:11010101011010010100101110010000
R[9]:00100100011111001100011001111010

L[10]:00100100011111001100011001111010
R[10]:10110111110101011101011110110010

L[11]:10110111110101011101011110110010
R[11]:11000101011110000011110001111000

L[12]:11000101011110000011110001111000
R[12]:01110101101111010001100001011000

L[13]:01110101101111010001100001011000
R[13]:00011000110000110001010101011010

L[14]:00011000110000110001010101011010
R[14]:11000010100011001001011000001101

L[15]:11000010100011001001011000001101
R[15]:01000011010000100011001000110100

L[16]:01000011010000100011001000110100
R[16]:00001010010011001101100110010101

 

4.逆转区位和最终变换

逆转区位：64位输出：R16+L16

然后将R16L16组合的64位数据进行最终置换，64位->64位

最终置换表：

40 8 48 16 56 24 64 32

39 7 47 15 55 23 63 31

38 6 46 14 54 22 62 30

37 5 45 13 53 21 61 29

36 4 44 12 52 20 60 28

35 3 43 11 51 19 59 27

34 2 42 10 50 18 58 26

33 1 41 9 49 17 57 25

比如，如果我们使用了上述方法得到了第16轮的左右两个区块：

1. L16 = 0100 0011 0100 0010 0011 0010 0011 0100
2. R16 = 0000 1010 0100 1100 1101 1001 1001 0101

我们将这两个区块调换位置，然后执行最终变换：