最长公共上升子序列

题目链接
题意：求最长公共上升子序列

题解：我们分析题意求最长公共上升子序列。

我们首先抛开模型：
本题有两个限制条件，一个上升，一个公共，所以我们用二维表示
（任何的dp题我们都要赋予他一个状态含义）
我们这时联系模型，上升的通常做法就是f[i]就是以a[i]为结尾然后求最大值，求公共的话就是二维数组，所以我们f[i][j]就表示以b[j]为结尾，然后上升公共子序列的最大值，每次求的时候就是包含a[i]和不包含a[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 (a[i] > b[k])
                    maxv = max(maxv, f[i - 1][k] + 1);
            f[i][j] = max(f[i][j], maxv);
        }
    }
}

下面是改进循环的代码（很简单，这里不再赘述）：

//f[i][j]表示以b[j]为结尾的最长公共上升子序列的最大值
//状态转移方程(包含a[i]，不包含a[i])
//f[i]][j]=f[i-1][j]//不包含a[i-1]
//包括的前提是s[i]==s[j]这时我们就要暴力了
//总结：无论做什么dp题，多么难或简单，都要首先赋予他一个含义

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=3010;
int f[N][N];
int a[N],b[N];
int 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++)
    {
        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(a[i]>b[j]) 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;
}

这里很多同学可能会疑问这里为什么求最长公共序列的时候没有用到那个f[i][j-1]，我们在这里可以分析一下，最长公共子序列的状态表达式f[i][j]表示以a[i]和b[j]为结尾的串里面含有的最大公共子序列。所以我们会分情况讨论，包含a[i]不包含b[j]（f[i][j - 1]）；不包含a[i]包含b[j]（f[i-1][j]）；两个都包含f[i-1][j-1]+1(a[i]==b[j])；这样分的话两层循环就可以，而因为本题呢，还需要求最长上升的子序列，所以我们必须确定以某一个串a[i]或者b[j]为最后一个来确定最长的公共子序列，这样做的后果就是多加了一维。

但是这样做的结果是正确的，我们可以试着用分析本题的方法去分析一下最长公共子序列：
f[i][j]表示以b[j]为结尾的最长公共子序列的个数， f[i][j]=f[i-1][j]（不包含a[i]）包含a[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);
            maxv=max(maxv,f[i-1][j]+1);//将此处去掉
        }
    }

可以找一道模板题稍做改动之后交一下代码，可以找一下这个
戳我
下面是AC代码：

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

我们从上面的例子可以得出，对一个f[i][j]赋予不同的含义，我们可以得到不同的转移方程，最终的答案是相同的。