密码学入门（一）：用Python实现对称加密算法

0.前言

最开始只是想整理一下密码学课程的做业，后面越写越多，就索性写成一篇入门的介绍。我会把本身对对称加密的理解和一些做业的代码串起来，力图清晰明白地展现出来，文中全部代码都放在个人Github上，若是有错误之处还请轻拍。

• 文章地址：https://gooong.cn/post/crypto...
• 代码地址：https://github.com/Gooong/Cry...

1.对称密码基础

加密是为了防止要传达的内容被别人知道。例如，你若是想在课堂上传小纸条給后位小红说：i love coding，但又怕在递纸条的过程当中被老师看到，知道了你的心思，因而将每一个字母变字母表中的后一个字母（如a变成b，i变成j，z变成a），获得密文：j mpwf dpejoh，这样即老师人拿到这纸条，也不知道你说的是什么。

这就是一个加密的过程，把本来的内容称为明文，通常用p表示；加密后获得的内容称为密文，通常用c表示；而加密的这个过程能够看作是一个加密函数E，即

$$c = E( p )$$

E是指Encrypt，函数输入是明文，输出是加密以后的密文。上面的例子中i love coding即是明文，j mpwf dpejoh即是密文，而把字母在字母表中向后移动一位的操做就是加密函数。

在小红获得小纸条后，能够根据你加密的方法，将每一个字母变成字母表中的前一个字母，就能够从你的密文小纸条获得你要说的内容i love coding，心照不宣，顺便还会怀疑一下你的脑壳……不管怎样，这个解密的过程就也能够看作是一个解密函数D，即

$$p = D( c )$$

D是指Decrypt，函数输入是密文，输出是解密以后的明文。

在这个过程这种，小红可以成功解密小纸条的前提是，你得和她在课前约定好你加密的时候移动的是1位，2位仍是几位，否则他就会和老师同样一脸懵逼，不知道你在说啥。大家提早约定好的这个“几位”，就是加密和解密的密钥k，你会根据这个秘钥来进行加密，小红会根据这个秘钥来进行解密。

因此你的传纸条的动做抽象成这个过程：

明文p---->加密函数E---->密文c---->传输---->密文c----->解密函数D---->明文p

或者用公式来表达是：

$$c = D_k( E_k( c ) )$$

用大白话说就是：明文用同一个密钥先加密再解密获得的仍是同一个明文（等于没说...）

从这里咱们能够总结出加密体质的五个要素：{明文p，密文c，密钥k，加密函数E，解密函数D}，对称解密的的意思就是说，加密和解密的密钥是同样的，上面的过程是否是正好很对称呢?

为了方便使用，不用每次本身手动掰手指数字符，你还写了Python程序：

# 移位密码

def _move_leter(letter, n):
    """
    把字母变为字母表后n位的字母,z后面接a
    :param letter: 小写字母
    :param n: 要移动的字母
    :return: 移动的结果
    """
    return chr((ord(letter) - ord('a') + n) % 26 + ord('a'))


def Encrypt(k, p):
    """
    移位密码加密函数E
    :param k: 秘钥k,每一个字母在字母表中移动k位
    :param p: 明文p
    :return: 密文c
    """
    letter_list = list(p.lower())
    c = ''.join([_move_leter(x, k) for x in letter_list])
    return c


def Decrypt(k, c):
    """
    移位密码解密函数D
    :param k: 秘钥k,每一个字母在字母表中移动k位
    :param c: 密文c
    :return: 明文p
    """
    letter_list = list(c.lower())
    p = ''.join([_move_leter(x, -k) for x in letter_list])
    return p


if __name__ == '__main__':
    p = 'ilovecoding'
    print('明文：' + p)
    print('密文：' + Encrypt(1, p))
    print('解密：' + Decrypt(1, Encrypt(1, p)))
    assert Decrypt(1, Encrypt(1, p)) == p

运行这段代码，就能够看到输出了：

明文：ilovecoding
密文：jmpwfdpejoh
解密：ilovecoding

