链接: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; }