bzoj4990 [Usaco2017 Feb]Why Did the Cow Cross the Road II(dp+树状数组)

dp[j]表示第二行选择的最后一个位置为j时，最多选择多少对。以第一行做到i划分阶段，则这次只会影响a[i]能连到的的位置j，
dp[j]=max{dp[k]| 0<=k< j}+1.这样复杂度是O(n2*9)的，过不了，我们可以用树状数组或线段树维护一下决策的前缀最大值，优化成O(9*nlogn)。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define N 100010
#define inf 0x3f3f3f3f
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,a[N],dp[N],c[N],dx[]={-4,-3,-2,-1,0,1,2,3,4};
vector<int>to[N];
inline void update(int x,int val){
    for(;x<=n;x+=x&(-x)) c[x]=max(c[x],val);
}
inline int ask(int x){
    int res=0;
    for(;x;x-=x&(-x)) res=max(res,c[x]);return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();for(int i=1;i<=n;++i) a[i]=read();
    for(int j=1;j<=n;++j){
        int x=read();
        for(int k=0;k<9;++k){
            int xx=x+dx[k];if(xx<1||xx>n) continue;to[xx].push_back(j);
        }
    }for(int i=1;i<=n;++i){
        for(int j=to[a[i]].size()-1;j>=0;--j){
            dp[to[a[i]][j]]=ask(to[a[i]][j]-1)+1;update(to[a[i]][j],dp[to[a[i]][j]]);
        }
    }int ans=0;for(int i=1;i<=n;++i) ans=max(ans,dp[i]);
    printf("%d\n",ans);
    return 0;
}