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

传送门

我已经连这种傻逼题都不会了orz

正常的dp是$O(n^2)$的，枚举第一个数组的$j$，然后第二个数组的$k$，如果相等，则$dp[i]=dp[j]+1$，否则$dp[i]=dp[j]$

然后发现可以用树状数组优化这个过程……

不知道讲清楚没有因为我自己都还有点懵

//minamoto
#include<iostream>
#include<cstdio>
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
const int N=1e5+5;
int c[N],f[N][5],cnt[N],n;
inline void add(int x,int y){
    for(;x<=n;x+=x&-x) cmax(c[x],y);
}
int query(int x){
    int res=0;for(;x;x-=x&-x) cmax(res,c[x]);return res;
}
int main(){
//    freopen("testdata.in","r",stdin);
    n=read()*5;
    for(int i=1,x;i<=n;++i)
    x=read(),f[x][cnt[x]++]=i;
    for(int i=1;i<=n;++i){
        int x=read();
        for(int j=4;j>=0;--j){
            int at=f[x][j];
            add(at,query(at-1)+1);
        }
    }
    printf("%d
",query(n));
    return 0;
}