注意!本帖不是算法介绍!只是贴代码(逃)

//嫌stdlib的rand太慢,手打了一个

 /*
Author: hotwords
*/ typedef unsigned int tkey;
class Random {
private:
tkey rs1,rs2,rs3;
public:
tkey operator()() {
rs1=(rs1<<)^(rs2>>)^rs3;
rs2=(rs2<<)^(rs3>>)^rs1;
rs3=(rs1<<)^(rs2>>)^rs2;
return rs1^rs2^rs3;
}
void seed(tkey x) {
rs1=x;
rs2=rs1^0x64db3c5e;
rs3=rs1^0xdb7a6f81;
}
Random() {
seed(0x86e241ca);
}
} random;
template<typename tp>
class Treap {
private:
class TNode {
private:
int p_size,p_csize;
public:
tp val;int cnt;
tkey key;
TNode *lc,*rc;
TNode(tp v) {
val=v;p_size=;
p_csize=cnt=;
key=random();
lc=rc=;
}
inline int size()const{return this?this->p_size:;}
inline int csize()const{return this?this->p_csize:;}
void combine() {
p_size=+lc->size()+rc->size();
p_csize=cnt+lc->csize()+rc->csize();
}
void split(int k,TNode*& l,TNode*& r) {
if(!this) {
l=r=;return;
}
if(k<=lc->size()) {
lc->split(k,l,r);
lc=r;r=this;
} else {
rc->split(k-lc->size()-,l,r);
rc=l;l=this;
}
combine();
}
tp min() const {
const TNode *cur=this;
while(cur->lc) cur=cur->lc;
return cur->val;
}
tp max() const {
const TNode *cur=this;
while(cur->rc) cur=cur->rc;
return cur->val;
}
} *root;
TNode* merge(TNode* l,TNode* r) {
if(!r) return l;
if(!l) return r;
if(l->key<r->key) {
l->rc=merge(l->rc,r);
l->combine();
return l;
} else {
r->lc=merge(l,r->lc);
r->combine();
return r;
}
}
public:
Treap(){root=;}
int size()const{return root->size();}
int count()const{return root->csize();}
int rank(tp v) const {
int ans=;
TNode *cur=root;
while(cur) {
if(v<cur->val) cur=cur->lc;
else {
ans+=cur->lc->csize();
if(v>cur->val) ans+=cur->cnt,cur=cur->rc;
else return ans+;
}
}
return ans;
}
int nrank(tp v) const {
int ans=;
TNode *cur=root;
while(cur) {
if(v<cur->val) cur=cur->lc;
else ans+=cur->lc->size()+,cur=cur->rc;
}
return ans;
}
void insert(tp v) {
int k=nrank(v);
TNode *l,*m,*r,*t;
root->split(k,t,r);
t->split(k-,l,m);
if(m&&m->val==v) ++m->cnt,m->combine();
else m=merge(m,new TNode(v));
root=merge(merge(l,m),r);
}
void erase(tp v) {
int k=nrank(v);
TNode *l,*m,*r,*t;
root->split(k,t,r);
t->split(k-,l,m);
if(!--m->cnt) {
delete []m;m=;
} else m->combine();
root=merge(merge(l,m),r);
}
tp get_kth(int k) const {
TNode *cur=root;
while(cur) {
if(k<=cur->lc->csize()) cur=cur->lc;
else {
k-=cur->lc->csize();
if(k<=cur->cnt) return cur->val;
k-=cur->cnt;cur=cur->rc;
}
}
}
tp prev(tp v) {
int k=nrank(v);
TNode *l,*m,*r,*t;
tp ans;
root->split(k,t,r);
t->split(k-,l,m);
ans=m->val<v?m->val:l->max();
root=merge(merge(l,m),r);
return ans;
}
tp next(tp v) {
int k=nrank(v);
TNode *l,*r;
tp ans;
root->split(k,l,r);
ans=r->min();
root=merge(l,r);
return ans;
}
};

[模板] FHQ Treap

