【Tyvj2133&BZOJ1146】网络管理Network(树套树,DFS序,树状数组,主席树,树上差分)
题意:有一棵N个点的树,每个点有一个点权a[i],要求在线实现以下操作:
1:将X号点的点权修改为Y
2:查询X到Y的路径上第K大的点权
n,q<=80000 a[i]<=10^8
思路:此题明显地体现了我对主席树理解不深
树上路径K大可以直接用树剖+二分答案+树做
但DFS序+主席树也可以
对于点U,它能影响DFS序上的区间(st[u],ed[u])
所以维护方法就是类似序列K大一样 s[st[u]]++ s[ed[u]+1]--
对于路径(x,y),信息为s[x]+s[y]-s[lca(x,y)]-s[fa[lca(x,y)]]
但因为是在树上的操作,所以每个点是从它在树上的父亲继承下来的,并不是DFS序上的前一个继承的
所以在DFS时就可以建立一棵没有修改过的原树
对于询问,再建立一棵树只做更改,询问时加上原树的信息即可
不开两棵树可以,不过需要主席树合并
var t:array[..]of record
l,r,s:longint;
end;
f:array[..,..]of longint;
head,vet,next,st,ed,dep,flag,hash,root1,root2,b:array[..]of longint;
d:array[..,..]of longint;
a:array[..]of longint;
q1,q2:array[..]of longint;
n,m,i,j,tot,time,up,la,lb,cnt,q,x,y,tmp,p,sum,z,f1,f2,f3,f4:longint; procedure add(a,b:longint);
begin
inc(tot);
next[tot]:=head[a];
vet[tot]:=b;
head[a]:=tot;
end; function max(x,y:longint):longint;
begin
if x>y then exit(x);
exit(y);
end; procedure swap(var x,y:longint);
var t:longint;
begin
t:=x; x:=y; y:=t;
end; function find(x:longint):longint;
var l,r,mid:longint;
begin
l:=; r:=up;
while l<=r do
begin
mid:=(l+r)>>;
if hash[mid]=x then exit(mid);
if hash[mid]<x then l:=mid+
else r:=mid-;
end;
end; procedure update(l,r:longint;var p:longint;v,x:longint);
var mid:longint;
begin
inc(cnt); t[cnt]:=t[p];
p:=cnt; t[p].s:=t[p].s+x;
if l=r then exit;
mid:=(l+r)>>;
if v<=mid then update(l,mid,t[p].l,v,x)
else update(mid+,r,t[p].r,v,x);
end; procedure dfs(u:longint);
var e,v,i,tmp:longint;
begin
for i:= to do
begin
if dep[u]<(<<i) then break;
f[u,i]:=f[f[u,i-],i-];
end;
tmp:=find(a[u]);
root1[u]:=root1[f[u,]];
update(,up,root1[u],tmp,);
inc(time); st[u]:=time; ed[u]:=time;
flag[u]:=;
e:=head[u];
while e<> do
begin
v:=vet[e];
if flag[v]= then
begin
f[v,]:=u;
dep[v]:=dep[u]+;
dfs(v);
ed[u]:=max(ed[u],ed[v]);
end;
e:=next[e];
end;
end; function lca(x,y:longint):longint;
var d,i:longint;
begin
if dep[x]<dep[y] then swap(x,y);
d:=dep[x]-dep[y];
for i:= to do
if d and (<<i)> then x:=f[x,i];
for i:= downto do
if f[x,i]<>f[y,i] then
begin
x:=f[x,i]; y:=f[y,i];
end;
if x=y then exit(x);
exit(f[x,]);
end; procedure qsort(l,r:longint);
var i,j,mid:longint;
begin
i:=l; j:=r; mid:=b[(l+r)>>];
repeat
while mid>b[i] do inc(i);
while mid<b[j] do dec(j);
if i<=j then
begin
swap(b[i],b[j]);
inc(i); dec(j);
end;
until i>j;
if l<j then qsort(l,j);
if i<r then qsort(i,r);
end; function lowbit(x:longint):longint;
begin
exit(x and (-x));
end; function query(l,r,k:longint):longint;
var s1,s2,i,mid:longint;
begin
if l=r then exit(l);
mid:=(l+r)>>;
s1:=; s2:=;
for i:= to la do s1:=s1+t[t[q1[i]].l].s;
for i:= to lb do s2:=s2+t[t[q2[i]].l].s;
s2:=s2+t[t[f1].l].s+t[t[f2].l].s;
s1:=s1+t[t[f3].l].s+t[t[f4].l].s;
if s2-s1>=k then
begin
f1:=t[f1].l; f2:=t[f2].l;
f3:=t[f3].l; f4:=t[f4].l;
for i:= to la do q1[i]:=t[q1[i]].l;
for i:= to lb do q2[i]:=t[q2[i]].l;
exit(query(l,mid,k));
end
else
begin
f1:=t[f1].r; f2:=t[f2].r;
f3:=t[f3].r; f4:=t[f4].r;
for i:= to la do q1[i]:=t[q1[i]].r;
for i:= to lb do q2[i]:=t[q2[i]].r;
exit(query(mid+,r,k-(s2-s1)));
end;
end; begin
assign(input,'data.in'); reset(input);
assign(output,'tyvj2133.out'); rewrite(output);
read(n,m);
for i:= to n do
begin
read(a[i]); b[i]:=a[i];
end;
q:=n;
for i:= to n- do
begin
read(x,y);
add(x,y);
add(y,x);
end; for i:= to m do
for j:= to do read(d[i,j]);
for i:= to m do
if d[i,]= then
begin
inc(q); b[q]:=d[i,];
end;
qsort(,q);
hash[]:=b[]; up:=;
for i:= to q do
if b[i]<>b[i-] then begin inc(up); hash[up]:=b[i]; end;
dfs(); for i:= to m do
if d[i,]= then
begin
tmp:=find(a[d[i,]]);
j:=st[d[i,]];
while j<=n do
begin
update(,up,root2[j],tmp,-);
j:=j+lowbit(j);
end;
j:=ed[d[i,]]+;
while j<=n do
begin
update(,up,root2[j],tmp,);
j:=j+lowbit(j);
end;
tmp:=find(d[i,]);
j:=st[d[i,]];
while j<=n do
begin
update(,up,root2[j],tmp,);
j:=j+lowbit(j);
end;
j:=ed[d[i,]]+;
while j<=n do
begin
update(,up,root2[j],tmp,-);
j:=j+lowbit(j);
end; a[d[i,]]:=d[i,];
end
else
begin
la:=; lb:=; p:=lca(d[i,],d[i,]);
j:=st[d[i,]];
while j> do
begin
inc(lb); q2[lb]:=root2[j];
j:=j-lowbit(j);
end; j:=st[d[i,]];
while j> do
begin
inc(lb); q2[lb]:=root2[j];
j:=j-lowbit(j);
end; j:=st[p];
while j> do
begin
inc(la); q1[la]:=root2[j];
j:=j-lowbit(j);
end; j:=st[f[p,]];
while j> do
begin
inc(la); q1[la]:=root2[j];
j:=j-lowbit(j);
end; f1:=root1[d[i,]]; f2:=root1[d[i,]];
f3:=root1[p]; f4:=root1[f[p,]];
sum:=dep[d[i,]]+dep[d[i,]]-dep[p]*+;
z:=sum-d[i,]+;
if sum<d[i,] then writeln('invalid request!')
else writeln(hash[query(,up,z)]);
end; close(input);
close(output);
end.
【Tyvj2133&BZOJ1146】网络管理Network(树套树,DFS序,树状数组,主席树,树上差分)的更多相关文章
- 【bzoj1146】[CTSC2008]网络管理Network 倍增LCA+dfs序+树状数组+主席树
题目描述 M公司是一个非常庞大的跨国公司,在许多国家都设有它的下属分支机构或部门.为了让分布在世界各地的N个部门之间协同工作,公司搭建了一个连接整个公司的通信网络.该网络的结构由N个路由器和N-1条高 ...
- zoj2112 树状数组+主席树 区间动第k大
Dynamic Rankings Time Limit: 10000MS Memory Limit: 32768KB 64bit IO Format: %lld & %llu Subm ...
- BZOJ_2120_数颜色_Set+树状数组+主席树
BZOJ_2120_数颜色_Set+树状数组+主席树 Description 墨墨购买了一套N支彩色画笔(其中有些颜色可能相同),摆成一排,你需要回答墨墨的提问.墨墨会像你发布如下指令: 1. Q L ...
- BZOJ_1901_Zju2112 Dynamic Rankings_树状数组+主席树
BZOJ_1901_Zju2112 Dynamic Rankings_树状数组+主席树 题意: 给定一个含有n个数的序列a[1],a[2],a[3]……a[n],程序必须回答这样的询问:对于给定的i, ...
- 【bzoj3744】Gty的妹子序列 分块+树状数组+主席树
题目描述 我早已习惯你不在身边, 人间四月天 寂寞断了弦. 回望身后蓝天, 跟再见说再见…… 某天,蒟蒻Autumn发现了从 Gty的妹子树(bzoj3720) 上掉落下来了许多妹子,他发现 她们排成 ...
- P1972 [SDOI2009]HH的项链[离线+树状数组/主席树/分块/模拟]
题目背景 无 题目描述 HH 有一串由各种漂亮的贝壳组成的项链.HH 相信不同的贝壳会带来好运,所以每次散步完后,他都会随意取出一段贝壳,思考它们所表达的含义.HH 不断地收集新的贝壳,因此,他的项链 ...
- 【BZOJ】1146: [CTSC2008]网络管理Network(树链剖分+线段树套平衡树+二分 / dfs序+树状数组+主席树)
http://www.lydsy.com/JudgeOnline/problem.php?id=1146 第一种做法(时间太感人): 第二种做法(rank5,好开心) ================ ...
- [BZOJ 1146] [CTSC2008]网络管理Network(树状数组+主席树)
题目描述 M公司是一个非常庞大的跨国公司,在许多国家都设有它的下属分支机构或部门.为了让分布在世界各地的N个部门之间协同工作,公司搭建了一个连接整个公司的通信网络.该网络的结构由N个路由器和N-1条高 ...
- [luogu2617][bzoj1901][Zju2112]Dynamic Rankings【树套树+树状数组+主席树】
题目网址 [传送门] 题目大意 请你设计一个数据结构,支持单点修改,区间查询排名k. 感想(以下省略脏话inf个字) 真的强力吹爆洛谷数据,一般的树套树还给我T了一般的点,加强的待修主席树还给我卡了几 ...
- hdu_5788_Level Up(树状数组+主席树)
题目链接:hdu_5788_Level Up 题意: 有一棵树,n个节点,每个节点有个能力值A[i],mid[i],mid的值为第i节点的子树的中位数(包括本身),现在让你将其中的一个节点的A值改为1 ...
随机推荐
- 使用lua做序列化和反序列化
-- lua对象序列化 function serialize(obj) local lua = "" local t = type(obj) if t == "numbe ...
- Virt-install用法:
#一般选项:指定虚拟机的名称.内存大小.VCPU个数及特性等 -n NAME, --name=NAME:虚拟机名称,需全局惟一: -r MEMORY, --ram=MEMORY:虚拟机内 ...
- vue中文本域限制字数的方法
用watch方法,来限制字数 <template> <div class="box"> <textarea v-model="title&q ...
- CSS基础(一)
一.CSS概述 CSS 指层叠样式表 (Cascading Style Sheets) 样式定义如何显示 HTML 元素 样式通常存储在样式表中 把样式添加到 HTML 4.0 中,是为了解决内容与表 ...
- LeetCode之Weekly Contest 102
第一题:905. 按奇偶校验排序数组 问题: 给定一个非负整数数组 A,返回一个由 A 的所有偶数元素组成的数组,后面跟 A 的所有奇数元素. 你可以返回满足此条件的任何数组作为答案. 示例: 输入: ...
- 【jquery】 form ajaxSubmit 问题
常见问题 这个插件跟哪些版本的jQuery兼容? 这个插件需要jQuery v1.0.3 或 以后的版本. 这个插件需要其它插件的支持吗? 不需要. 这个插件的运行效率高吗? 是的!请到 对比页面 查 ...
- bootmem_init_node
static unsigned long __init bootmem_init_node(int node, struct meminfo *mi) in arch/arm/mm/init.c 1. ...
- 求数组中两两相加等于20的组合(Python实现)
题目 求数组中两两相加等于20的组合. 例:给定一个数组[1, 7, 17, 2, 6, 3, 14],这个数组中满足条件的有两对:17+3=20, 6+14=20. 解析 分为两个步骤: 先采用堆排 ...
- Cinder配置多Ceph后端步骤
1. 检查cinder当前backend配置 使用cinder service-list,查看cinder-volume服务的Host字段格式. 旧版格式: 新版格式: 旧版中Host字段是cinde ...
- Apache不能启动: Unable to open logs
日志名称: Application来源: Apache Service日期: 2014/3/12 14:43:21事件 ID: ...