LCT 的基础 Splay简单食用指南

最新推荐文章于 2020-01-09 16:18:00 发布
置顶 最新推荐文章于 2020-01-09 16:18:00 发布
阅读量174 收藏 1
点赞数 2
分类专栏： ACM Splay
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/heucodesong/article/details/100995973
版权
ACM 同时被 2 个专栏收录
455 篇文章 0 订阅
订阅专栏
Splay
4 篇文章 0 订阅
订阅专栏

我是不会用Splay写平衡树的(真香
为了学LCT 记一份Splay板子

//普通平衡树
/*
    if you can't see the repay
    Why not just work step by step
    rubbish is relaxed
    to ljq
*/
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <cmath>
#include <map>
#include <stack>
#include <set>
#include <sstream>
#include <vector>
#include <stdlib.h>
#include <algorithm>
using namespace std;

#define dbg(x) cout<<#x<<" = "<< (x)<< endl
#define dbg2(x1,x2) cout<<#x1<<" = "<<x1<<" "<<#x2<<" = "<<x2<<endl
#define dbg3(x1,x2,x3) cout<<#x1<<" = "<<x1<<" "<<#x2<<" = "<<x2<<" "<<#x3<<" = "<<x3<<endl
#define max3(a,b,c) max(a,max(b,c))
#define min3(a,b,c) min(a,min(b,c))

typedef pair<int,int> pll;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int _inf = 0xc0c0c0c0;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const ll _INF = 0xc0c0c0c0c0c0c0c0;
const ll mod =  (int)1e9+7;

ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll ksm(ll a,ll b,ll mod){int ans=1;while(b){if(b&1) ans=(ans*a)%mod;a=(a*a)%mod;b>>=1;}return ans;}
ll inv2(ll a,ll mod){return ksm(a,mod-2,mod);}
void exgcd(ll a,ll b,ll &x,ll &y,ll &d){if(!b) {d = a;x = 1;y=0;}else{exgcd(b,a%b,y,x,d);y-=x*(a/b);}}//printf("%lld*a + %lld*b = %lld\n", x, y, d);

/*namespace sgt
{
    #define mid ((l+r)>>1)

    #undef mid
}*/
const int MAX_N = 100025;
int ch[MAX_N][2],fa[MAX_N],cnt[MAX_N],sz[MAX_N],val[MAX_N],rev[MAX_N],root,ncnt;
int chk(int x)
{
    return ch[fa[x]][1] == x;
}
void pushup(int x)
{
    sz[x] = sz[ch[x][0]] + sz[ch[x][1]] + cnt[x];
}
void Rotate(int x)
{
    int y = fa[x],z = fa[y],k = chk(x), w = ch[x][k^1];
    ch[y][k] = w;fa[w] = y;
    ch[z][chk(y)] = x;fa[x] = z;
    ch[x][k^1] = y;fa[y] = x;
    pushup(y);pushup(x);
}
void Splay(int x,int goal = 0)
{
    while(fa[x]!=goal)
    {
        int y = fa[x],z = fa[y];
        if(z!=goal)
        {
            if(chk(x)==chk(y)) Rotate(y);
            else Rotate(x);
        }
        Rotate(x);
    }
    if(!goal) root = x;
}
//辅助操作，将最大的小于等于 x 的数所在的节点splay到根。
void Find(int x)
{
    if(!root) return;
    int cur = root;
    while(ch[cur][x>val[cur]]&&x!=val[cur])
    {
        cur = ch[cur][x>val[cur]];
    }
    Splay(cur);
}
void Insert(int x)
{
    int cur = root,p = 0;
    while(cur&&val[cur]!=x)
    {
        p = cur;
        cur = ch[cur][x>val[cur]];
    }
    if(cur)
    {
        cnt[cur]++;
    }
    else
    {
        cur = ++ncnt;
        if(p) ch[p][x>val[p]] = cur;
        ch[cur][0] = ch[cur][1] = 0;
        val[cur] = x;fa[cur] = p;
        cnt[cur] = sz[cur] = 1;
    }
    Splay(cur);
}
void pushdown(int x)
{
    if(rev[x])
    {
        swap(ch[x][0],ch[x][1]);
        rev[ch[x][0]] ^= 1;
        rev[ch[x][1]] ^= 1;
        rev[x] = 0;
    }
}
int Kth(int k)
{
    int cur = root;
    while(1)
    {
        pushdown(cur);
        if(ch[cur][0]&&k<=sz[ch[cur][0]])
        {
            cur = ch[cur][0];
        }
        else if(k>sz[ch[cur][0]]+cnt[cur])
        {
            k -= sz[ch[cur][0]] + cnt[cur];
            cur = ch[cur][1];
        }
        else return cur;
    }
}
void Reverse(int l,int r)
{
    int x = Kth(l), y = Kth(r+2);
    Splay(x);Splay(y,x);
    rev[ch[y][0]] ^= 1;
}

int pre(int x)
{
    Find(x);
    if(val[root]<x) return root;
    int cur = ch[root][0];
    while(ch[cur][1])
    {
        cur = ch[cur][1];
    }
    return cur;
}
int suc(int x)
{
    Find(x);
    if(val[root]>x) return root;
    int cur = ch[root][1];
    while(ch[cur][0])
    {
        cur = ch[cur][0];
    }
    return cur;
}
void Remove(int x)
{
    int last = pre(x),next = suc(x);
    Splay(last);Splay(next,last);
    int del = ch[next][0];
    if(cnt[del]>1)
    {
        cnt[del]--;
        Splay(del);
    }
    else ch[next][0] = 0;
}
int main()
{
    //ios::sync_with_stdio(false);
    //freopen("a.txt","r",stdin);
    //freopen("b.txt","w",stdout);
    Insert(_inf);
    Insert(inf);
    int n,opt,x;
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d",&opt,&x);
        if(opt==1)
        {
            Insert(x);
        }
        else if(opt==2)
        {
            Remove(x);
        }
        else if(opt==3)
        {
            Find(x);
            printf("%d\n",sz[ch[root][0]]);
        }
        else if(opt==4)
        {
            printf("%d\n",val[Kth(x+1)]);
        }
        else if(opt==5)
        {
            printf("%d\n",val[pre(x)]);
        }
        else
        {
            printf("%d\n",val[suc(x)]);
        }
    }
    //fclose(stdin);
    //fclose(stdout);
    //cout << "time: " << (long long)clock() * 1000 / CLOCKS_PER_SEC << " ms" << endl;
    return 0;
}

