BZOJ3343 & 洛谷2801：教主的魔法—

http://www.lydsy.com/JudgeOnline/problem.php?id=3343

https://www.luogu.org/problemnew/show/2801

题目描述

教主最近学会了一种神奇的魔法，能够使人长高。于是他准备演示给XMYZ信息组每个英雄看。于是N个英雄们又一次聚集在了一起，这次他们排成了一列，被编号为1、2、……、N。

每个人的身高一开始都是不超过1000的正整数。教主的魔法每次可以把闭区间[L, R]（1≤L≤R≤N）内的英雄的身高全部加上一个整数W。（虽然L=R时并不符合区间的书写规范，但我们可以认为是单独增加第L（R）个英雄的身高）

CYZ、光哥和ZJQ等人不信教主的邪，于是他们有时候会问WD闭区间 [L, R] 内有多少英雄身高大于等于C，以验证教主的魔法是否真的有效。

WD巨懒，于是他把这个回答的任务交给了你。

输入输出格式

输入格式：

第1行为两个整数N、Q。Q为问题数与教主的施法数总和。

第2行有N个正整数，第i个数代表第i个英雄的身高。

第3到第Q+2行每行有一个操作：

（1）若第一个字母为“M”，则紧接着有三个数字L、R、W。表示对闭区间 [L, R] 内所有英雄的身高加上W。

（2）若第一个字母为“A”，则紧接着有三个数字L、R、C。询问闭区间 [L, R] 内有多少英雄的身高大于等于C。

输出格式：

对每个“A”询问输出一行，仅含一个整数，表示闭区间 [L, R] 内身高大于等于C的英雄数。

输入输出样例

输入样例#1：

输出样例#1：

2

3

————————————————————————————

带修改分块。

我们存两个数组，一个原数组，一个对每个块排好序的新数组。

用线段树lazy的思想对于每一块的修改存一个lazy进去，而对于非整块来说只要暴力改即可。

询问的时候非整块直接暴力，整块二分查新数组即可。

#include<cstdio>

#include<queue>

#include<cctype>

#include<cstring>

#include<cmath>

#include<iostream>

#include<algorithm>

using namespace std;

const int N=;

const int SQRTN=;

const int INF=;

inline int read(){

    int X=,w=;char ch=;

    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}

    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();

    return w?-X:X;

}

int n,m,lim,s,cnt,a[N],b[N],bl[SQRTN],br[SQRTN],lazy[SQRTN];

inline void intoblock(){

    for(int i=;i<=n;i++){

        if(i%s==){br[cnt]=i-;bl[++cnt]=i;}

    }

    br[cnt]=n;bl[cnt+]=n+;

    for(int i=;i<=cnt;i++)sort(b+bl[i],b+br[i]+);

    return;

}

inline void add(int l,int r,int w){

    int L=(l-)/s+,R=(r-)/s+;

    if(r-l+<=*s){

        for(int i=l;i<=r;i++)a[i]+=w;

    for(int i=L;i<=R;i++){

        for(int j=bl[i];j<=br[i];j++){

        b[j]=a[j];

        }

        sort(b+bl[i],b+br[i]+);

    }

        return;

    }

    for(int i=L+;i<=R-;i++)lazy[i]+=w;

    for(int i=l;i<=br[L];i++)a[i]+=w;

    for(int i=bl[L];i<=br[L];i++)b[i]=a[i];

    sort(b+bl[L],b+br[L]+);

    for(int i=bl[R];i<=r;i++)a[i]+=w;

    for(int i=bl[R];i<=br[R];i++)b[i]=a[i];

    sort(b+bl[R],b+br[R]+);

    return;

}

int find(int l,int r,int inc,int c){

    while(l<r){

    int mid=(l+r)>>;

    if(b[mid]+inc<c)l=mid+;

    else r=mid;

    }

    return r;

}

inline int query(int l,int r,int c){

    int ans=;

    if(r-l+<=*s){

    for(int i=l;i<=r;i++){

        if(a[i]+lazy[(i-)/s+]>=c)ans++;

    }

        return ans;

    }

    int L=(l-)/s+,R=(r-)/s+;

    for(int i=L+;i<=R-;i++){

    ans+=br[i]+-find(bl[i],br[i]+,lazy[i],c);

    }

    for(int i=l;i<=br[L];i++){

    if(a[i]+lazy[L]>=c)ans++;

    }

    for(int i=bl[R];i<=r;i++){

    if(a[i]+lazy[R]>=c)ans++;

    }

    return ans;

}

int main(){

    n=read();m=read();s=sqrt(n);

    for(int i=;i<=n;i++)a[i]=b[i]=read();

    intoblock();

    for(int i=;i<=m;i++){

    char op[];

    cin>>op;

    int l=read(),r=read(),w=read();

    if(op[]=='M')add(l,r,w);

    else printf("%d\n",query(l,r,w));

    }

    return ;

}