/*
按升序输出两个串所有的最长公共子序列
首先求出最长公共子序列,然后处理处两个数组,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)
最新推荐文章于 2022-02-17 15:08:40 发布