Description
两个字符串, 求两个字符串最长公共子序列的长度
字符串的长度不超过5000
Sample Input Copy
ABCBDAB
BACBBD
Sample Output Copy
4
公共子串为BCBD
两个字符串a和b。对比每一个字符串。
如果a[i] == b[j] ,dp[i][j] == 左上角+1
否者 dp[i][j] == 左边或者上边最大的
#include<bits/stdc++.h>
using namespace std;
char a[5009];
char b[5009];
int dp[5009][5009];
int main()
{
scanf("%s",&a);
scanf("%s",&b);
int L = strlen(a);
int R = strlen(b);
dp[0][0] = 0;
for(int i=1;i<=R;i++)
for(int j=1;j<=L;j++){
if(a[j-1]==b[i-1])
dp[i][j] = dp[i-1][j-1]+1;
else
dp[i][j] = max(dp[i-1][j],dp[i][j-1]);//左边与上面
}
cout<<dp[R][L]<<endl;
return 0;
}
若找到公共子串则需要记录一下方向,
不过一般会有多个比如样例中:
ABCBDAB
BACBBD
包含BCBB,BCDB,ABBD。
先输入a串,或者先输入b串公共子串会不同,当dp[i-1] [j] == dp[i] [j-1]时候,记录为左方向还是上方向也得到的公共子串也不同。图片是其中一种。
为1时向左上走,为2时向上走,为3时向左走。等于1时为子串其中一个字母。
#include<bits/stdc++.h>
using namespace std;
char a[5009];
char b[5009];
int dp[5009][5009];
int vis[5009][5009];
void LCS(int i,int j)
{
if(vis[i][j] == 0)
return;
if(vis[i][j]==1){
LCS(i-1,j-1);
cout<<b[i-1]<<" "; //a[j-1]
}
else if(vis[i][j]==2)
LCS(i-1,j);
else
LCS(i,j-1);
}
int main()
{
memset(dp,0,sizeof dp);
memset(vis,0,sizeof vis);
scanf("%s",&a);
scanf("%s",&b);
int L = strlen(a);
int R = strlen(b);
dp[0][0] = 0;
for(int i=1;i<=R;i++)
for(int j=1;j<=L;j++){
if(a[j-1]==b[i-1]){
dp[i][j] = dp[i-1][j-1]+1;
vis[i][j] = 1;
}
else if(dp[i-1][j]>dp[i][j-1]){
dp[i][j] = dp[i-1][j];
vis[i][j] = 2;
}
else {
dp[i][j] = dp[i][j-1];
vis[i][j] = 3;
}
}
cout<<dp[R][L]<<endl;
LCS(R,L);
return 0;
}