ZOJ 3765 Lights Splay Tree的几种基本操作

最新推荐文章于 2018-08-17 09:48:24 发布
最新推荐文章于 2018-08-17 09:48:24 发布
阅读量1k 收藏
点赞数
分类专栏： 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zmx354/article/details/27229937
版权

区间的查询，点修改，插入和删除。先姑且当作模板吧，略挫，慢慢补充，慢慢优化。

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <stack>
#include <map>

#pragma comment(linker, "/STACK:1024000000");
#define EPS (1e-8)
#define LL long long
#define ULL unsigned long long
#define _LL __int64
#define _INF 0x3f3f3f3f
#define Mod 9999991

using namespace std;

const int MAXN = 700100;

struct N
{
    //info
    int son[2],pre;

    //data
    int sta,ls,rs,s;
    int w,g[2];
    bool mark[2];

}st[MAXN];

struct W
{
    int sta,w;
}num[200100];

void Output(int root)
{
    if(root == -1)
        return ;
    Output(st[root].son[0]);
    printf("root = %d pre = %d data = %2d lg = %2d rg = %2d sta = %2d ls = %d rs = %d s = %d\n",root,st[root].pre,st[root].w,st[root].g[0],st[root].g[1],st[root].sta,st[root].ls,st[root].rs,st[root].s);
    Output(st[root].son[1]);
}

int Top;

int Gcd(int a,int b)
{
    int t;

    while(b)
    {
        t = a;
        a = b;
        b = t%b;
    }
    return a;
}

//更新st[root]的数据域
void Updata(int root)
{
    st[root].ls = 0,st[root].rs = 0;

    st[root].g[st[root].sta] = st[root].w;
    st[root].g[1^st[root].sta] = 1;
    st[root].mark[st[root].sta] = true;
    st[root].mark[1^st[root].sta] = false;

    if(st[root].son[0] != -1)
    {
        if(st[root].mark[0] == true)
            st[root].g[0] = (st[st[root].son[0]].mark[0] ?  Gcd(st[root].g[0],st[st[root].son[0]].g[0]) : st[root].g[0]);
        else
            st[root].g[0] = st[st[root].son[0]].g[0];
        if(st[root].mark[1] == true)
            st[root].g[1] = (st[st[root].son[0]].mark[1] ?  Gcd(st[root].g[1],st[st[root].son[0]].g[1]) : st[root].g[1]);
        else
            st[root].g[1] = st[st[root].son[0]].g[1];

        st[root].mark[0] = (st[root].mark[0] || st[st[root].son[0]].mark[0]);
        st[root].mark[1] = (st[root].mark[1] || st[st[root].son[0]].mark[1]);
        st[root].ls = st[st[root].son[0]].s;
    }

    if(st[root].son[1] != -1)
    {
        if(st[root].mark[0] == true)
            st[root].g[0] = (st[st[root].son[1]].mark[0] ?  Gcd(st[root].g[0],st[st[root].son[1]].g[0]) : st[root].g[0]);
        else
            st[root].g[0] = st[st[root].son[1]].g[0];
        if(st[root].mark[1] == true)
            st[root].g[1] = (st[st[root].son[1]].mark[1] ?  Gcd(st[root].g[1],st[st[root].son[1]].g[1]) : st[root].g[1]);
        else
            st[root].g[1] = st[st[root].son[1]].g[1];

        st[root].mark[0] = (st[root].mark[0] || st[st[root].son[1]].mark[0]);
        st[root].mark[1] = (st[root].mark[1] || st[st[root].son[1]].mark[1]);
        st[root].rs = st[st[root].son[1]].s;
    }

    st[root].s = st[root].ls + st[root].rs + 1;
}

void Init(int &root,int s,int e,int pre)
{
    if(s > e)
        return ;

    int mid = (s+e)>>1;

    root = Top++;
    st[root].w = num[mid].w;
    st[root].sta = num[mid].sta;

    st[root].pre = pre;
    st[root].son[0] = st[root].son[1] = -1;

    Init(st[root].son[0],s,mid-1,root);
    Init(st[root].son[1],mid+1,e,root);

    Updata(root);
}

