题目大意:给你一棵树,支持一下三种操作

1.获取某节点到根节点的路径上所有节点的权值和

2.更换某棵子树的父亲

3.某子树内所有节点的权值都增加一个值w

当时想到了splay维护dfs序,查完题解发现思路是对的,然后我就写了足足6个小时才A

st[x]代表入栈时间,ed[x]代表出栈时间

对于第一个操作,每个树上节点在splay中都有两个位置,分别对应入栈出栈序,然后把入栈的点权*1,出栈点权*-1,就可以避免其它子树干扰了

接着查询到根节点的路径权值和呢,splay把1~st[x]整个序列都扔到一个子树里,求一下子树权值和即可

第二个操作,可以找一下st[x]和ed[x]对应splay中节点的前驱和后继,把它们都搞到一个splay子树里,在把整棵树接到要求的父节点即可

第三个操作,和第二个操作类似,找前驱后继然后在子树打标记即可,此外,还要维护splay子树内的入栈节点数-出栈节点数来直接得到splay子树的权值和。

注意,一定要在所有能传标记的地方都传标记,不论是从下往上splay,还是从上到下找前驱or后继,而如果不在找前驱后继的过程中下传标记,上层的标记无法被传到下层,导致结果错误!就这个错误我调了2个小时!

 #include <cstdio>
#include <cstring>
#include <algorithm>
#define N 200100
#define il inline
#define ll long long
#define root d[0].ch[1]
using namespace std; char q[];
int n,m,cte,tot;
int st[N],ed[N],head[N],sq[N];
ll w[N];
struct EDGE{int to,nxt;}edge[N];
struct SPLAY{int fa,ch[],ss,fl;ll val,sum,tag;}d[N];
il int idf(int x){return d[d[x].fa].ch[]==x?:;}
il void con(int x,int ff,int p){d[ff].ch[p]=x;d[x].fa=ff;}
il void pushup(int x){
d[x].sum=d[d[x].ch[]].sum+d[d[x].ch[]].sum+d[x].val;
d[x].ss=d[d[x].ch[]].ss+d[d[x].ch[]].ss+d[x].fl;
}
il void pushdown(int x){
if(!d[x].tag) return;
int ls=d[x].ch[],rs=d[x].ch[];
d[ls].tag+=d[x].tag,d[rs].tag+=d[x].tag;
d[ls].val+=d[x].tag*(ll)d[ls].fl,d[rs].val+=d[x].tag*(ll)d[rs].fl;
d[ls].sum+=d[x].tag*(ll)d[ls].ss,d[rs].sum+=d[x].tag*(ll)d[rs].ss;
d[x].tag=;
}
il void rot(int x){
int y=d[x].fa;int ff=d[y].fa;int px=idf(x);int py=idf(y);
pushdown(y),pushdown(x);
con(d[x].ch[px^],y,px),con(y,x,px^),con(x,ff,py);
pushup(y),pushup(x);
}
il void splay(int x,int to){
to=d[to].fa;
while(d[x].fa!=to){
int y=d[x].fa;
if(d[y].fa==to) rot(x);
else if(idf(x)==idf(y)) rot(y),rot(x);
else rot(x),rot(x);}
}
int gc(){
int rett=,fh=;char p=getchar();
while(p<''||p>''){if(p=='-')fh=-;p=getchar();}
while(p>=''&&p<=''){rett=(rett<<)+(rett<<)+p-'';p=getchar();}
return rett*fh;
}
void ae(int u,int v){
cte++;edge[cte].to=v;
edge[cte].nxt=head[u];head[u]=cte;
}
void dfs1(int u,int ff)
{
st[u]=++tot,sq[tot]=u,d[tot].val=w[u],d[tot].fl=;
for(int j=head[u];j!=-;j=edge[j].nxt){
int v=edge[j].to;
if(v==ff) continue;
dfs1(v,u);}
ed[u]=++tot,sq[tot]=u,d[tot].val=-w[u],d[tot].fl=-;
}
int build(int l,int r)
{
if(l>r) return ;
int mid=(l+r)>>;
con(build(l,mid-),mid,);
con(build(mid+,r),mid,);
pushup(mid);
return mid;
}
int lower(int x)
{
pushdown(x),x=d[x].ch[];
while(x){
pushdown(x);
if(d[x].ch[]) x=d[x].ch[];
else break;}
return x;
}
int upper(int x)
{
pushdown(x),x=d[x].ch[];
while(x){
pushdown(x);
if(d[x].ch[]) x=d[x].ch[];
else break;}
return x;
}
ll query(int x)
{
splay(st[],root);
splay(st[x],d[root].ch[]);
int lw=upper(st[x]);
splay(lw,d[root].ch[]);
return d[d[lw].ch[]].sum;
}
void add(int x,ll ww)
{
int trl,trr,rt;
splay(st[x],root);
splay(ed[x],d[root].ch[]);
trl=lower(st[x]);
splay(trl,root);
trr=upper(ed[x]);
splay(trr,d[root].ch[]);
rt=d[trr].ch[];
d[rt].val+=ww*(ll)d[rt].fl;
d[rt].sum+=ww*(ll)d[rt].ss;
d[rt].tag+=ww;
}
void mov(int x,int y)
{
int trl,trr,frt,rtt;
splay(st[x],root);
splay(ed[x],d[root].ch[]);
trl=lower(st[x]);
splay(trl,root);
trr=upper(ed[x]);
splay(trr,d[root].ch[]);
pushdown(trl),pushdown(trr);
rtt=d[trr].ch[];
d[trr].ch[]=,d[rtt].fa=;
pushup(trr),pushup(trl); splay(st[y],root);
splay(ed[y],d[root].ch[]);
trl=lower(ed[y]);
if(trl) con(rtt,trl,);
else con(rtt,ed[y],);
splay(rtt,root);
}
void print(int x)
{
if(d[x].ch[]) print(d[x].ch[]);
printf("%d ",sq[x]);
if(d[x].ch[]) print(d[x].ch[]);
} int main()
{
scanf("%d",&n);
int x,y;
memset(head,-,sizeof(head));
for(int i=;i<=n;i++)
x=gc(),ae(x,i),ae(i,x);
for(int i=;i<=n;i++)
w[i]=gc();
st[]=++tot,dfs1(,-),ed[]=++tot;
root=build(,tot);
scanf("%d",&m);
for(int i=;i<=m;i++){
scanf("%s",q);
if(q[]=='Q'){
x=gc();
printf("%lld\n",query(x));
}else if(q[]=='F'){
x=gc(),y=gc();
add(x,(ll)y);
}else{
x=gc(),y=gc();
mov(x,y);
}
}
/*for(int i=1;i<=n;i++)
printf("%lld %lld %lld\n",query(i),d[st[i]].val,d[ed[i]].val);*/
return ;
}

