

inline void turn(int rt,int k){
int x=a[rt].f,y=a[x].f;
pushdown(x);pushdown(rt);
a[x].son[k^1]=a[rt].son[k];
if(a[rt].son[k])a[a[rt].son[k]].f=x;
a[rt].f=y;
if(y)a[y].son[a[y].son[1]==x]=rt;
a[x].f=rt;
a[rt].son[k]=x;
pushup(x);pushup(rt);
}
void spaly(int rt,int ancestry){
while(a[rt].f!=ancestry)turn(rt,a[a[rt].f].son[0]==rt);
if(ancestry==0)root=rt;
}
这样你就有20分暴力了。。。
然后,慢慢分析单旋的性质。
1. 插入操作:
我们会发现这样一个事实:
一个节点的插入所经过的$Spaly$路径,必然经过这个节点的$key$值在$Spaly$上的前驱与后继!
并且新插入的节点一定在前驱与后继中深度更深的节点的儿子节点上!
那就好办了,我们只要找到$key$的前驱与后继就可以了。
怎么找?
拿出了$Spaly$的原型——双旋的$Splay$!
$Splay$的复杂度是$O(\log_2n)$的,当然可以胜任。
于是只要$Splay$上的节点与$Spaly$的节点的编号一一对应就可以了。
为了方便,我们让它们相同。
找到$key$的前驱与后继了,就挂在前驱的右儿子,或者是后继的左儿子上,看那个深度更深就是了。
这个操作完美解决。
2. 单旋最小值:
我们知道,$Spaly$与$Splay$的最小值一定没有左儿子!
我们可以手动模拟一下$Spaly$上的最小值单旋到根的过程。(尽量把$Spaly$的深度画成3或者4,方便寻找规律)
我们又发现了一个事实:
当最小值单旋到根之后,最小值的深度为1(这不废话嘛)。
此时最小值仅有右子树。
但右子树即为除了最小值的$Spaly$!
并且形态只改变了一点:原来最小值的右子树 / 接到了 / 原来最小值的父亲的左子树 / 上!(“/”表停顿)
而且深度还没有变!
听得比较绕。。。
但是这个性质让我们可以用$Splay$维护某节点在$Spaly$上的深度!
把最小值的深度改为$1$不用多说。
之后呢?
先将最小值的右子树打上$-1$的标记(以后会抵消),再在$Splay$上旋转到根,然后对其右子树打上$+1$标记(抵消了前面的$-1$标记)。
那,$Spaly$的形态怎么维护呢?
把最小值与其父亲的父子关系断开,最小值与其右儿子的父子关系断开,将最小值的右儿子与最小值的父亲连上父子关系,最小值的右儿子挂在最小值的父亲的左儿子上。
再将最小值与原来$Spaly$的根节点连上父子关系,最小值为父亲。
然后每次旋转的时候不要忘了$pushdown$一下。。。
然后就没了。。。
3. 单旋最大值:
这个跟操作2差不多吧,省略。。
4. 单旋删除最小值:
题目都说了先执行操作2。。。
然后按照题目所说断开最小值(现在的根节点)与原来的根节点的父子关系,之后把根节点换成原来的根节点,并打上$+1$的标记。
然后就没了啊。。。
5. 单旋删除最大值:
同操作4,这里就省略了。
呼~终于解决了!
总复杂度应该是$Splay$的复杂度,即均摊$O(n\log_2n)$。
这题真的是毒瘤。。。
真心膜拜考场上写出来的那些大佬。。。
反正本蒟蒻顶多20分。。。
因为就算想出来了,也没有时间敲+调试了。。。
调试代码毁天灭地,鄙人把它们删掉了。
附上删掉毁天灭地的调试代码之后仍然毁天灭地的代码:
#include<iostream>
#include<algorithm>
#include<cstdio>
#define MAXN 100010
using namespace std;
int n,m=0,size=1,root=0,head;
struct Splay{
int f,s,son[2];
int v,c,key;
}a[MAXN];
struct Spaly{
int f,son[2];
}b[MAXN];
inline int read(){
int date=0,w=1;char c=0;
while(c<‘0‘||c>‘9‘){if(c==‘-‘)w=-1;c=getchar();}
while(c>=‘0‘&&c<=‘9‘){date=date*10+c-‘0‘;c=getchar();}
return date*w;
}
inline void clean(int rt){
a[rt].son[0]=a[rt].son[1]=a[rt].f=a[rt].s=a[rt].v=a[rt].c=0;
b[rt].f=b[rt].son[0]=b[rt].son[1]=0;
}
inline void pushup(int rt){
if(!rt)return;
a[rt].s=a[a[rt].son[0]].s+a[a[rt].son[1]].s+1;
}
inline void pushdown(int rt){
if(!rt||!a[rt].c)return;
if(a[rt].son[0]){
a[a[rt].son[0]].c+=a[rt].c;
a[a[rt].son[0]].v+=a[rt].c;
}
if(a[rt].son[1]){
a[a[rt].son[1]].c+=a[rt].c;
a[a[rt].son[1]].v+=a[rt].c;
}
a[rt].c=0;
}
inline void turn(int rt,int k){
int x=a[rt].f,y=a[x].f;
pushdown(x);pushdown(rt);
a[x].son[k^1]=a[rt].son[k];
if(a[rt].son[k])a[a[rt].son[k]].f=x;
a[rt].f=y;
if(y)a[y].son[a[y].son[1]==x]=rt;
a[x].f=rt;
a[rt].son[k]=x;
pushup(x);pushup(rt);
}
void splay(int rt,int ancestry){
while(a[rt].f!=ancestry){
int x=a[rt].f,y=a[x].f;
pushdown(y);pushdown(x);pushdown(rt);
if(y==ancestry)turn(rt,a[x].son[0]==rt);
else{
int k=a[y].son[0]==x;
if(a[x].son[k]==rt){turn(rt,k^1);turn(rt,k);}
else{turn(x,k);turn(rt,k);}
}
}
if(ancestry==0)root=rt;
}
inline int newnode(int fa,int key){
int rt=size++;
if(fa)a[fa].son[a[fa].key<key]=rt;
a[rt].son[0]=a[rt].son[1]=0;
a[rt].f=fa;a[rt].key=key;
a[rt].s=1;
return rt;
}
int insert(int rt,int key){
int fa=0;
while(rt){
fa=rt;
pushdown(rt);
rt=a[rt].son[a[rt].key<key];
}
rt=newnode(fa,key);
splay(rt,0);
return rt;
}
void find(int rt,int x){
if(!rt)return;
while(a[rt].son[a[rt].key<x])rt=a[rt].son[a[rt].key<x];
splay(rt,0);
}
int kth(int rt,int k){
while(1){
pushdown(rt);
int y=a[rt].son[0];
if(k>a[y].s+1){
k-=a[y].s+1;
rt=a[rt].son[1];
}
else if(k<=a[y].s)rt=y;
else return rt;
}
}
int front_next(int rt,int x,int k){
find(rt,x);
rt=root;
if((a[rt].key>x&&k)||(a[rt].key<x&&k^1))return rt;
rt=a[rt].son[k];
while(a[rt].son[k^1])rt=a[rt].son[k^1];
return rt;
}
int insert_spaly(int rt,int key){
rt=insert(rt,key);
int front=front_next(rt,key,0),next=front_next(rt,key,1),fa;
if(a[front].v>a[next].v)fa=front;
else fa=next;
a[rt].v=a[fa].v+1;
b[rt].f=fa;
if(fa)b[fa].son[a[fa].key<key]=rt;
return a[rt].v;
}
void splay_min(int rt,int k){
rt=kth(rt,k);
printf("%d\n",a[rt].v);
splay(rt,0);
if(b[rt].f)splay(b[rt].f,rt);
a[rt].v=1;
if(head!=rt){
if(b[rt].son[1])b[b[rt].son[1]].f=b[rt].f;
if(b[rt].f)b[b[rt].f].son[0]=b[rt].son[1];
b[head].f=rt;
b[rt].son[1]=head;
b[rt].f=0;head=rt;
if(!a[rt].son[1])return;
int front=a[rt].son[1];
a[front].v+=1;
if(a[front].son[1]){
a[a[front].son[1]].v+=1;
a[a[front].son[1]].c+=1;
}
}
}
void splay_max(int rt,int k){
rt=kth(rt,k);
printf("%d\n",a[rt].v);
splay(rt,0);
if(b[rt].f)splay(b[rt].f,rt);
a[rt].v=1;
if(head!=rt){
if(b[rt].son[0])b[b[rt].son[0]].f=b[rt].f;
if(b[rt].f)b[b[rt].f].son[1]=b[rt].son[0];
b[head].f=rt;
b[rt].son[0]=head;
b[rt].f=0;head=rt;
if(!a[rt].son[0])return;
int front=a[rt].son[0];
a[front].v+=1;
if(a[front].son[0]){
a[a[front].son[0]].v+=1;
a[a[front].son[0]].c+=1;
}
}
}
void work(){
int f,x;
n=read();
for(int cases=1;cases<=n;cases++){
f=read();
switch(f){
case 1:{
x=read();
if(m==0){
head=root=newnode(0,x);
a[root].v=1;
printf("1\n");
}
else printf("%d\n",insert_spaly(root,x));
m++;
break;
}
case 2:{
splay_min(root,1);
break;
}
case 3:{
splay_max(root,m);
break;
}
case 4:{
splay_min(root,1);
m--;
root=a[root].son[1];
head=b[head].son[1];
clean(a[root].f);
a[root].f=b[head].f=0;
if(root){
a[root].v-=1;
a[root].c-=1;
}
break;
}
case 5:{
splay_max(root,m);
m--;
root=a[root].son[0];
head=b[head].son[0];
clean(a[root].f);
a[root].f=b[head].f=0;
if(root){
a[root].v-=1;
a[root].c-=1;
}
break;
}
}
}
}
int main(){
work();
return 0;
}
原文:https://www.cnblogs.com/Yangrui-Blog/p/9463745.html