解题思路:
状态表示: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;
}