4553: [Tjoi2016&Heoi2016]序列

来源：互联网发布：sm3算法是对称算法吗? 编辑：程序博客网时间：2024/06/05 22:31

挺一般的思路，转变为三维偏序，然后cdq分治维护dp值就可以了
想了10min左右，写了20min左右，调了60min左右。。。
都是sb错误（只有1.5k啊，非得生成数据对拍才能看出吗？）
谨以此文，纪念debug能力越来越弱的自己

#include<bits/stdc++.h>#define rep(i,k,n) for(int i=k;i<=n;i++)#define rep2(i,k,n) for(int i=k;i>=n;i--)#define ls x<<1#define rs x<<1|1using namespace std;const int N=4e5+7;int n,a[N],mx[N],mi[N],m,dp[N],b[N],M[N],len=0;bool cmp_a(int x,int y){    return a[x]<a[y];}bool cmp_mx(int x,int y){    return mx[x]<mx[y];}void ch(int x,int l,int r,int pos,int d){    if(l==r){if(!d)M[x]=0;else M[x]=max(M[x],d);return;}    int mid=(l+r)>>1;    if(pos<=mid)ch(ls,l,mid,pos,d);    else ch(rs,mid+1,r,pos,d);    M[x]=max(M[ls],M[rs]);}int Q(int x,int l,int r,int pos){    if(r<=pos)return M[x];    int mid=(l+r)>>1;    if(pos<=mid)return Q(ls,l,mid,pos);    else return max(M[ls],Q(rs,mid+1,r,pos));}void solve(int l,int r){    if(l>=r)return;    int mid=(l+r)>>1;    solve(l,mid);    rep(i,l,r)b[i]=i;    sort(b+l,b+mid+1,cmp_mx);    sort(b+mid+1,b+r+1,cmp_a);    int pos=mid+1;    rep(i,l,mid){        ch(1,1,len,a[b[i]],dp[b[i]]);        for(;pos<=r && a[b[pos]]<mx[b[i]];pos++); ////        for(;pos<=r && a[b[pos]]>=mx[b[i]] && (a[b[pos]]<mx[b[i+1]] || i==mid);pos++){            dp[b[pos]]=max(dp[b[pos]],Q(1,1,len,mi[b[pos]])+1);        }    }    rep(i,l,mid)ch(1,1,len,a[b[i]],0);    solve(mid+1,r);}int main(){    scanf("%d%d",&n,&m);    rep(i,1,n)scanf("%d",&a[i]),mx[i]=mi[i]=a[i],dp[i]=1,len=max(len,a[i]);    int x,y;    rep(i,1,m){        scanf("%d%d",&x,&y);        mx[x]=max(mx[x],y);        mi[x]=min(mi[x],y);        len=max(len,y);    }    solve(1,n);    int ans=0;    rep(i,1,n)ans=max(ans,dp[i]);    printf("%d\n",ans);}

0 0