数列分块入门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 题目描述 给出 ...
随机推荐
- 简单实现人工智能：百度aip+tuling123
目录结构: app.py # -*- coding: utf-8 -*- # __author: ward # data: 2018/12/21 # @File: app from flask imp ...
- Wolfycz的娱乐赛题解
现在不会放题解的!比赛完了我会把题解放上来的 祝大家玩的愉快~ 等会,cnblogs不会显示更新时间?我禁赛我自己 UPD:2018.12.15 欢迎大家爆踩标程- painting 我们考虑转化题意 ...
- LightOj 1088 - Points in Segments (二分枚举)
题目链接: http://www.lightoj.com/volume_showproblem.php?problem=1088 题目描述: 给出一个n位数升序排列的数列,然后q个查询,每个查询问指定 ...
- [转]深入浅出WPF（7）——数据的绿色通道，Binding
本文转自:http://liutiemeng.blog.51cto.com/120361/95273 小序: 怎么直接从2蹦到7啦?!啊哦,实在是不好意思,最近实在是太忙了,忙的原因也非常简单——自己 ...
- CSS预处理less基本使用
中文API http://lesscss.cn 变量 @变量名:变量值 @maincolor:#aeeeee; @acolor:#ffffff; @ht200:200px; @ht50:50p ...
- vmware让虚拟机内外网络可互访
以下方法可使主机可以ping通虚拟机,虚拟机也可以ping通主机 首先对虚拟机设置 然后设置虚拟机,假设主机的ip是10.0.0.9,那虚拟机的ip应如下设置: 其中ip地址任意设置一个,但要求跟主机 ...
- mac下iterm2配置安装,通过expact实现保存账号,及通过跳板登陆配置
在参考了几款mac不错的ssh工具外,最终选择使用iterm2.本来打算用FinalShell,安装后发现其icon在访达中根本不现实,而且每次访问还需要输入管理员账号密码,强迫症根本受不了... 官 ...
- JS正则表达式匹配<div><style>标签
测试字符串: <style>v\:* { BEHAVIOR: url(#default#VML) } o\:* { BEHA ...
- vue2.0组件生命周期探讨
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- 遍历NSView下的子视图方法
如何遍历NSView下的子视图呢 for (NSView *aview in [SuperV subviews]) { if([aview isMemberOfClass:[NSButton clas ...