洛谷P10238 [yLCPC2024] F. PANDORA PARADOXXX

最新推荐文章于 2024-10-01 23:05:47 发布
最新推荐文章于 2024-10-01 23:05:47 发布
阅读量876 收藏 5
点赞数 25
分类专栏： 树的直径 LCA 文章标签： 算法 数据结构 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ZHUYINGYE_123456/article/details/141396048
版权

[Analysis] \color{blue}{\texttt{[Analysis]}} [Analysis]

删除，而且没有强制在线，一般就可以考虑时光倒流了。

离线，把删除改为加入。

在加入边前，整个图是由很多棵相互不联通的树组成的，两点间的最大长度就是所有树的直径长度的最大值。

求直径当然是 O ( n ) O(n) O(n) 的啦。

考虑加入一条边 ( u , v ) (u,v) (u,v) 对答案有什么影响。

假设加入边 ( u , v ) (u,v) (u,v) 前 u u u 所在子树的直径两端分别为 l u , r u l_{u},r_{u} lu​,ru​，同理， v v v 所在子树的直径两端分别为 l v , r v l_{v},r_{v} lv​,rv​。可以证明加入边 ( u , v ) (u,v) (u,v) 后这棵新树的直径两端一定是 l u , r u , l v , r v l_{u},r_{u},l_{v},r_{v} lu​,ru​,lv​,rv​ 四个点中的两个。

证明可以考虑分类讨论：

  1. 如果新树的直径不经过 ( u , v ) (u,v) (u,v)，那它的直径就是原来的两棵子树的直径较长的那一条；即端点要么是 l u , r u l_{u},r_{u} lu​,ru​，要么是 l v , r v l_{v},r_{v} lv​,rv​。
  2. 如果新树的直径经过 ( u , v ) (u,v) (u,v)，在原来 u u u 的子树中离 u u u 最远的点一定是 l u l_{u} lu​ 或者 r u r_{u} ru​，在原来 v v v 的子树中离 v v v 最远的点一定是 l v l_{v} lv​ 或者 r v r_{v} rv​。所以新树的直径一定是一端是 l u l_{u} lu​ 或者 r u r_{u} ru​，另一端是 l v l_{v} lv​ 或者 r v r_{v} rv​。

这样利用 LCA 可以 O ( log ⁡ n ) O(\log n) O(logn) 求出任意两点间的距离。即每次加入一条边的时间复杂度为 O ( log ⁡ n ) O(\log n) O(logn)。

总的时间复杂度为 O ( q log ⁡ n ) O(q \log n) O(qlogn)。

Code \color{blue}{\text{Code}} Code

struct UnionFind{
	int Fa[N];
	
	void clear(int n){
		for(int i=1;i<=n;i++) Fa[i]=i;
	}
	
	int getf(int x){
		if (Fa[x]==x) return x;
		return Fa[x]=getf(Fa[x]);
	}
	void merge(int x,int y){
		int u=getf(x),v=getf(y);
		Fa[v]=u;//注意次序 
	}
}Reach;//封装一个并查集

struct edge{
	int nxt,to,val;
	
	int operator = (int n){
		return nxt=to=val=n;
	}
}e[N<<1],E[N<<1];
int h[N],ecnt,H[N],Ecnt;
void add(int u,int v,int w){
	e[++ecnt]=(edge){h[u],v,w};h[u]=ecnt;
	e[++ecnt]=(edge){h[v],u,w};h[v]=ecnt;
}
void Add(int u,int v,int w){
	E[++Ecnt]=(edge){H[u],v,w};H[u]=Ecnt;
	E[++Ecnt]=(edge){H[v],u,w};H[v]=Ecnt;
}

int Fa[22][N],n,m,dep[N],Log[N];
ll len[N];

void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].to;
		if (v==fa) continue;
		
		len[v]=len[u]+e[i].val;
		Fa[0][v]=u;
		dfs(v,u);
	}
}
void dispose(){
	for(int i=1;i<=Log[n];i++)
		for(int u=1;u<=n;u++)
			Fa[i][u]=Fa[i-1][Fa[i-1][u]];
}
int LCA(int u,int v){
	if (dep[u]<dep[v]) swap(u,v);
	
	int k=dep[u]-dep[v];
	for(int i=Log[k];i>=0;i--)
		if (dep[u]-dep[v]>=(1<<i))
			u=Fa[i][u];
	
	if (u==v) return u;
	
	for(int i=Log[dep[u]];i>=0;i--)
		if (Fa[i][u]!=Fa[i][v]){
			u=Fa[i][u];v=Fa[i][v];
		}
	return Fa[0][u];
}
ll distance(int u,int v){
	return len[u]+len[v]-(len[LCA(u,v)]<<1);
}

