bzoj5312 冒险 线段树

Description

---

Kaiser终于成为冒险协会的一员，这次冒险协会派他去冒险，他来到一处古墓，却被大门上的守护神挡住了去路，守护神给出了一个问题，
只有答对了问题才能进入，守护神给出了一个自然数序列a,每次有一下三种操作。
1，给出l,r,x，将序列l,r之间的所有数都 and x
2，给出l,r,x，将序列l,r之间的所有数都 or x
3，给出l,r，询问l,r之间的最大值

第一行包含两个整数 n,m 接下来一行包含 n 个整数, 表示a序列，接下来 m 行, 每行描述了一个操作.
2<=n<=2e5 2<=q<=2e5,0<=ai<=2^20.

Solution

---

注意到普通线段树的lazy标记不管用了，因为一段区间中同时and或or上一个数字后最大值并不还是原来的最大值
我们考虑仍然用线段树，将一段连续的、and上一个数字后最大值不改变的区间成为同类区间，显然我们需要把tag打在这些同类区间上才是科学的，那么记录一下每个区间的and和，or和，以及两个操作的tag即可。注意tag下传时and和or的优先级，具体可以类比乘法和加法同时存在的tag
复杂度需要用到势能分析，实际操作可以检验这是不会T的（滑稽

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)

const int N=200005;

int max[N<<2],sum[N<<2][2],tag[N<<2][2];
int rev[1148576],a[N];

int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}

void mark1(int now,int x) {
    tag[now][0]&=x; tag[now][1]&=x; max[now]&=x;
    sum[now][0]&=x; sum[now][1]&=x;
}

void mark2(int now,int x) {
    tag[now][1]|=x; max[now]|=x;
    sum[now][0]|=x; sum[now][1]|=x;
}

void push_down(int now) {
    if (rev[tag[now][0]]) {
        mark1(now<<1,tag[now][0]); mark1(now<<1|1,tag[now][0]);
        tag[now][0]=rev[0];
    }
    if (tag[now][1]) {
        mark2(now<<1,tag[now][1]),mark2(now<<1|1,tag[now][1]);
        tag[now][1]=0;
    }
}

void push_up(int now) {
    max[now]=std:: max(max[now<<1],max[now<<1|1]);
    sum[now][0]=sum[now<<1][0]&sum[now<<1|1][0];
    sum[now][1]=sum[now<<1][1]|sum[now<<1|1][1];
}

void modify1(int now,int tl,int tr,int l,int r,int v) {
    if (l<=tl&&tr<=r&&(sum[now][0]&rev[v])==(sum[now][1]&rev[v])) {
        mark1(now,v);
        return ;
    }
    push_down(now);
    int mid=(tl+tr)>>1;
    if (l<=mid) modify1(now<<1,tl,mid,l,r,v);
    if (r>mid) modify1(now<<1|1,mid+1,tr,l,r,v);
    push_up(now);
}

void modify2(int now,int tl,int tr,int l,int r,int v) {
    if (l<=tl&&tr<=r&&(sum[now][0]&v)==(sum[now][1]&v)) {
        mark2(now,v);
        return ;
    }
    push_down(now);
    int mid=(tl+tr)>>1;
    if (l<=mid) modify2(now<<1,tl,mid,l,r,v);
    if (r>mid) modify2(now<<1|1,mid+1,tr,l,r,v);
    push_up(now);
}

int query(int now,int tl,int tr,int l,int r) {
    if (l<=tl&&tr<=r) return max[now];
    push_down(now);
    int mid=(tl+tr)>>1,qx=0,qy=0;
    if (l<=mid) qx=query(now<<1,tl,mid,l,r);
    if (r>mid) qy=query(now<<1|1,mid+1,tr,l,r);
    return std:: max(qx,qy);
}

void build_tree(int now,int tl,int tr) {
    tag[now][0]=rev[0];
    if (tl==tr) {
        sum[now][0]=sum[now][1]=max[now]=a[tl];
        return ;
    }
    int mid=(tl+tr)>>1;
    build_tree(now<<1,tl,mid);
    build_tree(now<<1|1,mid+1,tr);
    push_up(now);
}

int main(void) {
    int n=read(),m=read();
    rep(i,0,(1<<20)-1) rev[i]=((1<<20)-1)^i;
    rep(i,1,n) a[i]=read();
    build_tree(1,1,n);
    for (;m--;) {
        int opt=read(),l=read(),r=read();
        if (opt==1) {
            int x=read();
            modify1(1,1,n,l,r,x);
        } else if (opt==2) {
            int x=read();
            modify2(1,1,n,l,r,x);
        } else printf("%d\n", query(1,1,n,l,r));
    }
    return 0;
}