1264: [AHOI2006]基因匹配Match

题目链接

题目大意：给定n个数和两个长度为n*5的序列，两个序列中的数均由1..n组成，且1..n中每个数恰好出现5次，求两个序列的LCS

题解：考虑LCS的性质
LCS的决策+1的条件是a[i]==b[j] 于是我们记录a序列中每个数的5个位置

$f[i]表示以a[i]为末尾的最长公共子串长度(避免讨论)$

依次处理b[i]，找到5个位置转移，$f[k]=max(f[k],query(k-1)+1)$

用树状数组维护$f[]$数组的前缀最值

PS：需要倒序枚举k，因为k升序排列，如果从小到大枚举k，后来的k可能会使用之前的k转移，从大到小就可以避免这一点

我的收获：观察条件，减少枚举量，逆序避免干扰同层转移

#include <cstdio>  
#include <cstring>  
#include <iostream>  
#include <algorithm>  
using namespace std;
#define M 200200*5

int n,ans;
int a[M],b[M],c[M],f[M],pos[M][6];  

void updata(int x,int v){for(;x<=n;x+=x&-x) c[x]=max(c[x],v);} 
int query(int x){int ret=0;for(;x>0;x-=x&-x) ret=max(ret,c[x]);return ret;} 

void work()
{
    for(int i=1;i<=n;i++)  
    {  
        for(int j=5;j>=1;j--)  
        {  
            int k=pos[b[i]][j];  
            f[k]=max(f[k],query(k-1)+1);  
            updata(k,f[k]);
            ans=max(ans,f[k]);  
        }  
    }  
    cout<<ans<<endl; 
}

void init()
{
    cin>>n;n*=5;
    for(int i=1;i<=n;i++) scanf("%d",&pos[a[i]][++pos[a[i]][0]]=i;//省空间吼啊
    for(int i=1;i<=n;i++) scanf("%d",&b[i]);
}

int main()  
{  
    init();
    work();
    return 0; 
}