ll dis[N];
int dp(int u,int fa){
	int ret=u;
	for(int i=H[u];i;i=E[i].nxt){
		int v=E[i].to;
		if (v==fa) continue;
		
		dis[v]=dis[u]+E[i].val;
		int tmp=dp(v,u);
		if (dis[tmp]>dis[ret]) ret=tmp;
	}
	
	return ret;
}//求树中离 u 最远的点

ll ans,out[N];bool flag[N];
int del[N],l[N],r[N];

struct linkage{
	int u,v,val;
}Link[N];

void merge(int u,int v){
	u=Reach.getf(u);
	v=Reach.getf(v);
	
	ll d1=distance(l[u],r[u]);
	ll d2=distance(l[v],r[v]);
	ll d3=distance(l[u],r[v]);
	ll d4=distance(l[v],r[u]);
	ll d5=distance(l[u],l[v]);
	ll d6=distance(r[u],r[v]);
	ll D=max(d1,max(d2,max(d3,max(d4,max(d5,d6)))));
	
	ans=max(ans,D);
	if (d2==D){
	    l[u]=l[v];r[u]=r[v];
	}//这里别漏情况
	else if (d3==D){
		l[u]=l[u];r[u]=r[v];
	}
	else if (d4==D){
		l[u]=r[u];r[u]=l[v];
	}
	else if (d5==D){
		l[u]=l[u];r[u]=l[v];
	}
	else if (d6==D){
		l[u]=r[u];r[u]=r[v];
	}
	
	Reach.merge(u,v);//注意次序 
}

void initdata(){
	ecnt=Ecnt=ans=0;
	for(int i=1;i<=(n<<1);i++)
		E[i]=e[i]=0;
	for(int i=1;i<=n;i++){
		l[i]=r[i]=h[i]=H[i]=len[i]=dep[i]=0;
		flag[i]=true;
	}
	for(int i=0;i<=20;i++)
		for(int u=1;u<=n;u++)
			Fa[i][u]=0;
}//类似于 Codeforces 里的要求，不能直接用 memset 清零

int main(){
	Log[0]=-1;
	for(int i=1;i<=1e5;i++)
		Log[i]=Log[i>>1]+1;
	
	for(int T=read();T;T--){
		n=read();m=read();
		
		initdata();
		Reach.clear(n);
		
		for(int i=1;i<n;i++){
			Link[i].u=read();
			Link[i].v=read();
			Link[i].val=read();
			
			add(Link[i].u,Link[i].v,Link[i].val);
		}
		dfs(1,0);
		dispose();
		
		for(int i=1;i<=m;i++){
			del[i]=read();
			flag[del[i]]=false;
		}
		
		for(int i=1;i<n;i++)
			if (flag[i]){
				Add(Link[i].u,Link[i].v,Link[i].val);
				Reach.merge(Link[i].u,Link[i].v);
			}
		
		for(int i=1;i<=n;i++)
			if (Reach.getf(i)==i){
				dis[i]=0;l[i]=dp(i,0);
				dis[l[i]]=0;r[i]=dp(l[i],0);
				ans=max(ans,distance(l[i],r[i]));
			}
		
		for(int i=m;i>=1;i--){
			out[i]=ans;
			merge(Link[del[i]].u,Link[del[i]].v);
		}
		
		for(int i=1;i<=m;i++)
			printf("%lld\n",out[i]);
	}
	
	return 0;
}

ll 表示 long long

这份代码能过，不过显得很冗长。因为最开始求各棵子树的直径和最后加边时求新直径用了两种不同的算法。

我们可以使用相同的算法减小代码长度：最开始求子树直径的时候也不一定要 O ( n ) O(n) O(n) 的专门算法。开始时视为整个图由 n n n 棵只有 1 1 1 个节点的树组成，在求答案之前把所有不需要删除的边逐条加入图中也可以得到最开始时各棵子树的直径的最大值。