[模板] 无旋Treap (C++ class)的更多相关文章

  1. 模板 - 无旋Treap

    一般而言作为一棵平衡树只需要插入,删除,值求排名,排名求值,前驱,后继,六个接口. #include<bits/stdc++.h> using namespace std; typedef ...

  2. 模板——无旋Treap

    #include "bits/stdc++.h" using namespace std; inline int read(){ ,k=;char ch=getchar(); :, ...

  3. 无旋treap的简单思想以及模板

    因为学了treap,不想弃坑去学splay,终于理解了无旋treap... 好像普通treap没卵用...(再次大雾) 简单说一下思想免得以后忘记.普通treap因为带旋转操作似乎没卵用,而无旋tre ...

  4. 模板 - 数据结构 - 可持久化无旋Treap/PersistentFHQTreap

    有可能当树中有键值相同的节点时,貌似是要对Split和Merge均进行复制的,本人实测:只在Split的时候复制得到了一个WA,但只在Merge的时候复制还是AC,可能是恰好又躲过去了.有人说假如确保 ...

  5. 洛谷 - P3391 【模板】文艺平衡树(Splay) - 无旋Treap

    https://www.luogu.org/problem/P3391 使用无旋Treap维护序列,注意的是按顺序插入的序列,所以Insert实际上简化成直接root和Merge合并,但是假如要在序列 ...

  6. 浅谈无旋treap(fhq_treap)

    一.简介 无旋Treap(fhq_treap),是一种不用旋转的treap,其代码复杂度不高,应用范围广(能代替普通treap和splay的所有功能),是一种极其强大的平衡树. 无旋Treap是一个叫 ...

  7. [转载]无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    转自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182631.html 1500: [NOI2005]维修数列 Time Limit: 10 Sec  Mem ...

  8. [转载]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)

    转载自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182491.html 今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和t ...

  9. [您有新的未分配科技点]无旋treap:从好奇到入门(例题:bzoj3224 普通平衡树)

    今天我们来学习一种新的数据结构:无旋treap.它和splay一样支持区间操作,和treap一样简单易懂,同时还支持可持久化. 无旋treap的节点定义和treap一样,都要同时满足树性质和堆性质,我 ...

随机推荐

  1. 【编程语言】Kotlin之object关键字

    在一个体重秤项目里面使用Kotlin开发,考虑到项目比较小型轻量,所以和团队申请决定使用Kotlin开发,以此熟悉和尝试一下Kotlin. 首先使用Kotlin之后,发现能和Java很好的兼容一起,开 ...

  2. VS2017 MVC Spring net 环境配置

    首先打开管理NuGet程序包. 搜索 "spring web",安装Spring.Web,Spring.Web.Mvc5,Spring.Web.Extensions, 搜索Micr ...

  3. xlwt模块的使用

    前记:Python处理表格时会用到xlwt和xlrd模块 xlwt设置行高:row sheet.row(2).set_style(xlwt.easyxf('font:height 440;')) 13 ...

  4. Educational Codeforces Round 41 (Rated for Div. 2)F. k-substrings

    题意比较麻烦略 题解:枚举前缀的中点,二分最远能扩展的地方,lcp来check,然后线段树维护每个点最远被覆盖的地方,然后查询线段树即可 //#pragma GCC optimize(2) //#pr ...

  5. 设置Oracle数据库开机自启动-亲试ok

      1.oracle 用户下 修改$ORACLE_HOME/bin/dbstart [oracle@dev-oracle ~]$  vi  /opt/oracle/app/oracle/product ...

  6. C++标准模板库(STL)之Set

    1.Set的用法 Set:集合,一个内部自动有序而且不重复元素的容器.使用set,要加头文件#include<set>和using namespace std; 1.1.Set的定义 se ...

  7. python 两个 list 获取交集,并集,差集的函数

    1. 获取两个 list 的交集 a = [1, 2, 3, 4] b = [1, 2, 5] print(list(set(a).intersection(set(b)))) 2. 获取两个 lis ...

  8. EAP-MD5计算方法

    一.说明 领导要求确认设备进行802.1X认证时,是否直接将用户名密码明文传给交换机.配好端口镜像.搭好Radius环境后进行了抓包分析. 二.计算 2.1 802.1X认证过程 完整流程如下: 客户 ...

  9. Python数据类型——字符串

    概论 字符串顾名思义就是一串字符,由于Python中没有“字符”这种数据类型,所以单个的字符也依然是字符串类型的.字符串可以包含一切数据,无论是能从键盘上找到的,还是你根本都不认识的.与数一样,字符串 ...

  10. messageQ 消息队列

    之后就是对MessageQ的打开,关闭, 消息的发送和接受. MessageQ_create(),MessageQ_delete(): 对消息的创建和删除. MessageQ_open(); Mess ...