终于，如今你能和你的小红秘密地传达纸条内容了，迎来全班人羡慕的目光，今后走上人生巅峰，本文到此结束。

...Hey，醒醒...

2.密码分析

面对你俩日益频繁的纸条往来，老师终于坐不住了，他想知道你俩写的究竟是啥，因而在某次逮到你递纸条以后，决定下功夫破解你所使用的密码，也就是密码分析。

根据他的了解，以你的水平，最可能用的就是移位密码，但具体每次移动了几位，没法直接观察得出。不过他又一想，你移动的位数顶可能是25位，由于，移动26位的效果等于没移动，移27位的效果不就跟移动1位的效果是同样的嘛！这就是说，你的密码只能是0-25中的某一个数字，而不多是其余的，就这么二十几个秘钥，一个一个试就能知道你写的是啥！

老师果真聪明绝顶，关键是也还会Python，就索性写了一个程序，每次尝试用不一样的秘钥来进行解密，并观察解密出来的内容是否有意义：

def analyze(c):
    """
    移位密码分析
    :param c: 密文c
    :return:
    """
    for k in range(26):
        # 用不一样的秘钥k尝试解密
        print('秘钥%d：' % k + Decrypt(k, c))


if __name__ == '__main__':
    c = 'jmpwfdpejoh'
    analyze(c)

运行程序输出结果为：

秘钥0：jmpwfdpejoh
秘钥1：ilovecoding
秘钥2：hknudbnchmf
秘钥3：gjmtcambgle
...........

逐行观察输出结果，到第二行的时候就能看到原来的明文，也就知道了你要对小红说的内容以及大家所约定的秘钥。面对你冒着巨大风险在课堂上所传递的纸条内容，老师内心可能也是复杂的...

Anyway，你的小秘密已经被老师知道了，此时比较灰心，一直在想，到底是什么缘由导致纸条计划失败？其实缘由很明显，各位也看出来了，小明所使用的加密体制中，可用的秘钥太少，或者说秘钥空间过小，别人直接一一列举进行穷搜就能破解，这就提示咱们：一个好的加密体制，它的秘钥空间应该是足够大的。

其实，你这次所用的移位密码是古典的加密体制之一，听说凯撒打仗时就用这种方法与将军们联系，因此位移密码也叫凯撒密码（Caesar cipher）。相似的还有代换密码，仿设射密码等等，都是将单个字母替换成别的字母，来达到加密的目的。报纸上的猜谜游戏就常常用这些方法，通常根据字母频率进行破解，有兴趣能够进行进一步的了解。

因此到底要用什么样的加密方法，才能保证我和小红的秘密不被人偷窥呢？

2.1 密码分析情形

俗话说，知己知彼，百战不殆，了解破解者的密码分析方法，或许可以帮助咱们想出更安全的密码体制。能够在不一样的情形下考察密码体制的安全性，通常咱们都假设破解者知道咱们所使用的密码体制，也就是说，不把密码体制的安全性寄托在密码体制的保密性上，而是放在秘钥上。

破解者的目的就是找出所使用的秘钥，常见的有如下几种攻击情形：

• 惟密文攻击： 破解者拥有密文c。这就是老师破解纸条的情形。
• 已知明文攻击： 破解者拥有一些明文p及其对应的密文c。考虑到实际情形，这个假设是比较合理的，例如破解者得到一封邮件加密后的密文，能够猜想一个词极可能是'hi'或者'dear'，这样就可能找到一个明文--密文对。
• 选择明文攻击： 破解者可以指定一个明文p，得到其对应的密文c，较强的假设。
• 选择密文攻击： 破解者指定一个密文c，得到其对应的明文，较强的假设。

天啊，你不由惊呼，在这么强的假设下，真的会有密码体制可以存活吗？

答案是有，并且这种密码体制已经被普遍应用，甚至能够说无处不在，它就是AES（Advanced Encryption Standard）。

3.SPN网络

难道不是要介绍AES吗，怎么会变成SPN网络，这是啥？能够吃吗？

