[动态规划] LCS最长公共子序列问题

LCS问题简述：最长公共子序列

一个数列 S，如果分别是两个或多个已知数列的子序列，且是所有符合此条件序列中最长的，则S 称为已知序列的最长公共子序列。

子串（Substring）是串的一个连续的部分，子序列（Subsequence）则是从不改变序列的顺序，而从序列中去掉任意的元素而获得的新序列；更简略地说，前者（子串）的字符的位置必须连续，后者（子序列LCS）则不必。比如字符串acdfg同akdfc的最长公共子串为df，而他们的最长公共子序列是adf。

---

解决本题的核心思想是动态规划。我们先把问题分成几个小的子问题

dp[i][j]表示a[0]~a[i]与b[0]~b[j]这两个串的LCS

边界条件：当i=0或j=0，dp[i][j]=0;
当str[i]=str[j]时 说明这一个字符可以被加入LCS 所以dp[i][j]=dp[i-1][j-1]
当str[i]!=str[j]时 说明这个字符不能加入LCS，要继承上个状态的最优值（最大值）所以dp[i][j]=max(dp[i-1][j],dp[i][j-1])

---

附上代码

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define maxn 5000
using namespace std;
int a[maxn],b[maxn],c[maxn][maxn],mmax=0;
int n,m;
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        c[i][0]=0;
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&b[i]);
        c[0][i]=0;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(a[i]==b[j])
            {
                c[i][j]=c[i-1][j-1]+1;
            }
            else
            {
                c[i][j]=max(c[i-1][j],c[i][j-1]);
            }
            mmax=max(mmax,c[i][j]);
        }
    }
    printf("%d",mmax);
}