动态规划之最长公共子序列LCS

有两个数组a[ ],b[ ],若他们去掉某些元素之后，剩下的元素集合S全等，则S称为a[ ]和b[ ]的公共子序列，可能有多个集合S1,S2,...SN,称其中元素个数最长的集合为a[ ]和b[ ]的最长公共子序列。

我们定义dp[ i ][ j ]表示a[1~i ]和b[1~j ]的最长公共子序列的长度，则：

（1) 如果a[ i ]==a[ j ]，那么a[ i ]和b[ j ]是dp[ i ][ j ]的最后一个元素，则dp[ i ][ j ]=dp[ i-1 ][ j-1 ]+1；

（2) 如果a[ i ]!=a[ j ]，则dp[ i ][ j ]=max（dp[ i-1 ][ j ]，dp[ i ][ j-1 ]）。

为了更好地说明问题，我们举个例子：

a[ 6 ]={3 1 4 1 5 9}

b[ 4 ]={2 6 4 9}

因为 i=6,j=4, 且a[ i ]==b[ j ]=9,所以9是S的最后一个元素，
即dp[ 4 ][ 6 ]=dp[ 4-1 ][ 6-1 ]+1=dp[ 3 ][ 5 ]+1，
接着a[ 3 ]!=b[ 5 ]，（注意i和 j都已经递减了），
但是S的倒数第二个元素可能是a[ 3 ]或者是b[ 5 ]，即可能是5或者是4，

所以dp[ 3 ][ 5 ]=max（dp[ 3-1 ][ 5 ]，dp[ 3 ][ 5-1 ] = max（dp[ 2 ][ 5 ]，dp[ 3 ][ 4 ]）。

算法书上的式子是：

最后附上代码：

#include<string.h>
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int a[10001],b[10010],dp[10005][10005]={{0,0}};
int i, j, k;
int len1, len2;
int main()
{
    while(scanf("%d%d",&len1,&len2)!=EOF)
    {
        for(i=1;i<=len1;i++)
            cin>>a[i];
        for(j=1;j<=len2;j++)
            cin>>b[j];

        memset(dp,0,sizeof(dp));

        dp[0][0]=0;
        for(i=1;i<=len1;i++)
        {
             for(j=1;j<=len2;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]);
             }
        }

        cout<<dp[len1][len2]<<endl;
    }
    return 0;
}

代码结果：