动态规划 最长公共子序列

来源:互联网 发布:文艺复兴 知乎 编辑:程序博客网 时间:2024/06/06 08:33
描述
给定两个整数序列,写一个程序求它们的最长上升公共子序列。
当以下条件满足的时候,我们将长度为N的序列S1 , S2 , . . . , SN 称为长度为M的序列A1 , A2 , . . . , AM 的上升子序列:

存在 1 <= i1 < i2 < . . . < iN <= M ,使得对所有 1 <= j <=N,均有Sj = Aij,且对于所有的1 <= j < N,均有Sj < Sj+1

输入
每个序列用两行表示,第一行是长度M(1 <= M <= 500),第二行是该序列的M个整数Ai (-231 <= Ai < 231 )
输出
在第一行,输出两个序列的最长上升公共子序列的长度L。在第二行,输出该子序列。如果有不止一个符合条件的子序列,则输出任何一个即可。
样例输入
51 4 2 5 -124-12 1 2 4
样例输出

2

1 4


分析:
既要找最长子序列,还有一个限定条件是公共的,所以想到用结构体存放;

代码:
#include <bits/stdc++.h>
using namespace std;
const int maxn=510;
int m1,m2,a1[maxn],a2[maxn];
struct node{
    int len;
    vector<int> iv;
}ans[maxn],cur;
int main(){
    scanf("%d",&m1);
    for(int i=0;i<m1;i++)scanf("%d",&a1[i]);
    scanf("%d",&m2);
    for(int i=0;i<m2;i++)scanf("%d",&a2[i]);
    for(int i=0;i<maxn;i++) ans[i].len=0;
    for(int i2=0;i2<m2;i2++){
        cur.len=0; cur.iv.clear();
        for(int i1=0;i1<m1;i1++){
            if(a2[i2]>a1[i1]&&ans[i1].len>cur.len) cur=ans[i1];
            if(a2[i2]==a1[i1]){
                ans[i1]=cur; ans[i1].len++;
                ans[i1].iv.push_back(a1[i1]);
            }
        }
    }
    int p=0;
    for(int i=1;i<m1;i++){
        if(ans[p].len<ans[i].len) p=i;
    }
    printf("%d\n",ans[p].len);
    if(ans[p].iv.size()){
        printf("%d",ans[p].iv[0]);
        for(int i=1;i<ans[p].iv.size();i++)
            printf(" %d",ans[p].iv[i]);
    }
    printf("\n");
    return 0;
}
0 0
原创粉丝点击