BZOJ 3786: 星系探索 欧拉游览树
一个叫 Euler-Tour-Tree 的数据结构,说白了就是用 Splay_Tree 维护欧拉序
#include <cstring> #include <algorithm> #include <string> #include <cstdio> using namespace std; void setIO(string a){ freopen((a+".in").c_str(),"r",stdin); freopen((a+".out").c_str(),"w",stdout); } #define maxn 300000 #define ll long long int euler[maxn], w[maxn], cnt=1; int head[maxn],to[maxn],nex[maxn],edges,root; ll sumv[maxn], val[maxn], lazy[maxn]; void addedge(int u,int v){ nex[++edges]=head[u],head[u]=edges,to[edges]=v; } void dfs(int u){ euler[++cnt]=u*2, val[u*2]=w[u]; for(int v=head[u];v;v=nex[v]) dfs(to[v]); euler[++cnt]=u*2+1, val[u*2+1]=-w[u]; } struct Splay_Tree{ int f[maxn],ch[maxn][2],sta[maxn],siz[maxn]; int rson(int x) { return ch[x][1]; } int lson(int x) { return ch[x][0]; } int get(int x) { return ch[f[x]][1]==x; } void update(int x,int c) { val[x]+=(x%2==0?1:-1)*c; lazy[x]+=c; sumv[x]+=siz[x]*c; } void pushup(int x) { sumv[x]=sumv[lson(x)]+sumv[rson(x)]+val[x]; siz[x]=siz[lson(x)]+siz[rson(x)]; siz[x]+=(x>=200099 ? 0: (x%2==0?1:-1)); } void pushdown(int x) { if(lazy[x]) update(lson(x),lazy[x]), update(rson(x),lazy[x]),lazy[x]=0; } int pre(int x) { splay(x,root); x=lson(root); while(rson(x)) x=rson(x); return x; } int las(int x) { splay(x,root); x=rson(root); while(lson(x)) x=lson(x); return x; } void build(int l,int r,int &o,int fa) { if(l>r)return; int mid=(l+r)>>1; o=euler[mid], f[o]=fa; build(l,mid-1,ch[o][0],o); build(mid+1,r,ch[o][1],o); pushup(o); } void rotate(int x) { int old=f[x], oldf=f[old], which=get(x); ch[old][which]=ch[x][which^1], f[ch[old][which]]=old; ch[x][which^1]=old,f[old]=x,f[x]=oldf; if(oldf) ch[oldf][ch[oldf][1]==old]=x; pushup(old),pushup(x); } void splay(int x,int &tar) { int v=0,u=x,a=f[tar]; while(u!=a) sta[++v]=u,u=f[u]; while(v) pushdown(sta[v--]); for(int fa;(fa=f[x])!=a;rotate(x)) if(f[fa]!=a) rotate(get(fa)==get(x)?fa:x); tar=x; } void opt1(int a) { int x=las(a*2); splay(200100,root),splay(x,ch[root][1]); printf("%lld\n",sumv[ch[ch[root][1]][0]]); } void opt2(int a,int b) { int x=pre(a*2),y=las(a*2+1),tmp; splay(x,root),splay(y,ch[root][1]); tmp=ch[ch[root][1]][0], ch[ch[root][1]][0]=f[tmp]=0; pushup(ch[root][1]), pushup(root); x=las(b*2); splay(b*2,root), splay(x,ch[root][1]); ch[ch[root][1]][0]=tmp,f[tmp]=ch[root][1]; pushup(ch[root][1]),pushup(root); } void opt3(int a,int b) { int x=pre(a*2),y=las(a*2+1); splay(x,root),splay(y,ch[root][1]); update(ch[ch[root][1]][0],b); pushup(ch[root][1]),pushup(root); } }tree; int main(){ //setIO("input"); int n,x,m; scanf("%d",&n); for(int i=2;i<=n;++i) scanf("%d",&x),addedge(x,i); for(int i=1;i<=n;++i) scanf("%d",&w[i]); dfs(1); euler[1]=200100, euler[++cnt]=200101; tree.build(1,cnt,root,0); scanf("%d",&m); while(m--) { char opt[10]; int a,b; scanf("%s",opt); if(opt[0]=='Q') { scanf("%d",&a); tree.opt1(a); } if(opt[0]=='C') { scanf("%d%d",&a,&b); tree.opt2(a,b); } if(opt[0]=='F') { scanf("%d%d",&a,&b); tree.opt3(a,b); } } return 0; }
BZOJ 3786: 星系探索 欧拉游览树的更多相关文章
- 数据结构&图论:欧拉游览树
ETT可以称为欧拉游览树,它是一种和欧拉序有关的动态树(LCT是解决动态树问题的一种方案,这是另一种) dfs序和欧拉序是把树问题转化到区间问题上然后再用数据结构去维护的利器 通过借助这两种形式能够完 ...
- BZOJ 3786: 星系探索 解题报告
3786: 星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅 ...
- BZOJ 3786 星系探索
Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球.主星球 ...
- BZOJ 3786: 星系探索 [伪ETT]
传送门 数据,标程 题意: 一颗有根树,支持询问点到根路径权值和,子树加,换父亲 欧拉序列怎么求路径权值和? 一个点的权值只会给自己的子树中的点贡献,入栈权值正出栈权值负,求前缀和就行了! 和上题一样 ...
- bzoj 3786 星系探索 dfs+splay
[BZOJ3786]星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球 ...
- BZOJ 3786 星系探索 (splay+dfs序)
题目大意:给你一棵树,支持一下三种操作 1.获取某节点到根节点的路径上所有节点的权值和 2.更换某棵子树的父亲 3.某子树内所有节点的权值都增加一个值w 当时想到了splay维护dfs序,查完题解发现 ...
- BZOJ 3786: 星系探索 ETT
Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球.主星球 ...
- BZOJ 3786 星系探索 ——Splay
子树可以移动,唔. 还是用Splay维护DFS序即可. 子树的话直接截取出来就好了. 然后求前驱后继可能麻烦一些. 添加两个虚拟节点会比较好写. #include <map> #inclu ...
- BZOJ 2818 Gcd 线性欧拉
题意:链接 方法:线性欧拉 解析: 首先列一下表达式 gcd(x,y)=z(z是素数而且x,y<=n). 然后我们能够得到什么呢? gcd(x/z,y/z)=1; 最好还是令y>=x 则能 ...
随机推荐
- caffe中LetNet-5卷积神经网络模型文件lenet.prototxt理解
caffe在 .\examples\mnist文件夹下有一个 lenet.prototxt文件,这个文件定义了一个广义的LetNet-5模型,对这个模型文件逐段分解一下. name: "Le ...
- 寻找两个有序数组的中位数 C++实现leetcode系列(四)
给定两个大小为 m 和 n 的有序数组 nums1和 nums2. 请你找出这两个有序数组的中位数,并且要求算法的时间复杂度为 O(log(m + n)). 你可以假设 nums1 和 nums2 不 ...
- POJ 1950暴搜
思路: 暴力枚举好了..每回判断一下-- 用long long会超时 但是10^20会爆int... 不过仔细想一想 超过10^9的数肯定拼不回0啊-- 猥琐用int AC了 (当然可以打表 ) // ...
- C#压缩和解压缩字节(GZip)
标题:C#压缩和解压缩字节(GZip的使用) 作用:此类在 .NET Framework 2.0 版中是新增的.提供用于压缩和解压缩流的方法和属性.定义:表示 GZip 数据格式,它使用无损压缩和解压 ...
- vue路由跳转传参
this.$router.push({ path: '/message/result_notice', query: { id: id } }) // let type = this.$route.n ...
- How Javascript works (Javascript工作原理) (二) 引擎,运行时,如何在 V8 引擎中书写最优代码的 5 条小技巧
个人总结: 一个Javascript引擎由一个标准解释程序,或者即时编译器来实现. 解释器(Interpreter): 解释一行,执行一行. 编译器(Compiler): 全部编译成机器码,统一执行. ...
- 系统级脚本 rpcbind
[root@web02 ~]# vim /etc/init.d/rpcbind #! /bin/sh # # rpcbind Start/Stop RPCbind # # chkconfig: 234 ...
- python学习笔记:第五天
day05: 1.字符串格式化输出: 1.占位符:%s (字符串) %d(整型) %f (浮点型) 打印格式:print("字符串为%s" %s) 2.字符串:判断是否是 ...
- 紫书 习题8-6 UVa 1611 (构造法)
这道题和例题8-1相当的像. 例题8-1https://blog.csdn.net/qq_34416123/article/details/80112017 一开始我还以为用归并的思想, 用交换把局部 ...
- jQuery实现点击图标div循环放大缩小功能
很基础的一个功能,点击左下角的图标按钮,地图的整个div会变大,变大预览之后,再次点击图标按钮,地图的整个div会变小,恢复原样,两个图标在地图界面的放大和缩小时间不断的切换图标状态(箭头向里面,或者 ...