POJ 2127

 1 #include <iostream>
 2 #define MAXN 501
 3 using namespace std;
 4 
 5 int a[MAXN],b[MAXN],ans[MAXN];
 6 
 7 int GCIS(int l1, int *a, int l2, int *b);
 8 
 9 int main()
10 {
11     //freopen("acm.acm","r",stdin);
12     int l_1;
13     int l_2;
14     int i;
15     int j;
16     int ans_max;
17 //    cin>>l_1;
18     scanf("%d",&l_1);
19 
20     for(i = 0; i < l_1; ++ i)
21     {
22     //    cin>>a[i];
23         scanf("%d",&a[i]);
24     }
25 
26 //    cin>>l_2;
27     scanf("%d",&l_2);
28     for(i = 0; i < l_2; ++ i)
29     {
30         //cin>>b[i];
31         scanf("%d",&b[i]);
32     }
33 
34     cout<<(ans_max = GCIS(l_1,a,l_2,b) )<<endl;
35 
36     for(i = 0; i < ans_max; ++ i)
37     {
38         cout<<ans[i]<<" ";
39     }
40     cout<<endl;
41 
42 }
43 
44 
45 /////////////////////////////////////
46 //最长公共上升子序列～
47 /////////////////////////////////////
48 int GCIS(int l1, int *a, int l2, int *b)//ans[0...DP[max]-1]为序列,最长公共递增子序列！
49 {
50     int f[MAXN+1][MAXN+1];
51     int DP[MAXN+1];
52     int i,j,k,max; 
53     memset(f,0,sizeof(f));
54     memset(DP,0,sizeof(DP));
55     for (i=1;i<=l1;i++)
56     {
57         k=0;
58         for(int kk = 1;kk <= l2;++ kk)
59         {
60             f[i][kk] = f[i-1][kk];
61         }
62         for (j=1;j<=l2;j++)
63         {
64             if(b[j-1] < a[i-1] && DP[j] > DP[k])
65                 k=j;
66             if(b[j-1]==a[i-1]&&DP[k]+1>DP[j])
67             {
68                 DP[j]=DP[k]+1;
69                 f[i][j]=i*(l2+1)+k; 
70             }    
71         } 
72     }
73     max=0;
74     for(i=1;i<=l2;i++)
75     {
76         if (DP[i]>DP[max])
77         max=i;
78     }
79     i=l1*l2+l1+max; 
80     for(j = DP[max];j > 0;j --)
81     {
82         ans[j-1] = b[i%(l2+1)-1]; 
83         i=f[i/(l2+1)][i%(l2+1)];
84     } 
85     return DP[max];
86 }
87 ///////////////////////////////////////////////////
88 //返回值是子序列的容量，容器为ans[]
89 ///////////////////////////////////////////////////
关注我的公众号，当然，如果你对Java, Scala, Python等技术经验，以及编程日记，感兴趣的话。
技术网站地址： vmfor.com