洛谷 P4303 [AHOI2006]基因匹配（LCS+树状数组）

这题一看不是裸的lcs嘛。。
但看看数据范围，显然不支持传统$O(n^2)$算法。
注意到这样一个性质：每个数均只出现5次。

考虑这样的lcs做法：
设f[i][1~5]表示s2[1~i]与整个s1的lcs最大长度，且s2[i]必须在lcs中，且s2[i]在s1中第(1~5)次出现。
那么枚举s2[i]在s1中出现的位置pos[i][j]，f[i][j]=max{s2[k]在s1中位置位于pos[i][j]之前的f[k]}+1。
本题中每个i只用枚举5个pos，并且从后往前枚举以保证一个s2[i]不会对应多个pos；max值用树状数组维护一下就可以了。

复杂度$O(nlogn)$

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=100010;
int n,s1[maxn],s2[maxn],pos[maxn][5],cnt[maxn];
int arr[maxn];

inline int read(){
    int ret=0,sign=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') {if(ch=='-') sign=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();
    return ret*sign;
}

inline int lbt(int x){return x&-x;}
inline void add(int p,int x){while(p<=n*5) arr[p]=max(arr[p],x),p+=lbt(p);}
inline int query(int p){int ret=0; while(p>=1) ret=max(ret,arr[p]),p-=lbt(p); return ret;}

int main(){
    n=read();
    for(int i=1;i<=n*5;i++) s1[i]=read();
    for(int i=1;i<=n*5;i++) s2[i]=read();
    for(int i=1;i<=n*5;i++) pos[s1[i]][cnt[s1[i]]++]=i;
    for(int i=1;i<=n*5;i++)
        for(int j=4;j>=0;j--)
            add(pos[s2[i]][j],query(pos[s2[i]][j]-1)+1);
    printf("%d\n",query(n*5));
    return 0;
}