K大数查询（bzoj 3110）

Description

有N个位置，M个操作。操作有两种，每次操作如果是1 a b c的形式表示在第a个位置到第b个位置，每个位置加入一个数c
如果是2 a b c形式，表示询问从第a个位置到第b个位置，第C大的数是多少。

Input

第一行N，M
接下来M行，每行形如1 a b c或2 a b c

Output

输出每个询问的结果

Sample Input

2 5
1 1 2 1
1 1 2 2
2 1 1 2
2 1 1 1
2 1 2 3

Sample Output

1
2
1

HINT

【样例说明】

第一个操作 后位置 1 的数只有 1 ， 位置 2 的数也只有 1 。 第二个操作 后位置 1

的数有 1 、 2 ，位置 2 的数也有 1 、 2 。 第三次询问 位置 1 到位置 1 第 2 大的数 是

1 。 第四次询问 位置 1 到位置 1 第 1 大的数是 2 。 第五次询问 位置 1 到位置 2 第 3

大的数是 1 。‍

N,M<=50000,N,M<=50000

a<=b<=N

1操作中abs(c)<=N

2操作中c<=Maxlongint

/*
  线段树套线段树的经典题目。
  这道题乍一看我觉得应该是外层区间线段树，内层权值线段树，但是不会搞。
  题解写的都是外层是权值线段树，在询问的时候采用类似于二分答案的方法。
  PS：题解貌似有点问题，有种被黄学长坑了的即时感 
*/
#include<cstdio>
#include<iostream>
#define N 50010
using namespace std;
int root[N*4],ls[N*400],rs[N*400],tag[N*400],sum[N*400],tot,n,m;
void push_up(int k){
    sum[k]=sum[ls[k]]+sum[rs[k]];
}
void push_down(int k,int l,int r){
    if(!tag[k]||l==r) return;
    if(!ls[k])ls[k]=++tot;
    if(!rs[k])rs[k]=++tot;
    tag[ls[k]]+=tag[k];
    tag[rs[k]]+=tag[k];
    int mid=(l+r)>>1;
    sum[ls[k]]+=(mid-l+1)*tag[k];
    sum[rs[k]]+=(r-mid)*tag[k];
    tag[k]=0;
}
void modify(int &k,int l,int r,int a,int b){
    if(!k) k=++tot;
    if(l==a&&b==r){
        tag[k]++;
        sum[k]+=r-l+1;
        return;
    }
    push_down(k,l,r);
    int mid=l+r>>1;
    if(b<=mid) modify(ls[k],l,mid,a,b);
    else if(a>mid) modify(rs[k],mid+1,r,a,b);
    else modify(ls[k],l,mid,a,mid),modify(rs[k],mid+1,r,mid+1,b);
    push_up(k);
}
void add(int a,int b,int c){
    int l=1,r=n,k=1;
    while(l!=r){
        int mid=l+r>>1;
        modify(root[k],1,n,a,b);
        if(c<=mid)r=mid,k*=2;
        else l=mid+1,k=k*2+1;
    }
    modify(root[k],1,n,a,b);
}
int query(int k,int l,int r,int a,int b){
    if(!k)return 0;
    push_down(k,l,r);
    if(l==a&&r==b) return sum[k];
    int mid=l+r>>1;
    if(b<=mid) return query(ls[k],l,mid,a,b);
    else if(a>mid) return query(rs[k],mid+1,r,a,b);
    else return query(ls[k],l,mid,a,mid)+query(rs[k],mid+1,r,mid+1,b);
}
int solve(int a,int b,int c){
    int l=1,r=n,k=1;
    while(l!=r){
        int mid=l+r>>1;
        int t=query(root[k*2],1,n,a,b);
        if(t>=c) r=mid,k*=2;
        else l=mid+1,k=k*2+1,c-=t;
    }
    return l;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        int opt,a,b,c;
        scanf("%d%d%d%d",&opt,&a,&b,&c);
        if(opt==1){
            add(a,b,n-c+1);
        }
        else printf("%d
",n-solve(a,b,c)+1);
    }
    return 0;
}