BZOJ 3110 [ZJOI2013]K大数查询 (整体二分+线段树)

和dynamic rankings这道题的思想一样

只不过是把树状数组换成线段树区间修改，求第$K$大的而不是第$K$小的

这道题还有负数，需要离散

 #include <vector>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define N1 50500
 #define M1 100500
 #define ll long long
 #define dd double
 #define inf 233333333
 using namespace std;

 int gint()
 {
     int ret=,fh=;char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }
 ll gll()
 {
     ll ret=;int fh=; char c=getchar();
     while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
     while(c>=''&&c<=''){ret=ret*+c-'';c=getchar();}
     return ret*fh;
 }

 struct SEG{
 ll sum[N1<<],tag[N1<<];int use[N1<<],que[N1<<],tl;
 void pushdown(int l,int r,int rt)
 {
     if(!tag[rt]) return;
     int mid=(l+r)>>;
     sum[rt<<]+=tag[rt]*(mid-l+); tag[rt<<]+=tag[rt];
     sum[rt<<|]+=tag[rt]*(r-mid); tag[rt<<|]+=tag[rt];
     if(!use[rt<<]) use[rt<<]=,que[++tl]=rt<<;
     if(!use[rt<<|]) use[rt<<|]=,que[++tl]=rt<<|;
     tag[rt]=;
 }
 void pushup(int rt){ sum[rt]=sum[rt<<]+sum[rt<<|]; }
 void update(int L,int R,int l,int r,int rt,int w)
 {
     if(!use[rt]) use[rt]=,que[++tl]=rt;
     if(L<=l&&r<=R){ sum[rt]+=(r-l+)*w; tag[rt]+=w; return;}
     int mid=(l+r)>>; pushdown(l,r,rt);
     if(L<=mid) update(L,R,l,mid,rt<<,w);
     if(R>mid) update(L,R,mid+,r,rt<<|,w);
     pushup(rt);
 }
 ll query(int L,int R,int l,int r,int rt)
 {
     if(L<=l&&r<=R) return sum[rt];
     int mid=(l+r)>>; ll ans=; pushdown(l,r,rt);
     if(L<=mid) ans+=query(L,R,l,mid,rt<<);
     if(R>mid) ans+=query(L,R,mid+,r,rt<<|);
     return ans;
 }
 void clr(){ while(tl){ use[que[tl]]=; sum[que[tl]]=; tag[que[tl]]=; que[tl--]=; } }
 }s;
 struct node{int fl,l,r,ans,t;ll K;}a[N1],tmp[N1];
 int n,m,ma,nn;
 int f[N1],w[N1],use[N1];

 void alldic(int l,int r,int ql,int qr)
 {
     if(l>r||ql>qr) return;
     int mid=(l+r)>>,i,j,S=ql,E;ll sum;
     for(i=ql;i<=qr;i++)
     {
         if(a[i].fl==){
             if(a[i].K>=mid){ s.update(a[i].l,a[i].r,,n,,); use[i]=; }
             if(a[i].K>mid){ tmp[S++]=a[i]; }
         }else{
             sum=s.query(a[i].l,a[i].r,,n,);
             if(sum<a[i].K){ a[i].K-=sum; }
             else{ f[a[i].t]=w[mid]; tmp[S++]=a[i]; use[i]=; }
         }
     }
     for(i=ql,E=S;i<=qr;i++)
     {
         if(!use[i]){ tmp[E++]=a[i]; }
         else{ use[i]=; }
     }
     for(i=ql;i<=qr;i++) a[i]=tmp[i];
     s.clr();
     alldic(mid+,r,ql,S-); alldic(l,mid-,S,E-);
 } 

 int main()
 {
     scanf("%d%d",&n,&m);
     int i,j; memset(f,-,sizeof(f));
     for(i=;i<=m;i++){ a[i].fl=gint(); a[i].l=gint(); a[i].r=gint(); a[i].K=gll(); a[i].t=i; if(a[i].fl==) w[++ma]=a[i].K; }
     sort(w+,w+ma+); ma=unique(w+,w+ma+)-(w+);
     for(i=;i<=m;i++){ if(a[i].fl==) a[i].K=lower_bound(w+,w+ma+,a[i].K)-w; }
     alldic(,ma,,m);
     for(i=;i<=m;i++) if(f[i]!=-) printf("%d\n",f[i]);
     return ;
 }