bzoj3489: A simple rmq problem (主席树)
//==========================
蒟蒻Macaulish:http://www.cnblogs.com/Macaulish/ 转载要声明!
//==========================
说好的“因为是OJ上的题,就简单点好了。”呢?
一开始看不懂,不会写。
然后跪了一个晚上决定看云的题解&……似乎是主席树套主席树!吓傻,还开了40000000的数组。然后一交tle……
然后p是不可能玩常数的。
找不到其他做法。
然后找到了神牛dwjshift,好心地提供了题解(http://tieba.baidu.com/p/2947256742#47989538012l )之后还好心地提供了代码!简直业界良心!(dwj必拿AU)
一开始写的是主席树套配对堆(因为要个优先队列嘛)然后发现套个配对堆各种问题(就不说了,其实主要还是我太弱,这个里面很恶心。。。。。)
然后就跪了dwj的代码,用treap实现。
关于题解,正题开始:
主席树套主席树的做法详见帖子中吧主的话
然后发现没有那么必要。
主席树套优先队列详见帖子下面。
“给定一个序列,支持两个操作: 1.给区间[l,r]中塞进或去掉一个数 2.查询某个点的值”
然后类似标记永久化。就可以变成单点修改的线段树套优先队列了
(原来主席树也是可以区间修改区间查询?……只要套个标记永久化)
写的时候有个地方傻叉了,就是主席树修改的时候,如果区间被要修改的区间覆盖时就不用改他儿子,这时候注意不要忘记儿子们=旧的点的儿子们,要是一开始新的点的儿子都是旧的点的不就行了,这是个人风格问题!)
{$inline on}
const
maxn=;
maxm=; var
left,right,size,fix,max2,value:array[..maxm]of longint;
root,lson,rson,max:array[..maxm]of longint;
time,last,next,num:array[..maxn]of longint;
n,m,tot1,tot2,ll,rr:longint;
flag:boolean; procedure swap(var x,y:longint);inline;
var
i:longint;
begin
i:=x;
x:=y;
y:=i;
end; function mmax(x,y:longint):longint;inline;
begin
if x<y then exit(y);
exit(x);
end; procedure update(x:longint);inline;
begin
max2[x]:=mmax(max2[left[x]],max2[right[x]]);
if size[x]> then max2[x]:=mmax(max2[x],value[x]);
end; procedure leftr(var x:longint);inline;
var
k:longint;
begin
k:=right[x];
right[x]:=left[k];
left[k]:=x;
max2[k]:=max2[x];
update(x);
x:=k;
end; procedure rightr(var x:longint);inline;
var
k:longint;
begin
k:=left[x];
left[x]:=right[k];
right[k]:=x;
max2[k]:=max2[x];
update(x);
x:=k;
end; procedure insert(var x:longint;y:longint);inline;
begin
if x= then begin
inc(tot2);
x:=tot2;
size[x]:=;
value[x]:=y;
fix[x]:=random(maxn);
max2[x]:=y;
left[x]:=;
right[x]:=;
exit;
end;
if value[x]=y then begin
if flag then inc(size[x])
else dec(size[x]);
update(x);
end
else
if value[x]<y then begin
insert(right[x],y);
update(x);
if fix[right[x]]>fix[x] then leftr(x);
end
else begin
insert(left[x],y);
update(x);
if fix[left[x]]>fix[x] then rightr(x);
end;
end; procedure change(x,old,l,r,y:longint;var new:longint);inline;
var
mid:longint;
begin
inc(tot1);
new:=tot1;
if (ll<=l) and (r<=rr) then begin
insert(root[x],y);
max[new]:=max2[root[x]];
lson[new]:=lson[old];
rson[new]:=rson[old];
exit;
end;
max[new]:=max2[root[x]];
mid:=(l+r)>>;
if ll>mid then begin
change(x<<+,rson[old],mid+,r,y,rson[new]);
lson[new]:=lson[old];
end
else
if rr<=mid then begin
change(x<<,lson[old],l,mid,y,lson[new]);
rson[new]:=rson[old];
end
else begin
change(x<<,lson[old],l,mid,y,lson[new]);
change(x<<+,rson[old],mid+,r,y,rson[new]);
end;
end; function query(x,l,r,y:longint):longint;inline;
var
mid:longint;
begin
//writeln(x);
if l=r then exit(max[x]);
mid:=(l+r)>>;
if y>mid then exit(mmax(max[x],query(rson[x],mid+,r,y)))
else exit(mmax(max[x],query(lson[x],l,mid,y)))
end; procedure into;inline;
var
i,j:longint;
begin
readln(n,m);
for i:= to n do read(num[i]);
for i:= to n do last[i]:=n+;
for i:=n downto do begin
next[i]:=last[num[i]];
last[num[i]]:=i;
end;
for i:= to n do
if last[i]<>n+ then begin
ll:=last[i];
rr:=next[ll]-;
flag:=true;
change(,time[],,n,i,time[]);
end;
for i:= to n do begin
ll:=i-;
rr:=next[i-]-;
flag:=false;
change(,time[i-],,n,num[i-],time[i]);
if next[i-]<>n+ then begin
ll:=next[i-];
rr:=next[ll]-;
flag:=true;
change(,time[i],,n,num[i-],time[i]);
end;
end;
end; procedure work;inline;
var
lastans,l,r:longint;
begin
lastans:=;
while m> do begin
// writeln(m,'===========');
dec(m);
readln(l,r);
l:=(l+lastans) mod n+;
r:=(r+lastans) mod n+;
if l>r then swap(l,r);
lastans:=query(time[l],,n,r);
writeln(lastans);
end;
end; begin
randomize;
into;
work;
end.
bzoj3489: A simple rmq problem (主席树)的更多相关文章
- 【题解】BZOJ3489 A Hard RMQ problem(主席树套主席树)
[题解]A simple RMQ problem 占坑,免得咕咕咕了,争取在2h内写出代码 upd:由于博主太菜而且硬是要用指针写两个主席树,所以延后2hQAQ upd:由于博主太菜而且太懒所以他决定 ...
- bzoj 3489 A simple rmq problem——主席树套线段树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3489 题解:http://www.itdaan.com/blog/2017/11/24/9b ...
- BZOJ.3489.A simple rmq problem(主席树 Heap)
题目链接 当时没用markdown写,可能看起来比较难受...可以复制到别的地方看比如DevC++. \(Description\) 给定一个长为n的序列,多次询问[l,r]中最大的只出现一次的数.强 ...
- bzoj 3489 A simple rmq problem —— 主席树套线段树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3489 题解:http://www.itdaan.com/blog/2017/11/24/9b ...
- BZOJ3489 A simple rmq problem 【可持久化树套树】*
BZOJ3489 A simple rmq problem Description 因为是OJ上的题,就简单点好了.给出一个长度为n的序列,给出M个询问:在[l,r]之间找到一个在这个区间里只出现过一 ...
- bzoj 3489: A simple rmq problem k-d树思想大暴力
3489: A simple rmq problem Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 551 Solved: 170[Submit][ ...
- bzoj3489 A simple rmq problem 可持久化树套树
先预处理出两个个数组pre,next.pre[i]表示上一个与i位置数字相同的位置,若不存在则设为0:next[i]表示下一个与i位置数字相同的位置,若不存在则设为n+1.那么一个满足在区间[L,R] ...
- 【bzoj3489】 A simple rmq problem k-d树
由于某些原因,我先打了一个错误的树套树,后来打起了$k-d$.接着因不明原因在思路上被卡了很久,在今天中午蹲坑时恍然大悟...... 对于一个数字$a_i$,我们可以用一组三维坐标$(i,pre,nx ...
- bzoj 3489 A simple rmq problem - 线段树
Description 因为是OJ上的题,就简单点好了.给出一个长度为n的序列,给出M个询问:在[l,r]之间找到一个在这个区间里只出现过一次的数,并且要求找的这个数尽可能大.如果找不到这样的数,则直 ...
随机推荐
- UItraIso 制作ubentu 系统失败
设备忙,请退出所有正在运行的应用程序,按确定按钮重试. 解决方法: 不要使用UItraIso,不知道为什么一直不行.重启了电脑几次都不行.用Rufus吧 https://rufus.ie/ 注意: r ...
- springboot 读写excel
添加两个坐标: <dependency> <groupId>org.apache.poi</groupId> <artifactId>poi</a ...
- CakePHP 总结
1. 处理任何保存或获取数据的操作最好都写在模型(Model)中.这个概念通常被称为fat model. 2. 返回上一条插入记录的ID, $this->Example->getInser ...
- Ruby 基础教程1-9
异常 1.异常结构 [ begin] ... rescue [retry] ... [ensure] . ...
- mysql中的直方图采样逻辑
int handler::sample_next(uchar *buf) { // Temporary set inited to RND, since we are calling rnd_next ...
- hdu1175连连看(dfs+细节)
连连看 Time Limit: 20000/10000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submi ...
- Linux命令应用大词典-第18章 磁盘分区
18.1 fdisk:分区表管理 18.2 parted:分区维护程序 18.3 cfdisk:基于磁盘进行分区操作 18.4 partx:告诉内核关于磁盘上分区的号码 18.5 sfdisk:用于L ...
- Android 简介
一 Android起源 android: 机器人 android是google公司开发的基于Linux2.6的免费开源操作系统 2005 Google收购 Android Inc. 开始 Dalvik ...
- Shader-水流效果
效果图:(贴图类似于泥石流) 代码: Shader "CookbookShaders/Chapter02/ScrollingUVs" { Properties { _MainTin ...
- 聊一聊session
最近从上家公司离职了,到了一家新公司,这几天一直在了解他们的项目,所以我自己的那个小项目也暂时搁浅了.. 今天差不多把他们的项目了解了,来院子写写我在这里边遇到的问题,影响最深刻的是seesion的. ...