求最长公共子序列
最长公共子序列问题
(Longest-Common-Subsequence problem)给定两个序列 X X X = < x 1 x_1 x1, x 2 x_2 x2, … {\ldots} …, x m x_m xm> 和 Y Y Y= < y 1 y_1 y1, y 2 y_2 y2, … {\ldots} …, y n y_n yn> ,求 X X X和 Y Y Y长度最长的公共子序列。
一个字符串 s s s 被称作另一个字符串 S S S 的子序列,说明从序列 S S S
通过去除某些元素但不破坏余下元素的相对位置(在前或在后)可得到序列 s s s 。
例如:有两个字符串 S S S 和 T T T , S S S = “abcde”, T T T = “acbfe”,那么 S S S, T T T的最长公共子序列就是 “ace”,长度为3。
刻画最长公共子序列的特征
LCS的最优子结构:令 X X X = < x 1 x_1 x1, x 2 x_2 x2, … {\ldots} …, x m x_m xm> 和 Y Y Y=
< y 1 y_1 y1, y 2 y_2 y2, … {\ldots} …, y n y_n yn> 为两个序列, Z Z Z =
< z 1 z_1 z1, z 2 z_2 z2, … {\ldots} …, z k z_k zk> 为 X X X和 Y Y Y的任意 L C S LCS LCS。
- 如果 x m x_m xm = y n y_n yn ,则 z k z_k zk = x m x_m xm = y n y_n yn,且 Z k − 1 Z_{k-1} Zk−1是 X m − 1 X_{m-1} Xm−1和 Y n − 1 Y_{n-1} Yn−1的一个 L C S LCS LCS。
- 如果 x m x_m xm ≠ \neq = y n y_n yn,那么 z k z_k zk ≠ \neq = x m x_m xm意味着 Z Z Z是 X m − 1 X_{m-1} Xm−1 和 Y Y Y 的一个 L C S LCS LCS。
- 如果 x m x_m xm ≠ \neq = y n y_n yn,那么 z k z_k zk ≠ \neq = y n y_n yn意味着 Z Z Z是 X X X 和 Y n − 1 Y_{n-1} Yn−1 的一个 L C S LCS LCS。
递归公式
根据
L
C
S
LCS
LCS 问题的最优子结构性质,可以得到如下公式:
其中dp[ i ] [ j ]表示
S
1
i
S1_i
S1i 和
S
2
j
S2_j
S2j 的
L
C
S
LCS
LCS长度。
C++代码
int dp[1010][1010];
int main()
{
string s1,s2;
cin>>s1>>s2
int len1=s1.length();
int len2=s2.length();
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]);
}
cout<<dp[len1][len2]<<endl;
}
int dp[3010][3010];
int b[3010][3010];//记录路径
void print_LCS(string s,int i,int j)
{
if(i==0 || j==0) return;
if(b[i][j]==0)
{
print_LCS(s,i-1,j-1);
cout<<s[i-1];
}
else if(b[i][j]==-1)
{
print_LCS(s,i-1,j);
}
else print_LCS(s,i,j-1);
}
int main()
{
string s,t;
cin>>s>>t;
int len1=s.length();
int len2=t.length();
memset(dp,0,sizeof(dp));
for(int i=1;i<=len1;i++)
for(int j=1;j<=len2;j++)
{
if(s[i-1]==t[j-1])
{
dp[i][j]=dp[i-1][j-1]+1;
b[i][j]=0;//0表示指向左上方
}
else if(dp[i-1][j]>=dp[i][j-1])
{
dp[i][j]=dp[i-1][j];
b[i][j]=-1;//-1表示指向上方
}
else
{
dp[i][j]=dp[i][j-1];
b[i][j]=1;//1表示指向左方
}
}
print_LCS(s,len1,len2);
return 0;
}
模板题
求最长公共子序列的长度
打印最长公共子序列