Python实现替换密码频率分析与破解工具

1. 替换密码的基本原理与破解思路

替换密码是最古老的加密方式之一，其核心思想是将明文中的每个字母按照固定的规则替换为另一个字母。举个生活中的例子，就像我们把手机键盘上的字母对应到数字键上（2对应ABC，3对应DEF等），只不过替换密码的映射关系是人为设定的。

在密码学中，这种加密方式被称为单表替换密码（Monoalphabetic Cipher）。它的典型特点是：

• 每个明文字母唯一对应一个密文字母
• 字母频率分布特征保持不变
• 不改变字母的位置顺序

要破解这类密码，我们需要利用英文文本的统计特征。根据语言学分析，英文中字母的出现频率从高到低大致是：E, T, A, O, I, N... 这个规律就像DNA一样具有辨识度。当我们拿到足够长的密文时，统计其中各字母的出现频率，就能推测出可能的对应关系。

2. Python实现频率分析工具

2.1 基础频率统计

我们先构建一个基础的频率分析工具。以下代码可以统计文本中字母的出现频率：

python复制from collections import Counter
import string

def frequency_analysis(ciphertext):
    # 过滤非字母字符并转为大写
    filtered = [c.upper() for c in ciphertext if c.isalpha()]
    total = len(filtered)
    
    # 统计字母频率
    freq = Counter(filtered)
    
    # 计算百分比
    for char in freq:
        freq[char] = round(freq[char]/total * 100, 2)
    
    # 按频率排序
    sorted_freq = sorted(freq.items(), key=lambda x: x[1], reverse=True)
    
    return dict(sorted_freq)

# 示例用法
ciphertext = "GWC uivg qc q iuwt qh gwc bsoos..."
print(frequency_analysis(ciphertext))

这段代码会输出类似这样的结果：

code复制{'G': 15.38, 'W': 11.54, 'C': 11.54, 'U': 7.69, ...}

2.2 增强型频率分析

单纯的字母频率有时不够准确，我们还可以加入以下分析维度：

1. 双字母组合频率：统计常见的字母对（如TH, HE, IN等）
2. 单词首尾字母频率：单词开头常见字母（T, A, S等）和结尾常见字母（E, S, D等）
3. 单字母单词识别：英文中单字母单词通常是"A"或"I"

改进后的分析函数：

python复制def enhanced_analysis(ciphertext):
    words = ciphertext.split()
    single_letters = [w for w in words if len(w) == 1]
    
    # 首字母统计
    first_letters = [w[0] for w in words if w]
    first_freq = Counter(first_letters)
    
    # 双字母统计
    bigrams = [ciphertext[i:i+2] for i in range(len(ciphertext)-1)]
    bigram_freq = Counter(bigrams)
    
    return {
        'single_letters': Counter(single_letters),
        'first_letters': first_freq.most_common(5),
        'top_bigrams': bigram_freq.most_common(5)
    }

3. 交互式破解工具实现

3.1 密码映射管理

我们需要一个可以交互式修改字母映射的系统：

python复制class SubstitutionSolver:
    def __init__(self, ciphertext):
        self.ciphertext = ciphertext
        self.mapping = {c: None for c in string.ascii_uppercase}
        self.plaintext = []
        
    def update_mapping(self, cipher_char, plain_char):
        self.mapping[cipher_char.upper()] = plain_char.lower()
        
    def decrypt(self):
        result = []
        for char in self.ciphertext:
            if char.upper() in self.mapping and self.mapping[char.upper()]:
                result.append(self.mapping[char.upper()])
            else:
                result.append(char)
        return ''.join(result)

3.2 可视化界面

使用Python的curses库创建简单命令行界面：

python复制import curses

