最长公共子序列

最长公共子序列就是给定两个字符串，求这两个字符串中按严格不下降的下标求得他们最长的且一样的子序列。
eg： 字符串A=s a d w e，字符串B=d a w e t t。那A和B的最长公共子序列D=d w e。
那如何求最长公共子序列呢？
显然使用动态规划，我们第一个二维数组c[i][j]表示从字符串A开始到第i个字符和从字符串B开始到底j个字符所能达到的最长公共子序列的长度。
如果A[i]==B[j],那显然d[i][j]=d[i-1][j-1]+1。如果A[i]!=B[j]那我们就看究竟是c[i-1][j]大还是c[i][j-1]大，然后将最大值赋给c[i][j]。那如何求最长公共子序列的内容呢?
我们可以定义一个数组d[],lena,lenb和siz,d[]存储最长公共子序列的内容，lena是字符串A的长度，lenb是字符串B的长度,siz是e[]数组的长度。可以用for循环，当A[lena]!=B[lenb]时，可以比较c[lena-1][lenb]和c[lena][lenb-1]的大小，由此判断是将A[lena–]最优还是B[lenb–]最优。然后当A[lena]==B[lenb]时将A[lena]或B[lenb]存入d[siz++]。最后逆序输出就好了。
贴一下代码：
ps（代码中non就是siz）

#include<iostream>
#include<cstring>
using namespace std;
int c[1000][1000];
int lena,lenb;
char a[1000],b[1000];
char d[1000];
int non;
int main ()
{
    memset(c,0,sizeof(c));
    cin>>lena>>lenb;
    for(int i=1;i<=lena;i++)
    cin>>a[i];
    for(int i=1;i<=lenb;i++)
    cin>>b[i];
    for(int i=2;i<=lena+1;i++)
    {
        for(int j=2;j<=lenb+1;j++)
        {
            if(a[i-1]==b[j-1])
            {
                c[i][j]=c[i-1][j-1]+1;
            }
            else
            {
                c[i][j]=max(c[i-1][j],c[i][j-1]);
            }
        }
    }
    cout<<c[lena+1][lenb+1]<<endl;
    for(;;)
    {
        if(lena==0||lenb==0)
        break;
        else
        {
            if(a[lena]==b[lenb])
            {
                non++;
                d[non]=a[lena];
                lena--;
                lenb--;
            }
            else
            {
                if(c[lena-1][lenb]>c[lena][lenb-1])
                {
                    lena--;
                }
                else
                lenb--;
            }
        }
    }
    for(int i=non;i>=1;i--)
    cout<<d[i]<<" ";
    return 0;
}