P2801 教主的魔法

传送门

$N$ 太大了主席树过不了

考虑分块

对每个块内的元素排序，询问就对大块二分查找，对两边小的部分暴力枚举

修改时维护 $add[i]$ 标记，维护当前块内整块已经加的数

那么整块的就直接修改 $add$ ，两边小的部分就把那两个的块暴力修改然后重新排序

然后注意一下边界就完了

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=2e6+,M=2e3+,INF=1e9+;
int n,m;
int a[N];
int bel[N],L[M];//bel[i]存位置i所在的块，L[i]存第i个块的左端点
int blk[M][M],tmp[M],add[M];//blk就是分出来的块排序后的东西
inline void change(int l,int r,int k)//修改
{
    if(bel[l]==bel[r])//要特判只有小块的情况
    {
        int tot=,be=bel[l];
        for(int i=L[be];i<l;i++) tmp[++tot]=a[i];//把数直接扔到tmp里
        for(int i=l;i<=r;i++) a[i]+=k,tmp[++tot]=a[i];//有修改要修改a[i]
        for(int i=r+;i<L[be+];i++) tmp[++tot]=a[i];
        sort(tmp+,tmp+tot+);
        memcpy(blk[be],tmp,sizeof(tmp));//排序后扔给blk
    }
    int pl=bel[l-]+,pr=bel[r+]-,tot=;//找到整块
    for(int i=pl;i<=pr;i++) add[i]+=k;
    if(pl>bel[])//处理左边小块
    {
        for(int i=L[pl-];i<l;i++) tmp[++tot]=a[i];
        for(int i=l;i<L[pl];i++) a[i]+=k,tmp[++tot]=a[i];
        sort(tmp+,tmp+tot+);
        memcpy(blk[pl-],tmp,sizeof(tmp));
    }
    if(pr<bel[n])//右边
    {
        for(int i=L[pr+];i<=r;i++) a[i]+=k,tmp[++tot]=a[i];
        for(int i=r+;i<L[pr+];i++) tmp[++tot]=a[i];
        sort(tmp+,tmp+tot+);
        memcpy(blk[pr+],tmp,sizeof(tmp));
    }
}
inline int query(int l,int r,int k)//处理询问
{
    if(bel[l]==bel[r])//特判没有整块
    {
        int res=;
        for(int i=l;i<=r;i++) if(a[i]>=k-add[bel[i]]) res++;
        return res;
    }
    int pl=bel[l-]+,pr=bel[r+]-,res=;
    for(int i=pl;i<=pr;i++)//对大块的二分
        res+=(L[i+]-L[i]) -( lower_bound(blk[i]+ , blk[i]+(L[i+]-L[i])+ , k-add[i]) -blk[i])+;
    for(int i=l;i<L[pl];i++) if(a[i]>=k-add[bel[i]]) res++;//剩下的暴力计算
    for(int i=L[pr+];i<=r;i++) if(a[i]>=k-add[bel[i]]) res++;
    return res;
}
int main()
{
    n=read(),m=read(); int t=sqrt(n);
    for(int i=;i<=n;i++)
    {
        bel[i]=(i-)/t+;
        if(bel[i]!=bel[i-]) L[bel[i]]=i;
        else L[bel[i]]=L[bel[i-]];
        a[i]=read();
    }
    bel[n+]=bel[n]+; L[bel[n+]]=n+;//注意可能会访问到n+1
    for(int i=;i<=bel[n];i++)//预处理blk
    {
        int tot=;
        for(int j=L[i];j<L[i+];j++) tmp[++tot]=a[j];
        sort(tmp+,tmp+tot+);
        memcpy(blk[i],tmp,sizeof(tmp));
    }
    char s[]; int x,y,z;
    while(m--)
    {
        scanf("%s",s); x=read(),y=read(),z=read();
        if(s[]=='M') change(x,y,z);
        else printf("%d\n",query(x,y,z));
    }
    return ;
}