1120: 病毒+csuOJ+最长公共升子序列

1120: 病毒

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 259  Solved: 96
[Submit][Status][Web Board]

Description

你有一个日志文件，里面记录着各种系统事件的详细信息。自然的，事件的时间戳按照严格递增顺序排列（不会有两个事件在完全相同的时刻发生）。

遗憾的是，你的系统被病毒感染了，日志文件中混入了病毒生成的随机伪事件（但真实事件的相对顺序保持不变）。备份的日志文件也被感染了，但由于病毒采用的随机感染方法，主日志文件和备份日志文件在感染后可能会变得不一样。

给出被感染的主日志和备份日志，求真实事件序列的最长可能长度。

Input

输入第一行为数据组数T (T<=100)。每组数据包含两行，分别描述感染后的主日志和备份日志。

每个日志文件的格式相同，均为一个整数n (1<=n<=1000)（代表感染后的事件总数）和n 个不超过100,000的正整数（表示感染后各事件的时间戳）。

注意，感染后可能会出现时间戳完全相同的事件。

Output

对于每组数据，输出真实事件序列的最长可能长度。

Sample Input

1
9 1 4 2 6 3 8 5 9 1
6 2 7 6 3 5 1

Sample Output

3

解决方案：首先，在a[i]!=b[j]的时候有F[i][j]=F[i-1][j]。为什么呢？因为F[i][j]是以b[j]为结尾的LCIS，如果F[i][j]>0那么就说明a[1]..a[i]中必然有一个字符a[k]等于b[j]（如果F[i][j]等于0呢？那赋值与否都没有什么影响了）。因为a[k]!=a[i]那么a[i]对F[i][j]没有贡献，于是我们不考虑它照样能得出F[i][j]的最优值。所以在a[i]!=b[j]的情况下必然有F[i][j]=F[i-1][j]。这一点参考LCS的处理方法。  那如果a[i]==b[j]呢？首先，这个等于起码保证了长度为1的LCIS。然后我们还需要去找一个最长的且能让b[j]接在其末尾的LCIS。之前最长的LCIS在哪呢？首先我们要去找的F数组的第一维必然是i-1。
因为i已经拿去和b[j]配对去了，不能用了。并且也不能是i-2，因为i-1必然比i-2更优。第二维呢？那就需要枚举b[1]..b[j-1]了，因为你不知道这里面哪个最长且哪个小于b[j]。这里还有一个问题，可不可能不配对呢？也就是在a[i]==b[j]的情况下，需不需要考虑F[i][j]=F[i-1][j]的决策呢？答案是不需要。因为如果b[j]不和a[i]配对，那就是和之前的a[1]..a[j-1]配对（假设
F[i-1][j]>0，等于0不考虑），这样必然没有和a[i]配对优越。（为什么必然呢？因为b[j]和a[i]配对之后的转移是max(F[i-1][k])+1，而和之前的i`配对则是max(F[i`-1][k])+1。显然有F[i][j]>F[i`][j],i`>i） 于是我们得出了
状态转移方程： 
 a[i]!=b[j]:   F[i][j]=F[i-1][j] 
 a[i]==b[j]:   F[i][j]=max(F[i-1][k])+1 1<=k<=j-1&&b[j]>b[k] 
 不难看到，这是一个时间复杂度为O(n^3)的DP，离平方还有一段距离。

code：
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1004;
int num1[maxn];
int num2[maxn];
int dp[maxn][maxn];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        int n1,n2;
        scanf("%d",&n1);
        for(int i=1;i<=n1;i++){
            scanf("%d",&num1[i]);
        }
        scanf("%d",&n2);
        for(int i=1;i<=n2;i++){
            scanf("%d",&num2[i]);
        }
        memset(dp,0,sizeof(dp));
        int mmax=0;
        for(int i=1;i<=n1;i++){
            for(int j=1;j<=n2;j++){
                if(num1[i]!=num2[j]) dp[i][j]=dp[i-1][j];
                if(num1[i]==num2[j]){
                    dp[i][j]=1;
                    for(int k=1;k<j;k++){
                        if(num1[i]>num2[k]){
                            dp[i][j]=max(dp[i][j],dp[i-1][k]+1);
                        }
                    }
                }
                if(dp[i][j]>mmax) mmax=dp[i][j];
            }
        }
        printf("%d\n",mmax);
    }
    return 0;
}

但是，这个算法最关键的是，如果按照一个合理的递推顺序，max(F[i-1][k])的值我们可以在之前访问F[i][k]的时候通过维护更新一个max变量得到。怎么得到呢？首先递推的顺序必须是状态的第一维在外层循环，第二维在内层循环。也就是算好了F[1][len(b)]再去算F[2][1]。 如果按照这个递推顺序我们可以在每次外层循环的开始加上令一个max变量为0，然后开始内层循环。当a[i]>b[j]的时候令max=F[i-1][j]。如果循环到了a[i]==b[j]的时候，则令F[i][j]=max+1。  最后答案是F[len(a)][1]..F[len(a)][len(b)]的最大值。 

code：
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=1004;
int num1[maxn];
int num2[maxn];
int dp[maxn][maxn];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n1,n2;
        scanf("%d",&n1);
        for(int i=1; i<=n1; i++)
        {
            scanf("%d",&num1[i]);
        }
        scanf("%d",&n2);
        for(int i=1; i<=n2; i++)
        {
            scanf("%d",&num2[i]);
        }
        memset(dp,0,sizeof(dp));
        int mmax=0;
        int Max=0;
        for(int i=1; i<=n1; i++)
        {
            Max=0;
            for(int j=1; j<=n2; j++)
            {
                dp[i][j]=dp[i-1][j];
                if(num1[i]>num2[j]&&Max<dp[i-1][j]) Max=dp[i-1][j];
                if(num1[i]==num2[j]) dp[i][j]=Max+1;
                if(dp[i][j]>mmax) mmax=dp[i][j];
            }

        }
        printf("%d\n",mmax);
    }
    return 0;
}