poj 2127 zoj 2432 codeforces 10D 最长上升公共子序列

来源：互联网发布：渣浆泵选型软件编辑：程序博客网时间：2024/04/29 09:33

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