【BZOJ2989】数列(二进制分组,主席树)
【BZOJ2989】数列(二进制分组,主席树)
题面
BZOJ
权限题啊。。。
Description
给定一个长度为n的正整数数列a[i]。
定义2个位置的graze值为两者位置差与数值差的和,即graze(x,y)=|x-y|+|a[x]-a[y]|。
2种操作(k都是正整数):
1.Modify x k:将第x个数的值修改为k。
2.Query x k:询问有几个i满足graze(x,i)<=k。因为可持久化数据结构的流行,询问不仅要考虑当前数列,还要
考虑任意历史版本,即统计任意位置上出现过的任意数值与当前的a[x]的graze值<=k的对数。(某位置多次修改为
同样的数值,按多次统计)
Input
第1行两个整数n,q。分别表示数列长度和操作数。
第2行n个正整数,代表初始数列。
第3--q+2行每行一个操作。
Output
对于每次询问操作,输出一个非负整数表示答案
Sample Input
3 5
2 4 3
Query 2 2
Modify 1 3
Query 2 2
Modify 1 2
Query 1 1
Sample Output
2
3
3
HINT
N<=60000 修改操作数<=40000 询问<=50000 Max{a[i]}含修改<=100000
题解
其实这题并不需要支持强制在线,所以二进制分组被CDQ分治打爆了。
我们对于每个组开主席树,每次询问相当于在每个组内都查询一次求和,
好暴力啊。。。。。
然而复杂度\(O(nlog^2n)\),惨遭\(CDQ\)分治打爆。
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
#define ll long long
#define inf 1000000000
#define MAX 200200
#define MAXNODE MAX<<5
#define py 65000
#define N 200000
#define pi pair<int,int>
#define mp make_pair
#define pb push_back
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
char ch[20];
struct Node{int ls,rs,v;}t[MAXNODE];
int top,St[MAXNODE];bool vis[MAXNODE];
int NewNode(){vis[St[top]]=true;return St[top--];}
int Query(int A,int B,int l,int r,int L,int R)
{
if(L<=l&&r<=R)return t[B].v-t[A].v;
int mid=(l+r)>>1,ret=0;
if(L<=mid)ret+=Query(t[A].ls,t[B].ls,l,mid,L,R);
if(R>mid)ret+=Query(t[A].rs,t[B].rs,mid+1,r,L,R);
return ret;
}
void modify(int &x,int ff,int l,int r,int p)
{
t[x=NewNode()]=t[ff];t[x].v+=1;if(l==r)return;
int mid=(l+r)>>1;
if(p<=mid)modify(t[x].ls,t[ff].ls,l,mid,p);
else modify(t[x].rs,t[ff].rs,mid+1,r,p);
}
void Del(int &x)
{
if(!vis[x])return;
St[++top]=x;vis[x]=false;
Del(t[x].ls);Del(t[x].rs);
t[x].v=0;x=0;
}
vector<pi> v[20];
int rt[20][MAX],Top;
void insert(int x,int y)
{
v[++Top].pb(mp(x,y));
modify(rt[Top][1],rt[Top][0],1,N,y);
while(Top>1&&v[Top-1].size()==v[Top].size())
{
int t1=0,t2=0;vector<pi> t;
for(int i=1;i<=v[Top-1].size();++i)Del(rt[Top-1][i]);
for(int i=1;i<=v[Top].size();++i)Del(rt[Top][i]);
while(t1<v[Top-1].size()||t2<v[Top].size())
{
if(t1<v[Top-1].size()&&(t2==v[Top].size()||v[Top-1][t1]<v[Top][t2]))
t.pb(v[Top-1][t1++]);
else t.pb(v[Top][t2++]);
modify(rt[Top-1][t1+t2],rt[Top-1][t1+t2-1],1,N,t[t1+t2-1].second);
}
v[Top-1]=t;v[Top].clear();Top--;
}
}
int Query(int x,int y,int d)
{
int ret=0;
for(int i=1;i<=Top;++i)
{
int k1=lower_bound(v[i].begin(),v[i].end(),mp(x+d,inf))-v[i].begin();
int k2=lower_bound(v[i].begin(),v[i].end(),mp(x-d,000))-v[i].begin();
ret+=Query(rt[i][k2],rt[i][k1],1,N,max(1,y-d),min(y+d,N));
}
return ret;
}
int n,Q,a[MAX];
int main()
{
n=read();Q=read();
for(int i=1;i<MAXNODE;++i)St[++top]=i,vis[i]=false;
for(int i=1;i<=n;++i)a[i]=read(),insert(a[i]+i,a[i]-i+py);
while(Q--)
{
scanf("%s",ch);int x=read(),y=read();
if(ch[0]=='Q')printf("%d\n",Query(a[x]+x,a[x]-x+py,y));
else a[x]=y,insert(a[x]+x,a[x]-x+py);
}
return 0;
}
【BZOJ2989】数列(二进制分组,主席树)的更多相关文章
- [BZOJ 2989]数列(二进制分组+主席树)
[BZOJ 2989]数列(二进制分组+主席树) 题面 给定一个长度为n的正整数数列a[i]. 定义2个位置的graze值为两者位置差与数值差的和,即graze(x,y)=|x-y|+|a[x]-a[ ...
- bzoj2989&&4170数列——二进制分组+主席树
题意的转化挺巧妙的 可以联想到曼哈顿距离! 并且,所谓的修改还要查询历史版本,并且修改之间不动只算一次,不就是给平面上加一个点吗? 看成(x,a[x])的点 就是一个菱形区域 转切比雪夫距离,变成矩形 ...
- 2019.01.21 bzoj2989: 数列(二进制分组+主席树)
传送门 二进制分组入门题. 主席树写错调题2h+2h+2h+体验极差. 题意简述:给一堆点,支持加入一个点,询问有多少个点跟(x,y)(x,y)(x,y)曼哈顿距离不超过kkk. 思路:题目要求的是对 ...
- Dynamic Rankings ZOJ - 2112(主席树+树状数组)
The Company Dynamic Rankings has developed a new kind of computer that is no longer satisfied with t ...
- BZOJ2989 数列(二进制分组)
这题其实可以cdq分治做,但是如果强制在线的话,这里有个牛逼方法叫二进制分组. 它的基本思想是把修改操作按二进制分组,遇到修改就在尾部加一个,并与之前的合并,比如之前有23(16+4+2+1)个,加了 ...
- UOJ46. 【清华集训2014】玄学 [线段树,二进制分组]
UOJ 思路 模拟赛出了这题,结果我没学过二进制分组--一波主席树然后空间就爆炸了-- 用线段树维护时间序列,每个节点维护\(a_i\to x_i\times a_i+b_i,i\in [1,n]\) ...
- CDOJ 1104 求两个数列的子列的交集 查询区间小于A的数有多少个 主席树
求两个数列的子列的交集 Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.uestc.edu.cn/#/problem/show/1104 ...
- 【BZOJ3821/UOJ46】玄学(二进制分组,线段树)
[BZOJ3821/UOJ46]玄学(二进制分组,线段树) 题面 BZOJ UOJ 题解 呜,很好的题目啊QwQ. 离线做法大概可以线段树分治,或者直接点记录左右两次操作时的结果,两个除一下就可以直接 ...
- Codeforces 960 二进制构造子序列 完全二叉树shift模拟 主席树/MAP DP
A #include <bits/stdc++.h> #define PI acos(-1.0) #define mem(a,b) memset((a),b,sizeof(a)) #def ...
随机推荐
- react-native初体验(1) — hello world
没有简介,直接开始干活吧. 默认阅读本文的你已经安装好 nodejs, windows用户需要升级yarn到最新版本. 并且设置安装源为国内的淘宝源: npm config set registry ...
- Ubuntu下LimeSDR Mini使用说明
本文内容.开发板及配件仅限用于学校或科研院所开展科研实验! 淘宝店铺名称:开源SDR实验室 LimeSDR链接:https://item.taobao.com/item.htm?spm=a230r.1 ...
- RAID中条带的概念
raid把数据分成条带,一个条带横跨所有数据磁盘,每个磁盘上存储条带的一部分,称为sagment,也称为条带深度.一个条带包含的扇区或块的个数,称为条带长度. raid向操作系统提供的是卷,是连续的扇 ...
- python-两个筛子数据可视化(直方图)
""" 作者:zxj 功能:模拟掷骰子,两个筛子数据可视化 版本:3.0 日期:19/3/24 """ import random impo ...
- lscpu命令详解
基础命令学习目录首页 一.lscpu输出 使用lscpu查看的结果如下图,这里会显示很多信息,如下: 使用lscpu -p会详细的numa信息,如下: [root@localhost ~]# lscp ...
- python中数据分析常用函数整理
一. apply函数 作用:对 DataFrame 的某行/列应用函数之后,Apply 返回一些值.函数既可以使用默认的,也可以自定义.注意:在第二个输出中应用 head() 函数,因为它包含了很多行 ...
- 2018-2019-20172329 《Java软件结构与数据结构》第七周学习总结
2018-2019-20172329 <Java软件结构与数据结构>第七周学习总结 教材学习内容总结 <Java软件结构与数据结构>第十一章-二叉查找树 一.概述 1.什么是二 ...
- 如何知道一个App的包名呢
包名(Package name)是Android系统中判断一个APP的唯一标识 记录我获取包名的几种方式 方法一:通过cmd命令,打开你要获取包名的APP 1.adb shell 2.dumpsys ...
- 【每日scrum】第一次冲刺day2
和小伙伴一起找地图 ,学习了mapinfo地图格式的基本知识,数据和图像分开存储
- java微信开发之接口连接
个人学习注册订阅号就行,把资料尽可能完善,不然开发时权限不够,然后登陆,点击订阅号然后扫码就可以微信关注.