LibreOJ #6279. 数列分块入门 3 分块

https://loj.ac/problem/6279

思路：区间加法同上，因为要询问比$x$小的最大元素，$copy$原序列分块后对每一块进行排序，这样在完整的块内就可以利用二分来快速判断。因为修改操作是对原序列进行的，对完整的块没有影响，但是对不完整的块是有影响的，这时候暴力修改后应该对该块重新$copy$并排序。且在查询不完整的块的时候，用的是原序列的值而不是$copy$的序列的值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#define INF 0x3f3f3f3f
using namespace std;

int a[100005];//原始序列
int b[100005];//进行操作的序列
int l[2000];//第i个块的左区间端点
int r[2000];//第i个块的右区间端点
int pos[100005];//第i个点属于第几个块
int lazy[2000];

int read()
{
	int temp=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')
			f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		temp=temp*10+ch-48;
		ch=getchar();
	}
	return temp*f;
}

void add(int ll,int rr,int v)
{
	int posl=pos[ll];
	int posr=pos[rr];
	if(posl==posr)
	{
		for(int i=ll;i<=rr;i++)	//不完整的块 暴力
		{
			a[i]+=v;
			b[i]=a[i];
		}
		for(int i=l[posl];i<ll;i++)
			b[i]=a[i];
		for(int i=rr+1;i<=r[posl];i++)
			b[i]=a[i];
		sort(b+l[posl],b+r[posl]+1);
	}
	else
	{
		for(int i=ll;i<=r[posl];i++)	//不完整的块 暴力
		{
			a[i]+=v;
			b[i]=a[i];
		}
		for(int i=l[posl];i<ll;i++)
			b[i]=a[i];
		sort(b+l[posl],b+r[posl]+1);
		for(int i=pos[ll]+1;i<=pos[rr]-1;i++)//完整的块打上标记
			lazy[i]+=v;
		for(int i=l[posr];i<=rr;i++)//不完整的块 暴力
		{
			a[i]+=v;
			b[i]=a[i];
		}
		for(int i=rr+1;i<=r[posr];i++)
			b[i]=a[i];
		sort(b+l[posr],b+r[posr]+1);
	}
}

int query(int ll,int rr,int v)
{
	int MAX=-1;
	int posl=pos[ll];
	int posr=pos[rr];
	if(posl==posr)
	{
		for(int i=ll;i<=rr;i++)
			if(a[i]+lazy[posl]<v&&a[i]+lazy[posl]>MAX)
				MAX=a[i]+lazy[posl];
		return MAX;
	}
	for(int i=ll;i<=r[posl];i++)
			if(a[i]+lazy[posl]<v&&a[i]+lazy[posl]>MAX)
				MAX=a[i]+lazy[posl];
	for(int i=l[posr];i<=rr;i++)
		if(a[i]+lazy[posr]<v&&a[i]+lazy[posr]>MAX)
				MAX=a[i]+lazy[posr];
	int temp;
	for(int i=pos[ll]+1;i<=pos[rr]-1;i++)
	{
		temp=lower_bound(b+l[i],b+r[i]+1,v-lazy[i])-(b+l[i]);
		if(temp==0)//最小的元素都不满足题意
			continue;
		MAX=max(MAX,b[l[i]+temp-1]+lazy[i]);
	}
	return MAX;
}
int main()
{
	int n;
	n=read();
	int dis=sqrt(n);
	int num=ceil(n*1.0/dis);
	for(int i=1;i<=num;i++)
	{
		l[i]=(i-1)*dis+1;
		r[i]=i*dis;
	}
	r[num]=n;
	for(int i=1;i<=n;i++)
	{
		a[i]=read();
		pos[i]=(i-1)/dis+1;
		b[i]=a[i];
	}
	for(int i=1;i<=num;i++)
		sort(b+l[i],b+r[i]+1);
	int t1,t2,t3,t4;
	for(int i=0;i<n;i++)
	{
		t1=read();
		t2=read();
		t3=read();
		t4=read();
		if(t1==0)
			add(t2,t3,t4);
		else
			printf("%d\n",query(t2,t3,t4));
	}
	return 0;
}