【分块】教主的魔法 @洛谷P2801/upcexam3138
时间限制: 1 Sec 内存限制: 128 MB
题目描述
教主最近学会了一种神奇的魔法,能够使人长高。于是他准备演示给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
第一次写分块题
将数据分成sqrt(n)块,对每个块分别排序。
对于修改操作,区间内完整的块修改其add数组的值,不完整的块暴力修改在区间内的部分,然后对整个块排序。
对于查询操作,区间内完整的块二分查找C-add[i]在区间内的位置,然后用该位置减去块的首位置,不完整的块暴力统计。
#define FILE_PC() freopen("C:\\Users\\hz\\Desktop\\in.txt","r",stdin)
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
using namespace std;
const int maxn = 1000005;
const int maxdiv = 1005;
int a[maxn],b[maxn];
int l[maxdiv],r[maxdiv],add[maxdiv];
int lendiv,n,m;
int main() {
// FILE_PC();
scanf("%d%d",&n,&m);
for(int i=1; i<=n; i++) {
scanf("%d",a+i);
b[i] = a[i];
}
lendiv = (int)sqrt(n); //分块长度
int cntdiv = 0,flag = 1;
for(int i=1; i<=n; i++) {
if(cntdiv==0)l[++cntdiv] = i;//块的编号从1开始
if(!flag) {
l[++cntdiv] = i;
flag = 1;
}
if(i%lendiv==0) {
r[cntdiv] = i;
flag = 0;
}
}
r[cntdiv] = n; //最后一块的右区间
for(int i=1; i<=cntdiv; i++) {
sort(b+l[i],b+r[i]+1);
}//每个块里面的元素排序
for(int ca=0; ca<m; ca++) {
char s[20];
int ll,rr,cc;
scanf("%s%d%d%d",s,&ll,&rr,&cc);
if(s[0]=='M') {
int ld = (int)(lower_bound(l+1,l+cntdiv+1,ll)-l);
if(l[ld]!=ll) {
for(int i=ll; i<min(l[ld],rr+1); i++) {
a[i]+=cc;
}
for(int i=l[ld-1]; i<=r[ld-1]; i++) {
b[i] = a[i];
}
sort(b+l[ld-1],b+r[ld-1]+1);//重新对这个块排序
}
if(l[ld]<rr+1) {
int rd = (int)(lower_bound(r+1,r+cntdiv+1,rr)-r);
if(rr!=r[rd]) {
for(int i=max(ll,r[rd-1]+1); i<=rr; i++) {
a[i] += cc;
}
for(int i=l[rd]; i<=r[rd]; i++) {
b[i] = a[i];
}
sort(b+l[rd],b+r[rd]+1);//同上
}
for(int i=ld; i<=((rr==r[rd])?rd:rd-1); i++) { //完整的块
add[i]+=cc;
}
}
} else {
int ans = 0;
int ld = (int)(lower_bound(l+1,l+cntdiv+1,ll)-l);//给定的左区间落在哪一块
if(l[ld]!=ll) {
for(int i=ll; i<min(l[ld],rr+1); i++) { //暴力统计不完整的块
if(a[i]+add[ld-1]>=cc)ans++;
}
}
if(l[ld]<rr+1) {
int rd = (int)(lower_bound(r+1,r+cntdiv+1,rr)-r);//给定的右区间落在哪一块
if(rr!=r[rd]) {
for(int i=max(ll,r[rd-1]+1); i<=rr; i++) { //同上
if(a[i]+add[rd]>=cc)ans++;
}
}
for(int i=ld; i<=((rr==r[rd])?rd:rd-1); i++) { //统计每一个完整的块
ans += (int)(b+r[i]+1-lower_bound(b+l[i],b+r[i]+1,cc-add[i]));
}
}
printf("%d\n",ans);
}
}
return 0;
}
【分块】教主的魔法 @洛谷P2801/upcexam3138的更多相关文章
- BZOJ——3343: 教主的魔法 || 洛谷—— P2801 教主的魔法
http://www.lydsy.com/JudgeOnline/problem.php?id=3343 || https://www.luogu.org/problem/show?pid=280 ...
- 洛谷P2801 教主的魔法 [分块,二分答案]
题目传送门 教主的魔法 题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. ...
- 洛谷——P2801 教主的魔法(线段树or分块)
P2801 教主的魔法 (1) 若第一个字母为“M”,则紧接着有三个数字L.R.W.表示对闭区间 [L, R] 内所有英雄的身高加上W. (2) 若第一个字母为“A”,则紧接着有三个数字L.R.C.询 ...
- 洛谷 P2801 教主的魔法 解题报告
P2801 教主的魔法 题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.--.N. ...
- 洛谷 P2801 教主的魔法
题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. 每个人的身高一开始都是 ...
- 【算法】分块——教主的魔法&不勤劳的图书管理员
由不勤劳的图书管理员带入了分块的坑,深深地被其暴力与优雅所征服.分块的实质就是将暴力块状封装起来,一整块的部分就一整块处理,零碎的部分就怎么暴力怎么来.因为分块大小的原因,限制了零碎部分数据的数量级, ...
- 洛谷P2801 教主的魔法 分块
正解:分块 解题报告: 哇之前的坑还没填完就又写新博客? 不管不管,之前欠的两三篇题解大概圣诞节之前会再仔细想想然后重新写下题解趴,确实还挺难的感觉没有很好的理解呢QAQ还是太囫囵吞枣不求甚解了,这样 ...
- [洛谷P2801]教主的魔法
题目大意:有$n$个数,$q$个操作.两种操作: $M\;l\;r\;w:$把$[l,r]$所有数加上$w$ $A\;l\;r\;c:$查询$[l,r]$内大于等于$c$的元素的个数. 题解:分块,对 ...
- 洛谷 P2801 教主的魔法 题解
题面 刚看到这道题的时候用了个树状数组优化前缀和差分的常数优化竟然AC了?(这数据也太水了吧~) 本人做的第一道分块题,调试了好久好久,最后竟然没想到二分上还会出错!(一定要注意)仅此纪念: #inc ...
随机推荐
- 详解 Java 中的三种代理模式
代理模式 代理(Proxy)是一种设计模式,提供了对目标对象另外的访问方式;即通过代理对象访问目标对象.这样做的好处是:可以在目标对象实现的基础上,增强额外的功能操作,即扩展目标对象的功能. 这里使用 ...
- P1026 统计单词个数 区间dp
题目描述 给出一个长度不超过200200的由小写英文字母组成的字母串(约定;该字串以每行2020个字母的方式输入,且保证每行一定为2020个).要求将此字母串分成kk份(1<k \le 401& ...
- objectmapper使用
https://www.cnblogs.com/liam1994/p/5629356.html
- Maya mayapy.exe 安装 Cython,编译 pyd
Maya mayapy.exe 安装 Cython,编译 pyd 前言 在 Python 2.7 cython cythonize py 编译成 pyd 谈谈那些坑 中最后提到,使用 VCForPy ...
- HDU 1045 Fire Net 【二分图匹配】
<题目链接> 题目大意: 这题意思是给出一张图,图中'X'表示wall,'.'表示空地,可以放置炮台,同一条直线上只能有一个炮台,除非有'X'隔开,问在给出的图中最多能放置多少个炮台. 解 ...
- Python多重继承顺序---C3算法
什么是多重继承C3算法 MRO即 method resolution order (方法解释顺序),主要用于在多继承时判断属性的路径(来自于哪个类). 在python2.2版本中,算法基本思想是根据每 ...
- B. Divisiblity of Differences
B. Divisiblity of Differencestime limit per test1 secondmemory limit per test512 megabytesinputstand ...
- Spring:获取容器中的Bean
某些情况下我们要获取 IOC 容器中指定注解.类型.名字的 Bean 要获取 IOC 容器中指定条件的 Bean 可以通过 ApplicationContext 相应的方法 @Autowired pr ...
- 利用Webpack+React(antd)+ES6+python(flask)实现代码转换
之前的几篇博客是将flask 结合 antd本地化,但是这样使得antd无法按需加载(也不支持ES6的语法),而且在写的过程中还需要把每个组件都用antd对象,这样的做法虽然是实现了antd的本地化, ...
- DRF的视图
DRF的视图 APIView 我们django中写CBV的时候继承的是View,rest_framework继承的是APIView,那么他们两个有什么不同呢~~~ urlpatterns = [ ...