AES、DES等不少现代对称加密方法的核心就是SPN网络，它是代换-置换网络（Substitution-Permutation Network）的缩写，是现代对称加密方法设计的蓝本。能够说，了解SPN网络，就基本了解了AES。

很巧的是，这个网络正好是容易理解的。SPN网络的思想很简单：既然加密一次不够安全，那我就加密屡次，把第一次加密产生的密文再进行加密，解密的时候我连续进行两次解密就能够了，这样是否是就安全了一些呢？

对于密码体制\(S_1\)，其加密与解密函数为\(E_1\)与\(D_1\)，对于密码体制\(S_2\)，其加密与解密函数为\(E_2\)与\(D_2\)，我构造出一个新的密码体制\(S_3\)，其加密函数为：$$c = E_2( E_1( p ) )$$

解密函数为$$p=D_1( D_2( c ) )$$

记为$$S_3 = S_1 * S_2$$这样破解\(S_3\)就可能会困难些。这个想法是否是很直接呢？这个思想在1949年才被提出，而提出者，可能理科生都多少听过他的名字——香农（Shannon）。

注意，不是任何的加密体制均可以这样“乘”起来变得更强，例如对于你的移位密码，嵌套起来仍是移位密码（为何？），没有任何改善，即\(S_1*S_1=S_1\)，这样的密码体制被称为幂等的。

若是密码体制不是幂等的，那么屡次迭代就可能可以提升安全性，SPN就是使用这种思想，包含多轮的迭代，每轮的操做都是相同的。下面，介绍SPN单轮的操做：

3.1 SPN单轮操做

SPN网络是对必定长度的比特进行操做的，在本文中的SPN网络中，一次加密的长度为16个比特，即2字节，也就是说每次加密16比特的明文，输出16比特的密文。

一个SPN网络包含多轮迭代，每轮迭代的操做内容都同样是：异或运算-->分组代换-->单比特置换

3.1.1 第一步——异或运算

异或运算是比较常见的二元比特运算，用⊕表示，其规则就是“相同得0，不一样得1”：

0 ⊕ 0 = 0
1 ⊕ 1 = 0
1 ⊕ 0 = 1
0 ⊕ 1 = 1

对于比特串，直接按每一位对应进行计算便可以了：

0011 ⊕ 1010 = 1001

异或的有比较有意思的性质：一个比特串亦或另外一个比特串两遍，仍是等于他本身，即a ⊕ b ⊕ b = a，这是由于a ⊕ b ⊕ b = a ⊕ ( b ⊕ b ) =a ⊕ 0 = a，能够带入一些例子试试看。

SPN网络中，每一轮的第一步就是把输入的比特串w和秘钥k进行亦或：u = w ⊕ k，如：

0001110000100011 = 0010011010110111 ⊕ 0011101010010100

这一步的目的是根据秘钥对明文进行混淆。若是你只知道输出u而不知道秘钥k，那么你就猜不出实际输入的w是什么，它是什么均可能，并且是等几率的。例如对于1 = a ⊕ b，不告诉你b是0仍是1，你就不知道a是什么。而对于和操做，若是知道1 = a and b，那么就能肯定a与b都是1。

这就是第一步，是否是很简单呢？

3.1.2 第二步——分组代换

这一步也很简单，将第一步输出的16比特的串分为4组，每组4比特，即0001110000100011写成0001 1100 0010 0011。而后对于每组再根据事先所定的表进行代换，代换表长这样：

代换前	0	1	2	3	4	5	6	7	8	9	A	B	C	D	E	F
代换后	E	4	D	1	2	F	B	8	3	A	6	C	5	9	0	7

就拿第一列来讲，表的意思是：若是你是0（0000），那么我要把你换成E（1110），就是一个简单的映射操做。

原比特串长这样：0001 1100 0010 0011 <==> 1 C 2 3，再对每一个字母查表获得：4 5 D 1 <==> 0100 0101 1101 0001，这样就获得代换后的比特串0100 0101 1101 0001，完成了第二步。

