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

最长公共子序列问题
Time Limit: 1000 ms Memory Limit: 65536 KiB
Submit Statistic Discuss
Problem Description

给定两个序列 X={x1,x2,…,xm} 和 Y={y1,y2,…,yn}，找出X和Y的最长公共子序列。
Input

输入数据有多组，每组有两行 ，每行为一个长度不超过500的字符串（输入全是大写英文字母（A,Z）），表示序列X和Y。
Output

每组输出一行，表示所求得的最长公共子序列的长度，若不存在公共子序列，则输出0。
Sample Input

ABCBDAB
BDCABA
Sample Output

4

题解：
maxLen(i,j) 第一个长度为i子序列和第二个长度为j的子序列形成的最长公共子序列长度。

len1=s1.length();
len2=s2.length();

则本题要求为 maxLen(len1,len2);
考虑边界问题：

maxLen(n,0)=0;//(n=0、……、len1）
maxLen(0,n)=0;//(n=0、……、len2）

递推公式：

if(s1[i-1]==s2[j-1]) 
    maxLen(i,j)=maxLen(i-1,j-1)+1;
if(s1[i-1]!=s2[j-1])
    maxLen(i,j)=max{maxlen(i-1,j),maxlen(i,j-1)};

时间复杂度是O（mn）(m、n是两个子序列长度）。

代码：

#include<iostream>

#include<cstring>

#include<cstdio>

#include<algorithm>

#include<cmath>

using namespace std;

int maxLen[1005][1005];

int max(int a,int b)

{

    return a>b?a:b;

}

int main()

{

    string s1,s2;

    while(cin>>s1>>s2)

    {

        int len1=s1.length();

        int len2=s2.length();

        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;

    }

}