简单分块2

数列分块2

整体思想：对每一个块内重新排序，然后在√n个块内做二分，并且暴力2√n个元素。但小段修改后，可能会乱序，所以我们要重新排序。至于区间加，用上一种方法即可。

时间复杂度：快排复杂度：O(nlogn)。而现在每一块是√n，一共√n块，所以是O(√n√nlog√n) ==O(nlog√n)  n是5000左右，√n是500多，log√n差不多是8左右，复杂度可以了。

区间统一修改：大段维护（记在标记数组里），小段暴力（直接改），两小短暴力改完要重新排序，保持有序。

查询：左到右区间内，比x小的数有多少？仍然大段维护（做二分，lower_bound()之前的都是小于x的，但注意要减去标记）；两头小段，直接暴力去做(两短总共大小不会超过2√n)

#include <bits/stdc++.h>
using namespace std;
const int MAX=50000+5;
typedef long long ll;

int v[MAX],bl[MAX],tar[MAX];
int n,blo;
vector<ll> ve[250];//用向量，方便啊~

inline ll read()//快读
{
	ll x=0,f=1;
	char ch=getchar();
	while(!isdigit(ch))
	{
		if(ch=='-') f=-1;
		ch=getchar();
	}
	while(isdigit(ch))
	{
		x=x*10+ch-'0';
		ch=getchar();
	}
	
	return x*f;
}

inline void Sort(int x)//对被暴力改过的块，重新排序
{
	ve[x].clear();
	for(int i=(x-1)*blo+1;i<=min(x*blo,n);i++)
	 ve[i].push_back(v[i]);
	sort(ve[x].begin(),ve[x].end());
}
inline void add(int a,int b,int c)//区间统一修改
{
	for(int i=a;i<=min(bl[a]*blo,b);i++) //小段暴力
	 v[i]+=c;
	Sort(bl[a]);//重新排序
	
	if(bl[a]!=bl[b])
	{
		for(int i=(bl[b]-1)*blo+1;i<=b;i++)
	     v[i]+=c;
	     Sort(bl[b]);
	 } 
	
	for(int i=bl[a]+1;i<=bl[b]-1;i++)
	 tar[i]+=c;
}

inline int search(int a,int b,int c)//查询
{
	int ans=0;
	for(int i=a;i<=min(bl[a]*blo,b);i++)//暴力查询
	if(v[i]+tar[i]<c) ans++;
    
    if(bl[a]!=bl[b])
     for(int i=(bl[b]-1)*blo+1;i<=b;i++)
      if(v[i]+tar[i]<c) ans++;
    
    for(int i=bl[a]+1;i<=bl[b]-1;i++)
    {
    	int x=c-tar[i];
    	ans+=lower_bound(ve[i].begin(),ve[i].end(),x)-ve[i].begin();//二分
	}

   return ans;
}
int main()
{
	n=read();
	blo=sqrt(n);
	
	for(int i=1;i<=n;i++) v[i]=read();
	
	for(int i=1;i<=n;i++)
	{
		bl[i]=(i-1)/blo+1;//每个元素属于哪一块
		ve[bl[i]].push_back(v[i]);
	}
	
      for(int i=1;i<=bl[n];i++) 
        sort(ve[i].begin(),ve[i].end());

	for(int i=1;i<=n;i++)
	{
		int f=read(),a=read(),b=read(),c=read();
		if(f==0) add(a,b,c);
		if(f==1) cout<<search(a,b,c*c)<<endl;
	}
	
	return 0;
}

写作不易，求赞ing 呜呜呜~~~