【2020模拟赛day6】C. 寻找序列

 

这道题是LCS和LIS的结合，考虑设计状态仍然是f[i][j]表示a串前i位和b串前j位能够得到的最长公共上升子序列

然后呢考虑递推

当a[i]==b[j] 枚举一维k 如果a[i]>b[k] f[i][j]=max(f[i][j],f[i-1][k]+1);

当a[i]!=b[j] 直接f[i][j=f[i-1][j];

 

这样的算法是立方级的

 

观察发现，在计算i的时候，i-1行的k的大小已经确定，所以我们可以在计算时维护一个maxx去记录它

 

至于方案，也就是要求我们记录一下每个值得到的路径记录pre[i][j]表示这个位置是由哪里推出来的，每次f[i][j]有更优解的时候也更新一下pre[i][j]就好了

 

代码

#include<bits/stdc++.h>
using namespace std;
const int maxn=5005;
int n,m,a[maxn],b[maxn],f[maxn][maxn],pre[maxn][maxn],r[maxn];
int main()
{
	freopen("okarin.in","r",stdin);
	freopen("okarin.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	scanf("%d",&m);
	for(int j=1;j<=m;j++) scanf("%d",&b[j]);
	int maxx=0,pos;
	for(int i=1;i<=n;i++)
	{
		maxx=0; pos=0;
		for(int j=1;j<=m;j++)
		{
			f[i][j]=f[i-1][j];
			pre[i][j]=j;
			if(a[i]==b[j])
				if(f[i][j]<maxx+1)
				{
					f[i][j]=maxx+1;
					pre[i][j]=pos;	
				}
			if(a[i]>b[j])
				if(f[i-1][j]>maxx)
				{
					maxx=f[i-1][j];
					pos=j;
				}
		}	
	} 
	maxx=0; int last=0;
	for(int i=1;i<=m;i++)
		if(maxx<f[n][i])
		{
			maxx=f[n][i];
			last=i;
		}
	printf("%d\n",maxx);
	int tot=0; int i=n,j=last;
	if(maxx)
	{
		while(i||j)
		{
			if(pre[i][j]!=j)
				r[++tot]=b[j];
			j=pre[i][j];
			i--;
		}
		for(int i=tot;i>=1;i--)	printf("%d ",r[i]);
	}
	return 0;
}