Acwing 272. 最长公共上升子序列（线性DP）

解题思路：

状态表示：f[i][j]表示a[1~i]和b[1~j]中，以b[j]结尾的最长公共子序列的长度

状态转移：首先肯定要设置两重循环，分别循环i和j下标，每次循环到一个i和j时，会有两种情况出现，如果a[i]!=b[j],那么当前的f[i][j]=f[i-1][j]，因为a[i]和b[j]已经是不相等了，所以以b[j]结尾a[1~i],b[1~j]的最长公共子序列的最长长度就等于从a[1~i-1]，b[1~j]中选出的最长长度，第二种情况就是a[i]==b[j]，那么当前的a[i]和b[j]就构成了一个长度为1的公共子序列，那它可以从前面的哪些状态转移过来呢？我们再次设置一个变量k,1<=k<j,从众多的f[i-1][k]中选出一个最大的记为maxv，所以f[i][j]=max(f[i-1][j],maxv+1),当然了在选择k的时候是有条件的，也就是说并不是每一个k都可以供我们选择，因为我们要选的是最长的上升公共子序列，我们前面说的都是公共子序列，那上升这个性质就要在这里体现出来了，我们选的k必须满足，b[k]<b[j]才可以，代码如下：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=3010;
int f[N][N];
int n;
int a[N],b[N];
int main()
{
    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++)
        {
            f[i][j]=f[i-1][j];
            if(a[i]==b[j])
            {
               int maxv=1;
                for(int k=1;k<j;k++)
                if(b[k]<a[i])
                maxv=max(maxv,f[i-1][k]+1);
                f[i][j]=maxv;
            }
            
        }
    }
    int res=0;
    for(int i=1;i<=n;i++)
    res=max(res,f[n][i]);
    cout<<res<<endl;
    return 0;
}

我们分析一下这个代码，不难发现这是一个O(n^3)时间复杂度的做法，因为我们的n最大是3000，所以n^3一定是会超时的，所以现在我们的任务就是要对这个算法进行一个优化，仔细观察状态转移的部分，每当a[i]==b[j]的时候，我们都会把前面的f[i-1][k],1<=k<j都给遍历一遍，每次都会重复遍历到前面的状态，所以我们可以在每一步的j中，对maxv进行一个更新，一旦我们需要用到前面状态的最大值的时候直接拿过来用就行了，那如何在遍历j的过程中对maxv进行更新呢？每当我们发现b[j]<a[i]的时候，maxv=max(maxv,f[i-1][j]+1),因为如果后来可能会用到maxv的话必然会有a[i]==b[j1]这里的b[j]<a[i]，就相当于是b[j]<b[j1]了，优化后的代码如下：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=3010;
int f[N][N];
int n;
int a[N],b[N];
int main()
{
    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++)
    {
        int maxv=1;
        for(int j=1;j<=n;j++)
        {
            f[i][j]=f[i-1][j];
            if(a[i]==b[j])
            f[i][j]=max(f[i][j],maxv);
            if(b[j]<a[i])
            maxv=max(maxv,f[i-1][j]+1);
        }
    }
    int res=0;
    for(int i=1;i<=n;i++)
    res=max(res,f[n][i]);
    cout<<res<<endl;
    return 0;
}