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;
}