暑假acwing算法总结31：线性DP

• 一般的DP问题都没有特定的解法，当然除了背包问题，所以
  只能找典例，下面介绍一些经典的线性DP问题

• 1、数字三角形

• 代码如下

#include<iostream>
using namespace std;
const int INF=1e9;
int a[510][510];
int f[510][510];
int main()
{
    int n;cin>>n;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=i;j++)
            cin>>a[i][j];
    }
    for(int i=1;i<=n;i++)
    f[n][i]=a[n][i];
    for(int i=n-1;i>=1;i--)
    {
        for(int j=1;j<=i;j++)
        {
            f[i][j]=max(f[i+1][j]+a[i][j],f[i+1][j+1]+a[i][j]);
        }
    }
    cout<<f[1][1]<<endl;
    return 0;
}

• 2、最长上升子序列

• f[i]记录以当前为结尾的前面的最长上升序列的长度
• 代码

#include<iostream>
using namespace std;
const int N=1010;
int f[N],a[N];
int main()
{
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    for(int i=1;i<=n;i++)
    {
        f[i]=1;
        for(int j=1;j<i;j++)
        {
            if(a[j]<a[i])
            f[i]=max(f[i],f[j]+1);
            //1 2 4 2 5 3
            //1 2 3 2 4 3
        }
    }
    int res=0;
    for(int i=1;i<=n;i++)res=max(res,f[i]),cout<<f[i]<<' ';
    cout<<res<<endl;
    return 0;
}

• 打印路径的一般方法，每次更新的时候用per数组记录该点下标的前驱节点，然后知道有多少个点后，找到最大点的下标，依次静态链表倒序输出

    for(int i=1;i<=n;i++)
    {
        f[i]=1;
        per[i]=0;
        for(int j=1;j<i;j++)
        {
            if(a[j]<a[i])
            {
                if(f[i]<f[j]+1)
                {
                    f[i]=f[j]+1;
                    per[i]=j;
                }
            }
        }
    }
    int k=1;
    for(int i=1;i<=n;i++)
        if(f[k]<f[i])
        k=i;
    int len=f[k];
    while(len--)
    {
        cout<<a[k]<<' ';
        k=per[k];
    }

• 更大数据范围选择二分：遍历a[i]在p数列中找小于a[i]的最大值(二分出其下标也就是剩余有序数列的长度)，后面的全部删去，更新此时的长度，依次遍历所有值，最大值就是答案
• 代码

#include<iostream>
using namespace std;
const int N=100010;
int a[N],p[N];
int main()
{
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    int len=0;
    for(int i=1;i<=n;i++)
    {
        int l=0,r=len;
        while(l<r)
        {
            int mid=l+r+1>>1;
            if(p[mid]<a[i])l=mid;
            else r=mid-1;
        }
        len=max(len,r+1);
        p[r+1]=a[i];
    }
    cout<<len<<endl;
    return 0;
}

• 3、最长公共子序列

• 如果遍历到的字符不相等，则选择a的前i-1和b的前j或者a的前i和b的前j-1比较大小，否则更新f[i][j]等与f[i-1][j-1]+1
• 代码

#include<iostream>
using namespace std;
const int N=1010;
char a[N],b[N];
int f[N][N];
int main()
{
    int n,m;
    cin>>n>>m;
    scanf("%s%s",a+1,b+1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            f[i][j]=max(f[i-1][j],f[i][j-1]);
            if(a[i]==b[j])f[i][j]=max(f[i][j],f[i-1][j-1]+1);
        }
    }
    cout<<f[n][m]<<endl;
    return 0;
}

• 4、编辑距离

• 对于两个字符串要是前面都相等删除一个相当于该点更新为f[i][j]=min(f[i][j]，f[i-1][j]+1)
• 增加相当于B串减少f[i][j]=min(f[i][j]，f[i][j-1]+1)
• 合并之后相当于第一次更新
• f[i][j]=min(f[i-1][j]+1，f[i][j-1]+1)
• 如果该点两个字符相等，则更新f[i][j]=min(f[i][j]，f[i-1][j-1]+1)
• 所以常数更新过程为
• f[i][j]=min(f[i-1][j]+1，f[i][j-1]+1)
• f[i][j]=min(f[i][j]，f[i-1][j-1]+(a[i]!=b[j])
• 与字符串匹配问题的更新过程基本相似
• 先判断能不能增删，再判断能不能改
• 代码

#include<iostream>
using namespace std;
const int N=1010;
char a[N],b[N];
int f[N][N];
int n,m;
int main()
{
    cin>>n>>a+1>>m>>b+1;
    for(int i=1;i<=n;i++) f[i][0]=i;
    for(int i=1;i<=m;i++) f[0][i]=i;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            f[i][j]=min(f[i-1][j]+1,f[i][j-1]+1);
            if(a[i]==b[j])f[i][j]=min(f[i][j],f[i-1][j-1]);
            else f[i][j]=min(f[i][j],f[i-1][j-1]+1);
        }
    cout<<f[n][m]<<endl;
    return 0;
}

• 逐段匹配代码

#include<iostream>
#include<string.h>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=11,M=1010;
char str[M][N];
int f[N][N];
int n,m;
int edit_dis(char a[],char b[])
{
    int lena=strlen(a+1),lenb=strlen(b+1);
    //cout<<lena<<' '<<lenb<<endl;
    for(int i=1;i<=lena;i++) f[i][0]=i;
    for(int i=1;i<=lenb;i++) f[0][i]=i;
    for(int i=1;i<=lena;i++)
        for(int j=1;j<=lenb;j++)
        {
            f[i][j]=min(f[i][j-1]+1,f[i-1][j]+1);
            f[i][j]=min(f[i][j],f[i-1][j-1]+(a[i]!=b[j]));
        }
    return f[lena][lenb];
}
int main()
{
    cin>>n>>m;
    for(int i=0;i<n;i++) scanf("%s",str[i]+1);
    while(m--)
    {
        char s[N];
        int limit;
        scanf("%s%d",s+1,&limit);
        int res=0;
        for(int i=0;i<n;i++)
        {
            if(edit_dis(str[i],s)<=limit)
                res++;
        }
        printf("%d\n",res);
    }
    return 0;
}

• 2021-08-02写于商丘