BZOJ 3110 [ZJOI2013]K大数查询

题目在这里呀~

这题被卡常了qaq（ZJOI临近了我也不想在这种题上花太多时间…）
可以想到要用二分答案（只是以前做过一道类似的题啦
然后常规的，求出左子树的贡献，如果大于c，就往左子树找，否则往右子树找。
然后就是树套树了？
外层权值内层记区间和。
我没看懂他们说的标记永久化什么的，可听说是没负数的而且不会爆int？？
然后就调了一个晚上，后来看评论发现…天呐要开longlong qwq！

于是改成longlong？
然后洛谷T一个点？可能常数写的太大了吧哎呀不管了啦！

//Suplex
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 20000000+1000
#define ll long long
using namespace std;
int n,m,opt,a,b,c,cnt,p,l,r,root[N];
ll tag[N],sum[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    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 segment{
    int l,r;
}t[N];

inline void pushdown(int rt,int l,int r)
{
    if(tag[rt] && l<r){
        if(!t[rt].l) t[rt].l=++cnt;
        if(!t[rt].r) t[rt].r=++cnt;
        int mid=(l+r)>>1;
        tag[t[rt].l]+=tag[rt];
        tag[t[rt].r]+=tag[rt];
        sum[t[rt].l]+=tag[rt]*(mid-l+1);
        sum[t[rt].r]+=tag[rt]*(r-mid);
    }
    tag[rt]=0;
}

inline void modify(int &rt,int l,int r,int x,int y)
{
    if(!rt) rt=++cnt;
    pushdown(rt,l,r);
    if(x<=l && r<=y){
        sum[rt]+=(r-l+1);
        tag[rt]++;
        return;
    }
    int mid=(l+r)>>1;
    if(y<=mid) modify(t[rt].l,l,mid,x,y);
    else if(x>mid) modify(t[rt].r,mid+1,r,x,y);
    else modify(t[rt].l,l,mid,x,mid),modify(t[rt].r,mid+1,r,mid+1,y);
    sum[rt]=sum[t[rt].l]+sum[t[rt].r];
}

inline ll query(int rt,int l,int r,int x,int y)
{
    if(!rt) return 0;
    pushdown(rt,l,r);
    if(x<=l && r<=y) return sum[rt];
    int mid=(l+r)>>1;
    if(y<=mid) return query(t[rt].l,l,mid,x,y);
    else if(x>mid) return query(t[rt].r,mid+1,r,x,y);
    else return query(t[rt].l,l,mid,x,mid)+query(t[rt].r,mid+1,r,mid+1,y);
}

void insert()
{
    p=1;l=1;r=n;
    while(l<r){
        int mid=(l+r)>>1;
        modify(root[p],1,n,a,b);
        if(c<=mid) r=mid,p=p+p;
        else l=mid+1,p=p+p+1;
    }
    modify(root[p],1,n,a,b);
}

int solve()
{
    p=1;l=1;r=n;
    while(l<r){
        int mid=(l+r)>>1;
        ll lsum=query(root[p+p],1,n,a,b);
        if(lsum>=c) r=mid,p=p+p;
        else l=mid+1,p=p+p+1,c-=lsum;
    }
    return r;
}

int main()
{
    n=read();m=read();
    while(m--){
        opt=read();a=read();b=read();c=read();
        if(opt==1){
            c=n-c+1;
            insert();
        }else printf("%d\n",n-solve()+1);
    }
    return 0;
}