【bzoj3343】教主的魔法

来源:互联网 发布:怎么能订阅知乎周刊 编辑:程序博客网 时间:2024/03/28 23:26

*题目描述:
教主最近学会了一种神奇的魔法,能够使人长高。于是他准备演示给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的英雄数。
*样例输入:
5 3
1 2 3 4 5
A 1 5 4
M 3 5 1
A 1 5 4
*样例输出:
2
3
*提示:
【输入输出样例说明】
原先5个英雄身高为1、2、3、4、5,此时[1, 5]间有2个英雄的身高大于等于4。教主施法后变为1、2、4、5、6,此时[1, 5]间有3个英雄的身高大于等于4。
【数据范围】
对30%的数据,N≤1000,Q≤1000。
对100%的数据,N≤1000000,Q≤3000,1≤W≤1000,1≤C≤1,000,000,000。
*题解:
分块,块内是有序的,然后每次查询的话如果在同一块就暴力查,不在同一块就两边暴力,中间在块内二分(lower_bound大法好~)。区间修改的话和查询差不多,残缺的块就暴力加完重构,完整的块直接打上标记(标记永久化)。
*代码:

#include <cstdio>#include <cstring>#include <algorithm>#include <cmath>#ifdef WIN32    #define LL "%I64d"#else    #define LL "%lld"#endif#ifdef CT    #define debug(...) printf(__VA_ARGS__)    #define setfile() #else    #define debug(...)    #define filename ""    #define setfile() freopen(filename".in", "r", stdin); freopen(filename".out", "w", stdout);#endif#define R register#define getc() (S == T && (T = (S = B) + fread(B, 1, 1 << 15, stdin), S == T) ? EOF : *S++)#define dmax(_a, _b) ((_a) > (_b) ? (_a) : (_b))#define dmin(_a, _b) ((_a) < (_b) ? (_a) : (_b))#define cmax(_a, _b) (_a < (_b) ? _a = (_b) : 0)#define cmin(_a, _b) (_a > (_b) ? _a = (_b) : 0)char B[1 << 15], *S = B, *T = B;inline int FastIn(){    R char ch; R int cnt = 0; R bool minus = 0;    while (ch = getc(), (ch < '0' || ch > '9') && ch != '-') ;    ch == '-' ? minus = 1 : cnt = ch - '0';    while (ch = getc(), ch >= '0' && ch <= '9') cnt = cnt * 10 + ch - '0';    return minus ? -cnt : cnt;}#define maxn 1000010#define maxsize 4010int a[maxn], b[maxn], block, id[maxn], n, q, ans, w, tag[maxsize];inline void build(R int x){    R int l = x * block, r = (x + 1) * block - 1;    cmin(r, n);    for (R int i = l; i <= r; ++i) b[i] = a[i];    std::sort(b + l, b + r + 1);}inline void voi_query(R int l, R int r){    for (R int i = l; i <= r; ++i)        if (a[i] >= w) ++ans;}inline void voi_add(R int l, R int r){    for (R int i = l; i <= r; ++i)        a[i] += w;    build(id[l]);}inline void query(R int x){    R int l = x * block, r = (x + 1) * block - 1;    cmin(r, n);    ans += (r - l + 1) - (std::lower_bound(b + l, b + r + 1, w - tag[x]) - b - l);}int main(){//  setfile();    n = FastIn(), q = FastIn();    block = sqrt(n);//  printf("sizes = %d\n", block );    for (R int i = 1; i <= n; ++i) a[i] = FastIn();    for (R int i = 1; i <= n; ++i) id[i] = i / block;    for (R int i = 0; i <= n / block; ++i) build(i);    for (R int i = 1; i <= q; ++i)    {        R char opt = getc();        while (opt < 'A' || opt > 'Z') opt = getc();        R int l = FastIn(), r = FastIn(); w = FastIn();        if (opt == 'A')        {            ans = 0;            if (id[l] == id[r]) voi_query(l, r);            else            {                voi_query(l, (id[l] + 1) * block - 1);                voi_query(id[r] * block, r);                for (R int i = id[l] + 1; i < id[r]; ++i)                    query(i);            }            printf("%d\n", ans );        }        else        {            if (id[l] == id[r]) voi_add(l, r);            else            {                voi_add(l, (id[l] + 1) * block - 1);                voi_add(id[r] * block, r);                for (R int i = id[l] + 1; i < id[r]; ++i)                    tag[i] += w;            }        }    }    return 0;}
0 0
原创粉丝点击