今天主要是学了最长公共子串问题
最长公共子序列(LCS)和最长公共子串不是一个概念,需要区分
子序列:可以不连续
子串:必须连续
例:
EABABCFDEC 与 EABCE 来说
EABABCFDEC 与EABCE 这个为公共子串
EABABCFDEC 与EABCE 这个为公共子序列,EABABCFDEC 与EABCE 这个也可以是公共子序列,按顺序的交集就可以了
(1)如果s1与s2最后一个元素都相等
那么s1和s2的LCS就等于s1减去最后一个元素与s2减去一个元素的LCS再加上s1和s2相等的最后一个元素
(2)如果s1与s2最后一个元素不相等
那么s1与s2的LCS就等于 S1减去最后一个元素与s2的LCS,s2减去最后一个元素去s1的LCS 中取一个最大的
dp公式
s1[0->i]与s2[0->j]的最长公共子序列为:dp[i][j]
dp公式
s1[0->i]与s2[0->j]的最长公共子序列为:dp[i][j]
推导dp公式
(1)当s1与s2最后一个元素相同时
此时LCS为s1减去最后一个元素与s2减去一个元素的LCS再加上s1和s2相等的最后一个元素
我们假设减去最后一个元素LCS为5,那么最终的LCS=5+1=6
dp[i][j]=dp[i-1][j-1]+1
(2)当s1与s2最后一个元素不相同时
那么s1与s2的LCS就等于 S1减去最后一个元素与s2的LCS,s2减去最后一个元素去s1的LCS 中取一个最大的
dp[i][j]=max(dp[i-1][j],dp[i][j-1])
P1439 【模板】最长公共子序列
#include<iostream>
using namespace std;
const int N=1005;
int a[N],b[N],dp[N][N];
signed main()
{
int n;
cin>>n;
for(int i=1;i<=n;i++)
cin>>a[i];
for(int i=1;i<=n;i++)
cin>>b[i];
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(a[i]==b[j])
dp[i][j]=dp[i-1][j-1]+1;
else
dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
}
}
cout<<dp[n][n]<<endl;
}
然而数据太大,只过了一半🥲
之后学了优化再写一遍
P2758 编辑距离
对于删除还是插入,取决于dp数组相对应的变化的值,那个小我们就取哪一个
然后就是对于是否替换的判断,取决于我们s1与s2最后一个数是否相等,相等的话就是前面的哪一个dp表dp[i-1][j-1],不相等的话,就需要在前面的基础上+1,也就是dp[i-1][j-1]+1;
注意:我们要先对dp表进行初始化,因为变化了,全都需要再变化
在i=0时全部需要添加,所以dp[0][j]=j
j=0时全部都需要删除,所以dp[i][0]=i
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=2005;
char arr[maxn],brr[maxn];
int dp[maxn][maxn],len1,len2;
int min(int x,int y){
if(x<y)
return x;
else
return y;
}
int main(){
scanf("%s",arr+1);
len1=strlen(arr+1);
scanf("%s",brr+1);
len2=strlen(brr+1);
memset(dp,-1,sizeof(dp));
for(int i=0;i<=len1;i++)
dp[i][0]=i;
for(int i=0;i<=len2;i++)
dp[0][i]=i;
for(int i=1;i<=len1;i++){
for(int j=1;j<=len2;j++){
dp[i][j]=min(dp[i-1][j]+1,dp[i][j-1]+1);
if(arr[i]==brr[j])
dp[i][j]=min(dp[i][j],dp[i-1][j-1]);
else
dp[i][j]=min(dp[i][j],dp[i-1][j-1]+1);
}
}
cout<<dp[len1][len2]<<endl;
}