//文艺平衡树
您需要写一种数据结构（可参考题目标题），来维护一个有序数列，其中需要提供以下操作：翻转一个区间，例如原有序序列是5 4 3 2 1，翻转区间是[2,4]的话，结果是5 2 3 4 1
/*
    if you can't see the repay
    Why not just work step by step
    rubbish is relaxed
    to ljq
*/
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <cmath>
#include <map>
#include <stack>
#include <set>
#include <sstream>
#include <vector>
#include <stdlib.h>
#include <algorithm>
using namespace std;

#define dbg(x) cout<<#x<<" = "<< (x)<< endl
#define dbg2(x1,x2) cout<<#x1<<" = "<<x1<<" "<<#x2<<" = "<<x2<<endl
#define dbg3(x1,x2,x3) cout<<#x1<<" = "<<x1<<" "<<#x2<<" = "<<x2<<" "<<#x3<<" = "<<x3<<endl
#define max3(a,b,c) max(a,max(b,c))
#define min3(a,b,c) min(a,min(b,c))

typedef pair<int,int> pll;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int _inf = 0xc0c0c0c0;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const ll _INF = 0xc0c0c0c0c0c0c0c0;
const ll mod =  (int)1e9+7;

ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll ksm(ll a,ll b,ll mod){int ans=1;while(b){if(b&1) ans=(ans*a)%mod;a=(a*a)%mod;b>>=1;}return ans;}
ll inv2(ll a,ll mod){return ksm(a,mod-2,mod);}
void exgcd(ll a,ll b,ll &x,ll &y,ll &d){if(!b) {d = a;x = 1;y=0;}else{exgcd(b,a%b,y,x,d);y-=x*(a/b);}}//printf("%lld*a + %lld*b = %lld\n", x, y, d);

