【BZOJ3786】星系探索 DFS序+Splay
【BZOJ3786】星系探索
Description
物理学家小C的研究正遇到某个瓶颈。
他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球。主星球没有依赖星球。
我们定义依赖关系如下:若星球a的依赖星球是b,则有星球a依赖星球b.此外,依赖关系具有传递性,即若星球a依赖星球b,星球b依赖星球c,则有星球a依赖星球c.
对于这个神秘的星系中,小C初步探究了它的性质,发现星球之间的依赖关系是无环的。并且从星球a出发只能直接到达它的依赖星球b.
每个星球i都有一个能量系数wi.小C想进行若干次实验,第i次实验,他将从飞船上向星球di发射一个初始能量为0的能量收集器,能量收集器会从星球di开始前往主星球,并收集沿途每个星球的部分能量,收集能量的多少等于这个星球的能量系数。
但是星系的构成并不是一成不变的,某些时刻,星系可能由于某些复杂的原因发生变化。
有些时刻,某个星球能量激发,将使得所有依赖于它的星球以及他自己的能量系数均增加一个定值。还有可能在某些时刻,某个星球的依赖星球会发生变化,但变化后依然满足依赖关系是无环的。
现在小C已经测定了时刻0时每个星球的能量系数,以及每个星球(除了主星球之外)的依赖星球。接下来的m个时刻,每个时刻都会发生一些事件。其中小C可能会进行若干次实验,对于他的每一次实验,请你告诉他这一次实验能量收集器的最终能量是多少。
Input
第一行一个整数n,表示星系的星球数。
接下来n-1行每行一个整数,分别表示星球2-n的依赖星球编号。
接下来一行n个整数,表示每个星球在时刻0时的初始能量系数wi.
接下来一行一个整数m,表示事件的总数。
事件分为以下三种类型。
(1)"Q di"表示小C要开始一次实验,收集器的初始位置在星球di.
(2)"C xi yi"表示星球xi的依赖星球变为了星球yi.
(3)"F pi qi"表示星球pi能量激发,常数为qi.
Output
对于每一个事件类型为Q的事件,输出一行一个整数,表示此次实验的收集器最终能量。
Sample Input
1
1
4 5 7
5
Q 2
F 1 3
Q 2
C 2 3
Q 2
Sample Output
15
25
HINT
n<=100000,m<=300000,1<di,xi<=n,wi,qi<=100000.保证操作合法。
题解:参考http://blog.csdn.net/PoPoQQQ/article/details/41649197
如果没有C操作,我们完全可以树剖搞定,但树剖不能实现修改子树
于是我们想到在DFS序中,先求出入栈出栈序,那么子树就变成了一段连续的区间,这样我们就可以用Splay来进行区间移动,并维护区间和
那么在入栈出栈序中怎么求出到根节点的路径上的点权和呢?
我们可以给入栈出栈序上的每个点乘以一个系数,入栈为1,出栈为-1,这样我们只要求一个前缀和就可以啦。但我们还要维护Splay中区间的系数和,因为这在区间修改时会用到
好久没写Splay了,以前写的的Splay都是直来直去的区间,今天调试时各种WA,各种TLE,各种RE,渣电脑快被我搞炸了~
还有注意这题要开long long
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn=100010;
typedef long long ll;
struct node
{
ll fa,ch[2],size;
ll tag,sw,sv,w,v;
}p[maxn<<1];
ll n,m,tot,cnt,root;
ll f[maxn],to[maxn],next[maxn],head[maxn],q[maxn<<1];
char str[5];
ll readin()
{
ll ret=0; char gc;
while(gc<'0'||gc>'9') gc=getchar();
while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
return ret;
}
void dfs(ll x)
{
q[++tot]=x;
for(ll i=head[x];i!=-1;i=next[i]) dfs(to[i]);
q[++tot]=x+n;
}
void add(ll a,ll b)
{
to[cnt]=b;
next[cnt]=head[a];
head[a]=cnt++;
}
void pushup(ll x)
{
p[x].sw=p[p[x].ch[0]].sw+p[p[x].ch[1]].sw+p[x].w;
p[x].sv=p[p[x].ch[0]].sv+p[p[x].ch[1]].sv+p[x].v;
p[x].size=p[p[x].ch[0]].size+p[p[x].ch[1]].size+1;
}
void pushdown(ll x)
{
if(p[x].ch[0]) p[p[x].ch[0]].sw+=p[p[x].ch[0]].sv*p[x].tag,p[p[x].ch[0]].w+=p[p[x].ch[0]].v*p[x].tag,p[p[x].ch[0]].tag+=p[x].tag;
if(p[x].ch[1]) p[p[x].ch[1]].sw+=p[p[x].ch[1]].sv*p[x].tag,p[p[x].ch[1]].w+=p[p[x].ch[1]].v*p[x].tag,p[p[x].ch[1]].tag+=p[x].tag;
p[x].tag=0;
}
void build(ll l,ll r,ll last)
{
if(l>r) return ;
ll mid=l+r>>1;
p[q[mid]].fa=q[last];
if(last) p[q[last]].ch[mid>last]=q[mid];
build(l,mid-1,mid),build(mid+1,r,mid);
pushup(q[mid]);
}
void rotate(ll x,ll &k)
{
ll y=p[x].fa,z=p[y].fa,d=(x==p[y].ch[1]);
if(y==k) k=x;
else p[z].ch[y==p[z].ch[1]]=x;
p[x].fa=z,p[y].fa=x,p[y].ch[d]=p[x].ch[d^1];
if(p[x].ch[d^1]) p[p[x].ch[d^1]].fa=y;
p[x].ch[d^1]=y;
pushup(y),pushup(x);
}
ll qrank(ll x)
{
ll ret;
if(x==root) ret=p[p[x].ch[0]].size+1;
else if(x==p[p[x].fa].ch[0]) ret=qrank(p[x].fa)-p[p[x].ch[1]].size-1;
else ret=qrank(p[x].fa)+p[p[x].ch[0]].size+1;
pushdown(x);
return ret;
}
void splay(ll x,ll &k)
{
while(x!=k)
{
ll y=p[x].fa,z=p[y].fa;
if(y!=k)
{
if((p[y].ch[0]==x)^(p[z].ch[0]==y)) rotate(x,k);
else rotate(y,k);
}
rotate(x,k);
}
}
ll find(ll x,ll y)
{
pushdown(x);
if(p[p[x].ch[0]].size+1==y) return x;
if(y<=p[p[x].ch[0]].size) return find(p[x].ch[0],y);
else return find(p[x].ch[1],y-p[p[x].ch[0]].size-1);
}
void _Q()
{
ll x;
x=readin();
qrank(x);
splay(x,root);
printf("%lld\n",p[p[x].ch[0]].sw+p[x].w);
}
void _C()
{
ll x,y,u,t;
x=readin(),y=readin();
splay(find(root,qrank(x)-1),root);
splay(find(root,qrank(x+n)+1),p[root].ch[1]);
u=p[p[root].ch[1]].ch[0];
p[p[root].ch[1]].ch[0]=0;
pushup(p[root].ch[1]),pushup(root);
qrank(y),qrank(y+n);
splay(y,root);
splay(y+n,p[root].ch[1]);
t=p[root].ch[1];
while(p[t].ch[0]) t=p[t].ch[0];
qrank(t);
p[t].ch[0]=u,p[u].fa=t;
while(t) pushup(t),t=p[t].fa;
}
void _F()
{
ll x,y;
x=readin(),y=readin();
splay(find(root,qrank(x)-1),root);
splay(find(root,qrank(x+n)+1),p[root].ch[1]);
p[p[p[root].ch[1]].ch[0]].sw+=p[p[p[root].ch[1]].ch[0]].sv*y;
p[p[p[root].ch[1]].ch[0]].w+=p[p[p[root].ch[1]].ch[0]].v*y;
p[p[p[root].ch[1]].ch[0]].tag+=y;
}
int main()
{
n=readin();
ll i;
memset(head,-1,sizeof(head));
for(i=2;i<=n;i++) f[i]=readin(),add(f[i],i);
for(i=1;i<=n;i++) p[i].w=readin(),p[i+n].w=-p[i].w,p[i].v=1,p[i+n].v=-1;
tot=1;
dfs(1);
root=q[n+1],q[1]=2*n+1,q[2*n+2]=2*n+2;
build(1,2*n+2,0);
m=readin();
for(i=1;i<=m;i++)
{
scanf("%s",str);
switch(str[0])
{
case 'Q':_Q(); break;
case 'C':_C(); break;
case 'F':_F(); break;
}
}
return 0;
}
【BZOJ3786】星系探索 DFS序+Splay的更多相关文章
- BZOJ3786 星系探索 【Splay维护dfs序】*
BZOJ3786 星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均 ...
- [BZOJ3786]星系探索
[BZOJ3786]星系探索 试题描述 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个 ...
- BZOJ3786: 星系探索 Splay+DFS序
题目大意:给你一个树,支持三种操作,子树加,点到根的路径和,改变某一个点的父亲. 分析: 看起来像一个大LCT,但是很显然,LCT做子树加我不太会啊... 那么,考虑更换一个点的父亲这个操作很有意思, ...
- BZOJ3786星系探索——非旋转treap(平衡树动态维护dfs序)
题目描述 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球均有且仅有一个依赖星球.主星球没有依赖星球. ...
- bzoj 3786 星系探索 dfs+splay
[BZOJ3786]星系探索 Description 物理学家小C的研究正遇到某个瓶颈. 他正在研究的是一个星系,这个星系中有n个星球,其中有一个主星球(方便起见我们默认其为1号星球),其余的所有星球 ...
- bzoj3786星系探索 splay
3786: 星系探索 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1314 Solved: 425[Submit][Status][Discuss ...
- [BZOJ3786] 星系探索(括号序列+Splay)
3786: 星系探索 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 2191 Solved: 644[Submit][Status][Discuss ...
- [BZOJ3786]星系探索(伪ETT)
3786: 星系探索 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1638 Solved: 506[Submit][Status][Discuss ...
- BZOJ3786: 星系探索(伪ETT)
题面 传送门 题解 坑啊--我好像把\(Splay\)的东西全忘光了-- \(ETT\)(\(Euler\ Tour\ Tree\))是一种可以资瓷比\(LCT\)更多功能的数据结构,然而不管是功能还 ...
随机推荐
- Log4j,Log4j2,logback,slf4j日志学习
日志学习笔记 Log4j Log4j是Apache的一个开放源代码项目,通过使用Log4j,我们可以控制日志信息输送的目的地是控制台.文件.数据库等:我们也可以控制每一条日志的输出格式:通过定义每一条 ...
- java英文缩写
JSR Java Specification Request Java规范请求 POJO Plain Old Java Object 简单的Java对象
- 用 CSS 实现打印显示底色
上一篇有讲到如何在浏览器端实现打印功能.后面发现有个问题,就是表格表头有背景颜色,但是实际打印出来无背景颜色.网上的方法主要有以下几种实现方式: 1.把背景颜色写成行内样式,如下图所示: 但是发现这样 ...
- js彈出層或者js彈出引用url Frame 層
function Popup() { var _this = this; this.CssName = "layermask";//樣式 //遮蓋層 this.hiddLayer ...
- Translating between qplot and base graphics
Translating between qplot and base graphics Description There are two types of graphics functions in ...
- main函数位置
c语言中main函数的位置可以任意位置.在执行一个c语言编写的程序时,main函数就相当于是执行程序的入口.只要是没有语法和逻辑上的错误,main函数可以放在任意位置.
- C# 将RichTextBox中内容的文档以二进制形式存
private void button1_Click(object sender, EventArgs e) { System.IO.MemoryStream mstream = new ...
- zabbix-agent 自动注册
1. 概述 上一篇内容<zabbix自动发现配置>,大概内容是zabbix server去扫描一个网段,把在线的主机添加到Host列表中.我们本篇内容与上篇相反,这次是Active age ...
- 教你如何架设linux邮件服务器postfix
检查linux是否有安装postfix和dovecot 检查命令如下: Rpm -qa |grep postfix; Rpm –qa |grep dovecot; 如果没有显示任何数据,表明没有安 ...
- UVA - 11920 0 s, 1 s and ? Marks
Description 0 s, 1 s and ? Marks Given a string consisting of 0, 1 and ? only, change all the ? to ...