这样做的话，这一部分代码的时间复杂度是线性对数的，不如第一份代码；但是这样做可以减少代码长度，方便调试。

struct UnionFind{
	int Fa[N];
	
	void clear(int n){
		for(int i=1;i<=n;i++) Fa[i]=i;
	}
	
	int getf(int x){
		if (Fa[x]==x) return x;
		return Fa[x]=getf(Fa[x]);
	}
	void merge(int x,int y){
		int u=getf(x),v=getf(y);
		Fa[v]=u;//注意次序 
	}
}Reach;

struct edge{
	int nxt,to,val;
	
	int operator = (int n){
		return nxt=to=val=n;
	}
}e[N<<1];int h[N],ecnt;
void add(int u,int v,int w){
	e[++ecnt]=(edge){h[u],v,w};h[u]=ecnt;
	e[++ecnt]=(edge){h[v],u,w};h[v]=ecnt;
}

int Fa[22][N],n,m,dep[N],Log[N];
ll len[N];

void dfs(int u,int fa){
	dep[u]=dep[fa]+1;
	
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].to;
		if (v==fa) continue;
		
		len[v]=len[u]+e[i].val;
		Fa[0][v]=u;
		dfs(v,u);
	}
}
void dispose(){
	for(int i=1;i<=Log[n];i++)
		for(int u=1;u<=n;u++)
			Fa[i][u]=Fa[i-1][Fa[i-1][u]];
}
int LCA(int u,int v){
	if (dep[u]<dep[v]) swap(u,v);
	
	int k=dep[u]-dep[v];
	for(int i=Log[k];i>=0;i--)
		if (dep[u]-dep[v]>=(1<<i))
			u=Fa[i][u];
	
	if (u==v) return u;
	
	for(int i=Log[dep[u]];i>=0;i--)
		if (Fa[i][u]!=Fa[i][v]){
			u=Fa[i][u];v=Fa[i][v];
		}
	return Fa[0][u];
}
ll distance(int u,int v){
	return len[u]+len[v]-(len[LCA(u,v)]<<1);
}

ll ans,out[N];bool flag[N];
int del[N],l[N],r[N];

struct linkage{
	int u,v,val;
}Link[N];

void merge(int u,int v){
	u=Reach.getf(u);
	v=Reach.getf(v);
	
	ll d1=distance(l[u],r[u]);
	ll d2=distance(l[v],r[v]);
	ll d3=distance(l[u],r[v]);
	ll d4=distance(l[v],r[u]);
	ll d5=distance(l[u],l[v]);
	ll d6=distance(r[u],r[v]);
	ll D=max(d1,max(d2,max(d3,max(d4,max(d5,d6)))));
	
	ans=max(ans,D);
	if (d2==D){
		l[u]=l[v];r[u]=r[v];
	}
	else if (d3==D){
		l[u]=l[u];r[u]=r[v];
	}
	else if (d4==D){
		l[u]=r[u];r[u]=l[v];
	}
	else if (d5==D){
		l[u]=l[u];r[u]=l[v];
	}
	else if (d6==D){
		l[u]=r[u];r[u]=r[v];
	}
	
	Reach.merge(u,v);//注意次序 
}

void initdata(){
	ecnt=ans=0;
	for(int i=1;i<=(n<<1);i++) e[i]=0;
	for(int i=1;i<=n;i++){
		l[i]=r[i]=i;
		h[i]=len[i]=dep[i]=0;
		flag[i]=true;
	}
	for(int i=0;i<=20;i++)
		for(int u=1;u<=n;u++)
			Fa[i][u]=0;
}

