数列分块入门1~9 loj6277~6285
一
给出一个长为 \(n\) 的数列，以及 \(n\) 个操作，操作涉及区间加法，单点查值。
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n, a[50005], opt, uu, vv, ww, tag[305], blc, bel[50005];
void add(int uu, int vv, int ww){
int p=bel[uu], q=bel[vv];
if(p==q)
for(int i=uu; i<=vv; i++)
a[i] += ww;
else{
for(int i=p+1; i<=q-1; i++) tag[i] += ww;
for(int i=uu; i<=bel[uu]*blc; i++) a[i] += ww;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) a[i] += ww;
}
}
int main(){
cin>>n;
blc = sqrt(n);
for(int i=1; i<=n; i++) scanf("%d", &a[i]);
for(int i=1; i<=n; i++) bel[i] = (i - 1) / blc + 1;
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) add(uu, vv, ww);
else printf("%d\n", a[vv]+tag[bel[vv]]);
}
return 0;
}
二
给出一个长为 \(n\) 的数列，以及 \(n\) 个操作，操作涉及区间加法，询问区间内小于某个值 \(x\) 的元素个数。
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
using namespace std;
int n, blc, a[50005], tag[50005], bel[50005], opt, uu, vv, ww;
vector<int> vec[50005];
void qwq(int tat){
vec[tat].clear();
for(int i=(tat-1)*blc+1; i<=min(tat*blc, n); i++) vec[tat].push_back(a[i]);
sort(vec[tat].begin(), vec[tat].end());
}
void update(int uu, int vv, int ww){
if(bel[uu]==bel[vv]){
for(int i=uu; i<=vv; i++) a[i] += ww;
qwq(bel[uu]);
}
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++) tag[i] += ww;
for(int i=uu; i<=bel[uu]*blc; i++) a[i] += ww;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) a[i] += ww;
qwq(bel[uu]); qwq(bel[vv]);
}
}
int query(int uu, int vv, int ww){
int re=0;
if(bel[uu]==bel[vv]){
for(int i=uu; i<=vv; i++)
if(a[i]+tag[bel[i]]<ww)
re++;
}
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++)
re += lower_bound(vec[i].begin(), vec[i].end(), ww-tag[i]) - vec[i].begin();
for(int i=uu; i<=bel[uu]*blc; i++)
if(a[i]+tag[bel[i]]<ww)
re++;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++)
if(a[i]+tag[bel[i]]<ww)
re++;
}
return re;
}
int main(){
cin>>n;
blc = sqrt(n/200);
for(int i=1; i<=n; i++){
bel[i] = (i - 1) / blc + 1;
scanf("%d", &a[i]);
vec[bel[i]].push_back(a[i]);
}
for(int i=1; i<=n; i=bel[i]*blc+1)
sort(vec[bel[i]].begin(), vec[bel[i]].end());
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) update(uu, vv, ww);
else printf("%d\n", query(uu, vv, ww*ww));
}
return 0;
}
三
给出一个长为n的数列，以及n个操作，操作涉及区间加法，询问区间内小于某个值x的前驱（比其小的最大元素）。
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
using namespace std;
int n, blc, a[100005], bel[100005], tag[100005], opt, uu, vv, ww;
vector<int> vec[100005];
void qwq(int u){
vec[u].clear();
for(int i=(u-1)*blc+1; i<=min(n, blc*u); i++)
vec[u].push_back(a[i]);
sort(vec[u].begin(), vec[u].end());
}
void update(int uu, int vv, int ww){
if(bel[uu]==bel[vv]){
for(int i=uu; i<=vv; i++)
a[i] += ww;
qwq(bel[uu]);
}
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++) tag[i] += ww;
for(int i=uu; i<=bel[uu]*blc; i++) a[i] += ww;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) a[i] += ww;
qwq(bel[uu]); qwq(bel[vv]);
}
}
int query(int uu, int vv, int ww){
int re=0xffffffff;
if(bel[uu]==bel[vv]){
for(int i=uu; i<=vv; i++)
if(a[i]+tag[bel[i]]<ww)
re = max(re, a[i]+tag[bel[i]]);
}
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++){
int pos=lower_bound(vec[i].begin(), vec[i].end(), ww-tag[i])-vec[i].begin();
if(pos) re = max(re, vec[i][pos-1]+tag[i]);
}
for(int i=uu; i<=bel[uu]*blc; i++)
if(a[i]+tag[bel[i]]<ww)
re = max(re, a[i]+tag[bel[i]]);
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++)
if(a[i]+tag[bel[i]]<ww)
re = max(re, a[i]+tag[bel[i]]);
}
if(re==0xffffffff) return -1;
else return re;
}
int main(){
cin>>n;
blc = sqrt(n*log(n)/log(2));
for(int i=1; i<=n; i++){
scanf("%d", &a[i]);
bel[i] = (i - 1) / blc + 1;
vec[bel[i]].push_back(a[i]);
}
for(int i=1; i<=n; i=bel[i]*blc+1)
sort(vec[bel[i]].begin(), vec[bel[i]].end());
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) update(uu, vv, ww);
else printf("%d\n", query(uu, vv, ww));
}
return 0;
}
四
给出一个长为n的数列，以及n个操作，操作涉及区间加法，区间求和。
五
给出一个长为n的数列，以及n个操作，操作涉及区间开方，区间求和。
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n, blc, bel[50005], opt, uu, vv, ww, tag[50005], sum[50005], a[50005];
int qwq(int uu, int vv){
for(int i=uu; i<=vv; i++){
sum[bel[i]] -= a[i];
a[i] = sqrt(a[i]);
sum[bel[i]] += a[i];
}
return sum[bel[uu]];
}
void qaq(int u){
if(tag[u]) return ;
int re=qwq((u-1)*blc+1, u*blc);
if(re>blc) tag[u] = false;
else tag[u] = true;
}
void update(int uu, int vv){
if(bel[uu]==bel[vv])
qwq(uu, vv);
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++) qaq(i);
qwq(uu, bel[uu]*blc);
qwq((bel[vv]-1)*blc+1, vv);
}
}
int query(int uu, int vv){
int re=0;
if(bel[uu]==bel[vv])
for(int i=uu; i<=vv; i++)
re += a[i];
else{
for(int i=bel[uu]+1; i<=bel[vv]-1; i++) re += sum[i];
for(int i=uu; i<=bel[uu]*blc; i++) re += a[i];
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) re += a[i];
}
return re;
}
int main(){
cin>>n;
blc = sqrt(n);
for(int i=1; i<=n; i++){
scanf("%d", &a[i]);
bel[i] = (i - 1) / blc + 1;
sum[bel[i]] += a[i];
}
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) update(uu, vv);
else printf("%d\n", query(uu, vv));
}
return 0;
}
六
给出一个长为n的数列，以及n个操作，操作涉及单点插入，单点询问，数据随机生成。
#include <algorithm>
#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
using namespace std;
typedef pair<int,int> par;
int n, blc, opt, uu, vv, ww, sta[200005], din, m;
vector<int> vec[200005];
par query(int uu){
int x=1;
while(uu>vec[x].size()){
uu -= vec[x].size();
x++;
}
return make_pair(x, uu-1);
}
void rebuild(){
din = 0;
for(int i=1; i<=m; i++){
for(int j=0; j<vec[i].size(); j++)
sta[++din] = vec[i][j];
vec[i].clear();
}
blc = sqrt(din);
for(int i=1; i<=din; i++)
vec[(i-1)/blc+1].push_back(sta[i]);
}
void update(int uu, int vv){
par re=query(uu);
vec[re.first].insert(vec[re.first].begin()+re.second, vv);
if(vec[re.first].size()>blc+blc) rebuild();
}
int main(){
cin>>n;
blc = sqrt(n);
m = (n - 1) / blc + 1;
for(int i=1; i<=n; i++){
scanf("%d", &uu);
vec[(i-1)/blc+1].push_back(uu);
}
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) update(uu, vv);
else{
par re=query(vv);
printf("%d\n", vec[re.first][re.second]);
}
}
return 0;
}
七
给出一个长为n的数列，以及n个操作，操作涉及区间乘法，区间加法，单点询问。
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n, a[100005], bel[100005], add[100005], mul[100005], blc;
int opt, uu, vv, ww;
const int mod=10007;
void pushDown(int x){
for(int i=(x-1)*blc+1; i<=min(n, x*blc); i++)
a[i] = (a[i] * mul[x] + add[x]) % mod;
mul[x] = 1; add[x] = 0;
}
void updAdd(int uu, int vv, int ww){
pushDown(bel[uu]);
if(bel[uu]==bel[vv])
for(int i=uu; i<=vv; i++)
a[i] = (a[i] + ww) % mod;
else{
pushDown(bel[vv]);
for(int i=bel[uu]+1; i<=bel[vv]-1; i++) add[i] = (add[i] + ww) % mod;
for(int i=uu; i<=bel[uu]*blc; i++) a[i] = (a[i] + ww) % mod;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) a[i] = (a[i] + ww) % mod;
}
}
void updMul(int uu, int vv, int ww){
pushDown(bel[uu]);
if(bel[uu]==bel[vv])
for(int i=uu; i<=vv; i++)
a[i] = (a[i] * ww) % mod;
else{
pushDown(bel[vv]);
for(int i=bel[uu]+1; i<=bel[vv]-1; i++){
mul[i] = (mul[i] * ww) % mod;
add[i] = (add[i] * ww) % mod;
}
for(int i=uu; i<=bel[uu]*blc; i++) a[i] = (a[i] * ww) % mod;
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++) a[i] = (a[i] * ww) % mod;
}
}
int main(){
cin>>n;
blc = sqrt(n);
for(int i=1; i<=n; i++){
mul[i] = 1;
scanf("%d", &a[i]);
bel[i] = (i - 1) / blc + 1;
}
for(int i=1; i<=n; i++){
scanf("%d %d %d %d", &opt, &uu, &vv, &ww);
if(!opt) updAdd(uu, vv, ww);
else if(opt==1) updMul(uu, vv, ww);
else printf("%d\n", (a[vv]*mul[bel[vv]]+add[bel[vv]])%mod);
}
return 0;
}
八
给出一个长为n的数列，以及n个操作，操作涉及区间询问等于一个数c的元素，并将这个区间的所有元素改为c。
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int n, a[100005], tag[100005], bel[100005], blc, uu, vv, ww;
void pushDown(int x){
if(tag[x]==-1) return ;
for(int i=(x-1)*blc+1; i<=min(n, x*blc); i++)
a[i] = tag[x];
tag[x] = -1;
}
int update(int uu, int vv, int ww){
int re=0;
if(tag[uu]) pushDown(bel[uu]);
if(bel[uu]==bel[vv])
for(int i=uu; i<=vv; i++){
if(a[i]==ww)
re++;
a[i] = ww;
}
else{
if(tag[vv]) pushDown(bel[vv]);
for(int i=bel[uu]+1; i<=bel[vv]-1; i++){
if(tag[i]!=-1)
re += tag[i]==ww?blc:0;
else
for(int j=(i-1)*blc+1; j<=i*blc; j++)
re += a[j]==ww;
tag[i] = ww;
}
for(int i=uu; i<=bel[uu]*blc; i++){
re += a[i]==ww;
a[i] = ww;
}
for(int i=(bel[vv]-1)*blc+1; i<=vv; i++){
re += a[i]==ww;
a[i] = ww;
}
}
return re;
}
int main(){
cin>>n;
blc = sqrt(n);
for(int i=1; i<=n; i++){
tag[i] = -1;
scanf("%d", &a[i]);
bel[i] = (i - 1) / blc + 1;
}
for(int i=1; i<=n; i++){
scanf("%d %d %d", &uu, &vv, &ww);
printf("%d\n", update(uu, vv, ww));
}
return 0;
}
九
区间众数
参考luogu4168 蒲公英
数列分块入门1~9 loj6277~6285的更多相关文章
- LOJ6277~6285 数列分块入门
Portals 分块需注意的问题 数组大小应为,因为最后一个块可能会超出的范围. 当操作的区间在一个块内时,要特判成暴力修改. 要清楚什么时候应该+tag[t] 数列分块入门 1 给出一个长为的数列, ...
- LOJ #6285. 数列分块入门 9-分块(查询区间的最小众数)
#6285. 数列分块入门 9 内存限制:256 MiB时间限制:1500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 2 题目描述 给 ...
- 数列分块入门九题（三）：LOJ6283~6285
Preface 最后一题我一直觉得用莫队是最好的. 数列分块入门 7--区间乘法,区间加法,单点询问 还是很简单的吧,比起数列分块入门 7就多了个区间乘. 类似于线段树,由于乘法的优先级高于加法,因此 ...
- 数列分块入门九题（一）：LOJ6277~6279
Preface 分块,一个神奇的暴力算法.可以把很多\(O(n^2)\)的数据结构题的暴力优化到常数极小的\(O(n\sqrt n)\).当一些毒瘤题无法用线段树,主席树,平衡树,树状数组...... ...
- [Loj] 数列分块入门 1 - 9
数列分块入门 1 https://loj.ac/problem/6277 区间加 + 单点查询 #include <iostream> #include <cstdio> #i ...
- 数列分块入门九题（二）：LOJ6280~6282
Preface 个人感觉这中间的三题是最水的没有之一 数列分块入门 4--区间加法,区间求和 这个也是很多数据结构完爆的题目线段树入门题,但是练分块我们就要写吗 修改还是与之前类似,只不过我们要维护每 ...
- LOJ6285 数列分块入门9（分块）
昨天对着代码看了一晚上 然后今天终于在loj上过了 数列分块入门9题撒花★,°:.☆(￣▽￣)/$:.°★ . 然后相当玄学 块的大小调成\(\sqrt{n}\)会TLE,改成150就过了 啧 然后就 ...
- LOJ 6277：数列分块入门 1（分块入门）
#6277. 数列分块入门 1 内存限制:256 MiB时间限制:100 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计讨论 3 测试数据 题目描述 给出一 ...
- LOJ #6284. 数列分块入门 8-分块(区间查询等于一个数c的元素，并将这个区间的所有元素改为c)
#6284. 数列分块入门 8 内存限制:256 MiB时间限制:500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: hzwer 提交提交记录统计测试数据讨论 2 题目描述 给出 ...
随机推荐
- Qt文本读写之一：输入输出设备和文件操作
一.输入输出设备 QIODevice类是Qt中所有I/O设备的基础接口类,为诸如QFile.QBuffer和 QTcpSocket等支持读/写数据块的设备提供了一个抽象接口.QIODevice类是抽象 ...
- A Refining Company LightOJ - 1036
A Refining Company LightOJ - 1036 描述好长啊... 题意:在m*n的矩阵上,每一格摆一个向上或者向左的传送带(不能同时摆,只能摆一个).同时,每一格有两种物资Uran ...
- jmeter（十九）调试工具Debug Sampler
一.Debug Sampler介绍: 使用Jmeter开发脚本时,难免需要调试,这时可以使用Jmeter的Debug Sampler,它有三个选项:JMeter properties,JMeter v ...
- 机器学习概念之特征选择（Feature selection）之VectorSlicer算法介绍
不多说,直接上干货! VectorSlicer 算法介绍: VectorSlicer是一个转换器,输入特征向量,输出原始特征向量子集.VectorSlicer接收带有特定索引的向量列,通过对这些索引的 ...
- Smart 组件 vs Dumb 组件
大家已经知道,只会接受 props 并且渲染确定结果的组件我们把它叫做 Dumb 组件,这种组件只关心一件事情 —— 根据 props 进行渲染. Dumb 组件最好不要依赖除了 React.js 和 ...
- 【转】码云source tree 提交超过100m 为什么大文件推不上去
码云source tree 提交超过100m 为什么大文件推不上去 2017年01月12日 16:50:51 阅读数:7634 git -c diff.mnemonicprefix=false -c ...
- poj2184 Cow Exhibition
思路: dp+滚动数组. 类似01背包. 实现: #include <iostream> #include <cstdio> #include <algorithm> ...
- 使用vbScript 链接SQLserver数据库和基础操作
使用vbs链接SQLserver数据库 数据库的创建.设计使用 management studio完成 1.本地链接数据库 set oCon = server.createObject("a ...
- java 获取ip地址
1.使用WIFI 首先设置用户权限 <uses-permission android:name="android.permission.ACCESS_WIFI_STATE"& ...
- 毕业设计：HomeFragment（一）
一.主要思路 主要是通过ListView实现. 考虑到以后会添加长按修改功能,所以好几个地方都是用的FramLayout,而且CheckBox初始状态是被隐藏的.给ListView添加OnItemCl ...