ZJOI 2013 K大数查询

题目链接：bzoj点我:-) 洛谷点我:-)
题目描述：
有N个位置，M个操作。操作有两种，每次操作如果是1 a b c的形式表示在第a个位置到第b个位置，每个位置加入一个数c如果是2 a b c形式，表示询问从第a个位置到第b个位置，第C大的数是多少。

输入格式：
第一行N，M接下来M行，每行形如1 a b c或2 a b c

输出格式：
输出每个询问的结果

数据说明：
N,M<=50000,N,M<=50000
a<=b<=N
1操作中abs(c)<=N
2操作中c<=Maxlongint

思路：
本题有多法可解，比如权值线段树套区间线段树，线段树套splay，线段树套主席树，树状数组套线段树，整体二分，CDQ分治等
介绍最简单的权值线段树套区间线段树：
外层即以权值为线段树区间，内层以原区间为线段树区间，树套树。
1操作为外层线段树单点修改，内层线段树区间修改
2操作为外层线段树区间修改，内层线段树区间修改

比较卡空间，需把内层线段树动态开点。并且注意离散化和long long

感想：
第一次写树套树，怎么这么卡我空间。。。洛谷不让过，硬是讲我MLE，bzoj较为善良让我过了，
但是。。居然有17000+ms。。莫非我天生自带大常数？！！

代码：

//miaomiao 2017.2.11
#include<cstdio>
#include<algorithm>
#include<map>

using namespace std;

#define LL long long
#define mid ((L+R)>>1)
#define For(i, a, b) for(int i = (a); i <= (int)(b); i++)
#define N (50000+1)
#define M (12000000+1)

struct node{
    int op, a, b;
    LL c;
}qus[N];

int lc[M], rc[M], root[N<<4], xn, num[N], val[N], ql, qr, cnt, mn;
LL sum[M], add[M];
map<LL, int> ms;

void modify(int &o, int L, int R){
    if(!o) o = ++cnt;
    if(ql <= L && R <= qr){add[o]++; sum[o]+=(R-L+1); return;}

    if(ql <= mid) modify(lc[o], L, mid); if(qr > mid) modify(rc[o], mid+1, R);
    sum[o] = sum[lc[o]]+sum[rc[o]]+add[o]*(R-L+1);
}

LL query(int &o, int L, int R, LL Add){
    if(!o) o = ++cnt;
    if(ql <= L && qr >= R) return sum[o]+Add*(R-L+1);
    LL ret = 0;
    if(ql <= mid) ret += query(lc[o], L, mid, Add+add[o]);
    if(qr > mid) ret += query(rc[o], mid+1, R, Add+add[o]);
    return ret;
}

int main(){
    int n, m, now, L, R;
    LL c, t;
    scanf("%d%d", &n, &m);
    For(i, 1, m){
        scanf("%d%d%d%lld", &qus[i].op, &qus[i].a, &qus[i].b, &qus[i].c);
        if(qus[i].op == 1) num[++xn] = qus[i].c;
    }
    sort(num+1, num+xn+1);
    For(i, 1, xn)
        if(num[i] != num[i-1] || i==1) ms[num[i]] = ++mn, val[mn] = num[i];

    For(i, 1, m){
        ql = qus[i].a; qr = qus[i].b; c = qus[i].c;
        now = 1; L = 1; R = n;
        if(qus[i].op == 1){
            c = ms[c];
            while(true){
                modify(root[now], 1, n);
                if(L == R) break;
                if(c <= mid) now = now<<1, R = mid;
                else now = now<<1|1, L = mid+1;
            }
        }else{
            while(L < R){
                t = query(root[now<<1|1], 1, n, 0);
                if(t >= c) L = mid+1, now = now<<1|1;    
                else R = mid, now = now<<1, c -= t;
            }
            printf("%d
", val[L]);
        }
    }

    return 0;
}