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
,dp
,c
,dx[]={-4,-3,-2,-1,0,1,2,3,4};
vector<int>to
;
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;
}