题解 P2801 【教主的魔法】
分块入门题,不错的,建议大家做一做
开始学习
先看一下数列分块入门 2
这道题想让我们求区间[l,r]>=c的个数,然后我们可以看到“数列分块入门 2”是求区间[l,r]<c(忽略平方)的个数,即求c在区间[l,r]的排名。所以我们可以每一次查询c的排名,然后用区间长度减c的排名就可以达到答案了呢QAQ。
那么题目就被转移成了求区间[l,r]中c的排名
一看题目,咦,求[l,r]区间c的排名,马上就可以想到平衡树啦,可是平衡树这么难写,而且还不支持区间加,那怎么办? 分块!
首先,我们一个一个步骤分着来看:
(1).操作 操作和普通分块的操作一模一样,类个tag,我们先不理。
(2).查询 如果一般查询排名,我们可以将这一段数截出来,排序,然后二分查找(lower_bound操作)就ok了。那我们可不可以先将每一个块的数都排好序,然后查找每一个块的时候就直接二分呢? ofcause,当然可以!那零零散散的剩下的非整块的数就暴力找就可以了。
(3).排序操作 我们查询需要先排序,那么我们排序也成为了一个操作。我们另外开一个数组ve[i](ve[i]为vector,i表示第i个块)来存每个块排序的数据(这里我使用了vector)。到时候查询的时候直接二分ve[i]就可以了。那每一次零散修改操作都会修改数值(整块操作还是直接加tag),那么我们要对有修改数值的那个块重新进行块内排序。(额,好慢啊)
所以这道题的复杂度是O(msqrt(nlogn))
那么我们这个分块的基本思想就是这样了。还有很多小细节可以看看我的代码:
//P2801 教主的魔法
#include<map>
#include<set>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
int n,m,blo;
int v[1000005],bl[1000005],atag[1000005];
vector<int>ve[10001];
void reset(int x) //排序操作
{
ve[x].clear(); //先初始化
for(int i=(x-1)*blo+1;i<=min(x*blo,n);i++)
ve[x].push_back(v[i]); //将整个块的数值重新压进vector里
sort(ve[x].begin(),ve[x].end()); //块内排序
}
ll read() //快读
{
ll x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void add(int a,int b,int c) //修改操作
{
for(int i=a;i<=min(bl[a]*blo,b);i++) v[i]+=c; //零散操作
reset(bl[a]); //重新块内排序
if(bl[a]!=bl[b])
{
for(int i=(bl[b]-1)*blo+1;i<=b;i++) v[i]+=c;//零散操作
reset(bl[b]); //重新块内排序
}
for(int i=bl[a]+1;i<=bl[b]-1;i++)
atag[i]+=c;
}
int query(int a,int b,int c) //查询操作
{
int ans=0;//排名
for(int i=a;i<=min(bl[a]*blo,b);i++)//零散暴力查询
{
if(v[i]+atag[bl[a]]<c)
ans++;
}
if(bl[a]!=bl[b])
{
for(int i=(bl[b]-1)*blo+1;i<=b;i++) //零散暴力查询
{
if(v[i]+atag[bl[b]]<c)
ans++;
}
}
for(int i=bl[a]+1;i<=bl[b]-1;i++) //整块查询
{
int x=c-atag[i]; //注意:这里要先剪去这个块的atag,因为这整个块的所以数值都应该加了atag,所以我们二分找的数也要先剪atag
ans+=lower_bound(ve[i].begin(),ve[i].end(),x)-ve[i].begin(); //注意:这里lower_bound操作返回的是地址,所以要想知道块内排名,要减块头的地址
} //不会用lower_bound可以直接写二分
return ans;
}
int main()
{
n=read();m=read();blo=sqrt(n);
for(int i=1;i<=n;i++)v[i]=read();
for(int i=1;i<=n;i++)bl[i]=(i-1)/blo+1,ve[bl[i]].push_back(v[i]); //初始化ve[i]数组
for(int i=1;i<=bl[n];i++)
sort(ve[i].begin(),ve[i].end()); //排序
for(int i=1;i<=m;i++)
{
char f;cin>>f;
int a=read(),b=read(),c=read();
if(f=='M')add(a,b,c);
if(f=='A')printf("%d\n",b-a+1-query(a,b,c));
}
return 0;
}
广告
如果需要更系统的学习分块可以看我的blog,内容更详细,分解更到位qwq
分块学习blog
题解 P2801 【教主的魔法】的更多相关文章
- P2801 教主的魔法(分块)
P2801 教主的魔法 区间加法,区间查询 显然就是分块辣 维护一个按块排好序的数组. 每次修改依然是整块打标记,零散块暴力.蓝后对零散块重新排序. 询问时整块二分,零散块暴力就好辣 注意细节挺多和边 ...
- 洛谷 P2801 教主的魔法 解题报告
P2801 教主的魔法 题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.--.N. ...
- 洛谷——P2801 教主的魔法(线段树or分块)
P2801 教主的魔法 (1) 若第一个字母为“M”,则紧接着有三个数字L.R.W.表示对闭区间 [L, R] 内所有英雄的身高加上W. (2) 若第一个字母为“A”,则紧接着有三个数字L.R.C.询 ...
- P2801 教主的魔法 (线段树)
题目 P2801 教主的魔法 解析 成天做水题 线段树,第一问区间加很简单 第二问可以维护一个区间最大值和一个区间最小值,若C小于等于区间最小值,就加上区间长度,若C大于区间最大值,就加0 ps:求教 ...
- 洛谷P2801 教主的魔法 [分块,二分答案]
题目传送门 教主的魔法 题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. ...
- luogu P2801 教主的魔法
题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. 每个人的身高一开始都是 ...
- 洛谷 P2801 教主的魔法
题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. 每个人的身高一开始都是 ...
- BZOJ——3343: 教主的魔法 || 洛谷—— P2801 教主的魔法
http://www.lydsy.com/JudgeOnline/problem.php?id=3343 || https://www.luogu.org/problem/show?pid=280 ...
- 洛谷 P2801 教主的魔法 题解
题面 刚看到这道题的时候用了个树状数组优化前缀和差分的常数优化竟然AC了?(这数据也太水了吧~) 本人做的第一道分块题,调试了好久好久,最后竟然没想到二分上还会出错!(一定要注意)仅此纪念: #inc ...
- 洛谷P2801 教主的魔法 分块
正解:分块 解题报告: 哇之前的坑还没填完就又写新博客? 不管不管,之前欠的两三篇题解大概圣诞节之前会再仔细想想然后重新写下题解趴,确实还挺难的感觉没有很好的理解呢QAQ还是太囫囵吞枣不求甚解了,这样 ...
随机推荐
- idea maven项目使用过程中遇到的问题
1. Error:Cannot build Artifact :war exploded because it is included into a circular depency 参考: http ...
- 临时解决执行 Composer Install 返回 Killed 的问题
昨天在 Linux 服务器上部署 PHP 项目时遇到了一个问题,系统为 Centos 7 ,1 核 1G 的配置.通过 Git 拉取代码后,由于是基于 Laravel 框架的项目,所以需要使用 Com ...
- Python版GPA计算器
最近在网申投简历时遇到一个需要计算GPA的问题,想起自己在上学时写的Excel公式版GPA计算器略显low,而且操作也比较复杂,于是一时兴起,写了个Python版的,在此分享给大家,也算是跟上时代的浪 ...
- [读书]The Man Who Solved the Market
出乎个人意料的是,西蒙斯是从FICC类品种起步的,包括量化投资方法获得第一次重大突破也是在FICC品种上. FICC市场的深度不够,所以文艺复兴科技实现规模扩张是股票策略成功之后的事情,很靠后. 虽然 ...
- 使用maven搭建web项目
在pom.xml中添加java ee相关的三个依赖包:<scope> jar的有效范围 provided 表示编译期生效,不会打包发布到 tomcat 中 <properties&g ...
- Thread的join方法
一个线程在执行的过程中,可能调用另一个线程,前者可以称为调用线程,后者成为被调用线程. Thread.Join方法的使用场景:调用线程挂起,等待被调用线程执行完毕后,继续执行. 如下案列: 当NewT ...
- Navigating to current location ("/") is not allowed
main.js import Router from 'vue-router' // 这个是为了避免一个报错 const originalPush = Router.prototype.push; R ...
- day5-2正则表达式
正则表达式: 正则表达式对象的创建 1,构造函数 var pattern =new RegExp("正则表达式","修饰符") var pattern =new ...
- 二 Spring的IOC入门,环境搭建,Spring测试类
IOC:inversion of Control 控制反转,Spring框架的核心.削减计算机程序的耦合问题,把对象(例如JDBC)的创建权交给Spring. IOC的两种类型: 依赖注入: 依赖查 ...
- ES5中的this
参考资料:>>> this的指向 在 ES5 中,其实 this 的指向,始终坚持一个原理: this 永远指向最后调用它的那个对象 下面我们来看一个最简单的例子:(例子均来自参考资 ...