poj 1934 Trip (LCS)

/*
按升序输出两个串所有的最长公共子序列
首先求出最长公共子序列，然后处理处两个数组，f1['a'~'z'][j]表示该字母在第一个串的前j个字母中出现的最大下标，
f2['a'~'z'][j]表示该字母在第二个串的前j个字母中出现的最大下标。现在我们已经知道了最长公共子序列的长度，
我们从最后一位开始枚举每一位放什么字母，用dfs来实现，并用前面处理出来的数组进行可行性剪枝。
由于最后答案的串不超过1000个，所以暴搜的复杂度可以接受。
*/
# include <stdio.h>
# include <algorithm>
# include <iostream>
# include <string.h>
# include <vector>
# include <string>
using namespace std;
int dp[110][110];
int f1[30][110];
int f2[30][110];
vector<string>ans;
void dfs(int x,int y,string s,int L)
{
    if(x<0||y<0)
        return ;
    if(L<=0)//搜到了最底层，得到一个解
    {
        ans.push_back(s);
        return ;
    }
    char ch;
    for(int i=0; i<26; i++)
    {
        int x1=f1[i][x];
        int x2=f2[i][y];
        if(dp[x1][x2]==L)//第L位放（i+'a')字母可行
        {
            ch='a'+i;
            dfs(x1-1,x2-1,ch+s,L-1);
        }
    }
}
int main()
{
    string  a,b;
    int i,j,lena,lenb;
    while(cin>>a>>b)
    {

        ans.clear();
        lena=a.size();
        lenb=b.size();
        memset(dp,0,sizeof(dp));
        for(i=1; i<=lena; i++)
        {
            for(j=1; j<=lenb; j++)
            {
                if(a[i-1]==b[j-1])
                    dp[i][j]=dp[i-1][j-1]+1;
                else
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            }
        }
        for(i=0; i<26; i++)
        {
            f1[i][0]=0;
            for(j=1; j<=lena; j++)
            {
                if(a[j-1]==i+'a')
                    f1[i][j]=j;
                else
                    f1[i][j]=f1[i][j-1];
            }
        }
        for(i=0; i<26; i++)
        {
            f2[i][0]=0;
            for(j=1; j<=lenb; j++)
            {
                if(b[j-1]==i+'a')
                    f2[i][j]=j;
                else
                    f2[i][j]=f2[i][j-1];
            }
        }
     //   printf("%d\n",dp[lena][lenb]);
        dfs(lena,lenb,"",dp[lena][lenb]);
        sort(ans.begin(),ans.end());
        for(i=0; i<(int)ans.size(); i++)
            cout<<ans[i]<<endl;

    }
    return 0;
}