题目:
题解:
先说几句废话
分块作为一种常见的数据结构设计思路,很经常用于部分分,博主比较蒟,考试时可能只写些部分分,那就切入分块吧,反正博主已经被计算几何吓哭了/(ㄒ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));
}
}