NOJ1041——算法实验四——最长公共子序列

最新推荐文章于 2020-11-24 01:46:48 发布

Alfred young

最新推荐文章于 2020-11-24 01:46:48 发布

阅读量778

点赞数

分类专栏：算法设计与分析

本文链接：https://blog.csdn.net/qq_45757722/article/details/109552273

版权

算法设计与分析专栏收录该内容

7 篇文章 6 订阅

订阅专栏

最长公共子序列问题

问题描述

描述：

一个给定序列的子序列是在该序列中删去若干元素后得到的序列。确切地说，若给定序列X=<x1, x2,…, xm>，则另一序列Z=<z1, z2,…, zk>是X的子序列是指存在一个严格递增的下标序列 <i1, i2,…, ik>，使得对于所有j=1,2,…,k有：

Xij ＝ Zj

如果一个序列S即是A的子序列又是B的子序列，则称S是A、B的公共子序列。
求A、B所有公共子序列中最长的序列的长度。

输入：

输入共两行，每行一个由字母和数字组成的字符串，代表序列A、B。A、B的长度不超过200个字符。

输出：

一个整数，表示最长各个子序列的长度。
格式：printf("%d\n");

输入样例：

programming
contest

输出样例：

问题分析

设序列
$X=\{x_1,x_2,\cdots,x_m\}$

$Y=\{y_1,y_2,\cdots,y_n\}$

$Z=\{z_1,z_2,\cdots,z_k\}$

若 $x_m = y_n$ 则 $Z_{k-1}$ 是 $X_{m-1}$ 和 $Y_{n-1}$ 的最长公共子序列；
若 $x_m\neq y_n$ ，且 $x_m\neq z_k$ ，则 $Z$ 是 $X_{m-1}$ 和 $Y$ 的最长公共子序列。
若 $y_n\neq x_m$ ，且 $y_n\neq z_k$ ，则 $Z$ 是 $X$ 和 $Y_{n-1}$ 的最长公共子序列。

分析结论

当 $x_m = y_n$ 时，有一个子问题，即：找出 $X_{m - 1}$ 和 $Y_{n -1}$ 的最长公共子序列。
当 $x_m \neq y_n$ ，有两个子问题
- （1）找出 $X_{m-1}$ 和 $Y$ 的最长公共子序列，
- （2）找出 $X$ 和 $Y_{n-1}$ 的最长公共子序列。
这些子问题都包含了同一个子问题（ $X_{m-1},Y_{n-1}$ ）

动态规划的三个条件：

最优性原理：问题( $X, Y$ )包含的子问题（ $X - 1, Y - 1$ ）也是最优的，复合最优性原理。
无后效性：从前到后遍历字符串 $X$ 和 $Y$ 时，后面字符串不影响前面的最长公共序列。
有重叠子问题：问题（ $X, Y$ ）的子问题包含了同一个子问题（ $X_{m-1},Y_{n-1}$ ）

所以最长公共子序列可以用动态规划法求解，

状态转移方程为：
$\begin{cases} 0,\ i=0,j=0 \\ c[i-1][j-1]+1,x_i=y_i\\ max\{c[i][j-1],c[i-1][j]\}，x_i \neq y_j \end{cases}$

算法流程

源代码

#include <iostream>
#include <string>
#define MAX 100
#define max(x,y) ((x)>(y)?(x):(y))
using namespace std;
int m, n;
string a, b;
int dp[MAX][MAX];
void LCSlength()
{
    int i, j;
    for (i = 1; i <= m;i++)
    {
        for (j = 1; j <= n;j++)
        {
            if(a[i-1] == b[j-1])
            {
                dp[i][j] = dp[i - 1][j - 1] + 1;
            }
            else{
                dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]);
            }
        }
    }
}
int main()
{
    cin >> a >> b;
    m = a.length();
    n = b.length();
    LCSlength();
    cout << dp[m][n] << endl;
}