int main(){
	Log[0]=-1;
	for(int i=1;i<=1e5;i++)
		Log[i]=Log[i>>1]+1;
	
	for(int T=read();T;T--){
		n=read();m=read();
		
		initdata();
		Reach.clear(n);
		
		for(int i=1;i<n;i++){
			Link[i].u=read();
			Link[i].v=read();
			Link[i].val=read();
			
			add(Link[i].u,Link[i].v,Link[i].val);
		}
		dfs(1,0);
		dispose();
		
		for(int i=1;i<=m;i++){
			del[i]=read();
			flag[del[i]]=false;
		}
		
		for(int i=1;i<n;i++)//这里是改变的地方
			if (flag[i]) merge(Link[i].u,Link[i].v);
		
		for(int i=m;i>=1;i--){
			out[i]=ans;
			merge(Link[del[i]].u,Link[del[i]].v);
		}
		
		for(int i=1;i<=m;i++)
			printf("%lld\n",out[i]);
	}
	
	return 0;
}
Yingye Zhu(HPXXZYY)
关注
专栏目录
uni-app云打包 No enum constant com.pandora.pack.core.PackagePlatform.ANDROID_CLC
nba136667854的博客
11-15 2411
这个问题是原生app云打包是选择【云端证书】出现的错误 为了避开这个问题 可以选择 使用自有证书 （每个项目都要重新生成自有证书） 生成自由证书官方教程 https://ask.dcloud.net.cn/article/35777 下面是精简过程 第一步执行下面命令testalias 和test.keystore 更改为自定义别名 keytool -genkey -alias testalias -keyalg RSA -keysize 2048 -...
使用EDAS时，在启动消费者Ali-tomcat时启动失败，Pandora container stopped. HSF服务启动失败
起岸星辰
04-14 1295
我开发使用的IntelliJIdea, 找到当前问题所在是在其tomcat日志文件中发现的，因此若你还没去查看过该日志文件建议去查看。该日志文件路径在： C:\Users\tom\.IntelliJIdea2017.3\system\tomcat\项目名\logs 在使用阿里的EDAS时，开发中启动消费者的web项目无法启动，控制台打印如下信息： INFO: spas-client-initializer start JM.Log:INFO Init JM logger with Log4jLoggerF
[yLCPC2024] F. PANDORA PARADOXXX（LCA+树上差分）：经典复杂的删边问题=＞连边问题
m0_60738889的博客
06-11 649
扶苏所在的城市的机厅联合举办了 KING of PerforPandora！但是因为大雪封路，有些机厅不能到达。她想知道在能互相到达的机厅中距离最远为多少。
Maven 加载com.taobao.pandora.boot.PandoraBootstrap 失败
SJB2010007的博客
05-20 7165
<?xml version="1.0" encoding="UTF-8"?> <settings xmlns="http://maven.apache.org/SETTINGS/1.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation=...
K2P 潘多拉固件PandoraBox
热门推荐
路由器相关
09-27 1万+
k2p 潘多拉固件，来自L大，请自行研究。 解压密码 123解压密码 123 网盘下载：https://tc5.us/file/13909188-399669204
Spring——Alibaba-pandora boot实战
庄小焱
07-10 3964
阿里的Pandora Boot的核心是Pandora，因此在介绍Pandora Boot之前需要先介绍Pandora。在阿里集体内部，几乎所有的应用都用到了各式各样的中间件，比如HSF、TDDL、Diamond等等。本身中间件之间可能就有版本依赖的问题，比如你的应用HSF和Diamond分别依赖了同名jar包的不同版本，maven只会引入其中一个版本。同样的中间件和应用之间也存在同样的Jar包依赖的问题，出于要解决这些依赖冲突的问题，阿里就开发了Pandora。简单的来说就是一个类隔离容器，对外以taoba
Pandora：Android调试利器
人生只是一种体验，不必用来演绎完美。
07-03 1735
Pandora 是一款无需ROOT、可以直接在 应用内 查看和修改包括网络、数据库、UI等的Android工具箱，适合开发和测试阶段的各种问题的快速定位。 功能 查看每条网络请求的详细日志，例如headers、response等； 查看自身应用的内部存储系统； 查看所有数据库，支持直接进行增删改查操作； 查看并编辑所有Shared Preference； 预览当前页面的视图层级、查看/...
阿里云项目实践一 PandoraBoot与CommandLineRunner
taylor的博客
02-03 839
HmacSHA256和HmacSHA1签名实践 阿里云开放api import com.taobao.pandora.boot.PandoraBootstrap; import org.apache.commons.codec.binary.Base64; import org.springframework.boot.CommandLineRunner; import org.springframework.boot.SpringApplication; import org.springframew
pandora k2p
Who_really_care的博客
04-27 1458
有线连接 先不插网线，wifi无密码连接，进去路由器登录页面 账号root 密码admin 无线下 设置wifi密码2.4g 5g 修改接口 lan口静态ip改为dhcp动态 保存 插网线至任一lan口即可 无线桥接 无线下 设置查找要桥接的wifi名 高级设置填入密码，设置路由器密码 即可 ...
PandoraBox-mt7621-R8.1.11-k2p-squashfs-sysupgrade.bin
02-02
k2p rom包 潘多拉版本
jp.sblo.pandora.jota_89.apk
08-05
一款很好的文本编辑软件，支持编程语言语法高亮，包括c++、java、yaml、python等常用语言。
K2P路由器最新固件 荒野无灯、H大、openwrt、pandorabox最新固件合集及eeprom
11-18
标题中的“K2P路由器最新固件 荒野无灯、H大、openwrt、pandorabox最新固件合集及eeprom”揭示了这个压缩包内容的主体，涉及的是K2P路由器的不同定制固件版本，以及与之相关的eeprom数据。这些固件是由不同的开发者...
二叉树最小深度
2301_80161204的博客
09-27 465
25[0, 105]
使用Python实现图形学的纹理映射算法
qq_42568323的博客
09-27 1277
纹理映射的基本概念是将二维图像（纹理）应用于三维模型的表面，以便在渲染时提供更多的细节和视觉效果。纹理坐标生成：为每个顶点生成对应的纹理坐标，通常使用归一化的二维坐标系统。纹理采样：根据生成的纹理坐标从纹理图像中获取颜色值。光照计算：结合光照模型，计算最终颜色，产生更真实的效果。图像合成：将计算得到的颜色与其他渲染信息合成，生成最终图像。纹理映射不仅可以用于简单的颜色填充，还可以模拟表面细节、光照变化等，提升渲染效果。
904. 水果成篮
YQ20210216的博客
09-27 365
转化，找出一个最长的子数组长度，子数组中不超过两种类型的水果，要记录水果类型和水果数量，所以要使用unord_map，key来记录水果中类，valuse记录水果数量，当right开始增加的时候，将水果种类加入到哈希表中，当水果中类超过两种时，控制left使得水果种类回到两种。你想要尽可能多地收集水果。你正在探访一家农场，农场从左到右种植了一排果树。这些树用一个整数数组。，返回你可以收集的水果的。
排序算法之归并排序
最新发布
qq_64039411的博客
10-01 365
需要注意的是，在实际应用中，归并排序的空间复杂度较高，需要额外的空间来存储临时数组。如果对空间复杂度有要求，可以考虑使用其他排序算法。归并排序是一种分治算法，它将一个数组分成两个子数组，然后递归地对子数组进行排序，最后将两个有序的子数组合并成一个有序的数组。归并排序的时间复杂度是O(nlogn)，其中n是数组的长度。它是一种稳定的排序算法，适用于各种数据类型的排序。合并排序的关键在于合并操作。
＜使用生成式AI对四种冒泡排序实现形式分析解释的探讨整理＞
m0_57859086的博客
10-01 987
/冒泡排序实现1——不常用（这个版本两层for循环分别是i++、j++）// 声明三个整型变量i、j和temp用于循环控制和交换元素for (i = 0;i < n - 1;i++) { // 外层循环，从数组的第一个元素开始遍历到倒数第二个元素j < n;j++) { // 内层循环，从当前外层循环位置的下一个元素开始遍历到最后一个元素if (A[i] > A[j]) { // 如果当前元素大于下一个元素// 将当前元素存储到临时变量temp中// 将下一个元素赋值给当前元素。
DAY17||654.最大二叉树 |617.合并二叉树 |700.二叉搜索树中的搜索 |
2301_79865280的博客
09-29 1275
给定一个不含重复元素的整数数组。通过给定的数组构建最大二叉树，并且输出这个树的根节点。示例 ：构造二叉树一定要，因为要先找到中，才能分割左右子树。本题图解思路：本题做过上一题构建数组确实变得容易很多。代码优化方向有，不用每次找左右子树的时候构建新的数组的话效率就会高很多。（即直接用下标索引直接在原数组上操作）关于加不加if的问题，如果允许空节点进入递归，就不加。
com.taobao.pandora.boot.PandoraBootstrap
08-09
您提到的类 `com.taobao.pandora.boot.PandoraBootstrap` 是淘宝（Taobao）公司的 Pandora 框架中的一个引导类。Pandora 是一种基于 Java 的应用框架，它提供了一系列的工具和组件，用于简化应用程序的开发和管理。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值