这个表通常称为S盒（Substitution），这个过程能够用v = S(u)表示，u是第一步异或的结果，也是第二步分组代换的输入，v是第二步的输出。须要注意，S盒的输入和输出通常是非线性的关系。

3.1.3 第三步——单比特置换

单比特置换是将16比特中的每一比特，根据P盒（Permutation）移动挪位，这样说很不直观，直接上例子，P盒长这样：

置换前的位置	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16
置换后的位置	1	5	6	13	2	6	10	14	3	7	11	15	4	8	12	16

拿第二列来讲，表的意思是：第二个比特要挪到第五个比特的位置，举个好看的例子：

0100 0000 0000 0000 置换后为==> 0000 1000 0000 0000

这个例子里面第二个比特的1挪到了第五的位置，而其余位置的比特都是0，挪位置以后仍是0。

对于第二部输出的结果1100 1101 1100 0100，置换后的比特串为0010 1110 0000 0111，这样就完成了第三步。

这一步能够用W = S(v)表示，v是第二部的输出，也是第三步的输入，W是第三步的输出，P盒置换是一种线性的变换。

这三步放在一块儿结果以下，建议读者本身计算一遍：

w = 0010 0110 1011 0111
k = 0011 1010 1001 0100
第一步，异或运算：
u = w ⊕ k = 0001 1100 0010 0011
第二步，分组代换：
v = S(u) = 0100 0101 1101 0001
第三步，单比特置换：
W = P(v) = 0010 1110 0000 0111

能够写成：W = P( S(w ⊕ k) )，这样就完成了一轮迭代，里面用到的参数有k，S盒与P盒，如图（图片来自维基百科）：

3.2 SPN的多轮迭代

弄清楚一轮的流程，SPN总体就很容易明白了，就是一轮一轮的乘起来，上一轮的输出做为这一轮的输入：

w0 = x
w1 = P(S(w0 ⊕ k1))
w2 = P(S(w1 ⊕ k2))
w3 = P(S(w2 ⊕ k3))
w4 = P(S(w3 ⊕ k4))
y = w4

w0就是16比特的明文，w4是4轮操做后的16比特密文结果，是否是很简单？须要注意的是，每一轮迭代的秘钥k是不同的，通常是由一个基础秘钥经特定秘钥编排算法生成的，而使用的S盒P盒都是相同的，会提早肯定好，而且是公开的。
下图是一个三轮SPN网络的示意图（图片来自维基百科）：

注意在最后一轮去掉了代换操做，这样作可使加密算法稍微作一些调整就能够用来进行解密。

OK! SPN网络就是这些内容，你已经掌握了它，若是你还想和小红传纸条的话，能够试试用它加密，会比移位密码更安全一些。

什么？本身手动代换置换太麻烦？不用怕，贴心的我已经为你准备好了Python代码。

3.3 用Python实现SPN网络

我实现的是4轮迭代的SPN网络，以及加密和解密算法，其结构图以下（图片来自 Cryptography Theory and Practice ）：

每次加密输入16比特的明文，输出16比特的密文，代码以下：

# S盒参数
S_Box = [14, 4, 13, 1, 2, 15, 11, 8, 3, 10, 6, 12, 5, 9, 0, 7]

# P盒参数
P_Box = [1, 5, 9, 13, 2, 6, 10, 14, 3, 7, 11, 15, 4, 8, 12, 16]


def gen_K_list(K):
    """
    秘钥编排算法，由一个32比特秘钥生成5个16比特子秘钥
    :param K: 32比特秘钥
    :return: [k1,k2,k3,k4,k5]，五个16比特子秘钥
    """
    Ks = []
    for i in range(5, 0, -1):
        ki = K % (2 ** 16)
        Ks.insert(0, ki)
        K = K >> 4
    return Ks


def pi_s(s_box, ur):
    """
    分组代换操做
    :param s_box:S盒参数
    :param ur:输入比特串，16比特
    :return:输出比特串，16比特
    """
    vr = 0
    for i in range(4):
        uri = ur % (2 ** 4)
        vri = s_box[uri]
        vr = vr + (vri << (4 * i))
        ur = ur >> 4
    return vr


