zoukankan      html  css  js  c++  java
  • 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之间的最大值

    Input

    第一行包含两个整数 n,m 接下来一行包含 n 个整数, 表示a序列,接下来 m 行, 每行描述了一个操作.

    2<=n<=2e5 2<=q<=2e5,0<=ai<=2^20.

    Output

    对于每个第三类询问, 输出一个数字.

    Sample Input

    3 5
    9 19 0
    3 2 3
    2 3 3 18
    1 2 2 10
    3 1 2
    1 1 3 11

    Sample Output

    19
    9


    看到题目的时候相当僵硬,然后YY了一个算法然后僵硬了几个小时最后GG,我自己的错误算法还是不在这里说了。。说多了都是泪
    重大更新,我的代码终于调出来了,比正解更好理解!!!在正解后给予解释!!!!

    正解如下:
    首先我们可以发现,与和或的操作一个是有0变成0,一个是有1变成1
    那么如果一个数与、或上另一个数vl,只跟另一个数的二进制位上为0、1的位有关

    我们考虑在什么情境之下我们可以通过更新的vl值获得我们需要的最大值maxn,既然是区间操作,我们不难想到用线段树进行求解,但是单单是vl,我们并不能得到想要的maxn,所以我们需要维护其他的辅助变量

    我们发现&和|操作所关联的二进制位只有vl上的0或者1,所以我们可以考虑定义线段树上区间的sam,如果一个区间的所有数在第i个二进制位上的数码相等,sam的这个二进制位上的数为1,否则为0

    通过更新sam,我们是可以很方便的计算和更新maxn的

    但是当一个区间存在多个不同的sam怎么办?我们用vl值进行更新的时候依然不能很方便的进行计算,所以我们定义check函数来限制线段树修改的边界问题,显然,当满足所有(&vl且vl上为0的位sam上为1)或者(|vl且vl上为1的位sam上为1)我们只需要用vl&、|上当前区间最大值就好,不满足就向下递归问题

    现在我们考虑向上维护sam值,显然,对于一个二进制位,只有当左右两区间的sam在这个二进制上的值都为1且左右两区间的数在这两个二进制位上相等才行,可以表示成sam[t]=(sam[LD] & sam[RD]) & (INF ^ maxn[LD] ^ maxn[RD]),这里INF定义为二进制上的极大值((1<<20)-1),然后maxn直接左右区间选择max就好了

    那么如何向下更新呢?正解的思路真的比较神奇
    我们把|操作强行通过某种等价的方式转化成&操作,然后只用一个修改函数进行修改,这个有兴趣的照着代码枚举二进制情境验证一下吧。
    其次,还有一个比较玄学的是正解没有加lazy标记,直接用父亲节点t的maxn和sam值对儿子节点s的maxn和sam值进行更新,这个稍微讲一下,因为所有在sam[t]上出现的1一定会在sam[s]上出现,所以sam[s]|=sam[t]就可以维护,然后对于maxn,我们先把maxn[s]上左右和sam[t]有关的二进制位全部变成零,然后再或上sam[t]和maxn[t]均为1的二进制就好了


    #include<bits/stdc++.h>
    using namespace std;
    #define N 200010
    #define INF ((1<<20)-1)
    #define LD (t<<1)
    #define RD ((t<<1)|1)
    struct Segment_Tree{
    	//1->& 2->|
    	int num[N],l[N<<2],r[N<<2];
    	int maxn[N<<2],sam[N<<2];
    	void pushup(int t){
    		sam[t]=(sam[LD]&sam[RD])&(INF^(maxn[LD]^maxn[RD]));
    		maxn[t]=max(maxn[LD],maxn[RD]);
    	}
    	void pushnow(int t,int v1,int v2){
    		sam[t]|=v1;
    		maxn[t]=(maxn[t]&(INF^v1))|(v1&v2);
    	}
    	void pushdown(int t){
    		pushnow(LD,sam[t],maxn[t]);
    		pushnow(RD,sam[t],maxn[t]);
    	}
    	void build(int t,int ll,int rr){
    		if(ll>rr)return;
    		l[t]=ll;r[t]=rr;
    		if(ll==rr){
    			maxn[t]=num[ll];
    			sam[t]=INF;
    			return;
    		}
    		int mid=(ll+rr)>>1;
    		build(LD,ll,mid);
    		build(RD,mid+1,rr);
    		pushup(t);
    	}
    	void modify(int t,int ll,int rr,int v1,int v2){
    		if(l[t]>rr||r[t]<ll)return;
    		if(ll<=l[t]&&r[t]<=rr){
    			if(l[t]==r[t]||((v1&sam[t])==v1)){
    				maxn[t]=(maxn[t]&(INF^v1))|(v1&v2);
    				return;
    			}
    		}
    		pushdown(t);
    		modify(LD,ll,rr,v1,v2);
    		modify(RD,ll,rr,v1,v2);
    		pushup(t);
    	}
    	int query(int t,int ll,int rr){
    		if(l[t]>rr||r[t]<ll)return 0;
    		if(ll<=l[t]&&r[t]<=rr)return maxn[t];
    		pushdown(t);
    		return max(query(LD,ll,rr),query(RD,ll,rr));		
    	}
    }tree;
    int main(){
    	int n,m;scanf("%d%d",&n,&m);
    	for(int i=1;i<=n;i++)scanf("%d",&tree.num[i]);
    	tree.build(1,1,n);
    	for(int i=1;i<=m;i++){
    		int op,l,r,vl;
    		scanf("%d%d%d",&op,&l,&r);
    		if(op==1){
    			scanf("%d",&vl);
    			tree.modify(1,l,r,vl^INF,0);
    		}else if(op==2){
    			scanf("%d",&vl);
    			tree.modify(1,l,r,vl,vl);
    		}else printf("%d
    ",tree.query(1,l,r));
    	}
    	return 0;
    }
    

    更新:
    身为一个不顾一切强行刚的OIER,我肯定是会刚题刚到凌晨的,然后我就把它刚出来了,我的做法将&和|的操作分别进行,在modify过程中肯定是很好理解的,现在解释一下对&和|的边界设定:
    定义check1检查对&的modify:{
    我们发现当vl涉及的数码为0的所有二进制位在取件sam中都为1的时候,那么vl对于区间中的所有二进制的这几位的影响相同,不会改变大小顺序,所以可以直接进行maxn的更新,条件:((INF ^ vl)&sam[t]) == (INF^vl)
    }
    定义check2检查对|的modify:{
    我们发现当vl涉及的数码为1的所有二进制位在取件sam中都为1的时候,那么vl对于区间中的所有二进制的这几位的影响相同,不会改变大小顺序,所以可以直接进行maxn的更新,条件:(vl&sam[t]) == vl
    }
    然后注意这个时候需要在pushdown的时候进行边界条件判定,不然会出现奇奇怪怪的错误

    附上代码:

    #include<bits/stdc++.h>
    using namespace std;
    #define INF ((1<<20)-1)
    #define N 200010
    #define LD (t<<1)
    #define RD ((t<<1)|1)
    struct Segment_Tree{
    	//1->& 2->|
    	int num[N],sam[N<<2],maxn[N<<2],l[N<<2],r[N<<2];
    	void pushup(int t){
    		sam[t]=(sam[LD]&sam[RD])&(INF^(maxn[LD]^maxn[RD]));
    		maxn[t]=max(maxn[LD],maxn[RD]);
    	}
    	void pushnow(int t,int v1,int v2){
    		sam[t]|=v1; 
    		maxn[t]=(maxn[t]&(INF^v1))|(v1&v2);
    	}
    	void pushdown(int t){
    		if(l[t]==r[t])return;
    		pushnow(LD,sam[t],maxn[t]);
    		pushnow(RD,sam[t],maxn[t]);
    	}
    	void build(int t,int ll,int rr){
    		if(ll>rr)return;
    		l[t]=ll;r[t]=rr;
    		if(ll==rr){sam[t]=INF;maxn[t]=num[ll];return;}
    		int mid=(ll+rr)>>1;
    		build(LD,ll,mid);
    		build(RD,mid+1,rr);
    		pushup(t);
    	}
    	bool check1(int t,int vl){return ((INF^vl)&sam[t])==(INF^vl);}
    	bool check2(int t,int vl){return (vl&sam[t])==vl;}
    	void modify1(int t,int ql,int qr,int vl){
    		if(r[t]<ql||qr<l[t])return;
    		pushdown(t);
    		if(ql<=l[t]&&r[t]<=qr)
    			if(l[t]==r[t]||check1(t,vl)){maxn[t]&=vl;return;}
    		modify1(LD,ql,qr,vl);
    		modify1(RD,ql,qr,vl);
    		pushup(t);
    	}
    	void modify2(int t,int ql,int qr,int vl){
    		if(r[t]<ql||qr<l[t])return;
    		pushdown(t);
    		if(ql<=l[t]&&r[t]<=qr)
    			if(l[t]==r[t]||check2(t,vl)){maxn[t]|=vl;return;}
    		modify2(LD,ql,qr,vl);
    		modify2(RD,ql,qr,vl);
    		pushup(t);
    	}
    	int query(int t,int ql,int qr){
    		if(r[t]<ql||qr<l[t])return 0;
    		if(ql<=l[t]&&r[t]<=qr)return maxn[t];
    		pushdown(t);
    		return max(query(LD,ql,qr),query(RD,ql,qr));
    	}
    }tree;
    int main(){
    	int n,m;scanf("%d%d",&n,&m);
    	for(int i=1;i<=n;i++)scanf("%d",&tree.num[i]);
    	tree.build(1,1,n);
    	for(int i=1;i<=m;i++){
    		int op,l,r,vl;
    		scanf("%d%d%d",&op,&l,&r);
    		if(op==1){
    			scanf("%d",&vl);
    			tree.modify1(1,l,r,vl);
    		}else if(op==2){
    			scanf("%d",&vl);
    			tree.modify2(1,l,r,vl);
    		}else printf("%d
    ",tree.query(1,l,r));
    	}
    	return 0;
    }
    
  • 相关阅读:
    ngx-push-stream模块源码学习(四)——订阅
    nginx-push-stream模块源码学习(三)——发布
    nginx-push-stream模块源码学习(二)——模块初始化
    ngx-push-stream模块源码学习(一)——序言
    微信通过网页版定时推送消息脚本
    移动硬盘 或硬盘显示需先格式化问题解决
    SQL优化方法:
    sqlserver 电脑重启以后服务突然无法启动 报错
    sql server xp_cmdshell 过程中报错原因
    .net Core 2.0应用程序发布到IIS上注意事项
  • 原文地址:https://www.cnblogs.com/dream-maker-yk/p/9858132.html
Copyright © 2011-2022 走看看