算法模板——splay区间反转 1
时间:2022-05-07
本文章向大家介绍算法模板——splay区间反转 1,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
实现的功能:将序列区间反转,并维护
详见BZOJ3223
1 var
2 i,j,k,l,m,n,head,a1,a2:longint;
3 s1:ansistring;
4 a,b,c,d,fat,lef,rig:array[0..200000] of longint;
5 procedure swap(var x,y:longint);inline;
6 var z:longint;
7 begin
8 z:=x;x:=y;y:=z;
9 end;
10
11 procedure ext(x:longint);inline;
12 begin
13 if (x=0) then exit;
14 if c[x]=0 then exit;
15 swap(lef[x],rig[x]);
16 c[x]:=0;
17 c[lef[x]]:=1-c[lef[x]];
18 c[rig[x]]:=1-c[rig[x]];
19 c[0]:=0;
20 end;
21 procedure rt(x:longint);inline;
22 var f,l:longint;
23 begin
24 if x=0 then exit;
25 ext(x);
26 if lef[x]=0 then exit;
27 ext(lef[x]);
28 f:=x;l:=lef[x];
29 b[lef[x]]:=b[x];
30 b[x]:=b[rig[x]]+1+b[rig[l]];
31 lef[x]:=rig[l];
32 fat[rig[l]]:=x;
33 rig[l]:=x;
34 fat[l]:=fat[x];
35 fat[x]:=l;
36 if rig[fat[l]]=x then rig[fat[l]]:=l;
37 if lef[fat[l]]=x then lef[fat[l]]:=l;
38 fat[0]:=0;
39 end;
40 procedure lt(x:longint);inline;
41 var f,r:longint;
42 begin
43 if x=0 then exit;
44 ext(x);if rig[x]=0 then exit;
45 ext(rig[x]);
46 f:=x;r:=rig[x];
47 b[rig[x]]:=b[x];
48 b[x]:=1+b[lef[x]]+b[lef[r]];
49 rig[x]:=lef[r];
50 fat[lef[r]]:=x;
51 lef[r]:=x;
52 fat[r]:=fat[x];
53 fat[x]:=r;
54 if rig[fat[r]]=x then rig[fat[r]]:=r;
55 if lef[fat[r]]=x then lef[fat[r]]:=r;
56 fat[0]:=0;
57 end;
58 procedure ins(x,y:longint);inline;
59 begin
60 if a[y]<a[x] then
61 begin
62 if lef[x]=0 then
63 begin
64 lef[x]:=y;
65 fat[y]:=x;
66 end
67 else ins(lef[x],y);
68 end
69 else
70 begin
71 if rig[x]=0 then
72 begin
73 rig[x]:=y;
74 fat[y]:=x;
75 end
76 else ins(rig[x],y);
77 end;
78 b[x]:=1+b[lef[x]]+b[rig[x]];
79 end;
80 procedure up2(var x:longint);inline;
81 begin
82 if (fat[x]=0) or (x=0) then exit;
83 if lef[fat[x]]=x then
84 begin
85 if lef[fat[fat[x]]]=fat[x] then
86 begin
87 rt(fat[fat[x]]);
88 rt(fat[x]);
89 end
90 else
91 begin
92 rt(fat[x]);
93 lt(fat[x]);
94 end;
95 end
96 else
97 begin
98 if rig[fat[fat[x]]]=fat[x] then
99 begin
100 lt(fat[fat[x]]);
101 lt(fat[x]);
102 end
103 else
104 begin
105 lt(fat[x]);
106 rt(fat[x]);
107 end;
108 end;
109 end;
110 procedure up1(x:longint);inline;
111 begin
112 if (x=0) or (fat[x]=0) then exit;
113 if lef[fat[x]]=x then rt(fat[x]) else lt(fat[x]);
114 end;
115 procedure splay(x:longint);inline;
116 begin
117 if (x=0) or (fat[x]=0) then exit;
118 while fat[fat[x]]>0 do
119 up2(x);
120 if fat[x]>0 then up2(x);
121 head:=x;
122 end;
123 procedure splay2(x:longint);inline;
124 begin
125 if (x=0) or (fat[x]=0) then exit;
126 while fat[fat[fat[x]]]>0 do
127 up2(x);
128 if fat[fat[x]]>0 then up1(x);
129 end;
130 function getrank(x,y:longint):longint;inline;
131 begin
132 if (x=0) then exit(0);
133 ext(x);
134 if (b[lef[x]]+1)=y then exit(x);
135 if (b[lef[x]]+1)>y then exit(getrank(lef[x],y)) else exit(getrank(rig[x],y-1-b[lef[x]]));
136 end;
137 procedure turn(x,y:longint);inline;
138 var a1,a2:longint;
139 begin
140 if (x=1) and (y=n) then
141 c[head]:=1-c[head]
142 else
143 begin
144 if (x=1) then
145 begin
146 a1:=getrank(head,y+1);
147 splay(a1);
148 ext(a1);
149 c[lef[a1]]:=1-c[lef[a1]];
150 end
151 else
152 begin
153 if (y=n) then
154 begin
155 a2:=getrank(head,x-1);
156 splay(a2);
157 ext(a2);
158 c[rig[a2]]:=1-c[rig[a2]];
159 end
160 else
161 begin
162 a1:=getrank(head,x-1);
163 a2:=getrank(head,y+1);
164 splay(a2);splay2(a1);
165 ext(a2);ext(a1);
166 c[rig[a1]]:=1-c[rig[a1]];
167 end;
168 end;
169 end;
170 end;
171 function showoff(x:longint):ansistring;inline;
172 var s1:ansistring;
173 begin
174 if x=0 then exit('');
175 ext(x);
176 str(x,s1);
177 exit(showoff(lef[x])+s1+' '+showoff(rig[x]));
178 end;
179 begin
180 readln(n,m);
181 for i:=1 to n do
182 begin
183 a[i]:=i;c[i]:=0;b[i]:=1;
184 end;
185 head:=1;
186 for i:=2 to n do
187 begin
188 ins(head,i);
189 splay(random(i)+1);
190 end;
191 for i:=1 to m do
192 begin
193 readln(a1,a2);
194 turn(a1,a2);
195 end;
196 s1:=showoff(head);
197 writeln(s1);
198 readln;
199 end.
200
- 码农的瑞士军刀-脚本语言
- shell基础学习总结(一) (r3笔记第63天)
- 关于sysdba,sysoper,dba的区别(r3笔记第62天)
- 使用句柄实现特定场景的无备份恢复 (r3笔记第61天)
- 关于dual表的破坏性测试(r3笔记第60天)
- 哈希现金(Hashcash)与“工作量证明”
- 关于oracle中的sql数据类型(r3笔记第59天)
- 使用awr来分析sesson leak问题(r3笔记第78天)
- 弹窗层效果的实现(非jQuery实现)
- 简单的导航栏实现
- js实现css3的过渡,需要注意的一点(浏览器优化)
- 居中详解
- 通过pl/sql计算程序的运行时间(r3笔记第77天)
- 操作系统存储管理和oracle数据库(第一篇) (r3笔记第76天)
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- CentOS7如何执行PHP定时任务详解
- linux中ssh免密通信的实现
- Linux 启动停止SpringBoot jar 程序部署Shell 脚本的方法
- linux下安装boost库的完整步骤记录
- Linux下使用killall命令终止进程的8大用法实例详解
- 关于AIX挂载NFS写入效率低效的解决方法
- Linux查看进程的所有信息的方法示例
- 新版VPS主机管理面板WDCP安装及使用体验-国产简单易用型VPS面板
- linux中普通用户的定时任务详解
- 详解在Linux中清空或删除大文件内容的5种方法
- 教你一招实现Linux中的文本比对
- 怎么禁用 Ubuntu 服务器中终端欢迎消息中的广告
- Linux系统下部署项目的设置方法
- Linux中设置路由以及虚拟机联网图文详解
- 在Linux中如何一次重命名多个文件详解