算法学习之动态规划--最长公共子序列

题目：

给出两个字符串，求出这样的一个最长的公共子序列的长度：子序列中的每个字符都能在两个原串中找到，而且每个字符的先后顺序和原串中的先后顺序一致。

Sample Input ：
abcfbc abfcab
programming contest
abcd mnp
Sample Output ：
4
2
0

题目分析：

输入两个串s1,s2, 设MaxLen(i,j)表示: s1的左边i个字符形成的子串，与s2左边的j个字符形成的子串的最长公共子序列的长度(i,j从0开始算），MaxLen(i,j) 就是本题的“状态”
假定 len1 = strlen(s1), len2 = strlen(s2），那么题目就是要求 MaxLen(len1,len2)

显然有边界值：
MaxLen(n,0) = 0 ( n= 0…len1）
MaxLen(0,n) = 0 ( n=0…len2）
递推公式：
if ( s1[i-1] == s2[j-1] ) //s1的最左边字符是s1[0]
    MaxLen(i,j) = MaxLen(i-1,j-1) + 1;
else
    MaxLen(i,j) = Max(MaxLen(i,j-1),MaxLen(i-1,j) );
时间复杂度O(mn) m,n是两个字串长度

代码如下：

#include<iostream>
#include<cstring>
using namespace std;
const int maxn = 1001;
int maxLen[maxn][maxn];
char s1[maxn],s2[maxn];
int main(){
    while(cin>>s1>>s2){
        int len1 = strlen(s1), len2 = strlen(s2);
        for(int i=0; i<=len1; i++)                                  //边界条件赋值
            maxLen[i][0] = 0;
        for(int j=0; j<=len2; j++)                                  //边界条件赋值
            maxLen[0][j] = 0;
        for(int i=1; i<=len1; i++)
            for(int j=1; j<=len2; j++){
                if(s1[i-1]==s2[j-1])                                //状态转移方程
                    maxLen[i][j]=maxLen[i-1][j-1]+1;
                else
                    maxLen[i][j]=max(maxLen[i-1][j],maxLen[i][j-1]);
            }
        cout<<maxLen[len1][len2]<<endl;
    }
    return 0;
}