【封装】Splay
注意确保操作合法性,否则可能陷入死循环
以点权作为排序依据
struct Splay{
#define ls p[u].son[0]
#define rs p[u].son[1]
#define maxn 100000
int root, cnt;
struct Node{
int val, fa, size, sum;
int son[2];
}p[maxn];
inline void destroy(int u){
p[u].val = ls = rs = p[u].fa = p[u].sum = p[u].size = 0;
}
inline int identify(int u){
return p[p[u].fa].son[1] == u;
}
inline void update(int u){
if(u) p[u].sum = p[ls].sum + p[rs].sum + p[u].size;
}
void rotate(int u){
int f = p[u].fa, gf = p[f].fa, sta = identify(u), sta_f = identify(f);
p[f].son[sta] = p[u].son[sta ^ 1];
p[p[f].son[sta]].fa = f;
p[u].son[sta^1] = f, p[f].fa = u, p[u].fa = gf;
p[gf].son[sta_f] = u;
update(f);
}
void splay(int u, int goal){
for(int f; (f = p[u].fa) && (f != goal); rotate(u))
if(p[f].fa != goal) rotate(identify(u) == identify(f) ? f : u);
if(!goal) root = u;
update(u);
}
void insert(int u){ // 函数结束后权值为u的节点变为根节点
if(!root){
p[++cnt].val = u;
p[cnt].size = p[cnt].sum = 1;
root = cnt;
return ;
}
int now = root, f = 0;
while(true){
if(u == p[now].val){
++p[now].size;
splay(now, 0);
return ;
}
f = now, now = p[now].son[p[now].val < u];
if(!now){
p[++cnt].val = u;
p[cnt].size = p[cnt].sum = 1;
p[cnt].fa = f, p[f].son[p[f].val < u] = cnt;
++p[f].sum;
splay(cnt, 0);
return ;
}
}
}
int find_val(int rank){
int now = root;
while(true){
if(p[now].son[0] && rank <= p[p[now].son[0]].sum)
now = p[now].son[0];
else{
int temp = p[p[now].son[0]].sum + p[now].size;
if(rank <= temp) return p[now].val;
now = p[now].son[1], rank -= temp;
}
}
}
int find_rank(int u){ // 函数结束后权值为u的节点是根节点
int now = root, rank = 0;
while(true){
if(u < p[now].val) now = p[now].son[0];
else{
rank += p[p[now].son[0]].sum;
if(u == p[now].val){
splay(now, 0);
return rank + 1;
}
rank += p[now].size, now = p[now].son[1];
}
}
}
int find_pre(int x){ // 返回x前驱节点编号
insert(x);
int now = p[root].son[0];
while(p[now].son[1]) now = p[now].son[1];
delete_val(x);
return now;
}
int find_suffix(int x){ // 返回x后继节点编号
insert(x);
int now = p[root].son[1];
while(p[now].son[0]) now = p[now].son[0];
delete_val(x);
return now;
}
void delete_val(int u){
find_rank(u); // 将权值为u的节点旋转到根节点
if(p[root].size > 1){
--p[root].size, --p[root].sum;
return ;
}
if(!p[root].son[0] && !p[root].son[1]){
destroy(root), root = 0;
return ;
}
int old_root = root;
if(!p[root].son[0]){
root = p[root].son[1];
p[root].fa = 0;
destroy(old_root);
return ;
}
if(!p[root].son[1]){
root = p[root].son[0];
p[root].fa = 0;
destroy(old_root);
return ;
}
int left_max = find_pre(u);
splay(left_max, 0);
p[root].son[1] = p[old_root].son[1];
p[p[old_root].son[1]].fa = root;
destroy(old_root);
update(root);
}
}splay;
以序列下标作为排序依据,常用于需要增删序列的操作(区间平衡树)
struct Splay{
#define ls p[u].son[0]
#define rs p[u].son[1]
#define maxn N
static const int inf = 1e9;
// 因为有虚点直接初始化了
int root = 1, top = 0, temp = 5e5 + 45;
int id[N], c[N], cnt[N];
struct Node{
int fa, size, len, reset, sum; // size是子树大小,len是除去虚点的子树大小
int tag, val, l_max, r_max, mid_max;
int son[2];
Node() { reset = inf; }
}p[maxn];
inline int identify(int u){
return p[p[u].fa].son[1] == u;
}
// 空间回收
void destroy(int u){
if(!u) return ;
if(ls) destroy(ls);
if(rs) destroy(rs);
p[u] = p[temp];
id[++top] = u;
}
inline void update(int u){
p[u].size = p[ls].size + p[rs].size + 1;
p[u].len = p[ls].len + p[rs].len + (u > 2); // 判断u > 2是为了除去虚点的影响
p[u].sum = p[ls].sum + p[rs].sum + p[u].val;
p[u].l_max = max(p[ls].l_max, p[ls].sum + p[u].val + p[rs].l_max);
p[u].r_max = max(p[rs].r_max, p[rs].sum + p[u].val + p[ls].r_max);
p[u].mid_max = max(p[u].val + p[ls].r_max + p[rs].l_max, max(p[ls].mid_max, p[rs].mid_max));
}
void change(int u, int val){
p[u].val = p[u].reset = val;
p[u].sum = p[u].val * p[u].len;
p[u].l_max = p[u].r_max = max(0, p[u].sum);
p[u].mid_max = max(val, p[u].sum);
}
inline void pushdown(int u){
if(p[u].reset != inf){
if(ls) change(ls, p[u].reset);
if(rs) change(rs, p[u].reset);
p[u].reset = inf, p[u].tag = 0;
}
if(p[u].tag){
if(ls) p[ls].tag ^= 1, swap(p[ls].son[0], p[ls].son[1]), swap(p[ls].l_max, p[ls].r_max);
if(rs) p[rs].tag ^= 1, swap(p[rs].son[0], p[rs].son[1]), swap(p[rs].l_max, p[rs].r_max);
p[u].tag = 0;
}
}
void rotate(int u){
int f = p[u].fa, gf = p[f].fa, sta = identify(u), sta_f = identify(f);
p[f].son[sta] = p[u].son[sta ^ 1];
p[p[f].son[sta]].fa = f;
p[u].son[sta^1] = f, p[f].fa = u, p[u].fa = gf;
p[gf].son[sta_f] = u;
update(f);
}
void splay(int u, int goal){
for(int f; (f = p[u].fa) && (f != goal); rotate(u)){
if(p[f].fa != goal) rotate(identify(u) == identify(f) ? f : u);
}
if(!goal) root = u;
update(u);
}
int find_Kth(int k){
int u = root;
while(1){
pushdown(u);
if(p[ls].size + 1 == k) return u;
if(p[ls].size >= k) u = ls;
else k -= p[ls].size + 1, u = rs;
}
}
int build(int l, int r, int fa){
if(l > r) return 0;
int mid = (l + r) >> 1, now = cnt[mid];
if(l == r){
p[now].val = c[l];
p[now].fa = fa;
update(now);
p[now].mid_max = p[now].val;
return now;
}
p[now].fa = fa, p[now].val = c[mid];
p[now].son[0] = build(l, mid - 1, now);
p[now].son[1] = build(mid + 1, r, now);
update(now);
return now;
}
// 在第u个位置后面插入值插入tot个数
void insert(int u, int tot){
for(int i = 1; i <= tot; ++i) scanf("%d", &c[i]), cnt[i] = id[top--];
int rt = build(1, tot, 0);
int L = find_Kth(u), R = find_Kth(u + 1);
splay(L, 0), splay(R, L);
p[rt].fa = R, p[R].son[0] = rt;
splay(rt, 0);
}
// 区间删除
void delete_range(int pos, int tot){
int L = find_Kth(pos), R = find_Kth(pos + tot + 1);
splay(L, 0), splay(R, L);
destroy(p[R].son[0]);
p[R].son[0] = 0;
update(R), update(L);
}
// 区间修改
void modify_range(int pos, int tot, int set){
int L = find_Kth(pos), R = find_Kth(pos + tot + 1);
splay(L, 0), splay(R, L);
int u = p[R].son[0];
p[u].reset = p[u].val = set;
p[u].sum = p[u].len * set;
p[u].l_max = p[u].r_max = max(0, p[u].sum);
p[u].mid_max = max(p[u].sum, set);
update(R), update(L);
}
// 翻转区间
void reverse(int pos, int tot){
int L = find_Kth(pos), R = find_Kth(pos + tot + 1);
splay(L, 0), splay(R, L);
int u = p[R].son[0];
p[u].tag ^= 1;
swap(ls, rs);
swap(p[u].l_max, p[u].r_max);
update(R), update(L);
}
int get_sum(int pos, int tot){
int L = find_Kth(pos), R = find_Kth(pos + tot + 1);
splay(L, 0), splay(R, L);
return p[p[R].son[0]].sum;
}
// 非空子段最大值
int Max(){
return p[root].mid_max;
}
// 插入虚点,之后操作要注意下标
void init(){
for(int i = 3; i <= N - 45; ++i) id[++top] = i;
p[1].son[1] = 2, p[2].fa = 1;
p[1].size = 2, p[2].size = 1;
// 虚点赋值负无穷,消除影响,根据不同题而定
p[0].mid_max = p[1].mid_max = p[2].mid_max = p[0].val = p[1].val = p[2].val = -inf;
}
}splay;
【封装】Splay的更多相关文章
- [学习笔记] Splay Tree 从入门到放弃
前几天由于出行计划没有更博QwQ (其实是因为调试死活调不出来了TAT我好菜啊) 伸展树 伸展树(英语:Splay Tree)是一种二叉查找树,它能在O(log n)内完成插入.查找和删除操作.它是由 ...
- BZOJ 3196 Tyvj 1730 二逼平衡树:线段树套splay
传送门 题意 给你一个长度为 $ n $ 有序数列 $ a $ ,进行 $ m $ 次操作,操作有如下几种: 查询 $ k $ 在区间 $ [l,r] $ 内的排名 查询区间 $ [l,r] $ 内排 ...
- AVL树、splay树(伸展树)和红黑树比较
AVL树.splay树(伸展树)和红黑树比较 一.AVL树: 优点:查找.插入和删除,最坏复杂度均为O(logN).实现操作简单 如过是随机插入或者删除,其理论上可以得到O(logN)的复杂度,但是实 ...
- [C#] 简单的 Helper 封装 -- RegularExpressionHelper
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.T ...
- iOS开发之App间账号共享与SDK封装
上篇博客<iOS逆向工程之KeyChain与Snoop-it>中已经提到了,App间的数据共享可以使用KeyChian来实现.本篇博客就实战一下呢.开门见山,本篇博客会封装一个登录用的SD ...
- Ajax实现原理,代码封装
都知道实现页面的异步操作需要使用Ajax,那么Ajax到是怎么实现异步操作的呢? 首先需要认识一个对象 --> XMLHttpRequest 对象 --> Ajax的核心.它有许多的属性和 ...
- 用C语言封装OC对象(耐心阅读,非常重要)
用C语言封装OC对象(耐心阅读,非常重要) 本文的主要内容来自这里 前言 做iOS开发的朋友,对OC肯定非常了解,那么大家有没有想过OC中NSInteger,NSObject,NSString这些对象 ...
- 【知识必备】RxJava+Retrofit二次封装最佳结合体验,打造懒人封装框架~
一.写在前面 相信各位看官对retrofit和rxjava已经耳熟能详了,最近一直在学习retrofit+rxjava的各种封装姿势,也结合自己的理解,一步一步的做起来. 骚年,如果你还没有掌握ret ...
- 对百度WebUploader开源上传控件的二次封装,精简前端代码(两句代码搞定上传)
前言 首先声明一下,我这个是对WebUploader开源上传控件的二次封装,底层还是WebUploader实现的,只是为了更简洁的使用他而已. 下面先介绍一下WebUploader 简介: WebUp ...
- 封装集合(Encapsulate Collection)
封装就是将相关的方法或者属性抽象成为一个对象. 封装的意义: 对外隐藏内部实现,接口不变,内部实现自由修改. 只返回需要的数据和方法. 提供一种方式防止数据被修改. 更好的代码复用. 当一个类的属性类 ...
随机推荐
- VMware中的三种网络模式
1.桥接模式网络 通过桥接模式网络连接,虚拟机中的虚拟网络适配器可连接到主机中的物理网络适配器.虚拟机可通过主机网络适配器连接到主机系统所用的 LAN.桥接模式网络连接支持有线和无线主机网络适配器. ...
- 文盘Rust -- Mutex解决并发写文件乱序问题
在实际开发过程中,我们可能会遇到并发写文件的场景,如果处理不当很可能出现文件内容乱序问题.下面我们通过一个示例程序描述这一过程并给出解决该问题的方法. use std::{ fs::{self, Fi ...
- 基于C#的无边框窗体阴影绘制方案 - 开源研究系列文章
今天介绍无边框窗体阴影绘制的内容. 上次有介绍使用双窗体的方法来显示阴影,这次介绍使用API函数来进行绘制.这里使用的是Windows API函数,操作系统的窗体也是用的这个来进行的绘制. 1. 项目 ...
- CodeForces 1332D Walk on Matrix
题意 \(Bob\)想解决一个问题:一个\(n\cdot m\)的矩阵,从\((1,1)\)出发,只能走右和下,问从\((1,1)\)到\((n,m)\)的最大\(\&\)和 他的算法如下(\ ...
- [TSG开发日志](一)软件基础框架
目录 前言 说明 框架 TSG_Framework 一.底层信号机制 TSG_Caller 二.参数类型声明 TSG_Params 三.设备类声明 TSG_Device 四.设备配置文件控制 TSG_ ...
- DesignPattern-part3
title: "modern C++ DesignPattern-Part3" date: 2018-04-12T19:08:49+08:00 lastmod: 2018-04-1 ...
- 扩展ABP的Webhook功能,推送数据到第三方接口(企业微信群、钉钉群等)
前言 在上一篇文章[基于ASP.NET ZERO,开发SaaS版供应链管理系统]中有提到对Webhook功能的扩展改造,本文详细介绍一下具体过程. Webhook功能操作说明,请参见此文档链接:Web ...
- 开源.NetCore通用工具库Xmtool使用连载 - 图形验证码篇
[Github源码] <上一篇> 介绍了Xmtool工具库中的Web操作类库,今天我们继续为大家介绍其中的图形验证码类库. 图形验证码是为了抵御恶意攻击出现的一种设计:例如用户登录.修改密 ...
- 「hdu - 5780」gcd
link. 钦定 \(i>j\),研究得 \((x^i-1,x^j-1)\rightleftharpoons(x^i-x^j,x^j-1)\rightleftharpoons(x^j(x^{i- ...
- Cplex求解教程(基于OPL语言,可作为大规模运算输入参考)
最近导导让牛牛改篇论文,牛牛在她的指导下把非线性问题化成了线性.然鹅,化成线性后的模型决策变量和约束条件均达到上百甚至上千个,这让牛牛犯了难,以下方法或许能为这样大规模模型的变量和约束输入提供思路(๑ ...