首页 > 其他 > 详细

poj 1934 Trip (LCS)

时间:2015-02-02 23:18:13      阅读:345      评论:0      收藏:0      [点我收藏+]
/*
按升序输出两个串所有的最长公共子序列
首先求出最长公共子序列,然后处理处两个数组,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;
}

poj 1934 Trip (LCS)

原文:http://blog.csdn.net/lp_opai/article/details/43412173

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!