Time Limit: 1000MS | Memory Limit: 32768KB | 64bit IO Format: %I64d & %I64u |
Description
This is a problem from ZOJ 2432.To make it easyer,you just need output the length of the subsequence.
Input
Each sequence is described with M - its length (1 <= M <= 500) and M integer numbers Ai (-2^31 <= Ai < 2^31) - the sequence itself.
Output
output print L - the length of the greatest common increasing subsequence of both sequences.
Sample Input
1 5 1 4 2 5 -12 4 -12 1 2 4
Sample Output
2
认识到了最长公共上升子序列,我发现自己dp基础一塌糊涂,
后来仔细研究,终于明白它是怎么回事,下面介绍的代码
现在有两个字符串A[1.....n],B[1.....m]
dp[i][j]代表着A序列的前i个与以B[j]为结尾的最长公共上升子序列
这里我要说清楚,当时就是这里没有看明白,也许是自己文字理解能力不够,
这里就相当于是最长上升子序列一样,前i个不用将,以B[j]为结尾的意思是
就像最长上升子序列dp[i]=dp[j]+1,这里的意思是以i这个数结尾的最长上升子序列
上述的代表着A序列的前i个与以B[j]为结尾的最长公共上升子序列就是这个意思
然后就是不断更新了,用一个Max来记录比B[j]小的最长公共上升子序列
后来仔细研究,终于明白它是怎么回事,下面介绍的代码
现在有两个字符串A[1.....n],B[1.....m]
dp[i][j]代表着A序列的前i个与以B[j]为结尾的最长公共上升子序列
这里我要说清楚,当时就是这里没有看明白,也许是自己文字理解能力不够,
这里就相当于是最长上升子序列一样,前i个不用将,以B[j]为结尾的意思是
就像最长上升子序列dp[i]=dp[j]+1,这里的意思是以i这个数结尾的最长上升子序列
上述的代表着A序列的前i个与以B[j]为结尾的最长公共上升子序列就是这个意思
然后就是不断更新了,用一个Max来记录比B[j]小的最长公共上升子序列
下面是代码中的细节讲解
Max=0;这里的意思是最开始的最长公共上升子序列都为零
for(int j=1; j<=m; j++) {
dp[i][j]=dp[i-1][j];将前i-1位置的值直接复制
if(A[i]>B[j])Max=max(dp[i-1][j],Max);这里的意思是如果A[i]>B[j]就要更新Max
为什么,因为最长上升子序列,这是为下面这一行代码做准备的,因为往后走时,如果存在A[i]==B[j]
那么,我必须要有小于B[j]的最长公共上升子序列,这里我们可以知道如果A[i]==B[j]那么上面那一行的代码
求到的Max不正是我要的吗,因为A[i]==B[j],所以上面那一行相当于是变成了B[k(k是某一个,肯定不会与上面的j相等)]>B[j]
那么上面的代码单独提出来
:(其中两个j是不相等的)if(B[j]>B[j])Max=max(dp[i-1][j],Max);不正是求解我需要的小于B[j]的最长公共上升子序列
如此所有的代码都变得非常合理了。
if(A[i]==B[j])dp[i][j]=Max+1;
}
Max=0;这里的意思是最开始的最长公共上升子序列都为零
for(int j=1; j<=m; j++) {
dp[i][j]=dp[i-1][j];将前i-1位置的值直接复制
if(A[i]>B[j])Max=max(dp[i-1][j],Max);这里的意思是如果A[i]>B[j]就要更新Max
为什么,因为最长上升子序列,这是为下面这一行代码做准备的,因为往后走时,如果存在A[i]==B[j]
那么,我必须要有小于B[j]的最长公共上升子序列,这里我们可以知道如果A[i]==B[j]那么上面那一行的代码
求到的Max不正是我要的吗,因为A[i]==B[j],所以上面那一行相当于是变成了B[k(k是某一个,肯定不会与上面的j相等)]>B[j]
那么上面的代码单独提出来
:(其中两个j是不相等的)if(B[j]>B[j])Max=max(dp[i-1][j],Max);不正是求解我需要的小于B[j]的最长公共上升子序列
如此所有的代码都变得非常合理了。
if(A[i]==B[j])dp[i][j]=Max+1;
}
/*
Author: 2486
Memory: 2416 KB Time: 0 MS
Language: G++ Result: Accepted
VJ RunId: 4177842 Real RunId: 14209604
Public: No Yes
*/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=500+5;
int A[maxn],B[maxn];
int dp[maxn][maxn];
int t,n,m;
int main() {
scanf("%d",&t);
while(t--) {
scanf("%d",&n);
for(int i=1; i<=n; i++)scanf("%d",&A[i]);
scanf("%d",&m);
for(int i=1; i<=m; i++)scanf("%d",&B[i]);
int Max;
memset(dp,0,sizeof(dp));
for(int i=1; i<=n; i++) {
Max=0;
for(int j=1; j<=m; j++) {
dp[i][j]=dp[i-1][j];
if(A[i]>B[j])Max=max(dp[i-1][j],Max);
if(A[i]==B[j])dp[i][j]=Max+1;
}
}
Max=0;
for(int i=1; i<=m; i++) {
Max=max(Max,dp[n][i]);
}
printf("%d\n",Max);
if(t)printf("\n");
}
return 0;
}