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

//嫌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. [Linux]ubuntu安装基本流程

    ubuntu安装基本流程 1.设置分辨率2.设置语言环境3.设置服务器镜像源4.添加终端5.apt.apt-get更新和升级系统软件 sudo apt update/upgrade sudo apt- ...

  2. 712. Minimum ASCII Delete Sum for Two Strings

    题目: Given two strings s1, s2, find the lowest ASCII sum of deleted characters to make two strings eq ...

  3. helm安装redis+Sentinel集群搭建

    一.redis集群特点 数据 在多个Redis节点之间自动分片 sentinel特点: 它的主要功能有以下几点 不时地监控redis是否按照预期良好地运行; 如果发现某个redis节点运行出现状况,能 ...

  4. yml配置文件读取出错 Exception in thread "main" while scanning for the next token found character '\t(TAB)'

    这几天在学习springboot的微服务项目,在配置文件方面也想尝试下新的yml配置,就想把原来项目properties写的文件转换成yml试一下(老项目是之前检出在了eclipse里面),结果写好了 ...

  5. 后缀html和htm文件的区别

    后缀html和htm文件的区别: (1)如果一个网站有index.html和index.htm,默认情况下,优先访问.html. (2).htm后缀是为了兼容以前的dos系统8.3的命名规范

  6. Linux中文件夹访问权限不足

    经常操作Linux服务器,远程访问会遇到403问题,本篇讲述nginx下文件权限的修改. 1.命令修改该文件夹的权限 chmod -R 755 /usr/local/....   文件夹的路径 2.修 ...

  7. Win10安装.NetFamework3.5

    步骤1:装载Win10安装镜像 本人用的是"cn_windows_10_multiple_editions_x64_dvd_6848463.iso" 如图,我把镜像装载到H盘; 步 ...

  8. python -- 函数进阶

    一.函数参数-动态传参       1.形参:         *   在形参位置, 表示此参数为不定参数,接受的是位置参数            并且接收到的位置参数的动态传参都是元组 def fu ...

  9. 第二课 --- git的(管理修改和撤销修改、删除文件)

    读取文本内容: cat readme.txt 查看工作区与脚本库里面的区别: git diff HEAD -- readme.txt 丢弃工作区的修改内容: git checkout -- readm ...

  10. .net core 2.0 webapi部署iis操作

    1.安装 .net core 2.0 runtime, (dotnet-runtime-2.0.7-win-x64.exe) https://www.microsoft.com/net/downloa ...