3031 序列

2 篇文章 0 订阅

Task
序列有N个数,m种变化(a,x) 表示A[a]可以变成x,求最长子序列满足任意1个变化或都不变化,子序列都不降。
N<=1e5+3

Solution
子序列是由一个一个数构成的,且相邻两数之间满足
Mi[i]>=v[j],v[i]>=mx[i],i>j

由“最长子序列”和“不降”,我们想到了LIS,普通的LIS转移:
A[i]>A[j],i>j
只有2个限制条件,题目中有三个限制条件,是二维LIS。

多维度降维,复杂的问题化为简单。常见的消除维度的方法:排序,分治。
计算当前区间[L,R]的dp值,分成两个区间[L,mid],[mid+1,R],用左边的区间更新右边区间的dp值。除了左边的区间可以更新右边,右边也可以自己更新。
这里写图片描述
从左到右转移,下标关系一定符合,不用再维护。左边按照mx排序,右边按照v从小到大排序,可以归并转移,又消去了一维,再用树状数组维护第三维,权值BIT维护前缀信息。

const int M=1e5+3;
int n,m,an,bn,tot;
int val[M],mi[M],mx[M],A[M],B[M],dp[M];
struct Tree{//权值bit 
    int bit[M];//前缀信息 
    inline void add(int x,int v){
        while(x<=tot){
            bit[x]=max(bit[x],v);
            x+=lowbit(x);
        }
    }
    inline void clear(int x){
        while(x<=tot){
            bit[x]=0;
            x+=lowbit(x);
        }
    }
    inline int query(int x){
        int ans=0;
        while(x>0){
            ans=max(ans,bit[x]);
            x-=lowbit(x);
        }
        return ans;
    }
    inline bool check(){
        rep(i,1,M-1)
            if(bit[i]>0)return 1;
        return 0;
    }
}T;
inline bool cmpA(int x,int y){return mx[x]<mx[y];}
inline bool cmpB(int x,int y){return val[x]<val[y];}
inline void input(){
    int a,b;
    rd(n);rd(m);
    rep(i,1,n){
        rd(val[i]);
        mx[i]=mi[i]=val[i];
        tot=max(tot,val[i]);
    }
    rep(i,1,m){
        rd(a);rd(b);
        mx[a]=max(mx[a],b);
        mi[a]=min(mi[a],b);
        tot=max(tot,max(mx[a],mi[a]));
    }
}
inline void dfs(int l,int r){
    if(l==r){dp[l]=max(dp[l],1);return;}
    int mid=l+r>>1;
    dfs(l,mid);
    an=bn=0;
    rep(i,l,mid)A[++an]=i;
    rep(i,mid+1,r)B[++bn]=i;
    sort(A+1,A+1+an,cmpA);
    sort(B+1,B+1+bn,cmpB);

    int k=1;
    rep(i,1,bn){
        while(k<=an&&mx[A[k]]<=val[B[i]]){
            T.add(val[A[k]],dp[A[k]]);
            ++k;
        }
        dp[B[i]]=max(dp[B[i]],T.query(mi[B[i]])+1);
    }

    rep(i,1,an)
        T.clear(val[A[i]]);
    dfs(mid+1,r);
}
inline void print(){
    int ans=0;
    rep(i,1,n)ans=max(ans,dp[i]);
    sc(ans);
}
int main(){ 
    input();
    dfs(1,n);
    print();
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值