bzoj1264 [AHOI2006]基因匹配Match (树状数组优化DP)

bzoj1264 [AHOI2006]基因匹配Match

原题地址：http://www.lydsy.com/JudgeOnline/problem.php?id=1264

题意：
给定n，给定两个长为5*n的序列，两个序列中1~n这些数都刚好出现5次。
求两个序列的最长公共子序列。

数据范围
N<=20000

题解：

一般的求最长公共子序列时空都是O(n^2)
dp[i][j]表示a的前i位与b的前j位形成的最长公共子序列长度。
若a[i]==b[j] dp[i][j]=dp[i-1][j-1]+1
否则 dp[i][j]=max(dp[i-1][j],dp[i][j-1])
实际上只有当 a[i]==b[j]时，整个的值才会++

由 “ a[i]==b[j]时，整个的值才会++”这个关键点转换思路 ，
因为同一个值的pos只有5个 ，可以只取这些关键操作。

若我们让dp[i]表示 以i结尾的 a的前i位与b形成的最长公共子序列 长度
b[j]读取到一个x时，若a[i]==b[j]的位置为pos，dp[pos]都可以用max(dp[1~pos-1]+1)更新
（例如两个序列：
3 2 1 6 7 8 1
8 2 1 7 6 1 1
最长2161四位
3 2 1 6 7 8 1
8 2 1 7 6 1 1 (dp[6~7]用0+1=1更新)

3 2 1 6 7 8 1
8 2 1 7 6 1 1 (dp[2~7]用0+1=1更新)

3 2 1 6 7 8 1
8 2 1 7 6 1 1 (dp[3~7]用1+1=2更新,dp[7]用1+1=2更新)
…
这样一个过程。）

而这道题同一个值的pos只有5个 ，故可以用类似的操作

于是可以用树状数组维护dp值，每次就可以查询dp[1~pos-1]的最大值来更新pos，每次的pos只有五个。

（实际上当求两个n的排列的序列的最长公共子序列时，可以转化为LIS来做。由于此题的每个值也只出现五次，故可以用这个类似LIS的树状数组写法）

（传统问题的加强版：从不同点入手；思考原问题的关键点，是否可以只保留这些，优化冗余部分）

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int inf=0x3f3f3f3f;
const int N=100005;
const int M=20005;
int n,m,dp[N],a[N],b[N],c[N],pos[M][6],cnt[M];
void add(int x,int val)
{
    for(int i=x;i<=n;i=i+(i&(-i)))
    c[i]=max(c[i],val);
}
int query(int x)
{
    int ans=-inf;
    for(int i=x;i;i=i-(i&(-i)))
    ans=max(ans,c[i]);
    return ans;
}
void modify(int x)
{
    for(int i=1;i<=5;i++)
    dp[pos[x][i]]=max(dp[pos[x][i]],(pos[x][i]>1?query(pos[x][i]-1):0)+1);
    for(int i=1;i<=5;i++)
    add(pos[x][i],dp[pos[x][i]]);
}
int main()
{
    scanf("%d",&m);
    n=m*5;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        cnt[a[i]]++; pos[a[i]][cnt[a[i]]]=i;
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&b[i]);
        modify(b[i]);
    }
    printf("%d\n",*max_element(dp+1,dp+n+1));
    return 0;
}