LCS O（n+m）模板

LCS算法详解参见：动态规划解最长公共子序列问题

最长公共子序列求解：递归与动态规划方法

//LCS模板
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define eps 1e-9
#define INF 0x3f3f3f3f
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,m
#define rson m+1,r
const int MAXN=1e5+10;
const int maxn=510;
char a[maxn],b[maxn];
int num[maxn][maxn],flag[maxn][maxn];//num记录a[i]和b[j]的LSC的长度   flag记录num[i][j]是通过哪一个子问题的值求得的，并记录搜索方向
inline void LCS()//这是一个回溯的思想
{
    for(int i=1;i<=strlen(a);i++)//数组中从0开始存储，遍历从1开始
        for(int j=1;j<=strlen(b);j++)
        {
            if(a[i-1]==b[j-1])//注意下标是i-1和j-1
            {
                num[i][j]=num[i-1][j-1]+1;
                flag[i][j]=1;//斜向下标记
            }
            else if(num[i][j-1]>num[i-1][j])//这两个条件相当于num[i][j]=max(num[i-1][j],num[i][j-1])
            {
                num[i][j]=num[i][j-1];
                flag[i][j]=2;//向右标记
            }
            else
            {
                num[i][j]=num[i-1][j];
                flag[i][j]=3;//向下标记
            }
        }
}
inline void getLCS()
{
    char res[500];//存储最终结果，即两个数组中的公共部分
    int i=strlen(a),j=strlen(b);
    int k=0;//保存res数组的标志位
    while(i>0&&j>0)
    {
        if(flag[i][j]==1)//斜向下标记
        {
            res[k]=a[i-1];
            k++;
            i--;
            j--;
        }
        else if(flag[i][j]==2)//向右标记
                j--;
              else if(flag[i][j]==3)//向下标记
                    i--;
    }
    for(int i=k-1;i>=0;i--)//结果数组输出
        printf("%c",res[i]);
}
int main()
{
    std::ios::sync_with_stdio(false);
    while(cin>>a>>b)
    {
        LCS();
        getLCS();
        mem(a,0);
        mem(b,0);
    }
    return 0;
}