第二十六次总结

今天主要是学了最长公共子串问题

最长公共子序列（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;
}