[Tjoi2016&Heoi2016] 序列 CDQ分治

题意:

佳媛姐姐过生日的时候,她的小伙伴从某宝上买了一个有趣的玩具送给他。玩具上有一个数列,数列中某些项的值可能会变化,但同一个时刻最多只有一个值发生变化。

现在佳媛姐姐已经研究出了所有变化的可能性,她想请教你,能否选出一个子序列,使得在任意一种变化中,这个子序列都是不降的?请你告诉她这个子序列的最长长度即可。

注意:每种变化最多只有一个值发生变化。

在样例输入1中,所有的变化是

1 2 3
2 2 3
1 3 3
1 1 3
1 2 4

选择子序列为原序列,即在任意一种变化中均为不降子序列在样例输入2中,所有的变化是:

3 3 3

3 2 3

选择子序列为第一个元素和第三个元素,或者第二个元素和第三个元素,均可满足要求。

输入:

输入的第一行有两个正整数n, m,分别表示序列的长度和变化的个数。

接下来一行有n个数,表示这个数列原始的状态。

接下来m行,每行有2个数x, y,表示数列的第x项可以变化成y这个值。

1 ≤ x ≤ n

输出:

输出一个整数,表示对应的答案

思路:

看到最长不下降子序列首先想的是用树状数组优化\(dp\),但这道题相较于正常的\(LIS\)多了一个改变,实际上也就是多加了几个限制。

对于第\(i\)位和第\(j\)位如果选择\(j\)接在\(i\)后面,原来只用考虑\(i<j\)以及\(a_i<a_j\),这题要满足变化其中一个值后仍成立,所以需要考虑的有:
1、\(i<j\)
2、\(a_i<=Mi[j]\)
3、\(Mx[i]<=a_j\)
此时题目就变成了一个三维偏序问题,可以考虑用CDQ求解了。

即在外面维护条件1,CDQ中根据条件2排序,用树状数组直接维护条件3

但由于等号两边比较的东西不同,所以相较以前的写法,这里对于\([l,mid]\)和\([mid+1,r]\)的排序方式是不一样的。也就是对于左边区间直接用值进行升序排列,右区间根据Mi升序排列。

而树状数组中每次在Mx的位置上进行更新,查询时则根据val

应该没什么细节,只要不字母打错之类的

#include<bits/stdc++.h>
#define M 100005
#define lowbit(x) (x&-x)
using namespace std;
int n,m,ans[M];
struct node {
    int id,x,mx,mi;//最大可以变成的 和最小可以变成的值
    void add(int y) {
        mx=max(mx,y),mi=min(mi,y);
    }
} a[M],p[M];
struct Tree {
    int cnt[M];
    void add(int x,int y) {
        while(x<M)cnt[x]=max(cnt[x],y),x+=lowbit(x);
    }
    int Query(int x) {
        int res=0;
        while(x)res=max(res,cnt[x]),x-=lowbit(x);
        return res;
    }
    void del(int x) {
        while(x<M)cnt[x]=0,x+=lowbit(x);
    }
} T;
bool cmp_x(node A,node B) {
    if(A.x!=B.x)return A.x<B.x;
    return A.mx<B.mx;
}
bool cmp_mi(node A,node B) {
    if(A.mi!=B.mi)return A.mi<B.mi;
    return A.x<B.x;
}
bool cmp_id(node A,node B){
    return A.id<B.id;
} 
void CDQ(int l,int r) {
    while(l>=r)return;
    int mid=(l+r)>>1;
    CDQ(l,mid);
    sort(a+l,a+mid+1,cmp_x);
    sort(a+mid+1,a+r+1,cmp_mi);
    int j=l;
    for(int i=mid+1; i<=r; i++) {
        while(j<=mid&&a[j].x<=a[i].mi) {
            T.add(a[j].mx,ans[a[j].id]);
            j++;
        }
        int t=T.Query(a[i].x);
        ans[a[i].id]=max(ans[a[i].id],t+1);
    }
    for(int i=l; i<j; i++)T.del(a[i].mx);//撤销
//  sort(a+mid+1,a+r+1,cmp_id);
    for(int i=mid+1; i<=r; i++)a[i]=p[i];//还原成id升序排列
    CDQ(mid+1,r);
}
int main() {
//  freopen("data.in","r",stdin);
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++)scanf("%d",&a[i].x),a[i].mi=a[i].mx=a[i].x,a[i].id=i,ans[i]=1;//因为不管怎么样序列最短都有1
    for(int i=1; i<=m; i++) {
        int x,y;
        scanf("%d%d",&x,&y);
        a[x].add(y);
    }
    for(int i=1; i<=n; i++)p[i]=a[i];
    int Ans=0;
    CDQ(1,n);
    for(int i=1; i<=n; i++)Ans=max(Ans,ans[i]);
    printf("%d\n",Ans);
    return 0;
}
上一篇:Webview如何触发onReceivedLoginRequest;Webview怎么实现自动登录


下一篇:ubuntu下nodejs和npm的安装及升级