【分块】bzoj3343 教主的魔法

分块，每个块内部有序化，查询时二分。修改时直接修改覆盖在各个块上的delta。

块大小O(sqrt(n*log(n)))比较合适。

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
#define N 1000001
char op; int n,m,x,y,v,a[N],l[4500],r[4500],num[N],b[N],delta[4500];
int Num,CH[12],f,c;
inline void R(int &x){
    c=0;f=1;
    for(;c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
    for(x=0;c>='0'&&c<='9';c=getchar())(x*=10)+=(c-'0');
    x*=f;
}
inline void P(int x){
    if(x<10)putchar(x+'0');
    else{P(x/10);putchar(x%10+'0');}
}
void makeblock()
{
    for(int i=1;i<=n;++i) b[i]=a[i];
    int sz=sqrt((double)n*log2(n)),sum=1; if(!sz) sz=1;
    for(;sum*sz<n;++sum)
      {
          l[sum]=r[sum-1]+1; r[sum]=sum*sz;
          for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
          sort(b+l[sum],b+r[sum]+1);
      }
    l[sum]=r[sum-1]+1; r[sum]=n;
    for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
    sort(b+l[sum],b+r[sum]+1);
}
void work(const int &L,const int &R)
{
    for(int i=L;i<=R;++i) a[i]+=v;
    for(int i=l[num[L]];i<=r[num[L]];++i) b[i]=a[i];
    sort(b+l[num[L]],b+r[num[L]]+1);
}
inline void update()
{
    if(num[x]==num[y]) work(x,y);
    else
      {
          work(x,r[num[x]]); work(l[num[y]],y);
          for(int i=num[x]+1;i<num[y];++i) delta[i]+=v;
      }
}
inline void query()
{
    int res=0;
    if(num[x]==num[y]) {for(int i=x;i<=y;++i) if(a[i]+delta[num[x]]>=v) ++res;}
    else
      {
          for(int i=x;i<=r[num[x]];++i) if(a[i]+delta[num[x]]>=v) ++res;
          for(int i=l[num[y]];i<=y;++i) if(a[i]+delta[num[y]]>=v) ++res;
          for(int i=num[x]+1;i<num[y];++i)
          res+=(b+r[i]-lower_bound(b+l[i],b+r[i]+1,v-delta[i])+1);
      } P(res); puts("");
}
int main()
{
    R(n); R(m); for(int i=1;i<=n;++i) R(a[i]); makeblock();
    for(int i=1;i<=m;++i)
      {
          getchar(); op=getchar(); R(x); R(y); R(v);
          if(op=='M') update(); else query();
      }
    return 0;
}