链接:http://poj.org/problem?id=2127
也不知怎么的就一下子一起搜到了这三个题,结果被坑了,在CF上AC的代码拿到poj 和 zoj都不能A,然后开始疯狂的造数据,终于被我造出了一组,然后发现开一维数组搞不定这个问题(反正我是搞不定了),于是还是老老实实的开了二维的dp,递归输出路径
这种DP应该很经典也很简单吧,dp[i][j]表示a数组的第i个数与b数组的第j个数结尾的最长公共递增子序列的最大长度,然后转移的时候接在 前面比a[i]小的,且长度最长的序列的后面
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 510;
int dp[maxn][maxn];
pair<int,int> pre[maxn][maxn];
pair<int,pair<int,int> > lis[maxn];
int a[maxn],b[maxn];
void print(int x,int y){
if(x>0 && y>0 && dp[x][y]>0){
print(pre[x][y].first,pre[x][y].second);
printf("%d ",a[x]);
}
}
int main()
{
int m,n,x,y;
scanf("%d",&m);
for(int i=1;i<=m;i++) scanf("%d",&a[i]);
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&b[i]);
fill(lis,lis+n+1,make_pair(0,make_pair(-1,-1)));
x=y=0;
for(int i=1;i<=m;i++){
pair<int,pair<int,int> > k=make_pair(0,make_pair(-1,-1));
for(int j=1;j<=n;j++){
if(b[j]<a[i]){
k=max(k,lis[j]);
}
if(a[i]!=b[j]){
dp[i][j]=0;
pre[i][j]=make_pair(-1,-1);
}
else {
dp[i][j]=k.first+1;
pre[i][j]=k.second;
lis[j]=make_pair(dp[i][j],make_pair(i,j));
}
if(dp[i][j]>dp[x][y]) {
x=i;
y=j;
}
}
}
printf("%d\n",dp[x][y]);
print(x,y);
return 0;
}