[BZOJ1095][ZJOI2007]捉迷藏 Query on a tree IV(树链剖分)

首先,我们求出树的重链,然后对于每一条链,建一颗线段树
树大概长这样:

(其中用红边连起来的是一条条重链)
在线段树上,我们维护
Opt(u):经过 u节点代表的链的其中一段 的两个白点间的最长路径长度
MaxL(u):u节点代表的链的左端点到最远的白点的距离
MaxR(u):u节点代表的链的右端点到最远的白点的距离

怎么维护呢?
我们再定义一些辅助变量方便描述:
D(i):节点i到最远的白点的距离
D2(i):节点i到次远的白点的距离
Dist(x,y):节点x,y之间的距离
Lc:线段树上左儿子
Rc:线段树上右儿子

当 l==r 即链只有一个节点时:
若u为黑色:
MaxL(u)=MaxR(u)=D(L)
Opt(u)=D(L)+D2(L)
若u为白色:
MaxL(u)=MaxR(u)=Max{D(L),0}
Opt(u)=Max{D(L)+D2(L),D(L)}

然后考虑如何push_up:
MaxL(u)=Max{MaxL(Lc),Dist(L,mid+1)+MaxL(Rc)}
MaxR(u)=Max{MaxR(Rc),Dist(mid,R)+MaxR(Lc)}
Opt(u)=Max{Opt(Lc),Opt(Rc),MaxR(Lc)+MaxL(Rc)+Dist(mid,mid+1)}

这样我们就维护好线段树啦,接下来考虑怎么在树上修改一个点的颜色
在这里插入图片描述
假设我们要修改黄色点的颜色,那么我们肯定要修改被黄色笔框住的这条链
然后蓝色点肯定也被影响了,所以我们接下来修改被蓝色笔框住的这条链
再接下来绿色点也被影响力,所以我们修改被绿色笔框住的这条链
而下面已经没有点被影响了,修改结束
可见要修改一个点,我们只需要一层层链跳上去直到当前的链头无父节点为止
然后,我们还要考虑一下如何维护D(i)和D2(i) (节点到最远和次远白点的距离):
对此,我们对每个点维护一个大根堆,记录这个点到每个白点的距离
D(i)=s.top();
s.pop();
D2(i)=s.top();
s.push(D(i));

