当前位置: 代码迷 >> 综合 >> poj 1934 (LCS+输出路径 好题)
  详细解决方案

poj 1934 (LCS+输出路径 好题)

热度:55   发布时间:2024-01-30 04:16:32.0

题意:
输出LCS的所有不重路径

思路:
LCS很好求,但是怎么弄出路径呢,容易想到枚举最长下的每个位置的字符,这个可以不用去重。那么怎么枚举呢?设两个pos数组,pos[i][j]表示1到i的子串中字符j的最后一个的出现位置,为什么这么做呢?因为我们只需要这个字符存在即可,哪个位置我们不关心,所以我们尽量让他靠后,这样就可以使得同样状态下前面的枚举具有更优的可能。
这样乍看起来是O(26^80)次方,但题目保证不同的串的个数不超过一千个,而我们每次枚举的时候,只会往合法满足LCS的位置递归,所以其实是O(1000*80)而已,另外我们可以用set存储答案,保证有序

#include<iostream>
#include<cmath>
#include<algorithm>
#include<stdio.h>
#include<cstring>
#include<vector>
#include<map>
#include<set>using namespace std;char a[100], b[100], tmp[100];
int dp[100][100], pos1[100][100], pos2[100][100];
set<string> ans;
int lena, lenb;void solve(int len1, int len2, int len)
{if(len <= 0){string str=tmp + 1;ans.insert(str);return;}if(len1 > 0 && len2 > 0){for(int i = 0; i < 26; i++){int p1 = pos1[len1][i];int p2 = pos2[len2][i];if(dp[p1][p2] == len){tmp[len] = i + 'a';solve(p1 - 1, p2 - 1, len - 1);}}}
}int main()
{while(~scanf("%s%s", a + 1, b + 1)){lena = strlen(a + 1);lenb = strlen(b + 1);for(int i = 0; i <= lena; i++)dp[i][0] = 0;for(int j = 0; j <= lenb; j++)dp[0][j] = 0;for(int i = 1; i <= lena; i++){for(int j = 1; j <= lenb; j++){dp[i][j] = max(dp[i][j], max(dp[i-1][j],dp[i][j-1]));if(a[i] == b[j])dp[i][j] = dp[i - 1][j - 1] + 1;}}memset(pos1, -1, sizeof(pos1));memset(pos2, -1, sizeof(pos2));for(int i = 1; i <= lena; i++)for(int j = 0; j < 26; j++){if(a[i] == j +'a')pos1[i][j] = i;elsepos1[i][j] = pos1[i - 1][j];}for(int i = 1; i <= lenb; i++)for(int j = 0; j < 26; j++){if(b[i] == j + 'a')pos2[i][j] = i;elsepos2[i][j] = pos2[i - 1][j];}memset(tmp, 0, sizeof(tmp));solve(lena, lenb, dp[lena][lenb]);set<string>::iterator iter;for(iter = ans.begin(); iter != ans.end(); iter++){cout<<*iter<<"\n";}}return 0;
}