//dir == 0,root为pre的左儿子，dir == 1时为右儿子。
void Rotate(int root,int dir)
{
    //旋转操作
    st[st[root].pre].son[dir] = st[root].son[1^dir];
    st[root].son[1^dir] = st[root].pre;

    if(st[st[st[root].pre].pre].son[0] == st[root].pre)
        st[st[st[root].pre].pre].son[0] = root;
    else
        st[st[st[root].pre].pre].son[1] = root;
    int temp = st[root].pre;
    st[root].pre = st[st[root].pre].pre;
    st[temp].pre = root;

    if(st[temp].son[dir] != -1)
        st[st[temp].son[dir]].pre = temp;

    //更新数据域
    Updata(temp);
    Updata(root);
}

//将root旋转到goal下面,goal == 0时即表示将root旋转至根节点位置。
int Splay(int root,int goal)
{
    while(st[root].pre != goal)
    {
        Rotate(root,(st[st[root].pre].son[0] == root ? 0 : 1));
    }

    return root;
}

int Search_Site(int root,int site)
{
    while(st[root].ls + 1 != site)
    {
        if(st[root].ls + 1 < site)
        {
            site = site - st[root].ls - 1;
            root = st[root].son[1];
        }
        else
        {
            root = st[root].son[0];
        }
    }

    return root;
}

void Increase(int root,int site,int w,int sta,int &R)
{
    root = Search_Site(root,site);

    R = Splay(root,0);

    //在根节点的右子树的最左边插入一个点

    int tr = st[R].son[1];
    st[R].son[1] = -1;

    int New = Top++;

    st[New].w = w;
    st[New].sta = sta;

    st[New].pre = 0;
    st[0].son[0] = New;
    st[R].pre = New;
    st[New].son[0] = R;
    st[tr].pre = New;
    st[New].son[1] = tr;

    Updata(R);
    Updata(New);
    R = New;
}

void Link(int &root,int tr,int pre)
{
    if(root == -1)
    {
        root = tr;
        st[root].pre = pre;
    }
    else
    {
        Link(st[root].son[1],tr,root);
    }
    Updata(root);
}

//不存在删除所有节点的情况
void Decrease(int root,int site,int &R)
{
    root = Search_Site(root,site);
    //Output(R);
   // printf("D site = %d root = %d R = %d\n",site,root,R);

    R = Splay(root,0);
    //删除根节点

    R = st[root].son[0];
    st[R].pre = st[root].pre;
    st[0].son[0] = R;
    int tr = st[root].son[1];

    //将tr接到R的最右边。

    Link(R,tr,root);
    Updata(R);
}

void Reinit(int root,int site,int &R)
{
    root = Search_Site(root,site);

    R = Splay(root,0);

    st[R].sta ^= 1;
    Updata(R);
}

void Modify(int root,int site,int w,int &R)
{
    root = Search_Site(root,site);

    R = Splay(root,0);

    st[R].w = w;
    Updata(R);
}

int Query(int root,int l,int r,int sta,int &R)
{
    int lsite = Search_Site(root,l);
    int rsite = Search_Site(root,r);

    R = Splay(rsite,0);

    R = Splay(lsite,0);

    if(st[st[R].son[1]].son[0] == -1 || st[st[st[R].son[1]].son[0]].mark[sta] == false)
        return -1;
    return st[st[st[R].son[1]].son[0]].g[sta];
}