def interactive_solver(ciphertext):
    solver = SubstitutionSolver(ciphertext)
    
    # 初始化curses
    stdscr = curses.initscr()
    curses.noecho()
    curses.cbreak()
    stdscr.keypad(True)
    
    try:
        while True:
            stdscr.clear()
            
            # 显示当前解密结果
            stdscr.addstr(0, 0, "Current decryption:")
            stdscr.addstr(1, 0, solver.decrypt())
            
            # 显示频率分析
            freq = frequency_analysis(ciphertext)
            stdscr.addstr(3, 0, "Frequency analysis:")
            for i, (char, percent) in enumerate(freq.items()):
                stdscr.addstr(4+i, 0, f"{char}: {percent}%")
            
            # 显示当前映射
            stdscr.addstr(15, 0, "Current mapping:")
            for i, (cipher, plain) in enumerate(solver.mapping.items()):
                stdscr.addstr(16+i//5, (i%5)*15, f"{cipher}->{plain if plain else '_'}")
            
            # 获取用户输入
            stdscr.addstr(20, 0, "Enter cipher char then plain char (or q to quit):")
            input_str = stdscr.getstr(21, 0).decode('utf-8')
            
            if input_str.lower() == 'q':
                break
                
            if len(input_str) == 2:
                solver.update_mapping(input_str[0], input_str[1])
                
    finally:
        curses.nocbreak()
        stdscr.keypad(False)
        curses.echo()
        curses.endwin()
    
    return solver.decrypt()

4. 实战破解案例演示

让我们用以下密文演示破解过程：

code复制"QEB NRFZH YOLTK CLU GRJMP LSBO QEB IXWV ALD"

4.1 初始频率分析

运行频率分析函数得到：

code复制{'B': 11.11, 'E': 8.33, 'Q': 8.33, 'L': 8.33, ...}

4.2 逐步推理

1. 观察到单字母单词"Q"和"B"：

   • 英文中单字母单词通常是"A"或"I"
   • 假设Q->A，B->I

2. 三字母单词"QEB"：

   • 如果Q->A，B->I，则可能是"AND"或"THE"
   • 尝试E->H

3. 更新映射后部分解密：
   "A H_NFZH YOLTK CLU GRJMP LSIO A IXWV ALD"

4. 观察"ALD"：

   • 可能是"AND"，尝试L->N，D->D

5. 继续这个过程，最终得到完整映射

4.3 完整解密结果

最终破解出的明文是：

code复制"THE QUICK BROWN FOX JUMPS OVER THE LAZY DOG"

这是著名的包含所有英文字母的测试句子。

5. 算法优化与进阶技巧

5.1 自动化猜测

我们可以基于词频数据自动生成初始猜测：

python复制english_freq = ['E','T','A','O','I','N','S','H','R','D','L','C','U','M','W','F','G','Y','P','B','V','K','J','X','Q','Z']

def auto_initial_guess(ciphertext):
    cipher_freq = frequency_analysis(ciphertext)
    mapping = {}
    
    for (cipher_char, _), eng_char in zip(cipher_freq.items(), english_freq):
        mapping[cipher_char] = eng_char.lower()
    
    return mapping

5.2 字典验证

使用英语字典验证解密结果的合理性：

python复制import requests

def check_english_words(text, threshold=0.7):
    words = text.split()
    english_words = set(requests.get('https://raw.githubusercontent.com/dwyl/english-words/master/words.txt').text.splitlines())
    
    match_count = sum(1 for word in words if word.lower() in english_words)
    return match_count / len(words) > threshold

5.3 处理标点和大小写

改进解密函数保留原始格式：

python复制def decrypt_preserve_format(ciphertext, mapping):
    result = []
    for char in ciphertext:
        upper_char = char.upper()
        if upper_char in mapping and mapping[upper_char]:
            if char.isupper():
                result.append(mapping[upper_char].upper())
            else:
                result.append(mapping[upper_char].lower())
        else:
            result.append(char)
    return ''.join(result)

6. 密码学安全思考

虽然我们成功破解了简单替换密码，但需要明白：

1. 现代加密算法（如AES）使用更复杂的原理
2. 增加密钥空间（如多表替换的Vigenère密码）能提高安全性
3. 即使对于替换密码，缩短密文长度会增加破解难度

在实际应用中，应该使用经过验证的加密库，而不是自己实现加密算法。Python中推荐使用cryptography等专业库。

这个项目的主要价值在于：

• 理解基础密码学原理
• 学习频率分析技术
• 提升Python文本处理能力
• 培养系统性解决问题思维