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;
}