int main()
{
    //freopen("data.txt","r",stdin);
    //freopen("dataout.txt","w",stdout);
    int n,m;

    int i;

    int root;

    while(scanf("%d %d",&n,&m) != EOF)
    {
        for(i = 2;i <= n+1; ++i)
        {
            scanf("%d %d",&num[i].w,&num[i].sta);
        }

       // cout<<n<<" "<<m<<endl;

        num[1].w = 1,num[1].sta = 0;
        num[n+2].w = 1,num[n+2].sta = 0;

        //Top从1开始是为了将st[0]拿出来做
        //根节点的虚拟父节点，在缩短代码的同时防止RE。
        Top = 1,root = -1;

        Init(root,1,n+2,0);
      //  cout<<"dfadsf"<<endl;
        st[0].son[0] = root,st[0].son[1] = -1;

        char order[2];

        int site,sta,l,r,w;

        for(i = 1;i <= m; ++i)
        {
            scanf("%s",order);

            if(order[0] == 'I')
            {
                scanf("%d %d %d",&site,&w,&sta);
               // printf("i = %d site = %d\n",i+200000+1,site);
                Increase(root,site+1,w,sta,root);
            }
            else if(order[0] == 'D')
            {
                scanf("%d",&site);
               // printf("i = %d site = %d\n",i+200000+1,site);
                Decrease(root,site+1,root);
            }
            else if(order[0] == 'R')
            {
                scanf("%d",&site);
                Reinit(root,site+1,root);
            }
            else if(order[0] == 'M')
            {
                scanf("%d %d",&site,&w);
                Modify(root,site+1,w,root);
            }
            else if(order[0] == 'Q')
            {
                scanf("%d %d %d",&l,&r,&sta);
                printf("%d\n",Query(root,l,r+2,sta,root));
            }
            //cout<<"ans = "<<st[root].s<<endl;
        }

    }
    return 0;
}

