
在自然语言处理、生物信息学等众多领域中,我们常常需要比较两个字符串的相似程度。例如,在拼写检查时,需要判断输入的错误单词和正确单词之间的差异;在基因序列分析中,要比较不同生物的基因序列相似度。编辑距离(Edit Distance),也称为莱文斯坦距离(Levenshtein Distance),就是一种衡量两个字符串差异程度的有效方法。本文将深入探讨如何使用动态规划算法来计算字符串的编辑距离。
编辑距离是指两个字符串之间,由一个字符串转换成另一个字符串所需的最少编辑操作次数。允许的编辑操作包括:
例如,将字符串“kitten”转换为“sitting”,可以通过以下步骤实现:
动态规划(Dynamic Programming)是一种通过把原问题分解为相对简单的子问题,并保存子问题的解来避免重复计算,从而解决复杂问题的方法。计算编辑距离就可以通过动态规划算法来实现,具体步骤如下:
设 dp[i][j] 表示字符串 s1 的前 i 个字符和字符串 s2 的前 j 个字符之间的编辑距离。
i = 0 时,dp[0][j] = j,表示将空字符串转换为 s2 的前 j 个字符需要 j 次插入操作。j = 0 时,dp[i][0] = i,表示将 s1 的前 i 个字符转换为空字符串需要 i 次删除操作。对于 dp[i][j],有以下三种情况:
s1[i - 1] == s2[j - 1],则 dp[i][j] = dp[i - 1][j - 1],因为当前字符相同,不需要进行额外的编辑操作。s1[i - 1]!= s2[j - 1],则 dp[i][j] 取以下三种操作的最小值:dp[i][j - 1] + 1,表示在 s1 的前 i 个字符后面插入一个字符,使其与 s2 的前 j 个字符匹配。dp[i - 1][j] + 1,表示删除 s1 的第 i 个字符。dp[i - 1][j - 1] + 1,表示将 s1 的第 i 个字符替换为 s2 的第 j 个字符。用公式表示为:
[
dp[i][j] =
\begin{cases}
dp[i - 1][j - 1] & \text{if } s1[i - 1] == s2[j - 1] \
\min(dp[i][j - 1] + 1, dp[i - 1][j] + 1, dp[i - 1][j - 1] + 1) & \text{if } s1[i - 1]!= s2[j - 1]
\end{cases}
]
dp[m][n] 即为字符串 s1 和 s2 的编辑距离,其中 m 和 n 分别是 s1 和 s2 的长度。
def edit_distance(s1, s2):m, n = len(s1), len(s2)# 初始化 dp 数组dp = [[0] * (n + 1) for _ in range(m + 1)]# 初始化边界条件for i in range(m + 1):dp[i][0] = ifor j in range(n + 1):dp[0][j] = j# 填充 dp 数组for i in range(1, m + 1):for j in range(1, n + 1):if s1[i - 1] == s2[j - 1]:dp[i][j] = dp[i - 1][j - 1]else:dp[i][j] = min(dp[i][j - 1] + 1, # 插入dp[i - 1][j] + 1, # 删除dp[i - 1][j - 1] + 1) # 替换return dp[m][n]# 测试s1 = "kitten"s2 = "sitting"print(f"字符串 '{s1}' 和 '{s2}' 的编辑距离为: {edit_distance(s1, s2)}")
m 和 n 分别是两个字符串的长度。因为需要填充一个 $m+1$ 行 $n+1$ 列的二维数组。dp 数组。不过,可以通过优化将空间复杂度降低到 $O(min(m, n))$。| 要点 | 详情 |
|---|---|
| 编辑距离定义 | 两个字符串之间,由一个字符串转换成另一个字符串所需的最少编辑操作(插入、删除、替换)次数 |
| 动态规划思路 | 定义状态 dp[i][j] 表示 s1 前 i 个字符和 s2 前 j 个字符的编辑距离,通过初始化边界条件和状态转移方程求解 |
| 代码实现 | 使用 Python 语言,通过二维数组 dp 实现动态规划过程 |
| 复杂度 | 时间复杂度 $O(m n)$,空间复杂度 $O(m n)$ |
编辑距离是一种非常实用的字符串相似度度量方法,通过动态规划算法可以高效地计算出来。在实际应用中,我们可以根据编辑距离来进行拼写检查、文本相似度匹配等任务。希望本文能帮助你理解和掌握动态规划算法计算字符串编辑距离的方法。