背景:
有的时候树套树并不好打,所以我们需要一种优秀的做法。
题目传送门:
https://www.luogu.org/problem/P3834
题意:
求静态区间第
k
k
k小。
思路:
想打线段树的请移步。
整体二分是什么东西。
个人感觉就是二分+分治。
考虑离线。
二分答案
m
i
d
mid
mid。
对于当前的区间来说,我们只需要统计小于等于
m
i
d
mid
mid的有
x
x
x个,若
x
≤
k
x≤k
x≤k,显然
n
u
d
nud
nud可以更小,答案在
[
l
,
m
i
d
]
[l,mid]
[l,mid]里;否则,
m
i
d
mid
mid必须更大,答案在
(
m
i
d
,
r
]
(mid,r]
(mid,r]里。可以用树状数组维护。
最后将这些按照上面的顺序做好,二分+分治即可。
如果我们考虑修改操作,容易发现按照顺序做结果一定是对的,因此并不需要另外特判什么的。
详细可见:
luogu P2617 Dynamic Rankings
\text{luogu P2617 Dynamic Rankings}
luogu P2617 Dynamic Rankings。
代码:
这是支持修改的。
#include<cstdio>
#include<cstring>
#include<algorithm>
#define lowbit(x) ((x)&(-(x)))
using namespace std;
int n,m;
int 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()
{
int mi=(int)1e9,ma=0,x,y,z;
scanf("%d %d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
mi=min(mi,x);
ma=max(ma,x);
b[i]=(node){1,i,x,1,0};
}
for(int i=1;i<=m;i++)
{
scanf("%d %d %d",&x,&y,&z);
b[i+n]=(node){2,i,x,y,z};
}
dfs(mi,ma,1,n+m);
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
}