bzoj3110 [Zjoi2013]K大数查询（权值线段树套线段树）

外面一棵线段树维护权值，每个节点一颗线段树维护区间，要动态加点。
1操作相当于在外层线段树上单点修改，内层线段树区间加法。
最好维护永久标记，不要pushdown。
2操作在外层线段树上二分爬，内层线段树求区间和。
注意会爆int。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 50010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,owo,rt[N<<2],a[N],tot=0,L,R;
struct node{
    int lc,rc,cov;ll sum;//sum--不完全覆盖的和
}tr[N*400];
struct oper{
    int op,l,r,c;
}e[N];
inline ll qsum(int p,int l,int r,int x,int y){
    if(!p) return 0;ll res=(ll)tr[p].cov*(y-x+1);
    if(x==l&&r==y) return res+tr[p].sum;int mid=l+r>>1;
    if(y<=mid) return res+qsum(tr[p].lc,l,mid,x,y);
    if(x>mid) return res+qsum(tr[p].rc,mid+1,r,x,y);
    return res+qsum(tr[p].lc,l,mid,x,mid)+qsum(tr[p].rc,mid+1,r,mid+1,y);
}
inline void add(int &p,int l,int r,int x,int y){
    if(!p) p=++owo;
    if(x==l&&r==y){tr[p].cov++;return;}int mid=l+r>>1;tr[p].sum+=y-x+1;
    if(y<=mid) add(tr[p].lc,l,mid,x,y);
    else if(x>mid) add(tr[p].rc,mid+1,r,x,y);
    else add(tr[p].lc,l,mid,x,mid),add(tr[p].rc,mid+1,r,mid+1,y);
}
inline int ask(int p,int l,int r,int x){
    if(l==r) return l;
    ll sz=qsum(rt[p<<1|1],1,n,L,R);int mid=l+r>>1;
    if(x<=sz) return ask(p<<1|1,mid+1,r,x);
    return ask(p<<1,l,mid,x-sz);
}
inline void ins(int p,int l,int r,int x){
    add(rt[p],1,n,L,R);
    if(l==r) return;int mid=l+r>>1;
    if(x<=mid) ins(p<<1,l,mid,x);
    else ins(p<<1|1,mid+1,r,x);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=m;++i){
        e[i].op=read(),e[i].l=read(),e[i].r=read(),e[i].c=read();
        if(e[i].op==1) a[++tot]=e[i].c;
    }sort(a+1,a+tot+1);tot=unique(a+1,a+tot+1)-a-1;
    for(int i=1;i<=m;++i){
        L=e[i].l;R=e[i].r;
        if(e[i].op==2){printf("%d\n",a[ask(1,1,tot,e[i].c)]);continue;}
        e[i].c=lower_bound(a+1,a+tot+1,e[i].c)-a;ins(1,1,tot,e[i].c);
    }return 0;
}