BZOJ 3786 星系探索 (splay+dfs序)的更多相关文章

  1. 【BZOJ-3786】星系探索 Splay + DFS序

    3786: 星系探索 Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 647  Solved: 212[Submit][Status][Discuss] ...

  2. BZOJ3786: 星系探索 Splay+DFS序

    题目大意:给你一个树,支持三种操作,子树加,点到根的路径和,改变某一个点的父亲. 分析: 看起来像一个大LCT,但是很显然,LCT做子树加我不太会啊... 那么,考虑更换一个点的父亲这个操作很有意思, ...

  3. BZOJ 3786 星系探索 ——Splay

    子树可以移动,唔. 还是用Splay维护DFS序即可. 子树的话直接截取出来就好了. 然后求前驱后继可能麻烦一些. 添加两个虚拟节点会比较好写. #include <map> #inclu ...

  4. bzoj 3786 星系探索 dfs+splay

    [BZOJ3786]星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球 ...

  5. BZOJ 3786: 星系探索 解题报告

    3786: 星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅 ...

  6. BZOJ 3786 星系探索

    Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球.主星球 ...

  7. bzoj3786 星际探索 splay dfs序

    这道题 首先 因为他求的是当前点到根节点的路径和 我们可以将题目转换为括号序列的写法 将点拆为左括号以及右括号 左括号为正 右括号为负 这样题目就变为了求前缀和了 如果一个点是这个点的子树 那么他的左 ...

  8. BZOJ3786:星系探索(Splay,括号序)

    Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球.主星球 ...

  9. BZOJ 3786: 星系探索 欧拉游览树

    一个叫 Euler-Tour-Tree 的数据结构,说白了就是用 Splay_Tree 维护欧拉序 #include <cstring> #include <algorithm> ...

随机推荐

  1. 洛谷P1428 小鱼比可爱

    题目描述 人比人,气死人:鱼比鱼,难死鱼.小鱼最近参加了一个"比可爱"比赛,比的是每只鱼的可爱程度.参赛的鱼被从左到右排成一排,头都朝向左边,然后每只鱼会得到一个整数数值,表示这只 ...

  2. omap 移植qt4.7.0

    准备: 1.Qt源码包 qt-everywhere-opensource-src-4.7.0.tar.gz 2.交叉编译器 arm-eabi-4.4.0.tar.bz2 3.触摸屏校验工具:tslib ...

  3. 00068_Date类

    1.类Date表示特定的瞬间,精确到毫秒: 2.查阅Date类的描述,发现Date拥有多个构造函数,只是部分已经过时,但是其中有未过时的构造函数可以把毫秒值转成日期对象: //创建日期对象,把当前的毫 ...

  4. Shiro身份认证授权原理

    shiro在应用程序中的使用是用Subject为入口的, 最终subject委托给真正的管理者ShiroSecurityMannager Realm是Shiro获得身份认证信息和来源信息的地方(所以这 ...

  5. mybatis逆向工程不生成Example

    mybatis逆向生成映射文件时会生成一大堆example文件,没感觉有啥用,可以手动删除这些多余的东西,使项目变得好看许多 也可以通过配置达到目的: 原配置: <table tableName ...

  6. BA-siemens-insight-event builder使用

    event builder功能主要是用来给report使用的,作为一个独立的对象,这个对象的功能就是收集点位的信息,如果再使用report功能就可以显示或输出点位的信息.

  7. [jQuery]ajax同步请求在方法内问题

    在函数内写$.ajax() 的成功函数中return 的值,并没有return到外层方法返回至

  8. HDU 5168

    把边按权值排序后,就相当于求一个子序列以1开始和以n结束.由于边权递增,而且相差>=k,所以,边的顺序也必定是递增的.知道,当处理一条出边时,必定是从入边选择一条最优的边,考虑两个因素,入边的权 ...

  9. HDU 2815

    特判B不能大于等于C 高次同余方程 #include <iostream> #include <cstdio> #include <cstring> #includ ...

  10. [GraphQL] Fetch Server Data and Client-side State in One Query using React Apollo + GraphQL

    In this lesson we look at how the Apollo @client directive can be used to fetch client-side state al ...