codeforces 551E GukiZ and GukiZiana 分块

2024-06-12 11:32

本文主要是介绍codeforces 551E GukiZ and GukiZiana 分块,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意:

给你一个数列,有两种操作,第一种操作将l~r区间加上val,第二种操作是询问整个序列,是否存在val,若存在,则找出相差最远的两个val下标的差值。

没有则输出-1。


思路:

分块。

人生的第一道分块题,其实分块的概念顶容易理解的。

就是把序列分为多个长度为sqrt(n)的块。每次更新的区间更新的时候,如果是完整覆盖某个块,则直接对记录数组进行更新。对于两端的不完全覆盖的块,直接暴力更新。

可以参考这篇文章    http://www.cnblogs.com/sweetsc/archive/2012/08/15/2639395.html

这题一开始是用的set来搞,各种优化各种搞才跑8s过了。后来改用vector来搞,2s就过了。


说下vector的做法:

首先把序列分为每块长度为sqrt(n)块,添加到vec数组中。然后根据值大小还有下标排序。

操作1:对rec数组(记录数组)进行更新。对于不完整的块,则直接暴力更新,再重建对应的vec数组。(注意l,r在同一个块的情况)

操作2:运用upper_bound找所要的位置。


code:

#include <bits/stdc++.h>
using namespace std;const int N = 5e5+5;
const int INF = 0x3f3f3f3f;
const int INT = 1e9;
//typedef long long LL;struct PP {int val;int id;bool operator < (const PP &cmp) const {if(val == cmp.val) return id < cmp.id;return val < cmp.val;}
};int n, q;
int a[N], rec[N];
vector <PP> vec[1000];int preDeal() {int cnt = -1;int len = (int)sqrt(n);for(int i = 0;i < n; i++) {if(i%len==0) cnt++;vec[cnt].push_back((PP){a[i], i});}for(int i = 0;i <= cnt; i++)sort(vec[i].begin(), vec[i].end());return cnt;
}int main() {scanf("%d%d", &n, &q);for(int i = 0;i < n; i++)scanf("%d" ,&a[i]);int len = (int)sqrt(n);int cnt = preDeal()+1;while(q--) {int op;scanf("%d", &op);if(op == 1) {int l, r, val;scanf("%d%d%d", &l, &r, &val);if(val == 0) continue;l--, r--;int bl = l/len, br = r/len;for(int i = bl+1;i < br; i++) {if(rec[i] > INT) continue;rec[i] += val;}//left && rightvec[bl].clear(), vec[br].clear();if(bl == br) {for(int i = l;i <= r; i++) {if(a[i] > INT) continue;a[i] += val;}for(int i = bl*len; i < (bl+1)*len; i++)vec[bl].push_back((PP){a[i], i});sort(vec[bl].begin(), vec[bl].end());}else {//for(int i = l;i < (bl+1)*len;i++) {if(a[i] > INT) continue;a[i] += val;}for(int i = br*len; i <= r; i++) {if(a[i] > INT) continue;a[i] += val;}for(int i = bl*len; i < (bl+1)*len; i++)vec[bl].push_back((PP){a[i], i});sort(vec[bl].begin(), vec[bl].end());for(int i = br*len;i < (br+1)*len; i++)vec[br].push_back((PP){a[i], i});sort(vec[br].begin(), vec[br].end());}}else {int val;scanf("%d", &val);//lint l=-1, r;for(int i = 0;i < cnt; i++) {int tmp = val-rec[i];int idx = upper_bound(vec[i].begin(), vec[i].end(), (PP){tmp, -1})-vec[i].begin();if(idx == vec[i].size()) continue;if(vec[i][idx].val == tmp) {l = vec[i][idx].id;break;}}if(l == -1) {puts("-1");continue;}for(int i = cnt-1;i >= 0; i--) {int tmp = val-rec[i];int idx = upper_bound(vec[i].begin(), vec[i].end(), (PP){tmp, INF})-vec[i].begin();if(idx != 0) idx--;if(vec[i][idx].val == tmp) {r = vec[i][idx].id;break;}}//cout<<"l = "<<l<<" r = "<<r<<endl;printf("%d\n", r-l);}}return 0;
}



这篇关于codeforces 551E GukiZ and GukiZiana 分块的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1054093

相关文章

Codeforces Round #240 (Div. 2) E分治算法探究1

Codeforces Round #240 (Div. 2) E  http://codeforces.com/contest/415/problem/E 2^n个数,每次操作将其分成2^q份,对于每一份内部的数进行翻转(逆序),每次操作完后输出操作后新序列的逆序对数。 图一:  划分子问题。 图二: 分而治之,=>  合并 。 图三: 回溯:

Codeforces Round #261 (Div. 2)小记

A  XX注意最后输出满足条件,我也不知道为什么写的这么长。 #define X first#define Y secondvector<pair<int , int> > a ;int can(pair<int , int> c){return -1000 <= c.X && c.X <= 1000&& -1000 <= c.Y && c.Y <= 1000 ;}int m

Codeforces Beta Round #47 C凸包 (最终写法)

题意慢慢看。 typedef long long LL ;int cmp(double x){if(fabs(x) < 1e-8) return 0 ;return x > 0 ? 1 : -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}point op

Codeforces Round #113 (Div. 2) B 判断多边形是否在凸包内

题目点击打开链接 凸多边形A, 多边形B, 判断B是否严格在A内。  注意AB有重点 。  将A,B上的点合在一起求凸包,如果凸包上的点是B的某个点,则B肯定不在A内。 或者说B上的某点在凸包的边上则也说明B不严格在A里面。 这个处理有个巧妙的方法,只需在求凸包的时候, <=  改成< 也就是说凸包一条边上的所有点都重复点都记录在凸包里面了。 另外不能去重点。 int

Codeforces 482B 线段树

求是否存在这样的n个数; m次操作,每次操作就是三个数 l ,r,val          a[l] & a[l+1] &......&a[r] = val 就是区间l---r上的与的值为val 。 也就是意味着区间[L , R] 每个数要执行 | val 操作  最后判断  a[l] & a[l+1] &......&a[r] 是否= val import ja

Codeforces Round 971 (Div. 4) (A~G1)

A、B题太简单,不做解释 C 对于 x y 两个方向,每一个方向至少需要 x / k 向上取整的步数,取最大值。 由于 x 方向先移动,假如 x 方向需要的步数多于 y 方向的步数,那么最后 y 方向的那一步就不需要了,答案减 1 代码 #include <iostream>#include <algorithm>#include <vector>#include <string>

Codeforces#295(Div.2)A、B(模拟+BFS)

解题报告链接:点击打开链接 C. 题目链接:点击打开链接 解题思路: 对于给定的字符串,取出现次数最多的字母(可以同时有多个)。由这些字母组成长度为n的字符串,求有多少种组合。最后用数学知识即可。 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <climits>

Codeforces Round #281 (Div. 2)A(构造+暴力模拟)

题目链接:http://codeforces.com/problemset/problem/493/A 解题思路: 暴力的判断,分三种情况去判断即可。注意如果之前已经被罚下场后,那么在后面的罚下情况不应该算在输出结果内。 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <co

Codeforces Round #182 (Div. 2)A(水题)

题目链接:http://codeforces.com/contest/302/problem/A 解题思路: 只要通过重新排列使区间内和为0即是1,否则是0. 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <complex>#include <cstdio>#inc

Codeforces Round #233 (Div. 2)A(构造)

题目链接:http://codeforces.com/contest/399/problem/A 解题思路: 构造出来即可,考虑p-k和p+k两个边界分别于1和n作比较,对左右符号特殊处理。 完整代码: #include <algorithm>#include <iostream>#include <cstring>#include <complex>#include