数列分块入门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 题目描述 给出 ...
随机推荐
- 模板——扩展欧几里得算法（求ax+by=gcd的解）
Bryce1010模板 /**** *扩展欧几里得算法 *返回d=gcd(a,b),和对应等式ax+by=d中的x,y */ long long extend_gcd(long long a,long ...
- Subsequence HDU - 3530
Subsequence HDU - 3530 方法:单调队列区间最大最小 错误记录(本地写错)的原因:写成每次试着扩展右端点,却难以正确地处理"在多扩展右端点之后减去多扩展的部分" ...
- 双端队列 HDOJ 3530 Subsequence
题目传送门 题意:问最长子序列,满足区间最大值 - 最小值在[m, k]之间 分析:用双端队列维护最大值和最小值,保存的是位置.当满足条件时,更新最大值. /********************* ...
- DFS HDOJ 5348 Ponds
题目传送门 题意:有一张无向图,度数小于2的点会被去掉,直到全都大于等于2,问连通块顶点数为奇数的权值和为多少 分析:首先DFS把度数小于2的vis掉,第二次DFS把属于同一个连通块的vis掉,检查是 ...
- 1-13Object类之toString方法
Object中的toString方法 SUN在Object类中设计toString方法的目的:返回java对象的字符串表示形式. 在现实的开发过程中,Object中的toString方法就是要被重写的 ...
- javascript简单的表单验证
<html> <head> <title>用户登录</title> <script language="javascript" ...
- [未读]JavaScript高效图形编程
去年买来就一直搁置,因为是js游戏相关,暂时还用不到.
- asp.net 微信登录实现方式
之前我以为做微信登录跟微信公众号有关,后来发现是我想多了.原来微信还有一个叫开放平台的东西,见下图: 我的这个已经生成好了,没有的需要创建一个,https://open.weixin.qq.com/c ...
- AJPFX谈Java 性能优化之基本类型 vs 引用类型
★名词定义 先明确一下什么是“基本类型”,什么是“引用类型”. 简单地说,所谓基本类型就是 Java 语言中如下的8种内置类型: booleancharbyteshortintlongfloatdou ...
- android开发哪些坑需要注意
同一个应用的JNI代码,不要轻易换NDK编译的版本,否则会有很多问题(主要是一些方法实现不一样,并且高版本对代码的检测更严格),比如r8没有问题,但到r9就有问题了,这是个大坑: Android的JN ...