多表古典密码统计分析

一、前言

1. 多表古典密码分析

我们以Vigenere密码为例来说明多表古典密码的分析方法。确定密钥字长度的方法有Kasiski测试法（Kasiski Test）和重合指数法（index of coincidence).

2. Vigenere密码

• Vigenere（维吉尼亚）密码是使用一系列凯撒密码组成密码字母表的加密算法，属于多表密码的一种简单形式。
• 凯撒密码：一种替换加密的技术，明文中的所有字母都在字母表上向后（或向前）按照一个固定数目进行偏移后被替换成密文。例如，当偏移量是3的时候，所有的字母A将被替换成D，B变成E，以此类推。

二、Vigenere密码分析

1. Kasiski测试法

寻找密文中相同的片段对(一般长度大于三)，计算每对相同密文片段对之间的距离，不妨记为$d_1,d_2,\dots ,d_i$，若令密钥字的长度为$m$，则$m=gcd(d_1,d_2,\dots ,d_i)$

2. 重合指数法

设$X=x_1{x}_2\dots x_n$是一个长度为n的英文字母串，则x中任意选取两个字母相同的概率定义为重合指数，用 $I_c(x)$ 表示。
设英文字母$A，B,\dots ，Z$在$X$中出现的次数分别为：$f_0,f_1,\dots ,f_{25}$则从$X$中任意选取两个字母相同的概率为
$$I_c(X)=\frac{{\sum }_{i=0}^{25}{f}_i\left(f_i-1\right)}{n(n-1)}$$
下表为26个英文字母的出现频率：

已知每个英文字母出现的期望概率，分别记为$p_0,p_1,\dots ,p_{25}$，那么X中两个元素相同的概率计算为
$$I_c(X)\approx \sum _{i=0}^{25}{p}_i^2$$
结果约等于 0.065，若是均匀分布，则值约等于0.038。

三、实际密码分析

1. 问题描述

设某一段明文经过$Vigenere$密码加密后密文为
$CHREEVOAHMAERATBIAXXWTNXBEEOPHBSBQMQEQERBW$
$RVXUOAKXAOSXXWEAHBWGJMMQMNKGRFVGXWTRZXWIAK$
$LXFPSKAUTEMNDCMGTSXMXBTUIADNGMGPSRELXNJELX$
$VRVPRTULHDNQWTWDTYGBPHXTFALJHASVBFXNGLLCHR$
$ZBWELEKMSJIKNBHWRJGNMGJSGLXFEYPHAGNRBIEQJT$
$AMRVLCRREMNDGLXRRIMGNSNRWCHRQHAEYEVTAQEBBI$
$PEEWEVKAKOEWADREMXMTBJJCHRTKDNVRZCHRCLQOHP$
$WQAIIWXNRMGWOIIFKEE$
通过两种分析方法计算出此密码的密钥。

2. 问题解决

1. 密文片段CHR一共出现了5次，每次开始位置分别为1，166，236，276，286，第一次出现到各次出现的距离$d_i$分别为165，235，275，285，所以$m=gcd(165，235，275，285)=5$
2. 根据Kasiski测试法得到的$mv$，可以将密文$Yv$按照下列形式排列：
将$Y$排列成$m$行$n/m$列的形式，形如：
$$\begin{array}{ccccc}{k}_1& y_1& y_{m+1}& \cdots & y_{n-m+1}\\ k_2& y_2& y_{m+2}& \cdots & y_{n-m+2}\\ k_3& y_3& y_{m+3}& \cdots & y_{n-m+3}\\ ⋮& ⋮& ⋮& ⋮& ⋮\\ k_m& y_m& y_{2m}& \cdots & y_n\end{array}$$
设f0,f1,…,f25分别表示串yi中英文字母A，B,…，Z出现的次数，令n’=n/m表示串yi的长度，则26个英文字母在yi的概率分布为：
$$\frac{f_0}{n},\frac{f_1}{n},\frac{f_2}{n},\dots \dots ..\frac{f_{25}}{n}$$

考虑到$Vigenere$密码子串$y_i$是有对应的明文子集中的字母移动$K_i$个位置所得，因此，移动后的概率分布为：
$$\frac{f_{ki}}{n^{\prime }},\frac{f_{1+ki}}{n^{\prime }},\frac{f_{2+ki}}{n^{\prime }},\dots \dots \dots ,\frac{f_{25+ki}}{n^{\prime }}$$

$$M_g=\sum _{i=0}^{25}\frac{p_i{f}_{i+g}}{n^{\prime }}$$

每一行计算每一种字母的$M_g$值（$g$代表0-25即$a-z$),若$g$对应的$M_g$ 值与0.065接近，那么$g$则为密钥

四、编程

1. 代码

2. 运行结果

程序比较简单，但是代码比较长，我只把结果贴出来比对一下：