最长公共子序列（LCS）

描述：
给出两个序列X,Y，求解出这两个序列的最长公共子序列。（最长公共子序列是不要求连续的）；
例如 ： X=[1 ，，，m] ;  Y=[1，，，n];
利用穷举法，找出X的所有子序列，然后每个子序列对序列Y扫描一边，检查是否是序列Y的子序列；
X 的子序列的个数为2^m，与A进行比较所需时间为O(n) 。
则利用穷举法的时间复杂度为O[n*2^m]；有幂指数在，所以穷举法的时间复杂度太高了；


既然是要找到最长的公共子序列，那么我们就从长度入手
1）找出两个序列的最长子序列的长度；
2）返回寻找满足的子序列，并打印路径；


那么，
设A=“a0，a1，…，am-1”，B=“b0，b1，…，bm-1”，并Z=“z0，z1，…，zk-1”为它们的最长公共子序列。不难证明有以下性质：
（1） 如果am-1=bn-1，则zk-1=am-1=bn-1，且“z0，z1，…，zk-2”是“a0，a1，…，am-2”和“b0，b1，…，bn-2”的一个最长公共子序列；
（2） 如果am-1!=bn-1，则若zk-1!=am-1，蕴涵“z0，z1，…，zk-1”是“a0，a1，…，am-2”和“b0，b1，…，bn-1”的一个最长公共子序列；
（3） 如果am-1!=bn-1，则若zk-1!=bn-1，蕴涵“z0，z1，…，zk-1”是“a0，a1，…，am-1”和“b0，b1，…，bn-2”的一个最长公共子序列。
这样，在找A和B的公共子序列时，如有am-1=bn-1，则进一步解决一个子问题，找“a0，a1，…，am-2”和“b0，b1，…，bm-2”的一个最长公共子序列；如果am-1!=bn-1，则要解决两个子问题，找出“a0，a1，…，am-2”和“b0，b1，…，bn-1”的一个最长公共子序列和找出“a0，a1，…，am-1”和“b0，b1，…，bn-2”的一个最长公共子序列，再取两者中较长者作为A和B的最长公共子序列。
根据如上的三条性质，可以编写动态规划程序；
设一个数组c[ ][ ]记录序列X，Y的LCS的长度，b[ ][ ]数组用来记录c[ ][ ]是由哪一个子问题得来的；

问题的递归式写成：

回溯输出最长公共子序列过程：

由于每次调用至少向上或向左（或向上向左同时）移动一步，故最多调用（m+n）次就会遇到i==0 或 j==0的情况，此时开始返回。返回时与递归调用时方向相反，步数相同，故算法时间复杂度为O(m+n)。

代码：

/****LCS****/

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=1005;
int dp[N][N],b[N][N];
char s1[N],s2[N];
int len1,len2;

int LCS()
{
    memset(dp,0,sizeof(dp));
    for(int i=1; i<=len1; i++)
    {
        for(int j=1; j<=len2; j++)
        {
            if(s1[i-1]==s2[j-1])
            {
                dp[i][j]=dp[i-1][j-1]+1;
                b[i][j]=0;
            }
            else if(dp[i-1][j]>=dp[i][j-1])
            {
                dp[i][j]=dp[i-1][j];
                b[i][j]=1;
            }
            else
            {
                dp[i][j]=dp[i][j-1];
                b[i][j]=2;
            }
        }
    }
    return dp[len1][len2];
}

void print(int m,int n)
{
    if(m==0||n==0) return ;
    if(b[m][n]==0)
    {
        print(m-1,n-1);
        printf("%c ",s1[m-1]);
    }
    else if(b[m][n]==1)
        print(m-1,n);
    else
        print(m,n-1);
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        for(int i=0; i<n; i++)
        {
            scanf("%s%s",s1,s2);
            len1=strlen(s1),len2=strlen(s2);
            printf("%d\n",LCS());
        print(len1,len2);
            printf("\n");
        }
    }
    return 0;
}

注：参考转载这里；