整体二分学习笔记 & luogu P3834 【模板】可持久化线段树 1（主席树）

背景：

有的时候树套树并不好打，所以我们需要一种优秀的做法。

题目传送门：

https://www.luogu.org/problem/P3834

题意：

求静态区间第$k$小。

思路：

想打线段树的请移步。
整体二分是什么东西。
个人感觉就是二分+分治。
考虑离线。
二分答案$mid$。
对于当前的区间来说，我们只需要统计小于等于$mid$的有$x$个，若$x\le k$，显然$nud$可以更小，答案在$[l,mid]$里；否则，$mid$必须更大，答案在$(mid,r]$里。可以用树状数组维护。
最后将这些按照上面的顺序做好，二分+分治即可。

如果我们考虑修改操作，容易发现按照顺序做结果一定是对的，因此并不需要另外特判什么的。
详细可见：$\text{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]);
}