BZOJ 4942 NOI2017 整数 (压位+线段树)

题目大意：让你维护一个数x(x位数<=3*1e7)，要支持加/减a*2^b，以及查询x的第i位在二进制下是0还是1

作为一道noi的题，非常考验写代码综合能力，敲+调+借鉴神犇的代码 3个多小时才过...

思路并不难，题目里b<=30n暗示压位，每次压30位可过

先分析一下加法，加a*2^b相当于在第b-1位加a，如果进位了(即>)，那就在下一位+1，如果下一位还进位了，那就再下一位+1......

暴力进位显然不可取，那么用线段树维护它，维护连续的几大位之内是否全都是1

a<=1e9<=2^30，所以它最多被拆成两大位压进线段树里

那么每一个大位的加法操作可以变成 (方便描述，下面的数字是从左往右读)：

1.当前位+=val，拆分出的val<=2^30

2.如果进位了，就去线段树里找接下来连续是1的最大的位置，然后在它下一位+1 比如01011 11111 10110，就是在第三位+1(实际操作可以直接找它下一位)

3.把这两位之间(不包括两端)的所有的大位都改成，区间修改打标记

减法操作和加法非常类似，借位就是找最右边是

查询操作非常简单没什么好说的，别忘了先下推标记再查询...

最先我的代码写了1个多小时，维护的东西一大堆，感觉恶心得一匹而且巨难调

实在是没想到好的优化方法，就借鉴了神犇gxzlegend的一些神级操作(orzorz)：

a组(any)维护它子节点的&值，e数组(exist)维护它子节点的|值，这个维护方法很神，比我原来想的简单得多，而且不需要乱七八糟的打标记

find_inf和find_zero这两个函数也非常简洁明了

以及感谢LOJ的数据顺便吐槽一句辣鸡bzoj

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define ll long long
 #define inf (1<<30)-1
 #define il inline
 #define ls rt<<1
 #define rs rt<<1|1
 #define N 1000100
 using namespace std;
 //re
 int n,ma,t1,t2,t3;
 int a[N<<],e[N<<],tag[N<<];

 int gc()
 {
     int rett=,fh=;char p=getchar();
     while(p<''||p>''){if(p=='-')fh=-;p=getchar();}
     while(p>=''&&p<=''){rett=(rett<<)+(rett<<)+p-'';p=getchar();}
     return rett*fh;
 }
 il void pushup(int rt){
     e[rt]=e[rt<<]|e[rt<<|];
     a[rt]=a[rt<<]&a[rt<<|];
 }
 il void pushdown(int rt){
     if(tag[rt]==-) a[ls]=a[rs]=e[ls]=e[rs]=,tag[rt]=,tag[ls]=tag[rs]=-;
     if(tag[rt]==) a[ls]=a[rs]=e[ls]=e[rs]=inf,tag[rt]=,tag[ls]=tag[rs]=;
 }
 int find_inf(int x,int l,int r,int rt)
 {
     if(a[rt]==inf) return -;
     if(l==r) return l;
     pushdown(rt);
     int mid=(l+r)>>;
     if(x<=mid){
         int pos=find_inf(x,l,mid,rt<<);
         if(pos==-) return find_inf(x,mid+,r,rt<<|);
         else return pos;
     }
     return find_inf(x,mid+,r,rt<<|);
 }
 int find_zero(int x,int l,int r,int rt)
 {
     if(!e[rt]) return -;
     if(l==r) return l;
     pushdown(rt);
     int mid=(l+r)>>;
     if(x<=mid){
         int pos=find_zero(x,l,mid,rt<<);
         if(pos==-) return find_zero(x,mid+,r,rt<<|);
         else return pos;
     }
     return find_zero(x,mid+,r,rt<<|);
 }
 int upd1(int x,int l,int r,int rt,int p)
 {
     if(l==r)
     {
         a[rt]+=p,e[rt]+=p;
         if(a[rt]>inf){a[rt]-=(inf+);e[rt]-=(inf+);return ;}
         if(a[rt]<){a[rt]+=(inf+);e[rt]+=(inf+);return -;}
         return ;
     }
     pushdown(rt);
     int mid=(l+r)>>,ans=;
     if(x<=mid) ans=upd1(x,l,mid,rt<<,p);
     else ans=upd1(x,mid+,r,rt<<|,p);
     pushup(rt);return ans;
 }
 void upd2(int L,int R,int l,int r,int rt,int val)
 {
     if(L<=l&&r<=R){
         tag[rt]=(val)?:-;
         a[rt]=e[rt]=(val)?inf:;
         return;}
     pushdown(rt);
     int mid=(l+r)>>;
     if(L<=mid) upd2(L,R,l,mid,rt<<,val);
     if(R>mid) upd2(L,R,mid+,r,rt<<|,val);
     pushup(rt);
 }
 void add(int x,int y)
 {
     int p=upd1(x,,n,,y);
     if(!p) return;
     int pos=find_inf(x+,,n,);
     if(pos->=x+) upd2(x+,pos-,,n,,);
     upd1(pos,,n,,);
 }
 void del(int x,int y)
 {
     int p=upd1(x,,n,,-y);
     if(!p) return;
     int pos=find_zero(x+,,n,);
     if(pos->=x+) upd2(x+,pos-,,n,,);
     upd1(pos,,n,,-);
 }
 int query(int x,int l,int r,int rt,int p)
 {
     if(l==r) return ((<<p)&a[rt])?:;
     int mid=(l+r)>>;
     pushdown(rt);
     if(x<=mid) return query(x,l,mid,rt<<,p);
     else  return query(x,mid+,r,rt<<|,p);
 }
 int main()
 {
     scanf("%d%d%d%d",&n,&t1,&t2,&t3);
     int fl,x,y;
     for(int i=;i<=n;i++)
     {
         fl=gc();
         if(fl==){
             x=gc(),y=gc();
             if(x>=){
                 add(y/,(x&((<<(-y%))-))<<(y%));
                 add(y/+,x>>(-y%));
             }else{ x=-x;
                 del(y/,(x&((<<(-y%))-))<<(y%));
                 del(y/+,x>>(-y%));
             }
         }else{
             x=gc();
             printf("%d\n",query(x/,,n,,x%));
         }
     }
     return ;
 }