「BZOJ3343」教主的魔法(分块+二分查找）

题意：

给定一个数列，您需要支持以下两种操作：
给[l,r]同加一个数
询问[l,r]中有多少数字大于或等于v

（n<=1000000,m<=3000)

题解

块内排序二分查询
修改就用个数组存整块的修改值
不完整的部分都暴力修改和查询
时间复杂度大概O(Qsqrt(n)logn)

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 const int N=;
 int n,m,a[N],b[N],block[N],L[],R[],cnt[],ans,Block;
 int find(int x,int y){
     int l=L[x];int r=R[x];
     int tmp=;
     while(l<=r){
         int mid=(l+r)>>;
         if(b[mid]>=y){
             tmp=mid;
             r=mid-;
         }
         else l=mid+;
     }
     return max(,R[x]-tmp+);
 }
 int main(){
     scanf("%d%d",&n,&m);
     Block=sqrt(n);
     for(int i=;i<=n;i++){
         scanf("%d",&a[i]);
         b[i]=a[i];
         block[i]=(i-)/Block+;
         if(!L[block[i]])L[block[i]]=i;
         R[block[i]]=i;
     }
     for(int i=;i<=block[n];i++)sort(b+L[i],b++R[i]);
     char s[];
     while(m--){
         scanf("%s",s);
         if(s[]=='A'){
             int l,r,c;
             ans=;
             scanf("%d%d%d",&l,&r,&c);
             if(block[l]==block[r]){
                 for(int i=l;i<=r;i++)if(a[i]+cnt[block[i]]>=c)ans++;
                 printf("%d\n",ans);
                 continue;
             }
             for(int i=l;i<=R[block[l]];i++)if(a[i]+cnt[block[i]]>=c)ans++;
             for(int i=L[block[r]];i<=r;i++)if(a[i]+cnt[block[i]]>=c)ans++;
             if(block[l]+<block[r])
                 for(int i=block[l]+;i<=block[r]-;i++)ans+=find(i,c-cnt[i]);
             printf("%d\n",ans);
         }
         else{
             int l,r,c;
             scanf("%d%d%d",&l,&r,&c);
             if(block[l]+<block[r])
                 for(int i=block[l]+;i<=block[r]-;i++)cnt[i]+=c;
             for(int i=l;i<=R[block[l]];i++)a[i]+=c;
             for(int i=L[block[l]];i<=R[block[l]];i++)b[i]=a[i];
             sort(b+L[block[l]],b++R[block[l]]);
             for(int i=L[block[r]];i<=r;i++)a[i]+=c;
             for(int i=L[block[r]];i<=R[block[r]];i++)b[i]=a[i];
             sort(b+L[block[r]],b++R[block[r]]);
         }
     }
     return ;
 }