POJ 2127

#include <iostream>
#define MAXN 501
using namespace std;

int a[MAXN],b[MAXN],ans[MAXN];

int GCIS(int l1, int *a, int l2, int *b);

int main()
{
    //freopen("acm.acm","r",stdin);
    int l_1;
    int l_2;
    int i;
    int j;
    int ans_max;
//    cin>>l_1;
    scanf("%d",&l_1);

    for(i = 0; i < l_1; ++ i)
    {
    //    cin>>a[i];
        scanf("%d",&a[i]);
    }

//    cin>>l_2;
    scanf("%d",&l_2);
    for(i = 0; i < l_2; ++ i)
    {
        //cin>>b[i];
        scanf("%d",&b[i]);
    }

    cout<<(ans_max = GCIS(l_1,a,l_2,b) )<<endl;

    for(i = 0; i < ans_max; ++ i)
    {
        cout<<ans[i]<<" ";
    }
    cout<<endl;

}


/////////////////////////////////////
//最长公共上升子序列～
/////////////////////////////////////
int GCIS(int l1, int *a, int l2, int *b)//ans[0...DP[max]-1]为序列,最长公共递增子序列！
{
    int f[MAXN+1][MAXN+1];
    int DP[MAXN+1];
    int i,j,k,max; 
    memset(f,0,sizeof(f));
    memset(DP,0,sizeof(DP));
    for (i=1;i<=l1;i++)
    {
        k=0;
        for(int kk = 1;kk <= l2;++ kk)
        {
            f[i][kk] = f[i-1][kk];
        }
        for (j=1;j<=l2;j++)
        {
            if(b[j-1] < a[i-1] && DP[j] > DP[k])
                k=j;
            if(b[j-1]==a[i-1]&&DP[k]+1>DP[j])
            {
                DP[j]=DP[k]+1;
                f[i][j]=i*(l2+1)+k; 
            }    
        } 
    }
    max=0;
    for(i=1;i<=l2;i++)
    {
        if (DP[i]>DP[max])
        max=i;
    }
    i=l1*l2+l1+max; 
    for(j = DP[max];j > 0;j --)
    {
        ans[j-1] = b[i%(l2+1)-1]; 
        i=f[i/(l2+1)][i%(l2+1)];
    } 
    return DP[max];
}
///////////////////////////////////////////////////
//返回值是子序列的容量，容器为ans[]
///////////////////////////////////////////////////

关注我的公众号，当然，如果你对Java, Scala, Python等技术经验，以及编程日记，感兴趣的话。 

技术网站地址： vmfor.com