BZOJ.1901.Dynamic Rankings(整體二分)
阿新 • • 發佈:2018-06-27
gis esp #define dig strong efi 樹狀數組 fin pri
題目鏈接 BZOJ
洛谷
(以下是口胡)
對於多組的詢問、修改,我們可以發現:
假設有對p1,p2,p3...的詢問,在這之前有對p0的修改(比如+1),且p0<=p1,p2,p3...,那麽我們可以在修改完p0後對p1,p2,p3...這些詢問更改(比如需求-1),以後對於p1,p2,p3...都不需要管這個修改了。
這樣每次需要操作的序列長度只與當前區間有關。
而這個p0我們通過二分確定。不斷根據p0將操作劃分為兩部分(兩個隊列),當二分到底時(l==r),隊列裏所有詢問的答案就是l了。
復雜度\(O(nlog^2n)\)。
樹狀數組+主席樹見這.
線段樹套平衡樹見這.
#include <cstdio>
#include <cctype>
#define gc() getchar()
#define lb(x) (x&-x)
const int N=10005,M=N*3,INF=1e9;
int n,m,Q,A[N],Ans[N];
//char IN[MAXIN],*SS=IN,*TT=IN;
struct Operation
{
int K,l,r,pos;//K(opt)!=0:Query [l,r]:K; K(opt)=0:Modify 即將要改的數是l,貢獻值為r,下標為pos
Operation() {}
Operation(int K,int l,int r,int pos):K(K),l(l),r(r),pos(pos) {}
}q[M],q1[M],q2[M];
inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline bool Is_Query()
{
register char c=gc();
while (c!='Q'&&c!='C') c=gc();
return c=='Q';
}
namespace T
{
int n,t[N];
inline void Modify(int p,int v){
while(p<=n) t[p]+=v, p+=lb(p);
}
inline int Query(int p){
int res=0; while(p) res+=t[p], p-=lb(p);
return res;
}
inline void Clear(int p){
while(p<=n && t[p]) t[p]=0, p+=lb(p);
}
}
void Solve(int l,int r,int h,int t)
{
if(h>t /*|| l>r*/) return;
if(l==r){
for(int i=h; i<=t; ++i) if(q[i].K) Ans[q[i].pos]=l;
return;
}
int mid=l+r>>1, t1=0, t2=0;
for(int i=h; i<=t; ++i)
if(q[i].K)//Query
{
int tmp=T::Query(q[i].r)-T::Query(q[i].l-1);
if(tmp>=q[i].K) q1[t1++]=q[i];
else q[i].K-=tmp, q2[t2++]=q[i];
}
else
{
if(q[i].l<=mid) q1[t1++]=q[i], T::Modify(q[i].pos,q[i].r);
else q2[t2++]=q[i];
}
for(int i=0; i<t1; ++i) if(!q1[i].K) T::Modify(q1[i].pos,-q1[i].r);//T::Clear(q1[t].pos);
for(int i=0; i<t1; ++i) q[h+i]=q1[i];
for(int i=0; i<t2; ++i) q[h+t1+i]=q2[i];
Solve(l,mid,h,h+t1-1), Solve(mid+1,r,h+t1,t);
}
int main()
{
n=read(), m=read(), Q=n;
for(int i=1; i<=n; ++i) q[i]=Operation(0,A[i]=read(),1,i);
int tot=0;
for(int l,r,p,i=1; i<=m; ++i)
{
if(Is_Query()) l=read(), r=read(), q[++Q]=Operation(read(),l,r,++tot);
else p=read(), q[++Q]=Operation(0,A[p],-1,p), q[++Q]=Operation(0,A[p]=read(),1,p);
}
T::n=n, Solve(-INF,INF,1,Q);
for(int i=1; i<=tot; ++i) printf("%d\n",Ans[i]);
return 0;
}
BZOJ.1901.Dynamic Rankings(整體二分)