题目链接: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\n", val[L]);
}
}
return 0;
}