算法设计与分析 SCAU11084 多个串的最长公共子序列

11084 多个串的最长公共子序列

时间限制:1000MS 代码长度限制:10KB
提交次数:0 通过次数:0

题型: 编程题 语言: G++;GCC;VC;JAVA

Description
算法课本上有介绍“最长公共子序列”（书3.3节），那是对2个序列而言的公共子序列。
什么是子序列？子序列就是指原序列删除若干元素后形成的序列，或者说是原序列中挑选若干元素不改变它们的相对顺序形成的序列。
因此子序列中的元素并不要求在原序列中是连续的一段，完全允许跳着“挑选”。

现在要你来求n个非空串的最长公共子序列的长度（2<=n<=100）。如串13451，串5326582和串237564，
都含子序列35，因此长度为2，且这个长度已经是最大的了。

大家也不要把公共子序列简单理解为“公共的字符”，那是不对的，公共的字符，是不讲字符顺序的。
比如上面例子，公共字符为3和5，但若某个串有5在3的前面，公共字符个数依然为2，但公共子序列却形成不了“35”，
最长长度只能为1了，而不是2。

所以理解子序列的含义就是：子序列就是指原序列删除若干元素后形成的序列，或者说是原序列中挑选若干元素不改变
它们的相对顺序形成的序列。（注意：里面元素是讲顺序的）

为简化，n个串中的字符只会出现数字0_{9或小写英文字母a}z，串中不会出现空格等其他字符。每个串的串长
最多100个字符。这n个串的长度的乘积不会超过1000000。

输入格式

第一行整数n，表示有n个非空串。随后n行，每行一个串，串中无空格。

输出格式

一个数，表示这n个串的最长公共子序列的长度。

注意此题只要求出最长的公共子序列长度，而不需要记录公共子序列。


输入样例
3
13451
5326582
237564

输出样例
2

因为这三个串，都含有35子序列，长度为2，是最长的长度了。

解题思路

动态规划

以下解释中

• text1[0, i] 代表字符串 text1 的前 i + 1 个字符
• text2[0, j] 代表字符串 text2 的前 j + 1 个字符
  
  

状态定义：dp[i][j] 的值代表 text1 字符串前 i 个字符和 text2 字符串前 j 个字符的最长公共子序列的长度。

假设字符串 text1 和 text2 的长度分别为 m 和 n，创建 m+1 行 n+1 列的二维数组 dp，其中 dp[i][j] 表示 text1[0, i] 和 text2[0, j] 的最长公共子序列的长度。

例如：dp[3][4] 代表 text1 字符串前三个字符和 text2 字符串前四个字符的最长公共子序列的长度。

转移方程：

考虑动态规划的边界情况：

• 当 i = 0 时，text1[0, i] 为空，空字符串和任何字符串的最长公共子序列的长度都是 0，因此对任意 0 ≤ j ≤n，有 dp[0][j] = 0；

• 当 j=0 时，text2[0, j] 为空，同理可得，对任意 0 ≤ i ≤ m，有 dp[i][0]=0。

因此动态规划的边界情况是：

• 当 i = 0 或 j = 0 时，dp[i][j] = 0。

• 当 i > 0 且 j > 0 时，考虑 dp[i][j] 的计算：

  • 当 text1[i - 1] = text2[j - 1] 时，将这两个相同的字符称为公共字符，此时在 text1 [0, i − 1] 和 text2[0, j − 1] 的最长公共子序列基础上，再增加一个字符（即公共字符），即可得到 text1[0, i] 和 text2[0, j] 的最长公共子序列，因此 dp[i][j] = dp[i − 1][j − 1] + 1。

  • 当 text1[i − 1] != text2[j − 1] 时，考虑以下两项：

    • text1[0, i − 1] 和 text2[0, j] 的最长公共子序列；
    • text1[0, i] 和 text2[0, j − 1] 的最长公共子序列。

    要得到 text1[0, i] 和 text2[0, j] 的最长公共子序列，应取两项中的长度较大的一项，因此 dp[i][j] = max(dp[i − 1][j], dp[i][j − 1])。

由此可以得到如下状态转移方程：

最终计算得到 dp[m][n] 即为 text 1和 text 2的最长公共子序列的长度。

复杂度分析：

• 时间复杂度：O(mn)，其中 m 和 n 分别是字符串 text1 和 text2 的长度。二维数组 dp 有 m + 1 行和 n + 1 列，需要对 dp 中的每个元素进行计算。

• 空间复杂度：O(mn)，其中 m 和 n 分别是字符串 text 1和 text 2的长度。创建了 m + 1 行 n + 1 列的二维数组 dp。

输入输出处理

此题输入输出处理也比较繁琐，在上面实现了基本算法的基础上，还要好好思考一波这题的输入输出如何处理

由于此题的结果是要求这 n 个串的最长公共子序列的长度，所以需要用双 for 循环的形式来比较每两个字符串之间的最长公共序列，我这里选择用 minNum 来记录内循环每次比较时的最大公共长度，res 来记录每组外循环间的最大公共长度

例如

• 有4个字符串，外循环是 i，内循环是 j
• i = 0 时，minNum 记录的是 第1个字符串 和 第2，3，4个字符串的最大公共长度，那肯定要取最小那个（因为是多个串的公共，总不可能取你有，别人没有的那一部分吧）
• i = 1 时，minNum 记录的是 第2个字符串 和 第3，4个字符串的最大公共长度
• 则 res 记录的是 i = 0 和 i = 1 时小的那个

更多注释可查看代码中，有助于理解

代码如下

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <cstring>

using namespace std;

int maxCommonLength(char a[101], char b[101]) {
    int m = strlen(a), n = strlen(b);
    int dp[m + 1][n + 1];

    for(int i = 0; i <= m; i++) {
        for(int j = 0; j <= n; j++) {
            // 长度为0，公共肯定是0
            if(i == 0 || j == 0) {
                dp[i][j] = 0;
            } else {
                // 如果i索引的m串 和 j索引的n串字符相同，那就+1，同一行的前一个，因为是一行一行遍历的
                if(a[i - 1] == b[j - 1]) {
                    dp[i][j] = dp[i - 1][j - 1] + 1;
                } else {
                    // 为左边或上边更大的那个，可以理解成上一个状态中找更大的，由于是二维dp，所以有上一个状态有两种
                    dp[i][j] = max(dp[i][j - 1], dp[i - 1][j]);
                }
                //printf("i=%d j=%d，dp=%d\n", i, j, dp[i][j]);
            }
        }
    }
    return dp[m][n];
}

int main()
{
    int i, j, n;
    scanf("%d", &n);

    char s[n][101]; // 字符串数组，n行，每行即每个字符串最大长度是101

    for(i = 0; i < n; i++) {
        scanf("%s", s[i]);
    }

    int minNum; // 取最小那个，因为得多个串的公共长度
    int res = 101;

    // 最后一个串不用循环了
    for(i = 0; i < n - 1; i++) {
        minNum = 101;
        for(j = i + 1; j < n; j++) {
            minNum = min(minNum, maxCommonLength(s[i], s[j])); // 不断更新，记录该次 i 循环中彼此之间最小的，因为是多个串的公共，总不可能取你有，别人没有的那一部分吧
        }

        res = min(minNum, res);
        //printf("i=%d minNum=%d res=%d\n", i, minNum, res);
    }

    printf("%d", res);

    return 0;
}

最后

对我感兴趣的小伙伴可查看以下链接

• 我的掘金主页：https://juejin.cn/user/1302297507801358
• 博客主页：http://blog.zhangjiancong.top/
• 公众号：Smooth前端成长记录