def pi_p(p_box, vr):
    """
    单比特置换操做
    :param p_box:P盒参数
    :param vr:输入比特串，16比特
    :return:输出比特串，16比特
    """
    wr = 0
    for i in range(15, -1, -1):
        vri = vr % 2
        vr = vr >> 1
        wr = wr + (vri << (16 - p_box[i]))
    return wr


def reverse_Sbox(s_box):
    """
    求S盒的逆
    :param s_box:S盒参数
    :return:S盒的逆
    """
    re_box = [-1] * 16
    for i in range(16):
        re_box[s_box[i]] = i
    return re_box


def reverse_Pbox(p_box):
    """
    求P盒的逆
    :param s_box:P盒参数
    :return:P盒的逆
    """
    re_box = [-1] * 16
    for i in range(16):
        re_box[p_box[i] - 1] = i + 1
    return re_box


def do_SPN(x, s_box, p_box, Ks):
    """
    5轮的SPN网络，能够用来进行加密或解密
    :param x: 16比特输入
    :param s_box: S盒参数
    :param p_box: P盒参数
    :param Ks: [k1,k2,k3,k4,k5]，五个16比特子秘钥
    :return: 16比特输出
    """
    wr = x
    for r in range(3):
        ur = wr ^ Ks[r]  # 异或操做
        vr = pi_s(s_box, ur)  # 分组代换
        wr = pi_p(p_box, vr)  # 单比特置换

    ur = wr ^ Ks[3]
    vr = pi_s(s_box, ur)
    y = vr ^ Ks[4]
    return y


def encrypt(K, x):
    """
    根据秘钥K对16比特明文x进行加密
    :param K:32比特秘钥
    :param x:16比特明文
    :return:16比特密文
    """
    Ks = gen_K_list(K)
    return do_SPN(x, S_Box, P_Box, Ks)


def decrypt(K, y):
    """
    根据秘钥K对16比特密文y进行解密。
    :param K:32比特秘钥
    :param y:16比特密文
    :return:16比特明文
    """
    Ks = gen_K_list(K)
    Ks.reverse()  # 秘钥逆序编排
    # 秘钥置换
    Ks[1] = pi_p(P_Box, Ks[1])
    Ks[2] = pi_p(P_Box, Ks[2])
    Ks[3] = pi_p(P_Box, Ks[3])

    s_rbox = reverse_Sbox(S_Box)  # S盒求逆
    p_rbox = reverse_Pbox(P_Box)  # P盒求逆
    return do_SPN(y, s_rbox, p_rbox, Ks)


if __name__ == '__main__':
    x = 0b0010011010110111
    K = 0b00111010100101001101011000111111
    print('初始明文：', format(x, '016b'))
    print('加密密文：', format(encrypt(K, x), '016b'))
    print('解密结果：', format(decrypt(K, encrypt(K, x)), '016b'))
    assert decrypt(K, encrypt(K, x)) == x

能够直接看do_SPN函数，函数里面循环3次，对应3轮迭代，第4轮迭代没有置换操做。encrypt与decrypt函数调用do_SPN函数便可进行加密和解密操做（为何能够调用SPN进行解密？能够对照代码观察SPN的结构想想），运行程序输出为：

初始明文： 0010011010110111
加密密文： 1011110011010110
解密结果： 0010011010110111

至此，SPN网络已经彻底实现！那么它的安全性如何呢？

首先，咱们知道，这个SPN网络的秘钥是32位的，大约是有4百万的候选秘钥，这个数量的秘钥，手动穷搜是很难的，用计算机来穷搜就会比较容易了，不过咱们随时对它进行改造，增长秘钥长度，如256位，这时候机器穷搜也不行了。

其次，若是SPN层数较少，对其进行线性分析或者差分分析一般会比较容易成功，这些分析方法的代码我也已经写好了，不过今天的内容已经够多了，咱们下次再介绍它:)