最长公共子序列

最长公共子序列

动态规划的题目都是很灵活的，如果我们要解决问题，状态的设计是非常关键的

描述
我们称序列Z = < z1, z2, …, zk >是序列X = < x1, x2, …, xm >的子序列当且仅当存在 严格上升 的序列< i1, i2, …, ik >，使得对j = 1, 2, … ,k, 有xij = zj。比如Z = < a, b, f, c > 是X = < a, b, c, f, b, c >的子序列。

现在给出两个序列X和Y，你的任务是找到X和Y的最大公共子序列，也就是说要找到一个最长的序列Z，使得Z既是X的子序列也是Y的子序列。

输入
输入包括多组测试数据。每组数据包括一行，给出两个长度不超过200的字符串，表示两个序列。两个字符串之间由若干个空格隔开。
输出
对每组输入数据，输出一行，给出两个序列的最大公共子序列的长度。
样例输入

abcfbc         abfcab
programming    contest 
abcd           mnp

样例输出

4
2
0

如果小伙伴们想挑战一下百练POJ

我这里提供一种状态的设计
maxLen(i,j)表示s1字符前i个元素的最长公共序列，j就是s2的了，显然满足无后效性，最优子问题

//公共子序列
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <string>
#define len 1000
using namespace std;
char s1[len],s2[len];
int maxLen[len][len];//按照分析的思路来分解子问题

int main()
{
    while(cin>>s1>>s2)
    {
        int len1=strlen(s1);
        int 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][j-1],maxLen[i-1][j]);
                }
            }
        }
         cout<<maxLen[len1][len2]<<endl;
    }
    return 0;
}

maxLen[i][j]=max(maxLen[i][j-1],maxLen[i-1][j]);
  //这个情况可能比较难理解，大家可以证明max[i][j]绝对不会比他们两个都大，也不比任意一个小，画图看看吧

学会程序和算法，走遍天下都不怕

海南五指山