洛谷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[<<],*p1=buf,*p2=buf;
 template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,:;}
 int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 const int N=1e5+;
 int c[N],f[N][],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=;for(;x;x-=x&-x) cmax(res,c[x]);return res;
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     n=read()*;
     for(int i=,x;i<=n;++i)
     x=read(),f[x][cnt[x]++]=i;
     for(int i=;i<=n;++i){
         int x=read();
         for(int j=;j>=;--j){
             int at=f[x][j];
             add(at,query(at-)+);
         }
     }
     printf("%d\n",query(n));
     return ;
 }