LCS记录路径(26个字符以内，重复算一个)

声明：此文转载，方便自己以后分类查看

1）首先按照常规的方法求出最长公共子序列的长度
也就是用O(MN)的那个动态规划，结果放在二维数组dp里
dp[i][j] = { 字串a的1~i部分与字串b的1~j部分的最长公共子序列的长度 }
2）求辅助数组
last1[i][j] = { 到下标i为止，字符j在字串a中最后一次出现的下标 }
last2[i][j] = { 到下标i为止，字符j在字串b中最后一次出现的下标 }
3）枚举最长公共字串的每一个字符
从最后一个字符开始枚举
比如说现在枚举最后一个字符是'C'的情况。
那么 'CDCD' 与 'FUCKC' 这两个字串。
一共有 (0, 2) (0, 4)  (2, 2)  (2. 4) 这四种可能。
很明显前三个是可以舍弃的，因为第四个优于前三个，为后续的枚举提供了更大的空间。
last数组正好是用来做这个的。
4）排序输出
代码里用了stl的set。
注意，由于刚刚的枚举过程是针对每个字符，所以是不用判重的。

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <map>
#include <set>
#include <string>
#define frein freopen("D:\\cprogram\\acmprogram\\input", "r", stdin)
#define freout freopen("D:\\cprogram\\acmprogram\\ouput", "w", stdout)
#define mem(a) memset(a, 0, sizeof(a))
#define MAXN 100
#define MAXM 100
#define INF 1000000000
#define MOD 1000000000
using namespace std;
int dp[MAXM][MAXM];
int last1[MAXM][30],last2[MAXM][30];
int longest,len1,len2;
char s1[MAXM],s2[MAXM];
char temp[MAXM];
set<string>SET;

void find(int end1, int end2, int lens)
{
//    printf("%d %d %d\n",end1,end2,lens);
    if(lens==0)
    {
        SET.insert(temp+1);
        return;
    }

    for (int i = 0; i < 26; ++i)
    {
        int t1 = last1[end1][i];
        int t2 = last2[end2][i];
        if(dp[t1][t2]==lens)
        {
        //这个地方会把重复的去重,比如abcc abc只算了abcc最后个c，
        //如果都要输出就外层接2个循环,把所有结尾的都算上;
            temp[lens]=i+'a';
            find(t1-1,t2-1,lens-1);
        }
    }
}

void solve()
{
    len1 = strlen(s1+1);
    len2 = strlen(s2+1);
    for (int i = 1; i <= len1; ++i)
        for (int j = 1; j <= len2; ++j)
        {
            dp[i][j] = max(dp[i][j-1],dp[i][j]);
            dp[i][j] = max(dp[i][j],dp[i-1][j]);
            if(s1[i]==s2[j])
                dp[i][j] = dp[i-1][j-1]+1;
        }
    longest = dp[len1][len2];
    //    printf("%d\n", longest)；LCS结果

    //这里要保证last数组里边都是0;
        
    for (int i = 1; i <= len1; ++i)
        for (int j = 0; j < 26; ++j)
        {
            if(s1[i]==j+'a')
                last1[i][j] = i;
            else
            last1[i][j] = last1[i-1][j];
        }

    for (int i = 1; i <= len2; ++i)
        for (int j = 0; j < 26; ++j)
        {
            if(s2[i]==j+'a')
                last2[i][j] = i;
            else
            last2[i][j] = last2[i-1][j];
        }
//        printf("%d %d\n",last1[4][2],last1[3][2]);
        temp[longest+1]='\0';
        find(len1,len2,longest);
        set<string>::iterator it;
        for(it=SET.begin();it!=SET.end();it++)
            printf("%s\n",(*it).c_str());
}
int main(int argc, char const *argv[])
{
    scanf("%s%s",s1+1,s2+1);
    solve();
    return 0;
}