BZOJ 3110: [Zjoi2013]K大数查询

网上面貌似挺多人说cdq分治？（虽说和整体二分打起来是挺像的。）
可是这题就是整体二分吧？ 当然还有强大树套树 像我这种垃圾就不会打。
过程中会比int大 用unsigned int 或者 long long都可以

#include<bits/stdc++.h>
using namespace std;
const int N=50001;
inline int read()
{
    char ch=getchar(); int x=0,f=1;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
struct node{int l,r,c,s,id;}a[N],b[N],br[N];
int n,m,ans[N],lazy[N<<2]; unsigned s[N<<2];
bool v[N];
inline void pushdown(int x,int l,int r)
{
    int mid=(l+r)>>1,lc=x<<1,rc=lc|1,u=lazy[x];
    lazy[lc]+=u,lazy[rc]+=u,lazy[x]=0;
    s[lc]+=(mid-l+1)*u;
    s[rc]+=(r-mid)*u;
}
inline void add(int x,int l,int r,int ql,int qr,int u)
{
    if(ql==l && qr==r){s[x]+=u*(r-l+1),lazy[x]+=u; return;}
    pushdown(x,l,r);
    int mid=(l+r)>>1,lc=x<<1,rc=lc|1;
    if(qr<=mid)add(lc,l,mid,ql,qr,u);
    else if(ql>mid)add(rc,mid+1,r,ql,qr,u);
    else add(lc,l,mid,ql,mid,u),add(rc,mid+1,r,mid+1,qr,u);
    s[x]=s[lc]+s[rc];
}
inline unsigned query(int x,int l,int r,int ql,int qr)
{
    if(ql==l && qr==r)return s[x];
    pushdown(x,l,r);
    int mid=(l+r)>>1,lc=x<<1,rc=lc|1;
    if(qr<=mid) return query(lc,l,mid,ql,qr);
    else if(ql>mid)return query(rc,mid+1,r,ql,qr);
    return query(lc,l,mid,ql,mid)+query(rc,mid+1,r,mid+1,qr);
}
inline void solve(int l,int r,int lc,int rc)
{
    if(l>r)return;
    if(lc==rc)
    {
        for(int i=l;i<=r;i++) if(a[i].s) ans[a[i].id]=lc;
        return;
    }
    int mid=(lc+rc)>>1,bn=0,brn=0,i;
    for(i=l;i<=r;i++)
    {
        if(!a[i].s)
        {
            if(a[i].c>mid){br[++brn]=a[i]; add(1,1,n,a[i].l,a[i].r,1);}
            else b[++bn]=a[i];
        }
        else
        {
            unsigned u=query(1,1,n,a[i].l,a[i].r);
            if(u>=a[i].c)br[++brn]=a[i];
            else a[i].c-=u,b[++bn]=a[i];
        }
    }
    for(i=l;i<=r;i++)
        if(!a[i].s && a[i].c>mid)add(1,1,n,a[i].l,a[i].r,-1);
    for(i=1;i<=bn;i++)a[l+i-1]=b[i];
    for(i=1;i<=brn;i++)a[l+bn+i-1]=br[i];
    solve(l,l+bn-1,lc,mid);
    solve(l+bn,r,mid+1,rc);
}
int main()
{
    int i,x,l,r,c;
    n=read(),m=read();
    for(i=1;i<=m;i++)
    {
        x=read(),l=read(),r=read(),c=read();
        a[i]=(node){l,r,c,v[i]=x-1,i};
    }
    solve(1,m,-N,N);
    for(i=1;i<=m;i++)if(v[i])printf("%d\n",ans[i]);
    return 0;
}