bzoj3343: 教主的魔法 分块 标记

修改：两边暴力重构，中间打标记。复杂度：O(n^0.5)

查询：中间二分两边暴力。O(n^0.5logn^0.5)

总时间复杂度O(n*n^0.5logn^0.5)

空间复杂度是n级别的

标记不用下传因为标记不用下传也没时间下传，如果在访问时下传就造成更棘手的不整块，如果累加式下传不好记录，所以就带着就行。

细节很多......

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define MAXN 1010000
using namespace std;
typedef long long LL;
LL a[MAXN],b[][],mark[],n,pos[MAXN],lon,t;
inline LL Min(LL x,LL y)
{
   return x<y?x:y;
}
int comp(const int x,const int y)
{
   return x>y;
}
inline void work1()
{
   LL l,r,w;
   scanf("%lld%lld%lld",&l,&r,&w);
   LL z=pos[l]+,y=pos[r]-;
   for(LL i=z;i<=y;i++)
    mark[i]+=w;
   LL zzh=(pos[l]*lon);
   zzh=Min(zzh,r);
   for(LL i=l;i<=zzh;i++) a[i]+=w;
   zzh=pos[l]*lon;
   zzh=Min(zzh,n);
   memset(b[pos[r]],,sizeof(b[pos[r]]));
   for(LL i=zzh-lon+;i<=zzh;i++) b[pos[l]][i%lon+]=a[i];
   sort(b[pos[l]]+,b[pos[l]]+lon+);
   if(pos[l]!=pos[r])
   {
       for(LL i=(pos[r]-)*lon+;i<=r;i++) a[i]+=w;
       zzh=pos[r]*lon;
       zzh=Min(zzh,n);
       memset(b[pos[r]],,sizeof(b[pos[r]]));
       for(LL i=(pos[r]-)*lon+;i<=zzh;i++) b[pos[r]][i%lon+]=a[i];
       sort(b[pos[r]]+,b[pos[r]]+lon+);
   }
}
inline LL Count(LL p,LL c)
{
   LL ans=,l=,r=lon;
   while(l<=r)
   {
     LL mid=(l+r)>>;
     if(b[p][mid]+mark[p]<c)
      l=mid+;
     else
      ans+=(r-mid+),r=mid-;
   }
   return ans;
}
inline void work2()
{
   LL l,r,c,ans=;
   scanf("%lld%lld%lld",&l,&r,&c);
   LL z=pos[l]+,y=pos[r]-;
   for(LL i=z;i<=y;i++)
      ans+=Count(i,c);
   LL zzh=(pos[l]*lon);
   zzh=Min(zzh,r);
   for(LL i=l;i<=zzh;i++)
     if(a[i]+mark[pos[i]]>=c)
       ans++;
   if(pos[l]!=pos[r])
   {
      for(LL i=(pos[r]-)*lon+;i<=r;i++)
       if(a[i]+mark[pos[i]]>=c)
         ans++;
   }
   printf("%lld\n",ans);
}
int main()
{
   scanf("%lld",&n);
   LL Q;
   scanf("%lld",&Q);
   lon=(LL)sqrt(n+0.5);
   for(LL i=;i<=n;i++)
   {
      scanf("%lld",&a[i]);
      pos[i]=(i-)/lon+;
      b[pos[i]][i%lon+]=a[i];
   }
   t=pos[n];
   for(int i=;i<=t;i++)
    sort(b[i]+,b[i]+lon+);
   while(Q--)
   {
     char s[];
     scanf("%s",s);
     if(s[]=='M')work1();
     else work2();
   }
   return ;
}