[BZOJ3343]教主的魔法（分块）

题目：

我是超链接

题解：

先说几句废话
分块作为一种常见的数据结构设计思路，很经常用于部分分，博主比较蒟，考试时可能只写些部分分，那就切入分块吧，反正博主已经被计算几何吓哭了/(ㄒoㄒ)/~~

然而这道题就是正经不做作的分块了嘛
大概思想是这样的，如果修改覆盖了整个块，那么就在块的总标记上加上修改的数，如果是块的一部分，那么就直接在那个数值上加上修改的数，然后用一个b数组辅助处理，把b数组在每个块内按从小到大排序，用于二分查找第一个小于要求的数的位置。

代码：

#include <algorithm>
#include <cstdio>
#include <cmath>
using namespace std;
const int N=1000005;
int a[N],b[N],pos[N],cnt,n,add[N];
void insert(int x)//第x个块内的情况 
{
    int l=(x-1)*cnt+1,r=min(n,x*cnt);
    for (int i=l;i<=r;i++) b[i]=a[i];
    sort(b+l,b+r+1);
}
void change(int l,int r,int w)
{
    if (pos[l]==pos[r])
    {
        for (int i=l;i<=r;i++) a[i]+=w;
        insert(pos[l]);
    }else
    {
        int up=min(n,pos[l]*cnt);
        for (int i=l;i<=up;i++) a[i]+=w;
        insert(pos[l]);
        up=(pos[r]-1)*cnt+1;
        for (int i=up;i<=r;i++) a[i]+=w;
        insert(pos[r]);
        for (int i=pos[l]+1;i<pos[r];i++) add[i]+=w;
    }
}
int find(int x,int w)
{
    int l=(x-1)*cnt+1,r=min(n,x*cnt),R=r;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (b[mid]>=w) r=mid-1;
        else l=mid+1;
    }
    return R-l+1;
}
int qurry(int l,int r,int w)
{
    int ans=0;
    if (pos[l]==pos[r])
    {
        for (int i=l;i<=r;i++)
          if (a[i]+add[pos[l]]>=w) ans++; 
    }else
    {
        int up=min(n,pos[l]*cnt);
        for (int i=l;i<=up;i++) 
          if (a[i]+add[pos[l]]>=w) ans++; 
        up=(pos[r]-1)*cnt+1;
        for (int i=up;i<=r;i++)
          if (a[i]+add[pos[r]]>=w) ans++; 
        for (int i=pos[l]+1;i<pos[r];i++) ans+=find(i,w-add[i]);//因为add[i]并未被编入修改b 
    }
    return ans;
}
int main()
{
    int q,l,r,w,m;scanf("%d%d",&n,&q);
    cnt=sqrt(n);
    for (int i=1;i<=n;i++) pos[i]=(i-1)/cnt+1;
    for (int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i];
    if (n%cnt) m=n/cnt+1;else m=n/cnt;
    for (int i=1;i<=m;i++) insert(i);
    while (q--)
    {
        char st[5];scanf("%s",st);
        scanf("%d%d%d",&l,&r,&w);
        if (st[0]=='M') change(l,r,w);
        else printf("%d\n",qurry(l,r,w));
    }
}