
本教程深入探讨go语言中最长公共子序列(lcs)的回溯算法实现,重点解决在动态规划表格与字符串索引之间常见的偏移问题。文章详细分析了索引不匹配导致的问题,并提供了一个经过修正的go语言回溯函数,确保lcs能够被正确地重构。通过完整的代码示例和注意事项,帮助读者掌握lcs回溯的正确实践。
引言:最长公共子序列与回溯
最长公共子序列(Longest Common Subsequence, LCS)是计算机科学中一个经典的动态规划问题。给定两个序列,LCS的目标是找到一个最长的序列,该序列是这两个给定序列的子序列。LCS的长度计算通常通过构建一个二维动态规划(DP)表格来完成。然而,仅仅计算出LCS的长度是不够的,在许多实际应用中,我们还需要找出具体的LCS序列本身。这就需要用到“回溯”技术,即根据DP表格中的值逆向追踪,从而重构出LCS。
在go语言中实现LCS回溯时,一个常见的陷阱是处理DP表格索引与原始字符串索引之间的对应关系。如果处理不当,可能导致回溯结果不完整或错误。本教程将详细解析这一问题,并提供一个健壮且正确的Go语言回溯实现。
LCS长度计算的动态规划基础
在深入回溯之前,我们首先回顾LCS长度的动态规划计算。假设我们有两个字符串 str1 和 str2,长度分别为 m 和 n。我们通常会构建一个 (m+1) x (n+1) 的二维表格 dp。其中,dp[i][j] 表示 str1 的前 i 个字符与 str2 的前 j 个字符的最长公共子序列的长度。
表格填充规则如下:
立即学习“go语言免费学习笔记(深入)”;
- 初始化: dp[0][j] = 0 和 dp[i][0] = 0,表示空字符串与任何字符串的LCS长度为0。
- 递推:
- 如果 str1[i-1] == str2[j-1] (注意这里是字符串的0-based索引,对应DP表格的1-based索引),则 dp[i][j] = 1 + dp[i-1][j-1]。
- 如果 str1[i-1] != str2[j-1],则 dp[i][j] = max(dp[i-1][j], dp[i][j-1])。
最终,dp[m][n] 将包含 str1 和 str2 的LCS长度。
以下是一个Go语言实现LCS长度计算的示例:
func LCSLength(str1, str2 string) [][]int { m := len(str1) n := len(str2) // dp table is (m+1) x (n+1) dp := make([][]int, m+1) for i := range dp { dp[i] = make([]int, n+1) } for i := 1; i <= m; i++ { for j := 1; j <= n; j++ { if str1[i-1] == str2[j-1] { dp[i][j] = 1 + dp[i-1][j-1] } else { dp[i][j] = max(dp[i-1][j], dp[i][j-1]) } } } return dp } func max(a, b int) int { if a > b { return a } return b }
LCS回溯算法原理与常见陷阱
回溯LCS序列的核心思想是从DP表格的右下角(dp[m][n])开始,逆向追踪路径。
- 匹配字符: 如果 str1[i-1] == str2[j-1],这意味着当前字符是LCS的一部分。我们将 str1[i-1](或 str2[j-1])添加到LCS序列中,并向左上方移动 (i-1, j-1)。
- 不匹配字符: 如果 str1[i-1] != str2[j-1],这意味着当前字符不是LCS的一部分。我们需要查看 dp[i-1][j] 和 dp[i][j-1] 的值,选择其中较大的那个方向移动。
- 如果 dp[i-1][j] > dp[i][j-1],则向上移动 (i-1, j)。
- 如果 dp[i][j-1] >= dp[i-1][j],则向左移动 (i, j-1)。
- 终止条件: 当 i 或 j 达到0时,回溯结束。
常见陷阱:索引偏移问题
问题通常出现在 Back 函数中对字符串和DP表格索引的混淆。DP表格是 (m+1) x (n+1) 大小的,其索引 i 和 j 通常代表 str1 的前 i 个字符和 str2 的前 j 个字符。这意味着 dp[i][j] 对应的是 str1[0…i-1] 和 str2[0…j-1]。
如果回溯函数 Back(table, str1, str2, i, j) 中的 i 和 j 被设计为直接代表字符串的0-based索引(即 str1[i] 和 str2[j]),那么在访问DP表格时就需要进行调整。
原始代码的错误在于:
- 它将 i 和 j 视为字符串的0-based索引,因此 str1[i] 和 str2[j] 是合法的。
- 然而,它的基本情况 if i == 0 || j == 0 却将 0 视为字符串的起始边界,导致 str1[0] 和 str2[0] 可能被跳过。正确的字符串0-based索引的终止条件应该是 i == -1 || j == -1。
- 在不匹配的情况下,它使用了 table[i][j-1] 和 table[i-1][j]。如果 i 和 j 是字符串索引,那么对应的DP表格索引应该是 i+1 和 j+1。因此,要比较 str1 移除 str1[i] 后与 str2 的LCS长度(对应 table[i][j+1])和 str2 移除 str2[j] 后与 str1 的LCS长度(对应 table[i+1][j]),而不是 table[i][j-1] 和 table[i-1][j]。
正确的Go语言回溯实现
为了解决上述索引偏移问题,我们需要确保 Back 函数中的 i 和 j 始终表示字符串的0-based索引,并在访问DP表格时进行相应的调整。
以下是修正后的 Back 函数实现:
// Back 函数用于从LCS长度DP表格中回溯出最长公共子序列 // table: 预先计算好的LCS长度DP表格 (m+1)x(n+1) // str1, str2: 原始字符串 // i, j: 当前回溯点在字符串中的0-based索引 func Back(table [][]int, str1, str2 string, i, j int) string { // 基本情况:当任一字符串索引小于0时,表示已遍历完该字符串,返回空字符串 if i == -1 || j == -1 { return "" } // 如果当前字符匹配 (str1[i] == str2[j]) // 那么这个字符是LCS的一部分,将其添加到结果中,并向左上方移动 (i-1, j-1) if str1[i] == str2[j] { return Back(table, str1, str2, i-1, j-1) + string(str1[i]) } else { // 如果当前字符不匹配 // 比较移除 str1[i] 后的LCS长度 (table[i][j+1]) // 和移除 str2[j] 后的LCS长度 (table[i+1][j]) // 注意:table的索引比字符串索引大1 // table[i][j+1] 对应 str1[:i] 和 str2[:j+1] 的LCS长度 // table[i+1][j] 对应 str1[:i+1] 和 str2[:j] 的LCS长度 if table[i+1][j] > table[i][j+1] { // 如果移除 str2[j] 后的LCS更长,则向左移动 (即在str2中向前移动,j-1) return Back(table, str1, str2, i, j-1) } else { // 否则(移除 str1[i] 后的LCS更长或相等),则向上移动 (即在str1中向前移动,i-1) return Back(table, str1, str2, i-1, j) } } }
代码解释:
- i == -1 || j == -1: 这是正确的字符串0-based索引的终止条件。当 i 或 j 变为 -1 时,意味着我们已经处理了 str1 或 str2 的所有字符。
- str1[i] == str2[j]: 当字符匹配时,我们将 str1[i] 加入结果,并递归调用 Back(table, str1, str2, i-1, j-1)。这里的 i 和 j 仍然是字符串索引。
- table[i+1][j] > table[i][j+1]: 这是关键的修正点。
- i 和 j 是字符串的0-based索引。
- table[x][y] 存储的是 str1 的前 x 个字符和 str2 的前 y 个字符的LCS长度。
- 当比较 str1[i] 和 str2[j] 不匹配时,我们需要看:
- 忽略 str1[i]:LCS长度为 str1 的前 i 个字符与 str2 的前 j+1 个字符的LCS长度,这对应 table[i][j+1]。
- 忽略 str2[j]:LCS长度为 str1 的前 i+1 个字符与 str2 的前 j 个字符的LCS长度,这对应 table[i+1][j]。
- 因此,比较 table[i+1][j] 和 table[i][j+1] 是正确的。如果 table[i+1][j] 更大,说明不包含 str2[j] 的路径产生了更长的LCS,所以我们递归调用 Back(table, str1, str2, i, j-1)(在 str2 中向前移动)。反之,则递归调用 Back(table, str1, str2, i-1, j)(在 str1 中向前移动)。
完整示例代码
为了演示如何使用 LCSLength 和 Back 函数,我们提供一个完整的示例:
package main import ( "fmt" "strings" ) // max returns the maximum of two integers. func max(a, b int) int { if a > b { return a } return b } // LCSLength computes the length of the Longest Common Subsequence // and returns the DP table. func LCSLength(str1, str2 string) [][]int { m := len(str1) n := len(str2) // dp table is (m+1) x (n+1) dp := make([][]int, m+1) for i := range dp { dp[i] = make([]int, n+1) } for i := 1; i <= m; i++ { for j := 1; j <= n; j++ { if str1[i-1] == str2[j-1] { dp[i][j] = 1 + dp[i-1][j-1] } else {