Advanced Fruits
The company "21st Century Fruits" has specialized in creating new sorts of fruits by transferring genes from one fruit into the genome of another one. Most times this method doesn't work, but sometimes, in very rare cases, a new fruit emerges that tastes like a mixture between both of them.
A big topic of discussion inside the company is "How should the new creations be called?" A mixture between an apple and a pear could be called an apple-pear, of course, but this doesn't sound very interesting. The boss finally decides to use the shortest string that contains both names of the original fruits as sub-strings as the new name. For instance, "applear" contains "apple" and "pear" (APPLEar and apPlEAR), and there is no shorter string that has the same property.
A combination of a cranberry and a boysenberry would therefore be called a "boysecranberry" or a "craboysenberry", for example.
Your job is to write a program that computes such a shortest name for a combination of two given fruits. Your algorithm should be efficient, otherwise it is unlikely that it will execute in the alloted time for long fruit names.
A combination of a cranberry and a boysenberry would therefore be called a "boysecranberry" or a "craboysenberry", for example.
Your job is to write a program that computes such a shortest name for a combination of two given fruits. Your algorithm should be efficient, otherwise it is unlikely that it will execute in the alloted time for long fruit names.
Input is terminated by end of file.
apple peach ananas banana pear peach
appleach bananas pearch
解题思路:
本题不难看出是LCS最长公共子序列的一个升级版,将两个序列不是最长公共子序列的部分都输出,是最长公共子序列的部分只输出一次。回顾我们在求LCS的时候的算法。
if(s1[i-1] == s2[j-1])f[i][j] = f[i-1][j-1] + 1;
else f[i][j] = max(f[i-1][j],f[i][j-1]);
我们发现当s1[i-1]和s2[j-1]相等的时候(这里序列下标从0开始),f[i][j] = f[i-1][j-1] + 1。只是在原来的长度上增加1。而不等的时候取一个最大值便可。
通过观察我们发现f[i][j]只能有三种状态转移而来f[i-1][j-1],f[i-1][j],f[i][j-1]。而f[i][j]每次更新的路径其实就是我们要求的那个合并后的序列,所以我们再用一个数组记录路径,最后用一个递归输出即可。
代码如下;
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn = 100 + 5;
char s1[maxn],s2[maxn];
int dp[maxn][maxn],p[maxn][maxn];
void print(int ls1,int ls2)
{
if(ls1 == 0 && ls2 == 0)return;
if(p[ls1][ls2] == 0)print(ls1-1,ls2-1),printf("%c",s1[ls1-1]);
else if(p[ls1][ls2] == 1)print(ls1-1,ls2),printf("%c",s1[ls1-1]);
else print(ls1,ls2-1),printf("%c",s2[ls2-1]);
}
int main()
{
while(~scanf("%s%s",s1,s2))
{
int ls1 = strlen(s1);
int ls2 = strlen(s2);
memset(dp,0,sizeof(dp));
for(int i = 1;i <= ls1; i++)p[i][0] = 1;
for(int j = 1;j <= ls2; j++)p[0][j] = -1;
for(int i = 1;i <= ls1; i++)
{
for(int j = 1;j <= ls2; j++)
{
if(s1[i-1] == s2[j-1])dp[i][j] = dp[i-1][j-1] + 1,p[i][j] = 0;
else if(dp[i-1][j] >= dp[i][j-1])dp[i][j] = dp[i-1][j],p[i][j] = 1;
else dp[i][j] = dp[i][j-1],p[i][j] = -1;
}
}
print(ls1,ls2);
printf("\n");
}
return 0;
}