背景:
有的时候我们打树套树太麻烦了。
题目传送门:
https://www.luogu.org/problem/P2617
题意:
动态区间第
k
k
k小。
思路:
考虑支持修改的整体二分。
其实我们按照原顺序操作的答案一定是对的。
那么不妨按照顺序操作,如果是修改操作,直接用树状数组修改,即如果当前修改的数小于等于
m
i
d
mid
mid,直接按照要求加入柱状数组中。
在修改值的时候,其实就是将这个位置原来的数标记为
−
1
-1
−1,新的数标记为
1
1
1。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#define lowbit(x) ((x)&(-(x)))
using namespace std;
int n,m;
int a[200010],tr[200010],ans[200010];
struct node{int op,id,x,y,z;} b[400010],d1[400010],d2[400010];
void add(int x,int y)
{
for(;x<=n;x+=lowbit(x))
tr[x]+=y;
}
int findsum(int x)
{
int tot=0;
for(;x;x-=lowbit(x))
tot+=tr[x];
return tot;
}
int getsum(int l,int r)
{
return findsum(r)-findsum(l-1);
}
void dfs(int l,int r,int L,int R)
{
if(L>R) return;
if(l==r)
{
for(int i=L;i<=R;i++)
if(b[i].op==2) ans[b[i].id]=l;
return;
}
int mid=(l+r)>>1;
int t1=0,t2=0;
for(int i=L;i<=R;i++)
if(b[i].op==1)
{
if(b[i].x<=mid) add(b[i].id,b[i].y),d1[++t1]=b[i]; else d2[++t2]=b[i];
}
else
{
int op=getsum(b[i].x,b[i].y);
if(b[i].z<=op) d1[++t1]=b[i]; else b[i].z-=op,d2[++t2]=b[i];
}
for(int i=1;i<=t1;i++)
if(d1[i].op==1) add(d1[i].id,-d1[i].y);
for(int i=1;i<=t1;i++)
b[L+i-1]=d1[i];
for(int i=1;i<=t2;i++)
b[L+t1+i-1]=d2[i];
dfs(l,mid,L,L+t1-1);
dfs(mid+1,r,L+t1-1+1,R);
}
int main()
{
char s[5];
int op=0,cnt=0,mi=(int)1e9,ma=0,x,y,z;
scanf("%d %d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
mi=min(mi,a[i]);
ma=max(ma,a[i]);
b[++op]=(node){1,i,a[i],1,0};
}
for(int i=1;i<=m;i++)
{
scanf("%s",s+1);
if(s[1]=='Q')
{
scanf("%d %d %d",&x,&y,&z);
b[++op]=(node){2,++cnt,x,y,z};
}
else
{
scanf("%d %d",&x,&y);
b[++op]=(node){1,x,a[x],-1,0},b[++op]=(node){1,x,a[x]=y,1,0};
mi=min(mi,a[x]);
ma=max(ma,a[x]);
}
}
dfs(mi,ma,1,op);
for(int i=1;i<=cnt;i++)
printf("%d\n",ans[i]);
}