题目
题目描述
给出两个字符串,求两个字符串的最长公共字串。
输入
abcfbc abfcab
programming contest
abcd mnp
输出
4
2
0
题解
分析
可以看出:
F[i][j]=F[i-1][j-1]+1;(a[i]==b[j])
F[i][j]=max(F[i-1][j],F[i][j-1])(a[i]!=b[j]);
n由于F(i,j)只和F(i-1,j-1), F(i-1,j)和F(i,j-1)有关, 而在计算F(i,j)时, 只要选择一个合适的顺序, 就可以保证这三项都已经计算出来了, 这样就可以计算出F(i,j). 这样一直推到f(len(a),len(b))就得到所要求的解了.
代码
#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;
int dp[1100][1100];
char s1[1000],s2[1000];
int main()
{
while(~scanf("%s%s",s1,s2)){
int len1=strlen(s1);
int len2=strlen(s2);
memset(dp,0,sizeof(dp));
for(int i=1;i<=len1;++i){
for(int j=1;j<=len2;++j){
if(s1[i-1]==s2[j-1])
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
}
}
printf("%d\n",dp[len1][len2]);
}
return 0;
}
拓展 - 打印序列 - 求序列路径
给出两个字符串,求两个字符串的最长公共字串及打印公共子串
输入
abcfbc abfcab
programming contest
abcd mnp
输出
LCS=4
a b c b
LCS=2
o n
LCS=0
代码
#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;
int dp[110][110],path[110][110];
char s1[110],s2[110];
void LCS(int len1,int len2)
{
for(int i=1;i<=len1;i++){
for(int j=1;j<=len2;j++){
if(s1[i-1]==s2[j-1]){
dp[i][j]=dp[i-1][j-1]+1;
path[i][j]=1;
}else if(dp[i][j-1]<=dp[i-1][j]){
dp[i][j]=dp[i-1][j];
path[i][j]=2;
}else{
dp[i][j]=dp[i][j-1];
path[i][j]=3;
}
}
}
printf("LCS=%d\n",dp[len1][len2]);
}
void print(int x,int y)
{
if(x==0||y==0) return;//注意边界
if(path[x][y]==1){
print(x-1,y-1);
printf("%c ",s1[x-1]);
}else if(path[x][y]==2){
print(x-1,y);
}else if(path[x][y]==3){
print(x,y-1);
}
}
int main()
{
string s;
while(~scanf("%s%s",s1,s2)){
memset(dp,0,sizeof(dp));
int len1=strlen(s1),len2=strlen(s2);
LCS(len1,len2);
print(len1,len2);
cout<<endl;
}
return 0;
}
拓展 - 打印序列 - 求序列路径
给定两个字符串,输出两个字符串公共不重复部分。
输入
apple peach
ananas banana
pear peach
输出
appleach
bananas
pearch
代码
将两个字符串结合起来,他们的公共子串只输出一次。
根据LCS的原理,将每个字符都进行标记,看两个字符串中对应的字符究竟处于什么状态,然后输出,其标记为公共子串的字符只输出一次即可。
注意数组的初始化。
#include<stdio.h>
#include<string.h>
#include<math.h>
#include <iostream>
#include <algorithm>
using namespace std;
int dp[110][110],path[110][110];
char s1[110],s2[110];
void LCS(int len1,int len2)
{
for(int i=0;i<=len1;i++)
path[i][0]=2;
for(int i=0;i<=len2;i++)
path[0][i]=3;
for(int i=1;i<=len1;i++){
for(int j=1;j<=len2;j++){
if(s1[i-1]==s2[j-1]){
dp[i][j]=dp[i-1][j-1]+1;
path[i][j]=1;
}else if(dp[i][j-1]<=dp[i-1][j]){
dp[i][j]=dp[i-1][j];
path[i][j]=2;
}else{
dp[i][j]=dp[i][j-1];
path[i][j]=3;
}
}
}
/*for(int i=1;i<=len1;i++){
for(int j=1;j<=len2;j++){
printf("%d ",path[i][j]);
}
printf("\n");
}
printf("LCS=%d\n",dp[len1][len2]);
*/
}
void print(int x,int y)
{
if(x==0&&y==0) return;//注意边界
if(path[x][y]==1){
print(x-1,y-1);
printf("%c",s1[x-1]);
}else if(path[x][y]==2){
print(x-1,y);
printf("%c",s1[x-1]);
}else if(path[x][y]==3){
print(x,y-1);
printf("%c",s2[y-1]);
}
}
int main()
{
string s;
while(~scanf("%s%s",s1,s2)){
memset(dp,0,sizeof(dp));
int len1=strlen(s1),len2=strlen(s2);
LCS(len1,len2);
print(len1,len2);
cout<<endl;
}
return 0;
}