最长公共子序列 LCS 【长度及打印】

给出两个字符串A B，求A与B的最长公共子序列（子序列不要求是连续的）。
比如两个串为：

abcicba
abdkscab

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。
Input
第1行：字符串A 第2行：字符串B (A,B的长度 <= 1000)
Output
输出最长的子序列，如果有多个，随意输出1个。
Sample Input
abcicba
abdkscab
Sample Output
abca

DP 模型
代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<set>
#define CLR(a,b) memset((a),(b),sizeof(a))
#define inf 0x3f3f3f3f
#define mod 100009
#define LL long long
#define M  1000+10
#define ll o<<1
#define rr o<<1|1
#define lson o<<1,l,mid
#define rson o<<1|1,mid+1,r
using namespace std;
char a[M],b[M];
char an[M];
int dp[M][M];
int main()
{   
    gets(a+1);gets(b+1);
    int lena=strlen(a+1);
    int lenb=strlen(b+1);
    int i,j;
    for(i=1;i<=lena;i++)
    {
        for(j=1;j<=lenb;j++)
        {
            if(a[i]==b[j]) dp[i][j]=dp[i-1][j-1]+1;
            else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
        }
    }
    i--;j--;int c=dp[lena][lenb];
    for(c;c>0;  )
    {
        if(a[i]==b[j])
        {
            an[c]=a[i];
            i--;j--;c--;
         } 
        else {
            if(dp[i-1][j]==dp[i][j]) i--;
            else j--;
        }
    }
    puts(an+1);
    return 0;
}