【ZJOI2013】【BZOJ3110】K大数查询

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

a≤b≤N

1操作中abs(c)≤N

2操作中abs(c)≤Maxlongint

裸题.

树套树/主席树/整体二分

拿来练整体二分再好不过233

//AC code by CreationAugust
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 50010
#define MAXINT 0x7fffffff
#define lchild rt<<1,l,mid
#define rchild rt<<1|1,mid+1,r
#define ln rt<<1
#define rn rt<<1|1
using namespace std;
int op,a,b,c;
int n,m;
struct query
{
    int x;
    int op,a,b,c;
    int fin;//操作是否得到解决
    bool operator < (const query& a)const{
        if (fin!=a.fin) return fin<a.fin;
        return x<a.x;
    }
}ques[MAXN];//询问序列
int ans[MAXN];//答案序列,如果对应询问是插入,ans为-MAXINT
struct seg
{
    int l,r;
    int mem;//清空标记
    int flag;//操作标记
    int add;//加量标记
    int sum;
}tree[MAXN<<4];
bool comp(query a,query b)
{
    if (a.fin<b.fin) return 1;
    if (a.fin==b.fin) a.x<b.x;
    return 0;
}
void build(int rt=1,int l=1,int r=n)
{
    int mid=(l+r)>>1;
    tree[rt].l=l;tree[rt].r=r;
    if (l==r) return;
    build(lchild);build(rchild);
}
void push_up(int rt)
{
    tree[rt].sum=tree[ln].sum+tree[rn].sum;
}
void push_down(int rt)
{
    if (tree[rt].mem)
    {
        tree[ln].add=tree[rn].add=tree[ln].sum=tree[rn].sum=0;
        tree[ln].mem=tree[rn].mem=1;
        tree[rt].mem=0;
    }
    if (tree[rt].add)
    {
        tree[ln].add+=tree[rt].add;tree[rn].add+=tree[rt].add;
        tree[ln].sum+=(tree[ln].r-tree[ln].l+1)*tree[rt].add;
        tree[rn].sum+=(tree[rn].r-tree[rn].l+1)*tree[rt].add;
        tree[rt].add=0;
    }
}
void modify(int rt,int l,int r)
{
    push_down(rt);
    int L=tree[rt].l,R=tree[rt].r;
    if (L==l&&R==r)
    {
        tree[rt].add++;
        tree[rt].sum+=(R-L+1);
        return;
    }
    int mid=(L+R)>>1;
    if (r<=mid) modify(ln,l,r);
    else
    if (l>mid) modify(rn,l,r);
    else modify(ln,l,mid),modify(rn,mid+1,r);
    push_up(rt);
}
int query(int rt,int l,int r)
{
    push_down(rt);
    int L=tree[rt].l,R=tree[rt].r;
    if (L==l&&R==r) return tree[rt].sum;
    int mid=(L+R)>>1;
    if (r<=mid) return query(ln,l,r);
    else
    if (l>mid) return query(rn,l,r);
    else return query(ln,l,mid)+query(rn,mid+1,r);
}
void solve(int L,int R,int l,int r)//L,R是数据序列,l,r是询问序列
{
    if (l>r) return;
    if (L==R)//答案get
    {
        for (int i=l;i<=r;i++)
            if (ques[i].op==2) ans[ques[i].x]=L;
        return;
    }
    tree[1].add=tree[1].sum=0;tree[1].mem=1;
    int Mid=l-1;//用于确定下一次划分询问区间的端点
    int mid=(L+R)>>1;
    for (int i=l;i<=r;i++)
    {
        if (ques[i].op==1)
        {
            if (ques[i].c>mid)
                modify(1,ques[i].a,ques[i].b),ques[i].fin=1;
            else
                ques[i].fin=0,Mid++;
        }
        else
        {
            int temp=query(1,ques[i].a,ques[i].b);
            if (ques[i].c<=temp) ques[i].fin=1;
            else ques[i].fin=0,Mid++,ques[i].c-=temp;
        }
    }
    sort(ques+l,ques+r+1);
    solve(L,mid,l,Mid);
    solve(mid+1,R,Mid+1,r);
}
int main()
{
    scanf("%d%d",&n,&m);
    build();
    for (int i=1;i<=MAXN;i++) ans[i]=-MAXINT;
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&op,&a,&b,&c);
        ques[i].x=i;ques[i].op=op;ques[i].a=a;ques[i].b=b;ques[i].c=c;
    }
    solve(0,n,1,m);
    for (int i=1;i<=m;i++)
        if (ans[i]!=-MAXINT) printf("%d\n",ans[i]);
}