4553: [Tjoi2016&Heoi2016]序列

挺一般的思路,转变为三维偏序,然后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|1
using 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);
}
阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/limboman/article/details/53982725
文章标签: 调试 cdq
相关热词: 4553大写
上一篇3925: [Zjoi2015]地震后的幻想乡
下一篇4521: [Cqoi2016]手机号码
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