信息学奥赛一本通：1306：最长公共子上升序列

1306：最长公共子上升序列

【题目描述】

给定两个整数序列，写一个程序求它们的最长上升公共子序列。

当以下条件满足的时候，我们将长度NN的序列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。在第二行，输出该子序列。如果有不止一个符合条件的子序列，则输出任何一个即可。

【输入样例】

5
1 4 2 5 -12
4
-12 1 2 4

【输出样例】

2
1 4

//输入两个数组a和b，以保存输入的数，f是a的第i个数与b的j个数组成的最长值 
//本题用公共序列和最长上升序列来做
//先求出在a中的哪几个数，和b中的数相同，放在c数组中，再把c中的数取上升的最大值
//有一个相同时，只取第一个，后面有一个相同的不取了，这样保证后面有相同的一个的，不在c中； 
#include<bits/stdc++.h>
using namespace std;
int n,m;
int a[505],b[505],f[505][505],c[505];
int main()
{
    int sum=1,t=1;
    scanf("%d",&n);for(int i=1;i<=n;i++)scanf("%d",&a[i]);//输入a 
    scanf("%d",&m);for(int j=1;j<=m;j++)scanf("%d",&b[j]);//输入b 
    for(int i=1;i<=n;i++)//求最大公共序列 把a中的和b公共序列找出来，然后再弄个最大上升序列就可了 
        for(int j=1;j<=m;j++)
        {
            if(a[i]==b[j])//关键是找到从第一个相同的数，如例子中的1，然后怎么把4和2都能放在c中 
            {
            f[i][j]=f[i-1][j-1]+1;
            if(sum==1)//第一个相同的 a,b都是从第一位开始，这个相同的必用到 
            c[sum++]=a[i];//把a的值放在c中 
            if(f[i-1][j]==f[i][j-1]&&f[i-1][j-1])//是第二个及以上相同的数时，都保存下来 
            c[sum++]=a[i]; 
             }
            else
            f[i][j]=max(f[i-1][j],f[i][j-1]);
        }

    //f[n][m]是公共序列，但不一定是上升的最大序列 

    //用例子中，每个多加一个值在最后，如弄个3
    //    6
   //1 4 2 5 -12 3
   //5
   //-12 1 2 4 3
   //反复运算，得到上面的公式 
//下面是上升序列的最大值和分别是哪几个数 如1259：【例9.3】求最长不下降序列 
    for(int i=1;i<=sum;i++) b[i]=1;//用b来做计算的 
    //下面是公式：上升或是下降序列的 
    for(int i=2;i<sum;i++)
        for(int j=1;j<i;j++)
        if(c[i]>c[j]&&b[i]<b[j]+1)
        {
            b[i]=b[j]+1;
            t=max(t,b[i]);
            a[i]=j;//mark it; 
        }
        for(int i=1;i<sum;i++)
           if(t==b[i]){//找到最大值时 
                sum=i;//保存是几号的 
                break;//出来 
            }
        for(int i=t;i>=1;i--)//这是有几个要找出来的，上升的序列数 
        {
            b[i]=c[sum];//反向存，放在b中 
            sum=a[sum];//a中是标记的指向的数 
        }
/*    for(int i=1;i<=n;i++)//看看表上的数 
    {
        
        for(int j=0;j<=m;j++)
        cout<<f[i][j]<<" ";
        cout<<endl;
    } */
    if(f[n][m]==0)//如果没有公共序列，就输出这个，over了。 
    {
        cout<<0;
        return 0;
    }
    cout<<t<<endl;//有公共序列，就从这开始输出 
    for(int i=1;i<=t;i++)//从小到大的顺序 
    cout<<b[i]<<" ";
    return 0;
 }