这样就可以求出D(i)和D2(i)啦
最后,在全局用个堆维护每条链的Opt,就可以直接查询了
这是Query on a tree IV的代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int maxn=100005;
const int inf=1e9;
struct work{priority_queue <int> f,g;inline void ins(int v) { if(v != -inf) f.push(v); } inline void era(int v) { if(v != -inf) g.push(v); }inline int top() {while(1){if(f.empty()) return -inf;if(g.empty()) return f.top();if(f.top() == g.top()) f.pop(), g.pop();else return f.top();}}
}h[maxn], ans;
struct Edge{int u,v,w,next;
}edge[maxn<<1];
struct seg{int l,r,v,ls,rs;
}st[maxn<<2];
int wn,n,m,head[maxn],cnt,col[maxn],rt[maxn],onp;
int sz[maxn],fa[maxn],dep[maxn],son[maxn];
int tid[maxn],ord[maxn],ind,top[maxn],len[maxn];
inline int read() {int p=0,w=1; char ch=getchar();while(ch>'9'||ch<'0') {if(ch=='-') w = -1; ch=getchar();}while(ch>='0'&&ch<='9') p=p*10+ch-'0',ch=getchar();return p*w;
}
void add(int u,int v,int w){edge[cnt].u=u;edge[cnt].v=v;edge[cnt].w=w;edge[cnt].next=head[u];head[u]=cnt++;
}
#define dis(x) dep[ord[x]]
void push(int u,int l,int r){int ls=st[u].ls,rs=st[u].rs,mid=(l+r)>>1;st[u].l=max(st[ls].l,st[rs].l+dis(mid+1)-dis(l));//MaxL(u)=Max{MaxL(Lc),Dist(L,mid+1)+MaxL(Rc)}st[u].r=max(st[rs].r,st[ls].r+dis(r)-dis(mid));//MaxR(u)=Max{MaxR(Rc),Dist(mid,R)+MaxR(Lc)}st[u].v=max(max(st[ls].v,st[rs].v),st[ls].r+st[rs].l+dis(mid+1)-dis(mid));//Opt(u)=Max{Opt(Lc),Opt(Rc),MaxR(Lc)+MaxL(Rc)+Dist(mid,mid+1)}
}
void build(int &u,int l,int r){if(!u) u=++onp;if(l==r){int x=ord[l];for(int i=head[x];i!=-1;i=edge[i].next){int v=edge[i].v;if(v==fa[x]||v==son[x]) continue;h[x].ins(st[rt[v]].l+dep[v]-dep[x]);}int d1=h[x].top(); h[x].era(d1); int d2=h[x].top(); h[x].ins(d1);st[u].l=st[u].r=max(d1,0); st[u].v=max(0,max(d1,d1+d2));//Opt(u)=Max{D(L)+D2(L),D(L)} 白色 return;}int mid =(l+r)>>1;build(st[u].ls,l,mid);build(st[u].rs,mid+1,r);push(u,l,r);
}
void update(int u,int l,int r,int v,int s){if(l==r){if(v==s){int d1=h[v].top();h[v].era(d1); int d2 = h[v].top(); h[v].ins(d1);if(col[v]) st[u].l=st[u].r=d1,st[u].v=d1+d2;//Opt=D(L)+D2(L) 黑色 else st[u].l=st[u].r=max(d1,0),st[u].v=max(0,max(d1,d1+d2));//Opt(u)=Max{D(L)+D2(L),D(L)} 白色 }else{h[v].ins(st[rt[s]].l+dep[s]-dep[v]);int d1 = h[v].top(); h[v].era(d1);int d2 = h[v].top(); h[v].ins(d1);if(col[v]) st[u].l=st[u].r=d1,st[u].v=d1+d2;//Opt=D(L)+D2(L) 黑色 else st[u].l=st[u].r=max(d1,0),st[u].v=max(0,max(d1,d1+d2));//Opt(u)=Max{D(L)+D2(L),D(L)} 白色 }return;}int mid=(l+r)>>1;if(tid[v]<=mid) update(st[u].ls,l,mid,v,s);else update(st[u].rs,mid+1,r,v,s);push(u,l,r);
}
void dfs1(int u){sz[u]=1;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;int w=edge[i].w;if(v==fa[u]) continue;fa[v]=u; dep[v]=dep[u]+w;dfs1(v);sz[u]+=sz[v];if(sz[v]>sz[son[u]]) son[u]=v;}
}
void dfs2(int u,int tp){tid[u]=++ind;top[u]=tp;ord[ind]=u;len[tp]++;if(!son[u]) return;dfs2(son[u],tp);for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(v==fa[u]||v==son[u]) continue;dfs2(v,v);}
} 
int main(){ios::sync_with_stdio(false);memset(head,-1,sizeof(head));wn=n=read();int u,v,w;for(int i=1;i<n;i++){u=read();v=read();w=read();add(u,v,w);add(v,u,w);}dfs1(1);dfs2(1,1);ans.ins(0);for(int i=n;i;i--){int u=ord[i]; if(u!=top[u]) continue; build(rt[u],tid[u],tid[u]+len[u]-1);  //每一条链建一线段树 ans.ins(st[rt[u]].v);}m=read();char ch;for(int i=1;i<=m;i++){ch=getchar();while(ch!='C'&&ch!='A') ch=getchar();if(ch=='C'){int x=read(); col[x]^= 1;if(col[x]==0) wn++; else wn --;for(int u=x,p=u;u;u=fa[u]){int tp=top[u];int p1=st[rt[tp]].v,d1=st[rt[tp]].l;if(fa[tp]) h[fa[tp]].era(st[rt[tp]].l+dep[tp]-dep[fa[tp]]);update(rt[tp],tid[tp],tid[tp]+len[tp]-1,u,p);int p2=st[rt[tp]].v,d2=st[rt[tp]].l;if(p1!=p2) ans.era(p1),ans.ins(p2);p=u=tp;}}else {if(wn==0) printf("They have disappeared.\n");else printf("%d\n", ans.top());}}return 0;
}

捉迷藏的代码只需要改一下输入输出,就不贴了

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320333.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ASP.NET Core Web API + Identity Server 4 + Angular 6 实战小项目视频

今天开始尝试录制ASP.NET Core Web API的教学视频. 这是一个小项目的实战视频, 该项目采用了:ASP.NET Core 2.1 做APIIdentity Server 4Angular 6Angular Material这个项目比较简单, 适合ASP.NET Core Web API 和 Angular 初学者. 项目最终完成的大致效果如图:声明: 以前虽然经…

[USACO]Sprinklers 2: Return of the Alfalfa P(网格DP)

思路&#xff1a; 由题目易得&#xff1a;网格内种植的两种植物形成了一条 左上—右下 的分界线&#xff0c;考虑将问题转化成 DP求出有多少条合法的分界线 我们注意到&#xff1a; 分界线上的点都必须放洒水器&#xff0c;且所放洒水器有唯一选择&#xff1b; 其他的可以放洒…

.net core grpc consul 实现服务注册 服务发现 负载均衡(二)

在上一篇 .net core grpc 实现通信(一) 中&#xff0c;我们实现的grpc通信在.net core中的可行性&#xff0c;但要在微服务中真正使用&#xff0c;还缺少 服务注册&#xff0c;服务发现及负载均衡等&#xff0c;本篇我们将在 .net core grpc 通信 的基础上加上 服务注册&#x…

【BZOJ3218】a+b problem (最小割 + 主席树)

传送门 继续优化&#xff1a;把a[ ]离散化 #include<bits/stdc.h> using namespace std; const int inf1000000007; const int N200010; const int M1000010;struct Edge{int u,v,f,next; }edge[M]; int head[N],cnt; int s,t,flow,level[N];struct Node{int x,id; }e[…

微软技术直通车(第三期)现场实录

微软技术直通车&#xff08;第三期&#xff09;于本月18日14时&#xff0c;在中国微软总部圆满举行。本次活动得到了微软最佳有价值专家&#xff08;Microsoft MVP&#xff09;中国区项目组的鼎力支持。同时&#xff0c;Microsoft MVP中国区项目组负责人也亲临现场致辞&#xf…

[NOI2019] 序列(模拟费用流)

原先自己想的建图&#xff1a; 正确建图&#xff1a; 但是 n 太大了&#xff0c;所以考虑模拟费用流&#xff1a; 注意&#xff1a; 在1中&#xff0c; 若选的两个位置相同&#xff0c;则为情况2&#xff0c;不用减 f&#xff1b; 若选的位置在另一序列中已被选&#xff0c;…

【结论】游戏(jzoj 5536)

游戏 jzoj 5536 题目大意&#xff1a; 给一个范围和一个矩形&#xff0c;让你在这个范围内放若干个这样的矩形&#xff08;不能改变方向&#xff0c;不能重叠&#xff09;&#xff0c;让你求出最少放多个矩形可以使范围内无法再放矩形 输入样例#1 11 4 3 2输入样例#2 10 …

LCA总结

文章目录LCA介绍解决方法概括&#xff1a;倍增法&#xff1a;Tarjan欧拉序树剖解法&#xff1a;看了很多关于LCA的文章&#xff0c;这里是一个总结我们通过这个题洛谷P3379 【模板】最近公共祖先来讲LCA LCA介绍 lca是啥&#xff1f;最近公共祖先 就是&#xff1a;两个点在这…

ASP.NET Core Web API + Ng6 实战视频 Day 2

第一天课程&#xff1a;ASP.NET Core Web API Identity Server 4 Angular 6 实战小项目视频Day 2 第一部分:Day 2 第二部分:视频专辑持续更新中....地址请点击原文链接.原文地址: http://v.qq.com/vplus/4cfb00af75c16eb8d198c58fb86eb4dc/foldervideos/8hk0029019k2fft.NET…

浮沉子

浮沉子制作及其原理 浮沉子的下沉过程&#xff1a;

.NET Core 项目指定SDK版本

一. 版本里的坑自从 .NET Core 2.1.0版本发布以后&#xff0c;近几个月微软又进行了几次小版本的发布&#xff0c;可见 .NET Core 是一门生命力非常活跃的技术。经过一段时间的实践&#xff0c;目前做 ASP.NET Core 开发时&#xff0c;使用的 Nuget 包&#xff0c;比如 Microso…

牛客网【每日一题】4月17日题目精讲 华华给月月准备礼物

文章目录题目描述题解&#xff1a;代码&#xff1a;推荐例题&#xff1a;试题链接时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld 题目描述 二月中旬虐狗节前夕&#xff0c;华华决定给…

Bing.com在.NET Core 2.1上运行

相关知识请参考.netCore开发团队博客 https://blogs.msdn.microsoft.com/dotnet/ Bing.com是一种云服务&#xff0c;运行在遍布全球许多数据中心的数千台服务器上。Bing服务器每秒处理来自全球消费者的数千个用户查询&#xff0c;通过他们的浏览器&#xff0c;使用Microsoft认…

c#中关于协变性和逆变性(又叫抗变)帮助理解

今天回忆了之前看的《深入理解C#》这本书中的泛型章节&#xff0c;其中对泛型的可变性的理解。泛型可变性分两种&#xff1a;协变和逆变。逆变也又称为抗变。怎么理解这两个名词的意思&#xff1a;①&#xff1a;协变即为在泛型接口类型中使用out标识的类型参数。协变的字面意思…

C# 7.0中的解构功能---Deconstruct

解构元组C#7.0新增了诸多功能&#xff0c;其中有一项是新元组(ValueTuple)&#xff0c;它允许我们可以返回多个值&#xff0c;并且配合解构能更加方便的进行工作&#xff0c;如下面例子可以看到解构元组可以写出优雅的代码&#xff0c;并且可以使用类型推断&#xff0c;但在这里…

csp初赛复习(往年真题+解析)

排序算法 前缀/后缀表达式 二进制补码、反码 最短路 图片/音频/视频文件格式 前序/中序/后序遍历 以比较作为基本运算&#xff0c;在 N 个数中找最小数的最少运算次数为&#xff08; &#xff09;。 A. NNN B. N−1N-1N−1 C. N2N^2N2 D. logNlogNlogN 正确答案&#xff1a; B …

ASP.NET Core 2.1中基于角色的授权

授权是来描述用户能够做什么的过程。例如&#xff0c;只允许管理员用户可以在电脑上进行软件的安装以及卸载。而非管理员用户只能使用软件而不能进行软件的安装以及卸载。它是独立的而又与验证配合使用&#xff0c;需要身份验证机制。对于应用程序来说&#xff0c;首先需要进行…

CF449B Jzzhu and Cities(Dijkstra)

设每个点到1的距离为dis[x]dis[x]dis[x],特殊边为(1,vi,wi)(1,v_i,w_i)(1,vi​,wi​) 1、wi>dis[vi]w_i>dis[v_i]wi​>dis[vi​]的特殊边可以删除 2、widis[vi]且num[vi]>1w_idis[v_i]且num[v_i]>1wi​dis[vi​]且num[vi​]>1,特殊边可以删掉 思路上的偏差&…

分布式事务一致性解决方案

一、从数据一致性谈起↑一致性问题&#xff0c;“万恶之源”是数据冗余和分布并通过网络交互网络异常是常态。1、数据一致性的情形主库、从库和缓存数据一致性&#xff0c;相同数据冗余&#xff0c;关系数据库&#xff0c;为保证关据库的高可用和高性能&#xff0c;一般会采用主…

.net core 生成二维码

其实生成二维码的组件有很多种&#xff0c;如&#xff1a;QrcodeNet&#xff0c;ZKWeb.Fork.QRCoder&#xff0c;QRCoder等我选QRCoder&#xff0c;是因为小而易用、支持大并发生成请求、不依赖任何库和网络服务。既然是.net core 那当然要用依赖注入&#xff0c;通过构造函数注…