BZOJ4942 & UOJ314：[NOI2017]整数——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=4942

http://uoj.ac/problem/314

https://www.luogu.org/problemnew/show/P3822

题面是markdown形式的所以我传不上……

UPD：18.5.11改成对参考代码的理解失误。

参考：http://www.cnblogs.com/RabbitHu/p/UOJ314.html仔细思考之后发现lazy标记可以不下传，因为区间修改都是改0/INF，单点修改就直接改就行了……

（晚上睡前突然发现这个sb问题然后太晚了只好现在改。）

暴力修改显然是不行的。

考虑让你改的是a*2^b，已经明示了想让你直接修改二进制位，于是把a用O(log)时间拆成二进制然后log次加减，期间的进位退位问题用参考博客的方法可以O(log)实现。

但问题是位数一共有3e7……显然会T的。

于是压60位，这样时间复杂度就有保证了。

（我记得以前有人批评松松松就是因为他WC出神题然后NOI出签到题……woc我怕不是送分题都做不出来我退役吧。）

#include<cstdio>
#include<queue>
#include<cctype>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=5e5+;
const int B=;
const ll INF=(1LL<<B)-;
inline ll read(){
    ll X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
ll lazy[N*+],val[N*+];
int pos[N*+];
bool tag[N*+][];
inline void upt(int a){
    tag[a][]=tag[a<<][]&tag[a<<|][];
    tag[a][]=tag[a<<][]&tag[a<<|][];
}
inline void mdy(int a,ll x){
    if(pos[a]!=-)val[pos[a]]=x;
    if(x==)tag[a][]=,tag[a][]=;
    else if(x==INF)tag[a][]=,tag[a][]=;
    else tag[a][]=tag[a][]=;
    lazy[a]=x;
}
inline void push(int a){
    if(lazy[a]==-)return;
    mdy(a<<,lazy[a]);
    mdy(a<<|,lazy[a]);
    lazy[a]=-;
}
void build(int a,int l,int r){
    tag[a][]=;lazy[a]=-;pos[a]=-;
    if(l==r){
        pos[a]=l;
        return;
    }
    int mid=(l+r)>>;
    build(a<<,l,mid);build(a<<|,mid+,r);
}
void modify(int a,int l,int r,int l1,int r1,ll x){
    if(r<l1||r1<l||l1>r1)return;
    if(l1<=l&&r<=r1){
        mdy(a,x);
        return;
    }
    push(a);
    int mid=(l+r)>>;
    modify(a<<,l,mid,l1,r1,x);
    modify(a<<|,mid+,r,l1,r1,x);
    upt(a);
}
ll query(int a,int l,int r,int p){
    if(l==r)return val[l];
    push(a);
    int mid=(l+r)>>;
    if(p<=mid)return query(a<<,l,mid,p);
    else return query(a<<|,mid+,r,p);
}
int find(int a,int l,int r,int p,int o){
    if(tag[a][!o])return -;
    if(l==r)return l;
    push(a);
    int mid=(l+r)>>,tmp;
    if(p<=mid&&(tmp=find(a<<,l,mid,p,o))!=-)return tmp;
    return find(a<<|,mid+,r,p,o);
}
void add(int p,ll x){
    ll tmp=query(,,N,p);
    modify(,,N,p,p,(tmp+x)&INF);
    if(tmp+x>INF){
        int r=find(,,N,p+,);
        modify(,,N,r,r,val[r]+);
        modify(,,N,p+,r-,);
    }
}
void del(int p,ll x){
    ll tmp=query(,,N,p);
    modify(,,N,p,p,(tmp-x)&INF);
    if(tmp-x<){
        int r=find(,,N,p+,);
        modify(,,N,r,r,val[r]-);
        modify(,,N,p+,r-,INF);
    }
}
int main(){
    int n=read();read(),read(),read();
    build(,,N);
    for(int i=;i<=n;i++){
        int op=read();
        if(op==){
            ll a=read(),b=read();
            ll p=b/B,m=b%B;
            if(a>){
                ll x=(ll)a<<m&INF;
                if(x)add(p,x);
                p++;a>>=(B-m);
                if(a)add(p,a);
            }else{
                a=-a;
                ll x=(ll)a<<m&INF;
                if(x)del(p,x);
                p++;a>>=(B-m);
                if(a)del(p,a);
            }
        }else{
            ll k=read();
            printf("%lld\n",query(,,N,k/B)>>(k%B)&);
        }
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/ +

+++++++++++++++++++++++++++++++++++++++++++