Longest Common Subsequence

Longest Common Subsequence

最長共通部分列

最長共通部分列問題 (Longest Common Subsequence problem: LCS)は、２つの与えられた列 $X = \{x_1, x_2, ..., x_m\}$ と $Y = \{y_1, y_2, ..., y_n\}$ の最長共通部分列を求める問題です。

ある列 $Z$ が $X$ と $Y$ 両方の部分列であるとき、$Z$ を $X$ と$Y$ の共通部分列と言います。例えば、$X = \{a,b,c,b,d,a,b\}$, $Y = \{b,d,c,a,b,a\}$ とすると、列 $\{b,c,a\}$ は $X$ と $Y$ の共通部分列です。一方、列 $\{b,c,a\}$ は $X$ と $Y$ の最長共通部分列ではありません。なぜなら、その長さは 3 であり、長さ 4 の共通部分列 $\{b,c,b,a\}$ が存在するからです。長さが 5 以上の共通部分列が存在しないので、列 $\{b,c,b,a\}$ は $X$ と $Y$ の最長共通部分列の１つです。

与えられた２つの文字列 $X$、$Y$に対して、最長共通部分列 $Z$ の長さを出力するプログラムを作成してください。与えられる文字列は英文字のみで構成されています。
入力

複数のデータセットが与えられます。最初の行にデータセットの数 $q$ が与えられます。続く $2 \times q$ 行にデータセットが与えられます。各データセットでは２つの文字列 $X$, $Y$ がそれぞれ１行に与えられます。

出力

各データセットについて $X$, $Y$ の最長共通部分列 $Z$ の長さを１行に出力してください。

制約

$1 \leq q \leq 150$
$1 \leq X, Yの長さ \leq 1,000$
$X$ または $Y$ の長さが 100 を超えるデータセットが含まれる場合、$q$ は 20 以下である。

入力例 1

3
abcbdab
bdcaba
abc
abc
abc
bc

出力例 1

4
3
2

参考文献

Introduction to Algorithms, Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, and Clifford Stein. The MIT Press.

解题思路：
本题是一个裸的LCS，动态转移方程如下

当$s1[i-1] == s2[j-1]$ 时 $dp[i][j] = dp[i-1][j-1] + 1$
否则$dp[i][j] = max(dp[i-1][j],dp[i][j-1])$
i和j表示序列一的前i个字符和序列二的前j个字符的最长公共子序列的长度为dp[i][j]

代码如下：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;

const int maxn = 1050;

char s1[maxn],s2[maxn];
int f[maxn][maxn];

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s%s",s1,s2);
        int len1 = strlen(s1);
        int len2 = strlen(s2);
        memset(f,0,sizeof(f));
        for(int i = 0;i <= len1; i++)
        {
            for(int j = 0;j <= len2; j++)
            {
                if(i == 0 || j == 0)f[i][j] = 0;
                else
                {
                    if(s1[i-1] == s2[j-1])f[i][j] = f[i-1][j-1] + 1;
                    else f[i][j] = max(f[i-1][j],f[i][j-1]);
                }
            }
        }
        printf("%d\n",f[len1][len2]);
    }
}