• splay区间翻转


    原题P3391 【模板】文艺平衡树(Splay)

    题目背景
    这是一道经典的Splay模板题——文艺平衡树。
    
    题目描述
    您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作:翻转一个区间,例如原有序序列是5 4 3 2 1,翻转区间是[2,4]的话,结果是5 2 3 4 1
    
    输入输出格式
    输入格式:
    
     
    
    第一行为n,m n表示初始序列有n个数,这个序列依次是(1,2,⋯n−1,n) m表示翻转操作次数
    
    接下来m行每行两个数 [l,r] 数据保证1≤l≤r≤n
    
     
    
    输出格式:
    
     
    
    输出一行n个数字,表示原始序列经过m次变换后的结果
    
     
    
    输入输出样例
    输入样例#1: 复制
    5 3
    1 3
    1 3
    1 4
    输出样例#1: 复制
    4 3 2 1 5
    说明
    n,m≤100000
    查看题面

    题解

    这道题目要求支持区间翻转,这是splay的典型操作。

    基本思想:

    • 只要将一棵BST每一个节点的左右子树交换,就相当于翻转次树的中序遍历

    实现方法:

    1. 用节点中序遍历的序号表示当前节点的下标,用key存节点的序号,经过一系列操作后中序遍历即为所求

    2. 每次翻转区间【L,R】时只要将排名为L-1的点splay到根,将排名R+1的点splay到根的右节点,根的右子树的左子树的中序遍历就是【L,R】

    如何将一个点splay到指定位置?

    只要引入y标记,当前节点的父亲为y时跳出,就可以将x刚好splay到y的儿子

    inline void splay(int x,int y){
        for(int fa;(fa=f[x])!=y;rotate(x)){
            if(f[fa]!=y){
                rotate(get(fa)==get(x)?fa:x);
            }
        }
        if(y==0){
            root=x;
        }
    }

    1.标记是在每一次访问到一个新的节点是就要pushdown的(改变树的结构会破坏标记区间,所以先一步下传标记) 

    inline void push_down(int x){
        if(x&&tag[x]){
            tag[son[x][0]]^=1;
            tag[son[x][1]]^=1;
            swap(son[x][0],son[x][1]);
            tag[x]=0;
        }
    }

    2.区分一个节点的排名和这个节点的值:这个节点的排名是它是当前数组中的第几个,用左儿子的size+1表示;这个节点的值是题目中输入的数字,在本题中是1~N

    3.增加数字为1和N+2的两个哨兵节点+INF和-INF,因此第i个节点的下标为i+1。

    4.难道交换左右子树不会破坏BST的性质吗?这就是容易混淆的一点,我们的区间操作是根据下标翻转的,用数组实现时下标就是数组地址,子树交换时,只是存储内容的改变,下标位置(树的形状)只会在旋转时改变,保证BST性质。

    5.此题要求一次性插入N个节点,所以可以用线段树建树的方法O(N)建树

    int init(int l,int r,int fa){
        if(l>r){
            return 0;
        }
        int mid=(l+r)>>1,cur=++sz;
        key[cur]=val[mid];
        f[cur]=fa;
        son[cur][0]=init(l,mid-1,cur);
        son[cur][1]=init(mid+1,r,cur);
        upd(cur);
        return cur;
    }

    完整代码

      1 #include<bits/stdc++.h>
      2 using namespace std;
      3 typedef long long LL;
      4 const int INF=1e9+7,MAXN=1e5+5,MAXM=1e5+5;
      5 int N,M,root,sz;
      6 int key[MAXN],val[MAXN],f[MAXN],son[MAXN][2],tag[MAXN],siz[MAXN];
      7 inline void DEBUG(){
      8     printf("root=%d siz=%d
    ",root,sz);
      9     for(int i=1;i<=sz;i++){
     10         printf("(idx=%d	key=%d	lc=%d	rc=%d)
    ",i,key[i],son[i][0],son[i][1]);
     11     }
     12     puts("-------------------------------------------------------");
     13 }
     14 inline void upd(int x){
     15     if(x){
     16         siz[x]=1;
     17         if(son[x][0]){
     18             siz[x]+=siz[son[x][0]];
     19         }
     20         if(son[x][1]){
     21             siz[x]+=siz[son[x][1]];
     22         }
     23     }
     24 }
     25 inline void push_down(int x){
     26     if(x&&tag[x]){
     27         tag[son[x][0]]^=1;
     28         tag[son[x][1]]^=1;
     29         swap(son[x][0],son[x][1]);
     30         tag[x]=0;
     31     }
     32 }
     33 inline int get(int x){
     34     return x==son[f[x]][1];
     35 }
     36 inline void rotate(int x){
     37     int fa=f[x],gf=f[fa],which=get(x);
     38     son[fa][which]=son[x][which^1];
     39     f[son[fa][which]]=fa;
     40     son[x][which^1]=fa;
     41     f[fa]=x;
     42     f[x]=gf;
     43     if(gf){
     44         son[gf][son[gf][1]==fa]=x;
     45     }
     46     upd(fa);
     47     upd(x);
     48 }
     49 inline void splay(int x,int y){
     50     for(int fa;(fa=f[x])!=y;rotate(x)){
     51         if(f[fa]!=y){
     52             rotate(get(fa)==get(x)?fa:x);
     53         }
     54     }
     55     if(y==0){
     56         root=x;
     57     }
     58 }
     59 int init(int l,int r,int fa){
     60     if(l>r){
     61         return 0;
     62     }
     63     int mid=(l+r)>>1,cur=++sz;
     64     key[cur]=val[mid];
     65     f[cur]=fa;
     66     son[cur][0]=init(l,mid-1,cur);
     67     son[cur][1]=init(mid+1,r,cur);
     68     upd(cur);
     69     return cur;
     70 }
     71 inline int findx(int x){//return the index
     72     int cur=root;
     73     while(1){
     74         push_down(cur);
     75         if(son[cur][0]&&x<=siz[son[cur][0]]){
     76             cur=son[cur][0];
     77         }else{
     78             x-=siz[son[cur][0]]+1;
     79             if(!x){
     80                 return cur;
     81             }
     82             cur=son[cur][1];
     83         }
     84     }
     85 }
     86 void dfs(int x){
     87     push_down(x);
     88     if(son[x][0]){
     89         dfs(son[x][0]);
     90     }
     91     if(key[x]!=INF&&key[x]!=-INF){
     92         printf("%d ",key[x]);
     93     }
     94     if(son[x][1]){
     95         dfs(son[x][1]);
     96     }
     97 }
     98 int main(){
     99     scanf("%d%d",&N,&M);
    100     for(int i=2;i<=N+1;i++){
    101         val[i]=i-1;
    102     }
    103     val[1]=-INF;
    104     val[N+2]=INF;
    105     root=init(1,N+2,0);
    106     for(int i=1;i<=M;i++){
    107         int ii,jj;
    108         scanf("%d%d",&ii,&jj);
    109         int l=findx(ii),r=findx(jj+2);
    110         splay(l,0);
    111         splay(r,l);
    112         tag[son[son[root][1]][0]]^=1;
    113     }
    114     dfs(root);
    115     return 0;
    116 }
  • 相关阅读:
    微服务-Nacos
    微服务RPC框架-Feign
    读书笔记之《大型分布式网站架构设计与实践》
    读书笔记之《Java并发编程的艺术》—— 四
    读书笔记之《Java并发编程的艺术》—— 三
    读书笔记之《Java并发编程的艺术》—— 二
    读书笔记之《Java并发编程的艺术》—— 一
    noip2021 训练4 做题记录
    noip 训练 (线段树专项)
    UVA11671 矩阵中的符号 Sign of Matrix 题解
  • 原文地址:https://www.cnblogs.com/guoshaoyang/p/10631589.html
Copyright © 2020-2023  润新知