自在_飞花
关注
专栏目录
ZOJ-3765 Lights（Splay树）
十分残念的博客
07-24 296
Lights Time Limit: 8 Seconds      Memory Limit: 131072 KB Now you have N lights in a line. Don't worry - the lights don't have color. The only status they have is on and off. And, each ligh
zoj3765
a5199519的博客
01-03 102
题解： splay维护 注意是gcd 代码： #include<bits/stdc++.h> using namespace std; #define Key_value ch[ch[root][1]][0] const int N=300010; char op[100]; int pre[N],ch[N][2],key[N],size[N],root...
Lights ZOJ - 3765
sunyutian1998的博客
08-17 171
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3765 模板题 熟悉基本操作 每次收口袋之后 记得把口袋里的区间splay到根节点 否则慢的出奇 #include &lt;bits/stdc++.h&gt; using namespace std; int getgcd(int a,int b); struct ...
ZOJ 3765 splay
u013268685的专栏(有一种幸福叫AC，有一种期待叫AK)
12-01 1154
B - Lights Time Limit:8000MS     Memory Limit:131072KB     64bit IO Format:%lld & %llu Submit Status Description Now you have N lights in a line. Don't worry - the lights don't have co
【splay tree】 ZOJ 3765 Lights
yysys
08-02 553
有五种操作，
zoj 1002_zoj1002_
09-29
在深入学习和研究这个"zoj 1002"程序时，我们需要关注以下几个方面： 1. **输入输出格式**：ACM题目通常有严格的输入输出格式，程序需要正确处理标准输入和输出，或者读取文件数据。 2. **算法设计**：程序的核心...
zoj.gz_ ZOJ_ZOJ 1016_max flow_zoj 1045_zoj.rar
09-24
标题中的"ZOJ.gz_ ZOJ_ZOJ 1016_max flow_zoj 1045_zoj.rar" 提到了两个ZOJ（Zhejiang Online Judge）的题目，分别是1016和1045，这两个数字通常代表在线编程竞赛中的题目编号。这些题目通常涉及到算法和数据结构的...
zoj.rar_zoj_zoj4041
09-23
《ZOJ 4041问题的正确解法与程序分析》 ZOJ（Zhejiang Online Judge）是一个知名的在线编程竞赛平台，其中的题目编号为4041的题目吸引了众多程序员的关注。本篇文章将深入探讨ZOJ 4041的正确解法，并对提供的源代码...
ZOJ1014.zip_zoj code_zoj1004
09-23
标题“ZOJ1014.zip_zoj code_zoj1004”表明这是一个与ZOJ（ZeroJudge）在线判题系统相关的代码压缩包，其中可能包含了解决ZOJ问题1004的源代码。ZOJ是面向编程爱好者和学生的一个在线编程竞赛平台，它提供了各种算法...
ZOJ.zip_Jugs A_ZOJ NTA_zoj acm_zoj acm 1216_zoj code
09-20
【ZOJ.zip】是一个压缩包，里面包含了与ZOJ（Zhejiang Online Judge）相关的ACM（International Collegiate Programming Contest）题解。ZOJ是一个在线编程竞赛平台，它为参赛者提供了一系列算法题目进行练习，以...
ZOJ 3765 —— Lights（伸展树）
hongrock的专栏
03-03 1052
题目：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3765 昨天比赛的时候，别人都把这题KO我却只能干瞪眼，虽然感觉是伸展树，可是之前一直没真正去学这东西。 事后学了下这个东西马上回去ZOJ把这题QJ了！ 伸展树（Splay Tree）这东西太神奇啦！ 我是看cxlove大牛的文章学习的：http://blog.
ZOJ 3765 Lights (伸展树)
JayYe的专栏
03-03 1787
伸展树裸题，做了几题后已经不需要调了，像前几次经常调个好久。。听说可持久化treap更好用，调起来方便，下次学一下。 code： #include #include using namespace std; #define lson x->ch[0] #define rson x->ch[1] #define ket (root->ch[1]->ch[0]) con
伸展树zoj3765
龙崎的专栏
12-02 482
Lights Time Limit: 8 Seconds      Memory Limit: 131072 KB Now you have N lights in a line. Don't worry - the lights don't have color. The only status they have is on and off. And, each ligh
zoj 3765 Lights（Splay）
ok_again的专栏
03-05 1061
由于数组开小了，一直TLE。。。检查了一下午，都要哭了，题目没啥好说的，就是一个入门级别的Splay Tree。只有插入，删除，修改操作。 #include #include #include #include #include #define LL long long #define CLR(a, b) memset(a, b, memset(a)) using namespace std;
zoj 3765 Lights(伸展树)
绝望的乐园
03-02 1778
题意：最开始有N个灯，有Q个操作。操作有一下几种类型：①查询区间[L,R]中灯的状态为status的所有灯的gcd②在第i位后插入一个值为val，状态为status的灯③删除第i个灯④第i个灯是开的那么关掉，否则打开⑤修改第i个灯的值。 思路：用伸展树暴力点写就行了。用数组记录当前区间状态为0和1的gcd的值。对于0结点，这些值赋值为-1就行了，在gcd中修改一下，对-1进行处理就行了。剩下的都
ZOJ 3765 Lights （SPLAY）
Jinx Jinx again
04-15 861
gcd[x] =GCD( GCD( gcd[ls],gcd[rs]),val[x]) ); 注意判断没有亮灯状态的情况 #include #include #define inf 0x3f3f3f3f #define maxn 222222 #define keyTree (ch[ch[root][1]][0]) using namespace std; type
浅谈LCA的几种算法
自在飞花
01-10 9312
LCA,Lowest Common Ancestor,最近的公共祖先。在一棵树中对于两个节点u , v找出节点T，使得T同时为u，v的祖先。显然这样的T点肯定存在且有可能有多个，其中深度最大的那个点肯定为即为u，v两点的LCA。关于LCA的解法有很多种，暴力枚举，事先需要知道所有询问的离线的tarjan算法和基于RMQ的在线算法，下面说一下自己对这种几种算法的理解。   ⒈最容易想到的暴力搜索
浅谈二维中的树状数组与线段树
自在飞花
06-17 1447
一般来说，树状数组可以实现的东西线段树均可胜任，实际应用中也是如此。但是在二维中，线段树的操作变得太过复杂，更新子矩阵时第一维的lazy标记更是麻烦到不行。 但是树状数组在某些询问中又无法胜任，如最值等不符合区间减法的询问。此时就需要根据线段树与树状数组的优缺点来选择了。 做一下基本操作的对比，如下图。 因为线段树为自上向下更新，从而可以使用lazy标记使得矩阵的更新变的高校起来，几个不
二维树状数组的区间加减及查询 tyvj 1716 上帝造题的七分钟
自在飞花
06-17 1390
具体解释见小结。 #include #include #include #include #include #include #include #include #include #pragma comment(linker, "/STACK:1024000000"); #define EPS (1e-8) #define LL long long #define ULL u
Zoj4041题解及运行程序代码分析
资源摘要信息:"zoj.rar_zoj_zoj4041" 知识点: 1. ZOJ简介：ZOJ是浙江大学在线评测系统的缩写，它是一个为程序设计竞赛提供在线自动评测的系统，广泛应用于学校内或者各编程竞赛中。用户可以通过提交自己的代码来...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值