[bzoj1264]基因匹配

首先朴素dp的方程，即$f[i][j]=max(f[i][j-1],f[i-1][j],(a[i]==b[j])*(f[i-1][j-1]+1))$，这中间特殊的转移只在a[i]=b[j]时，而在这道题中与a[i]相等的b[j]只有5个，考虑利用这一性质
先预处理每一种数的位置，在枚举i的同时，用线段树来维护f[i]这个数组，然后对于特殊的5个位置从上一轮转移，并对后面的区间max即可

#include<bits/stdc++.h>
using namespace std;
#define N 100005
#define L (k<<1)
#define R (L+1)
#define mid (l+r>>1)
int n,x,a[N],b[N][11],f[N<<2];
void down(int k){
    f[L]=max(f[L],f[k]);
    f[R]=max(f[R],f[k]);
    f[k]=0;
}
void update(int k,int l,int r,int x,int y,int z){
    if ((l>y)||(x>r))return;
    if ((x<=l)&&(r<=y)){
        f[k]=max(f[k],z);
        return;
    }
    down(k);
    update(L,l,mid,x,y,z);
    update(R,mid+1,r,x,y,z);
}
int query(int k,int l,int r,int x){
    if (!x)return 0;
    if (l==r)return f[k];
    down(k);
    if (x<=mid)return query(L,l,mid,x);
    return query(R,mid+1,r,x);
}
int main(){
    scanf("%d",&n);
    n*=5;
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=n;i++){
        scanf("%d",&x);
        b[x][++b[x][0]]=i;
    }
    for(int i=1;i<=n/5;i++)b[i][6]=n+1;
    for(int i=1;i<=n;i++)
        for(int j=5;j;j--)
            update(1,1,n,b[a[i]][j],b[a[i]][j+1]-1,query(1,1,n,b[a[i]][j]-1)+1);
    printf("%d",query(1,1,n,n));
}