当前位置：首页 > wzjs >正文

网站建站和项目部署一样吗微信公众号怎么制作网页

wzjs 2025/9/5 19:11:59

网站建站和项目部署一样吗,微信公众号怎么制作网页,表格制作方法,资源网站推荐1. 理论原理推导假设有两个序列： 序列序列子问题定义： 令表示序列和的最长公共子序列的长度。递推公式： 对于任意和 ，我们考虑两个位置上的元素和 ： 如果 ，说明这两个元素可以同时…

1. 理论原理推导

假设有两个序列：

序列 $X = x_1, x_2, \dots, x_m$
序列 $Y = y_1, y_2, \dots, y_n$

子问题定义：
令 $dp[i][j]$ 表示序列 $X[1 \dots i]$ 和 $Y[1 \dots j]$ 的最长公共子序列的长度。

递推公式：
对于任意 $1 \le i \le m$ 和 $1 \le j \le n$ ，我们考虑两个位置上的元素 $x_i$ 和 $y_j$ ：

如果 $x_i = y_j$ ，说明这两个元素可以同时出现在公共子序列中，此时：
$dp[i][j] = dp[i-1][j-1] + 1$
如果 $x_i \neq y_j$ ，那么最长公共子序列不可能同时包含这两个元素，只能在两种可能中取最大值：
$dp[i][j] = \max\{ dp[i-1][j],\; dp[i][j-1] \}$

边界条件：

当 i=0 或 j=0 时，意味着其中一个序列为空，此时 $dp[i][j]=0$ 。

这一推导基于最优子结构性质：最长公共子序列问题的最优解可以由其子问题的最优解构造得到。

2. 时间复杂度推导

状态个数：
动态规划表 dp 的大小为 (m+1)×(n+1) 。
状态转移：
每个状态 $dp[i][j]$ 的计算只涉及常数时间的比较和加法操作。

因此，总的时间复杂度为：O(m×n)

空间复杂度也为 O(m×n) ；若通过空间优化（例如仅保留当前行和上一行），可降为 $O(\min(m, n))$ 。

3. 算法步骤

下面是基于上述理论推导的动态规划求解 LCS 的具体步骤：

初始化：
- 构建一个大小为 (m+1)×(n+1)的二维数组 dp，所有元素初始化为 0。
- 行和列的索引从 0 开始，令 $dp[0][j]=0$ （对于所有 j）和 $dp[i][0]=0$ （对于所有 i）。
状态转移：
- 对于 i 从 1 到 m：
  - 对于 j 从 1 到 n：
    - 若 $x_i = y_j$ ，则令：
      $dp[i][j] = dp[i-1][j-1] + 1$
    - 否则：
      $dp[i][j] = \max\{ dp[i-1][j],\; dp[i][j-1] \}$
结果输出：
- 最终， $dp[m][n]$ 即为序列 X 和 Y 的最长公共子序列的长度。
（可选）路径回溯：
- 若需要重构出一个最长公共子序列，可以从 $dp[m][n]$ 开始，逆向回溯：
  - 当 $x_i = y_j$ 时，该元素属于 LCS，记录下来，并同时减小 i 和 j；
  - 当 $dp[i-1][j] \ge dp[i][j-1]$ 时，移动到 $(i-1, j)$ ；否则移动到 $(i, j-1)$ 。

4. Python 代码示例

下面提供一个 Python 代码示例，既计算 LCS 的长度，也重构出一个 LCS：

def lcs_length(X, Y):m, n = len(X), len(Y)# 创建 (m+1) x (n+1) 的二维 dp 数组，并初始化为 0dp = [[0] * (n + 1) for _ in range(m + 1)]# 填充 dp 表for i in range(1, m + 1):for j in range(1, n + 1):if X[i - 1] == Y[j - 1]:dp[i][j] = dp[i - 1][j - 1] + 1else:dp[i][j] = max(dp[i - 1][j], dp[i][j - 1])return dpdef construct_lcs(dp, X, Y):"""根据 dp 表从 X 和 Y 中回溯构造最长公共子序列"""i, j = len(X), len(Y)lcs = []while i > 0 and j > 0:# 如果两个元素相等，加入 LCSif X[i - 1] == Y[j - 1]:lcs.append(X[i - 1])i -= 1j -= 1else:# 否则，沿着较大值的方向回溯if dp[i - 1][j] >= dp[i][j - 1]:i -= 1else:j -= 1lcs.reverse()  # 反转得到正确顺序return lcsdef main():# 示例序列，可以根据需要替换成任意两个序列X = "ABCBDAB"Y = "BDCAB"# 计算 dp 表dp = lcs_length(X, Y)print("最长公共子序列的长度为:", dp[len(X)][len(Y)])# 重构最长公共子序列lcs_seq = construct_lcs(dp, X, Y)print("最长公共子序列为:", ''.join(lcs_seq))if __name__ == '__main__':main()