/*namespace sgt
{
    #define mid ((l+r)>>1)

    #undef mid
}*/
const int MAX_N = 100025;
int ch[MAX_N][2],fa[MAX_N],cnt[MAX_N],sz[MAX_N],val[MAX_N],rev[MAX_N],root,ncnt,n;
int chk(int x)
{
    return ch[fa[x]][1] == x;
}
void pushup(int x)
{
    sz[x] = sz[ch[x][0]] + sz[ch[x][1]] + cnt[x];
}
void Rotate(int x)
{
    int y = fa[x],z = fa[y],k = chk(x), w = ch[x][k^1];
    ch[y][k] = w;fa[w] = y;
    ch[z][chk(y)] = x;fa[x] = z;
    ch[x][k^1] = y;fa[y] = x;
    pushup(y);pushup(x);
}
void Splay(int x,int goal = 0)
{
    while(fa[x]!=goal)
    {
        int y = fa[x],z = fa[y];
        if(z!=goal)
        {
            if(chk(x)==chk(y)) Rotate(y);
            else Rotate(x);
        }
        Rotate(x);
    }
    if(!goal) root = x;
}
//辅助操作，将最大的小于等于 x 的数所在的节点splay到根。
void Find(int x)
{
    if(!root) return;
    int cur = root;
    while(ch[cur][x>val[cur]]&&x!=val[cur])
    {
        cur = ch[cur][x>val[cur]];
    }
    Splay(cur);
}
void Insert(int x)
{
    int cur = root,p = 0;
    while(cur&&val[cur]!=x)
    {
        p = cur;
        cur = ch[cur][x>val[cur]];
    }
    if(cur)
    {
        cnt[cur]++;
    }
    else
    {
        cur = ++ncnt;
        if(p) ch[p][x>val[p]] = cur;
        ch[cur][0] = ch[cur][1] = 0;
        val[cur] = x;fa[cur] = p;
        cnt[cur] = sz[cur] = 1;
    }
    Splay(cur);
}
void pushdown(int x)
{
    if(rev[x])
    {
        swap(ch[x][0],ch[x][1]);
        rev[ch[x][0]] ^= 1;
        rev[ch[x][1]] ^= 1;
        rev[x] = 0;
    }
}
int Kth(int k)
{
    int cur = root;
    while(1)
    {
        pushdown(cur);
        if(ch[cur][0]&&k<=sz[ch[cur][0]])
        {
            cur = ch[cur][0];
        }
        else if(k>sz[ch[cur][0]]+cnt[cur])
        {
            k -= sz[ch[cur][0]] + cnt[cur];
            cur = ch[cur][1];
        }
        else return cur;
    }
}
void Reverse(int l,int r)
{
    int x = Kth(l), y = Kth(r+2);
    Splay(x);Splay(y,x);
    rev[ch[y][0]] ^= 1;
}

