LCS 最长公共子序列并输出该序列

#include <string.h>
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <string>
#include <vector>
using namespace std;
string s1,s2;
string lcs2(string s1,string s2)
{
    if(s1==""||s2=="")
        return "";
    int m=s1.size()+1;
    int n=s2.size()+1;
    printf("%d %d\n",m,n);
    int lcs[m][n];
    memset(lcs,0,sizeof(lcs));
    for(int i=1;i<m;i++)
        for(int j=1;j<n;j++)
        {
            if(s1[i-1]==s2[j-1])
                lcs[i][j]=lcs[i-1][j-1]+1;
            else
                lcs[i][j]=lcs[i-1][j]>=lcs[i][j-1]?lcs[i-1][j]:lcs[i][j-1];//取上侧或左侧的最大值
        }
        int i=m-2;
        int j=n-2;
        string ss="";
        while(i!=-1&&j!=-1)
        {
            if(s1[i]==s2[j])
            {
                //printf("%c\n",s1[i]);
                ss+=s1[i];
                i--;
                j--;
            }
            else
            {
                if(lcs[i+1][j+1]==lcs[i][j])
                {
                    i--;
                    j--;
                }
                else
                {
                    if(lcs[i][j+1]>=lcs[i+1][j])
                        i--;
                    else
                        j--;
                }
            }
        }
        reverse(ss.begin(),ss.end());//将字符串倒置
        return ss;
}
int main()
{
    while(cin>>s1>>s2)
    {
        string s=lcs2(s1,s2);
        cout << s<<endl;
    }
    return 0;
}

参考详解：http://hzzy-010.blog.163.com/blog/static/79692381200872024242126/

方法二：

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <string>
using namespace std;
char s1[1005],s2[1005],s[1005];
int a[1005][1005];
void LCS()
{
    int n=strlen(s1);
    int m=strlen(s2);
    memset(a,0,sizeof(a));
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++)
        {
            if(s1[i-1]==s2[j-1])
                a[i][j]=a[i-1][j-1]+1;
            else
                a[i][j]=a[i-1][j]>=a[i][j-1]?a[i-1][j]:a[i][j-1];
        }
    int i=n-1;
    int j=m-1;
    int k=0;
    while(i!=-1&&j!=-1)
    {
        if(s1[i]==s2[j])
        {
            s[k++]=s1[i];
            //printf("%c\n",s1[i]);
            j--;
            i--;
        }
        else
        {
            if(a[i+1][j+1]==a[i][j])
                i--,j--;
            else
            {
                if(a[i][j+1]>=a[i+1][j])
                    i--;
                else
                    j--;
            }
        }
    }
    for(int i=0;i<k/2;i++)
        swap(s[i],s[k-1-i]);
}
int main()
{
    memset(s1,0,sizeof(s1));
    memset(s2,0,sizeof(s2));
    while(~scanf("%s%s",s1,s2))
    {
        memset(s,0,sizeof(s));
        LCS();
        printf("%s\n",s);
        memset(s1,0,sizeof(s1));
        memset(s2,0,sizeof(s2));
    }
    return 0;
}