树板
突然发现没存过平衡树板子,趁这几天在复习平衡树,存一下还热乎的未封装版本。
update 2022.2.24 上午更新了Splay板子,但是依然没封装。下午突然意识到Splay不是时刻平衡,所以每次操作都要Splay,但是上午的版本kth没有,所以又改了一下。
update 2022.2.27 顺手再存一个LCT板子,将平衡树板子
更名为树板
。
$\rm FHQ\ Treap$
//12252024832524
#include <cstdio>
#include <cstring>
#include <algorithm>
#define TT template<typename T>
using namespace std;
typedef long long LL;
const int MAXN = 100005;
const int MOD = 998244353;
int n;
LL Read()
{
LL x = 0,f = 1;char c = getchar();
while(c > '9' || c < '0'){if(c == '-')f = -1;c = getchar();}
while(c >= '0' && c <= '9'){x = (x*10) + (c^48);c = getchar();}
return x * f;
}
TT void Put1(T x)
{
if(x > 9) Put1(x/10);
putchar(x%10^48);
}
TT void Put(T x,char c = -1)
{
if(x < 0) putchar('-'),x = -x;
Put1(x); if(c >= 0) putchar(c);
}
TT T Max(T x,T y){return x > y ? x : y;}
TT T Min(T x,T y){return x < y ? x : y;}
TT T Abs(T x){return x < 0 ? -x : x;}
int tot,rt;
struct FHQ_Treap
{
int l,r,val,rd,siz;
}t[MAXN];
void up(int x){t[x].siz = t[t[x].l].siz + t[t[x].r].siz + 1;}
void split(int now,int val,int &x,int &y)
{
if(!now) {x = y = 0;return;}
else
{
if(t[now].val <= val) x = now,split(t[now].r,val,t[now].r,y);
else y = now,split(t[now].l,val,x,t[now].l);
up(now);
}
}
int mge(int x,int y)
{
if(!x || !y) return x|y;
if(t[x].rd < t[y].rd) {t[x].r = mge(t[x].r,y);up(x);return x;}
else {t[y].l = mge(x,t[y].l);up(y);return y;}
}
int newnode(int val)
{
t[++tot] = {0,0,val,rand(),1};
return tot;
}
int x,y,z;
void ins(int val)
{
// printf("ins ins ins\n");
split(rt,val-1,x,y);
// printf("ins %d %d\n",x,y);
rt = mge(x,mge(newnode(val),y));
}
void del(int val)
{
split(rt,val,x,z);
split(x,val-1,x,y);
rt = mge(mge(x,mge(t[y].l,t[y].r)),z);
}
int qrk(int val)
{
split(rt,val-1,x,y);
int ret = t[x].siz+1;
rt = mge(x,y);
return ret;
}
int qval(int rk)
{
x = rt;
while(1)
{
if(t[t[x].l].siz+1 == rk) return t[x].val;
else if(t[t[x].l].siz >= rk) x = t[x].l;
else rk -= t[t[x].l].siz+1,x = t[x].r;
}
}
int pre(int val)
{
split(rt,val-1,x,y);
int now = x;
while(t[now].r) now = t[now].r;
rt = mge(x,y);
return t[now].val;
}
int suf(int val)
{
split(rt,val,x,y);
int now = y;
while(t[now].l) now = t[now].l;
rt = mge(x,y);
return t[now].val;
}
void dfs(int now)
{
if(!now) return;
// printf("dfs %d %d %d %d\n",now,t[now].val,t[now].l,t[now].r);
dfs(t[now].l);
dfs(t[now].r);
}
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
srand(42523);
for(int T = Read(); T ;-- T)
{
int opt = Read();
if(opt == 1) ins(Read());
else if(opt == 2) del(Read());
else if(opt == 3) Put(qrk(Read()),'\n');
else if(opt == 4) Put(qval(Read()),'\n');
else if(opt == 5) Put(pre(Read()),'\n');
else if(opt == 6) Put(suf(Read()),'\n');
}
return 0;
}
$\rm Splay$
//12252024832524
#include <bits/stdc++.h>
#define TT template<typename T>
using namespace std;
typedef long long LL;
const int MAXN = 100005;
const int INF = 0x3f3f3f3f;
int n;
LL Read()
{
LL x = 0,f = 1;char c = getchar();
while(c > '9' || c < '0'){if(c == '-')f = -1;c = getchar();}
while(c >= '0' && c <= '9'){x = (x*10) + (c^48);c = getchar();}
return x * f;
}
TT void Put1(T x)
{
if(x > 9) Put1(x/10);
putchar(x%10^48);
}
TT void Put(T x,char c = -1)
{
if(x < 0) putchar('-'),x = -x;
Put1(x); if(c >= 0) putchar(c);
}
TT T Max(T x,T y){return x > y ? x : y;}
TT T Min(T x,T y){return x < y ? x : y;}
TT T Abs(T x){return x < 0 ? -x : x;}
#define lc t[x].ch[0]
#define rc t[x].ch[1]
int tot,rt;
struct node{
int ch[2],f,siz,cnt,val;
}t[MAXN];
int newnode(int val,int fa){t[++tot] = node{0,0,fa,1,1,val};if(fa) t[fa].ch[val > t[fa].val] = tot;return tot;}
void up(int x){t[x].siz = t[lc].siz + t[rc].siz + t[x].cnt;}
bool ck(int x){return x == t[t[x].f].ch[1];}
void lk(int x,int fa,bool d){t[x].f = fa; if(fa) t[fa].ch[d] = x;}
void rotate(int x){
int fa = t[x].f,d = ck(x);
lk(t[x].ch[d^1],fa,d);
lk(x,t[fa].f,ck(fa));
lk(fa,x,d^1);
up(fa); up(x);
}
void splay(int x,int to = 0){
while(t[x].f ^ to){
if(t[t[x].f].f ^ to){
if(ck(x) == ck(t[x].f)) rotate(t[x].f);
else rotate(x);
}
rotate(x);
}
if(!to) rt = x;//attention!
}
void ins(int val){
int x = rt,fa = 0;
while(x && t[x].val != val) fa = x,x = t[x].ch[val > t[x].val];
if(x) ++t[x].cnt,++t[x].siz;
else x = newnode(val,fa);
splay(x);
}
void tort(int val){
int x = rt;
while(t[x].ch[val > t[x].val] && t[x].val != val) x = t[x].ch[val > t[x].val];
splay(x);
}
int qrk(int val){
tort(val);
if(t[rt].val >= val) return t[t[rt].ch[0]].siz;
else return t[t[rt].ch[0]].siz+t[rt].cnt;
}
int kth(int rk){
int x = rt; ++rk;
while(1){
if(t[lc].siz >= rk) x = lc;
else if(t[lc].siz + t[x].cnt >= rk) {splay(x);return x;}
else rk -= t[lc].siz+t[x].cnt,x = rc;
}
}
int pre(int val){return kth(qrk(val)-1);}
int suf(int val){return kth(qrk(val+1));}
void del(int val){
int p = pre(val),s = suf(val);
splay(p); splay(s,p);
if(t[t[s].ch[0]].cnt > 1) --t[t[s].ch[0]].cnt,--t[t[s].ch[0]].siz;
else t[s].ch[0] = 0;
up(p);//new
}
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
ins(-INF); ins(INF);
for(int T = Read(); T ;-- T){
int opt = Read();
if(opt == 1) ins(Read());
else if(opt == 2) del(Read());
else if(opt == 3) Put(qrk(Read()),'\n');
else if(opt == 4) Put(t[kth(Read())].val,'\n');
else if(opt == 5) Put(t[pre(Read())].val,'\n');
else if(opt == 6) Put(t[suf(Read())].val,'\n');
}
return 0;
}
$\rm LCT$
//12252024832524
#include <bits/stdc++.h>
#define TT template<typename T>
using namespace std;
typedef long long LL;
const int MAXN = 100005;
int n,m;
int val[MAXN];
LL Read()
{
LL x = 0,f = 1;char c = getchar();
while(c > '9' || c < '0'){if(c == '-')f = -1;c = getchar();}
while(c >= '0' && c <= '9'){x = (x*10) + (c^48);c = getchar();}
return x * f;
}
TT void Put1(T x)
{
if(x > 9) Put1(x/10);
putchar(x%10^48);
}
TT void Put(T x,char c = -1)
{
if(x < 0) putchar('-'),x = -x;
Put1(x); if(c >= 0) putchar(c);
}
TT T Max(T x,T y){return x > y ? x : y;}
TT T Min(T x,T y){return x < y ? x : y;}
TT T Abs(T x){return x < 0 ? -x : x;}
#define lc (t[x].ch[0])
#define rc (t[x].ch[1])
struct node{
int ch[2],f,val;
bool rev;
}t[MAXN];
bool ck(int x){return t[t[x].f].ch[1] == x;}
void lk(int x,int fa,bool d){t[x].f = fa; if(fa) t[fa].ch[d] = x;}
bool isroot(int x){return t[t[x].f].ch[0] != x && t[t[x].f].ch[1] != x;}
void calcrev(int x){if(!x) return; swap(lc,rc); t[x].rev ^= 1;}
void up(int x){t[x].val = t[lc].val ^ val[x] ^ t[rc].val;}
void down(int x){if(!t[x].rev) return; calcrev(lc); calcrev(rc); t[x].rev = 0;}
void rotate(int x){//这里没有down的原因是LCT的down是自上而下,因此Splay中的rotate执行前就已经down完了
int fa = t[x].f,d = ck(x);
if(isroot(fa)) t[x].f = t[fa].f;
else lk(x,t[fa].f,ck(fa));
lk(t[x].ch[d^1],fa,d);
lk(fa,x,d^1);
up(fa); up(x);
}
void downtag(int x){
if(!isroot(x)) downtag(t[x].f);
down(x);
}
void splay(int x){
downtag(x);
for(int y = t[x].f;!isroot(x);rotate(x),y = t[x].f)
if(!isroot(y)) rotate(ck(x) == ck(y) ? y : x);
// up(x);
}
void access(int x){
for(int lst = 0; x ;lst = x,x = t[x].f)
splay(x),rc = lst,up(x);
}
void makeroot(int x){
access(x); splay(x); calcrev(x);
}
int findroot(int x){
access(x); splay(x);
while(lc) down(x),x = lc;
splay(x); return x;
}
void split(int x,int y){
makeroot(x); access(y); splay(y);
}
void link(int x,int y){
if(findroot(x) == findroot(y)) return;
makeroot(x); t[x].f = y;
}
void cut(int x,int y){
if(findroot(x) ^ findroot(y)) return;
makeroot(x); access(y); splay(x);
if(t[y].f == x && !t[y].ch[0]) t[y].f = rc = 0,up(x);
}
int main()
{
// freopen(".in","r",stdin);
// freopen(".out","w",stdout);
n = Read(); m = Read();
for(int i = 1;i <= n;++ i) t[i].val = val[i] = Read();
while(m --> 0){
int opt = Read(),u,v;
if(!opt) u = Read(),v = Read(),split(u,v),Put(t[v].val,'\n');
else if(opt == 1) link(Read(),Read());
else if(opt == 2) cut(Read(),Read());
else if(opt == 3) u = Read(),splay(u),val[u] = Read(),up(u);
}
return 0;
}