int pre(int x)
{
    Find(x);
    if(val[root]<x) return root;
    int cur = ch[root][0];
    while(ch[cur][1])
    {
        cur = ch[cur][1];
    }
    return cur;
}
int suc(int x)
{
    Find(x);
    if(val[root]>x) return root;
    int cur = ch[root][1];
    while(ch[cur][0])
    {
        cur = ch[cur][0];
    }
    return cur;
}
void Remove(int x)
{
    int last = pre(x),next = suc(x);
    Splay(last);Splay(next,last);
    int del = ch[next][0];
    if(cnt[del]>1)
    {
        cnt[del]--;
        Splay(del);
    }
    else ch[next][0] = 0;
}
void Output(int x)
{
    pushdown(x);
    if(ch[x][0]) Output(ch[x][0]);
    if(val[x]&&val[x]<=n&&val[x]>=1) printf("%d ",val[x]);
    if(ch[x][1]) Output(ch[x][1]);
}
int main()
{
    //ios::sync_with_stdio(false);
    //freopen("a.txt","r",stdin);
    //freopen("b.txt","w",stdout);
    Insert(_inf);
    Insert(inf);
    int opt,x,Q,y;
    scanf("%d%d",&n,&Q);
    for(int i = 1;i<=n;++i) Insert(i);
    while(Q--)
    {
        scanf("%d%d",&x,&y);
        Reverse(x,y);
    }
    Output(root);
    //fclose(stdin);
    //fclose(stdout);
    //cout << "time: " << (long long)clock() * 1000 / CLOCKS_PER_SEC << " ms" << endl;
    return 0;
}
ljq真的不能再吃了
关注
  • 2
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
splay&LCT学习笔记
White_gugu的博客
01-25 250
伸展树(splay)和Link-Cut Tree
SPLAY，LCT学习笔记（一）
lleozhang的博客
09-15 383
写了两周数据结构，感觉要死掉了，赶紧总结一下，要不都没学明白。 SPLAY专题： 例：NOI2005 维修数列 典型的SPLAY问题，而且综合了SPLAY常见的所有操作，特别适合新手入门学习（比如我这种蒟蒻） 题目要求很多，我们一步一步来分析 首先，区间翻转是SPLAY一个很基础的操作，我们以他为基础分析这个SPLAY 例：luogu文艺平衡树 题目：读入一个序列，进行多次区间翻转操...
SPLAY，LCT学习笔记（四）
lleozhang的博客
09-18 343
前三篇好像变成了SPLAY专题... 这一篇正式开始LCT！ 其实LCT就是基于SPLAY的伸展操作维护树（森林）连通性的一个数据结构 核心操作有很多，我们以一道题为例： 例：bzoj 2049 洞穴勘测 要求：加边和删边，询问连通性 其实如果没有删边，裸跑并查集似乎就可以搞定 但由于存在删边，并查集思想受阻，我们要考虑更高级的数据结构 于是LCT横空出世！ LCT的核心思想：多...
Splay与LCT学习笔记
hanjinbo的博客
04-14 264
Splay 1:rotate1:rotate1:rotate操作无论是左边旋到右边还是右边旋到左边，节点的相对高度都上升了。 2:2:2:单旋splaysplaysplay可能会被卡掉，一条链单旋过后还是一条链。 3:3:3:双旋分一下几种情况：(x−&amp;amp;amp;gt;to)(x-&amp;amp;amp;gt;to)(x−&amp;amp;gt;to) tototo是xxx的fatherfatherfather，此时直接把x...
Splay和LCT的复杂度分析
weixin_30336061的博客
12-08 318
\(Splay\)的复杂度分析 不论插入，删除还是访问，我们可以发现它们的复杂度都和\(splay\)操作的复杂度同阶，只是一点常数的区别 我们不妨假设有\(n\)个点的\(splay\)，进行了\(m\)次\(splay\)操作 采用势能分析 我们记\(w(x) = \left \lceil \log_2 (size(x)) \right \rceil\)，注意以\(2\)为底和上取整 我们定...
[BZOJ]3159 决战 LCT&Splay
Le Petit Prince
11-28 700
3159: 决战 Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 420  Solved: 183 [Submit][Status][Discuss] Description Input 第一行有三个整数N、M和R，分别表示树的节点数、指令和询问总数，以及X国的据点。 接下来N-1行，每行两个整数X和Y，表示Kat
BZOJ 3159: 决战 LCT+Splay
EM LGH
01-09 99
挺厉害的一道大数据结构题. 由于 LCT 是维护树的形态的，所以说不支持翻转操作. 而在维护序列时 splay 是支持区间翻转的. 所以，我们对于 LCT 中每一个重链都维护一个 splay(这个不同于 LCT 中的 splay) 由于重链是一个序列，所以是支持序列的区间翻转的. 那么我们的翻转，链加和，链求和的操作就都在这个重链对应的 splay 上进行...
lct.zip_LCT
09-19
【标题】"lct.zip_LCT" 涉及到的是一个名为 "lct.zip" 的压缩文件，其中包含了与"LCT"相关的资料。LCT通常代表“逻辑控制技术”（Logical Control Technology）或者“线性电流变换器”（Linear Current Transformer...
lct.zip_lct文件
09-24
"流程图"文件则为使用者提供了直观的操作指南，帮助他们更好地理解和利用“lct”库。为了充分利用这个库，开发者需要了解其API，参考流程图来正确地导入和使用库中的功能，从而提升开发效率和代码质量。在实际应用中...
lct_LCT_
09-29
标题中的“lct_LCT_”很可能代表“Link-Cut Tree”，这是一项在计算机科学领域，特别是算法设计中广泛使用的数据结构。Link-Cut Tree（链接切割树）是一种自顶向下的动态树数据结构，主要用于解决与树相关的操作，如...
U2000 Web LCT 用户指南-(V100R005C00_01).
05-23
U2000 Web LCT 用户指南-(V100R005C00_01).
U2000 LCT 用户指南-(V100R005C00_01).
05-23
U2000 LCT 用户指南-(V100R005C00_01).
GDAL-2.4.1-cp36-cp36m-win_amd64.whl
最新发布
07-17
GDAL-2.4.1-cp36-cp36m-win_amd64.whl
注射用七叶皂苷钠说明书.doc
07-17
说明书
docker配置资料.txt
07-17
docker安装
2024年欧洲氯化聚氯乙烯市场主要企业市场占有率及排名.docx
07-17
2024年欧洲氯化聚氯乙烯市场主要企业市场占有率及排名.docx
H公司-职业与任职资格管理.pdf
07-17
H公司-职业与任职资格管理.pdf
平衡计分卡与战略地图.docx
07-17
平衡计分卡与战略地图.docx
帮我写一个c++的LCT
05-26
好的，下面是一个简单的 C++ LCT 实现，包含了 LCT 基本操作： ```cpp #include <bits/stdc++.h> using namespace std; const int N = 1e5 + 5; // 节点结构体 struct Node { int ch[2], fa, rev; int val, sum; } tr[N]; // 判断节点x是否是其父节点的左儿子 bool is_left(int x) { return tr[tr[x].fa].ch[0] == x; } // 反转标记 void reverse(int x) { swap(tr[x].ch[0], tr[x].ch[1]); tr[x].rev ^= 1; } // 维护反转标记 void pushdown(int x) { if (tr[x].rev) { reverse(tr[x].ch[0]); reverse(tr[x].ch[1]); tr[x].rev = 0; } } // 维护sum值 void pushup(int x) { tr[x].sum = tr[tr[x].ch[0]].sum + tr[tr[x].ch[1]].sum + tr[x].val; } // 旋转操作 void rotate(int x) { int y = tr[x].fa, z = tr[y].fa; int k = is_left(x), w = tr[x].ch[k ^ 1]; tr[y].ch[k] = w; tr[w].fa = y; tr[z].ch[is_left(y)] = x; tr[x].fa = z; tr[x].ch[k ^ 1] = y; tr[y].fa = x; pushup(y); pushup(x); } // 递归下传反转标记 void splay(int x) { static int stk[N]; int top = 0, y = x; stk[++top] = y; while (!stk[top]) { y = tr[y].fa; stk[++top] = y; } while (top) pushdown(stk[top--]); while (is_left(x)) { int y = tr[x].fa, z = tr[y].fa; if (is_left(y)) rotate(y); rotate(x); } while (!is_left(x)) { int y = tr[x].fa, z = tr[y].fa; if (!is_left(y)) rotate(y); rotate(x); } } // 将x节点到根节点的路径变为一条链 void access(int x) { for (int y = 0; x; y = x, x = tr[x].fa) { splay(x); tr[x].ch[1] = y; pushup(x); } } // 将x节点所在的树变为splay树的根节点 void make_root(int x) { access(x); splay(x); reverse(x); } // 查询x节点所在的树的根节点 int find_root(int x) { access(x); splay(x); while (tr[x].ch[0]) { pushdown(x); x = tr[x].ch[0]; } splay(x); return x; } // 连接x和y两个节点 void link(int x, int y) { make_root(x); tr[x].fa = y; } // 断开x和y两个节点 void cut(int x, int y) { make_root(x); if (find_root(y) == x && tr[y].fa == x && !tr[y].ch[0]) { tr[y].fa = tr[x].ch[1] = 0; pushup(x); } } // 修改x节点的值为val void modify(int x, int val) { splay(x); tr[x].val = val; pushup(x); } // 查询x节点到y节点路径上的所有节点的权值和 int query(int x, int y) { make_root(x); access(y); splay(y); return tr[y].sum; } int main() { int n, m; scanf("%d%d", &n, &m); for (int i = 1; i <= n; ++i) { scanf("%d", &tr[i].val); tr[i].sum = tr[i].val; } while (m--) { int opt, x, y; scanf("%d%d%d", &opt, &x, &y); if (opt == 0) printf("%d\n", query(x, y)); else if (opt == 1) link(x, y); else if (opt == 2) cut(x, y); else if (opt == 3) modify(x, y); } return 0; } ``` 这里只是一个简单的 LCT 实现，对于一些细节和优化并没有处理。如果需要了解更多关于 LCT 的内容，可以参考一些经典题目，如 P3834、P3836、P4219 等。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值