BZOJ3110:[ZJOI2013]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,N,M<=50000

a<=b<=N

1操作中abs(c)<=N

2操作中c<=Maxlongint

Solution

整体二分板子题哪来的题解，笔记又懒得写这样子。

线段树写错不做人了。爆$int$了不做人了。懒得一个一个改$long~long$直接$define$了……慢点就慢点吧……

Code

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #define N (50009)
 #define int long long
 using namespace std;

 struct Que{int opt,a,b,c,pos;}Q[N],q1[N],q2[N];
 struct Sgt{int val,add;}Segt[N<<];
 int n,m,cnt,ans[N];

 inline int read()
 {
     int x=,w=; char c=getchar();
     while (c<'' || c>'') {if (c=='-') w=-; c=getchar();}
     while (c>='' && c<='') x=x*+c-'', c=getchar();
     return x*w;
 }

 void Pushdown(int now,int l,int r)
 {
     if (Segt[now].add)
     {
         int mid=(l+r)>>;
         Segt[now<<].add+=Segt[now].add;
         Segt[now<<|].add+=Segt[now].add;
         Segt[now<<].val+=Segt[now].add*(mid-l+);
         Segt[now<<|].val+=Segt[now].add*(r-mid);
         Segt[now].add=;
     }
 }

 void Update(int now,int l,int r,int l1,int r1,int k)
 {
     if (l>r1 || r<l1) return;
     if (l1<=l && r<=r1)
     {
         Segt[now].add+=k;
         Segt[now].val+=(r-l+)*k;
         return;
     }
     int mid=(l+r)>>; Pushdown(now,l,r);
     Update(now<<,l,mid,l1,r1,k); Update(now<<|,mid+,r,l1,r1,k);
     Segt[now].val=Segt[now<<].val+Segt[now<<|].val;
 }

 int Query(int now,int l,int r,int l1,int r1)
 {
     if (l>r1 || r<l1) return ;
     if (l1<=l && r<=r1) return Segt[now].val;
     int mid=(l+r)>>; Pushdown(now,l,r);
     return Query(now<<,l,mid,l1,r1)+Query(now<<|,mid+,r,l1,r1);
 }

 void Solve(int l,int r,int L,int R)
 {
     if (l>r || L>R) return;
     if (l==r)
     {
         for (int i=L; i<=R; ++i) if (Q[i].opt==) ans[Q[i].pos]=l;
         return;
     }
     int mid=(l+r)>>,cnt1=,cnt2=;
     for (int i=L; i<=R; ++i)
         if (Q[i].opt==)
         {
             if (Q[i].c>mid) Update(,,n,Q[i].a,Q[i].b,), q2[++cnt2]=Q[i];
             else q1[++cnt1]=Q[i];
         }
         else
         {
             int now=Query(,,n,Q[i].a,Q[i].b);
             if (now>=Q[i].c) q2[++cnt2]=Q[i];
             else Q[i].c-=now, q1[++cnt1]=Q[i];
         }
     for (int i=L; i<=R; ++i) if (Q[i].opt== && Q[i].c>mid) Update(,,n,Q[i].a,Q[i].b,-);
     for (int i=; i<=cnt1; ++i) Q[L+i-]=q1[i];
     for (int i=; i<=cnt2; ++i) Q[L+cnt1+i-]=q2[i];
     Solve(l,mid,L,L+cnt1-); Solve(mid+,r,L+cnt1,R);
 }
 #undef int
 int main()
 #define int long long
 {
     n=read(); m=read();
     for (int i=; i<=m; ++i)
     {
         int opt=read(),a=read(),b=read(),c=read();
         if (opt==) ++cnt;
         Q[i]=(Que){opt,a,b,c,cnt};
     }
     Solve(,n,,m);
     for (int i=; i<=cnt; ++i) printf("%lld\n",ans[i]);
 }