最近公共祖先(LCA)算法详解与应用场景
1. 什么是最近公共祖先LCA最近公共祖先Lowest Common Ancestor简称LCA是图论中树结构的一个重要概念。给定一棵有根树和树中的两个节点它们的最近公共祖先就是这两个节点的所有公共祖先中距离它们最近的节点。换句话说LCA就是这两个节点到根节点的路径上最后一个相交的节点。举个例子假设我们有一棵家族树A是B和C的父母B是D和E的父母。那么D和E的LCA就是BD和C的LCA就是A。这个概念在计算机科学中有广泛应用比如在编译器优化、网络路由算法、生物信息学等领域都会用到。注意LCA问题通常假设树结构是有根的即有一个明确的根节点。对于无根树我们需要先选择一个根节点将其转化为有根树。2. 求解LCA的常见算法2.1 朴素算法最简单的LCA求解方法是朴素算法步骤如下从第一个节点开始向上遍历到根节点记录路径上的所有节点从第二个节点开始向上遍历到根节点记录路径上的所有节点比较两条路径找到最后一个相同的节点这种方法的时间复杂度是O(n)其中n是树的深度。在最坏情况下比如树退化成链表时间复杂度会达到O(n)。def findLCA(root, p, q): path_p getPath(root, p) path_q getPath(root, q) lca None for i in range(min(len(path_p), len(path_q))): if path_p[i] path_q[i]: lca path_p[i] else: break return lca def getPath(root, node): path [] while node ! root: path.append(node) node node.parent path.append(root) return path[::-1]2.2 倍增法Binary Lifting倍增法是求解LCA的高效算法时间复杂度为O(nlogn)预处理O(logn)查询。它的核心思想是通过预处理每个节点的2^i级祖先使得我们可以快速跳跃式地查找祖先。实现步骤预处理阶段计算每个节点的深度预处理每个节点的2^i级祖先查询阶段将两个节点调整到同一深度从最大可能的i开始尝试跳跃直到找到LCAclass LCA: def __init__(self, root, n): self.up [[-1]*(n1) for _ in range(20)] self.depth [0]*(n1) self.preprocess(root) def preprocess(self, root): queue [root] self.up[0][root] -1 # 假设根节点的父节点是-1 while queue: u queue.pop(0) for v in children[u]: self.depth[v] self.depth[u] 1 self.up[0][v] u queue.append(v) for k in range(1, 20): for v in range(1, n1): if self.up[k-1][v] ! -1: self.up[k][v] self.up[k-1][self.up[k-1][v]] def query(self, u, v): if self.depth[u] self.depth[v]: u, v v, u # 将u提升到与v同一深度 for k in range(19, -1, -1): if self.depth[u] - (1 k) self.depth[v]: u self.up[k][u] if u v: return u # 现在u和v在同一深度 for k in range(19, -1, -1): if self.up[k][u] ! -1 and self.up[k][u] ! self.up[k][v]: u self.up[k][u] v self.up[k][v] return self.up[0][u]2.3 Tarjan离线算法Tarjan算法是一种离线算法可以一次性处理多个LCA查询。它基于并查集和深度优先搜索时间复杂度为O(n q)其中n是节点数q是查询数。算法步骤对树进行DFS遍历当访问一个节点时将其与父节点合并处理与该节点相关的所有查询如果一个查询的另一个节点已经被访问过那么它们的LCA就是另一个节点所在集合的代表元素def tarjanOLCA(root, queries): parent [i for i in range(n1)] visited [False]*(n1) ancestor [0]*(n1) result {} def find(u): while parent[u] ! u: parent[u] parent[parent[u]] u parent[u] return u def union(u, v): root_u find(u) root_v find(v) if root_u ! root_v: parent[root_v] root_u def dfs(u): visited[u] True ancestor[u] u for v in children[u]: if not visited[v]: dfs(v) union(u, v) ancestor[find(u)] u for v in queries[u]: if visited[v]: result[(u, v)] ancestor[find(v)] dfs(root) return result3. LCA算法的应用场景3.1 树中两点间距离计算利用LCA可以高效计算树中任意两点间的距离。计算公式为 distance(u, v) depth[u] depth[v] - 2 * depth[LCA(u, v)]def distance(u, v, lca_obj): lca lca_obj.query(u, v) return depth[u] depth[v] - 2 * depth[lca]3.2 子树统计问题在某些子树统计问题中我们需要知道两个节点是否在同一个子树中或者需要统计某个子树的信息。LCA可以帮助我们快速判断节点间的关系。3.3 网络路由优化在计算机网络中LCA算法可以用于优化路由选择找到两个节点间的最短路径或者最优转发节点。3.4 基因序列分析在生物信息学中LCA用于分析基因序列的进化关系确定不同物种在进化树上的最近共同祖先。4. 算法选择与优化建议4.1 不同场景下的算法选择单次查询朴素算法足够多次查询但树结构不变倍增法或Tarjan离线算法动态树结构节点可能增加需要使用更高级的数据结构如Link-Cut Tree4.2 倍增法的优化技巧预处理时可以按需计算2^i级祖先而不是全部预计算对于深度很大的树可以考虑使用哈希表存储部分节点的祖先信息在实际实现中可以根据树的平均深度调整预处理的最大层级4.3 常见错误与调试技巧根节点处理不当确保根节点的父节点正确处理通常设为-1或自身深度计算错误在调整节点深度时注意比较和跳跃的顺序预处理不完整确保所有节点的2^i级祖先都被正确计算边界条件处理两个节点相同的情况或者一个节点是另一个节点的祖先的情况提示在实现倍增法时建议先实现朴素算法作为验证基准确保复杂算法的正确性。5. 实际案例分析5.1 LeetCode例题236. 二叉树的最近公共祖先题目描述给定一个二叉树找到两个节点的最近公共祖先。解决方案def lowestCommonAncestor(root, p, q): if not root or root p or root q: return root left lowestCommonAncestor(root.left, p, q) right lowestCommonAncestor(root.right, p, q) if left and right: return root return left if left else right这个解法利用了递归的思想时间复杂度O(n)空间复杂度O(n)递归栈。5.2 扩展问题带权树的LCA对于带权树我们可能需要计算两点间路径上的某些统计信息如最大边权、边权和等。这时可以在预处理阶段同时存储这些信息。class WeightedLCA: def __init__(self, root, n): self.up [[-1]*(n1) for _ in range(20)] self.max_edge [[0]*(n1) for _ in range(20)] self.depth [0]*(n1) self.preprocess(root) def preprocess(self, root): # 类似前面的预处理但需要同时处理边权信息 pass def query_max_edge(self, u, v): max_e 0 # 将u和v调整到同一深度同时记录最大边权 # 类似LCA查询但在跳跃时更新max_e return max_e6. 进阶话题与扩展阅读6.1 动态树上的LCA对于动态变化的树结构节点可能增加或删除我们需要更高级的数据结构来维护LCA信息Link-Cut Trees支持动态连接和断开树的边Euler Tour Trees基于欧拉序的表示方法Heavy-Light Decomposition轻重链剖分方法6.2 区间最小值查询RMQ与LCA的等价性LCA问题可以转化为RMQ问题反之亦然。这种转化使得我们可以使用高效的RMQ算法如稀疏表来解决LCA问题。转化步骤对树进行深度优先搜索记录访问顺序和每个节点的深度LCA(u, v)对应于欧拉序列中u和v首次出现位置之间的深度最小的节点6.3 并行算法对于大规模树结构可以考虑并行化的LCA算法并行DFS预处理使用MapReduce框架处理批量查询GPU加速的倍增法实现在实际工程实现中我发现在处理超大规模树结构时如社交网络图基于分布式计算的LCA算法往往能获得更好的性能。特别是在预处理阶段可以将树分割成多个子树并行处理。

相关新闻

社区自制游戏资源逆向工程:从文件解构到环境集成的通用方法论

社区自制游戏资源逆向工程:从文件解构到环境集成的通用方法论

最近在整理一些老游戏资源时,偶然翻到一个名为“DxS《blue》Rhythm hive 极困”的压缩包。文件名看起来像某种音游的谱面或资源包,但“DxS”、“blue”、“Rhythm hive”这些词组合在一起,加上“极困”这个后缀,让它在硬盘里显得格…

2026/8/3 1:47:40 阅读更多 →
Kairos05:全局门控线性注意力(GLA)数学推导

Kairos05:全局门控线性注意力(GLA)数学推导

从注意力到动态记忆:高中生也能看懂的 Kairos 门控线性注意力(GLA)数学推导 一、先用一句话理解 GLA Kairos 使用 门控线性注意力(Gated Linear Attention, GLA) 来维护长期信息。 它可以被理解为一块固定大小的“动态记忆板”: 每当新信息到来,模型先查看记忆中原来…

2026/8/3 1:47:40 阅读更多 →
OpenRGB终极指南:一个免费软件统一控制所有RGB设备,告别厂商软件束缚

OpenRGB终极指南:一个免费软件统一控制所有RGB设备,告别厂商软件束缚

OpenRGB终极指南:一个免费软件统一控制所有RGB设备,告别厂商软件束缚 【免费下载链接】OpenRGB Open source RGB lighting control that doesnt depend on manufacturer software. Supports Windows, Linux, MacOS. Mirror of https://gitlab.com/CalcPr…

2026/8/3 1:47:40 阅读更多 →

最新新闻

RF Explorer软件:手持频谱分析仪从入门到实战应用指南

RF Explorer软件:手持频谱分析仪从入门到实战应用指南

1. 从频谱仪到口袋里的“信号侦探”:RF Explorer 是什么?如果你玩过无线电、调试过无线模块,或者只是对身边的电磁波世界感到好奇,那你大概率听说过频谱分析仪。传统上,这玩意儿是实验室里的大家伙,价格动辄…

2026/8/3 2:32:00 阅读更多 →
Cocos Creator TypeScript编码规范:提升团队协作与代码质量

Cocos Creator TypeScript编码规范:提升团队协作与代码质量

1. 项目概述:为什么我们需要一份专属的Cocos编码规范?如果你正在用Cocos Creator开发游戏,无论是TypeScript还是JavaScript,大概率都遇到过这样的场景:项目初期,代码写得飞快,一切看起来井井有条…

2026/8/3 2:30:59 阅读更多 →
LaTeX文档嵌入视频:media9宏包实战指南与避坑技巧

LaTeX文档嵌入视频:media9宏包实战指南与避坑技巧

1. 从Flash到现代视频:LaTeX文档嵌入视频的困境与破局如果你还在为如何在PDF里优雅地嵌入一个视频而烦恼,或者你曾尝试过用movie15宏包,结果发现生成的PDF依赖早已被淘汰的Adobe Flash Player,那么这篇文章就是为你准备的。在制作…

2026/8/3 2:30:59 阅读更多 →
字符串处理实战:从大小写转换看编程基本功与数据清洗

字符串处理实战:从大小写转换看编程基本功与数据清洗

1. 项目概述:从一道字符串处理题看编程基本功的锤炼最近在带学生刷《信息学奥赛一本通》和OpenJudge的题目,又碰到了这道经典的“整理药名”(题目编号1139,对应OpenJudge NOI 1.7 15)。这道题表面上看,就是…

2026/8/3 2:30:59 阅读更多 →
【无人机巡航】基于卡尔曼状态估计、自适应混合控制、多入侵机处理及全姿态跟踪(横滚、俯仰、偏航)的3D无人机避障系统在MATLAB中实现

【无人机巡航】基于卡尔曼状态估计、自适应混合控制、多入侵机处理及全姿态跟踪(横滚、俯仰、偏航)的3D无人机避障系统在MATLAB中实现

✅作者简介:热爱科研的Matlab仿真开发者,擅长毕业设计辅导、数学建模、数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。🍎 往期回顾关注个人主页:Matlab科研工作室👇 关注我领取海量matlab电子书和…

2026/8/3 2:30:59 阅读更多 →
XIAO ESP32-C5 Zigbee开发实战:从环境搭建到双核通信全解析

XIAO ESP32-C5 Zigbee开发实战:从环境搭建到双核通信全解析

1. 从ESP32-C5到Zigbee:为什么是它?如果你最近在关注物联网开发板,尤其是那些主打无线连接和低功耗的,那么Seeed Studio的XIAO ESP32-C5这个名字大概率已经出现在你的视野里了。它最吸引人的地方,就是把一颗支持Wi-Fi …

2026/8/3 2:30:59 阅读更多 →

日新闻

3个让你工作效率翻倍的Umi-OCR实战技巧:免费离线文字识别完全指南

3个让你工作效率翻倍的Umi-OCR实战技巧:免费离线文字识别完全指南

3个让你工作效率翻倍的Umi-OCR实战技巧:免费离线文字识别完全指南 【免费下载链接】Umi-OCR OCR software, free and offline. 开源、免费的离线OCR软件。支持截屏/批量导入图片,PDF文档识别,排除水印/页眉页脚,扫描/生成二维码。…

2026/8/3 0:00:47 阅读更多 →
[具身智能-181]:PC+服务器+具身机器人:构建具身智能从仿真到量产的闭环迭代混合架构

[具身智能-181]:PC+服务器+具身机器人:构建具身智能从仿真到量产的闭环迭代混合架构

PC服务器具身机器人:构建具身智能从仿真到量产的闭环迭代混合架构一、前言:具身智能需要“混合算力闭环系统”传统人工智能依赖云端静态数据集训练,不具备物理交互能力,无法适应真实世界的不确定性。具身智能(Embodied…

2026/8/3 0:00:47 阅读更多 →
[具身智能-181]:大分布式通信模型对比:看懂为什么 DDS 是 ROS2 底层通信最优解

[具身智能-181]:大分布式通信模型对比:看懂为什么 DDS 是 ROS2 底层通信最优解

前言构建机器人、具身智能这类分布式实时系统,通信底座直接决定整套系统的实时性、容错性、组网能力。分布式领域长期存在 4 类经典通信架构:点对点模式、Broker 中间代理模式、广播模式、以数据为中心(DDS)模式。很多开发者疑惑&…

2026/8/3 0:00:47 阅读更多 →

周新闻

最大流算法详解:从水管网络到Ford-Fulkerson与Dinic实战

最大流算法详解:从水管网络到Ford-Fulkerson与Dinic实战

1. 从水管网络到最大流:一个核心问题的诞生想象一下,你是一个城市供水系统的总工程师。你的城市有多个水源(水库),需要通过一个复杂的地下管道网络,将水输送到各个居民区。每条管道都有其最大通水能力&…

2026/8/2 0:00:38 阅读更多 →
基于Springboot的企业门户网站(源码+LW+调试文档+讲解)

基于Springboot的企业门户网站(源码+LW+调试文档+讲解)

温馨提示:本人主页置顶文章(点我)开头有 CSDN 平台官方提供的学长联系方式的名片! 温馨提示:本人主页置顶文章(点我)开头有 CSDN 平台官方提供的学长联系方式的名片! 温馨提示:本人主页置顶文章(点我)开头有 CSDN 平台…

2026/8/3 1:53:31 阅读更多 →
MATLAB xcorr函数详解:从互相关原理到四大实战应用

MATLAB xcorr函数详解:从互相关原理到四大实战应用

1. 从一次信号“找茬”说起:为什么我们需要互相关几年前,我在处理一组声学传感器数据时遇到了一个棘手的问题。我有两个麦克风记录了一段相同的音频信号,理论上它们接收到的声音波形应该非常相似,只是由于麦克风位置不同&#xff…

2026/8/2 0:00:38 阅读更多 →

月新闻

免费解锁百度网盘SVIP加速:macOS用户必备的下载提速终极指南

免费解锁百度网盘SVIP加速:macOS用户必备的下载提速终极指南

免费解锁百度网盘SVIP加速:macOS用户必备的下载提速终极指南 【免费下载链接】BaiduNetdiskPlugin-macOS For macOS.百度网盘 破解SVIP、下载速度限制~ 项目地址: https://gitcode.com/gh_mirrors/ba/BaiduNetdiskPlugin-macOS 还在为百度网盘macOS版的龟速下…

2026/8/2 6:34:16 阅读更多 →
终极ncmdump指南:3分钟实现网易云NCM音乐解密与格式转换

终极ncmdump指南:3分钟实现网易云NCM音乐解密与格式转换

终极ncmdump指南:3分钟实现网易云NCM音乐解密与格式转换 【免费下载链接】ncmdump 项目地址: https://gitcode.com/gh_mirrors/ncmd/ncmdump 还在为网易云音乐下载的NCM格式文件无法在其他播放器播放而烦恼吗?ncmdump解密工具帮你轻松解决这个困…

2026/8/2 2:47:48 阅读更多 →
HarmonyOS 应用开发《掌上英语》第81篇: 智能体卡片:为英语学习 App 打造桌面级学习助手

HarmonyOS 应用开发《掌上英语》第81篇: 智能体卡片:为英语学习 App 打造桌面级学习助手

AgentCard 智能体卡片:为英语学习 App 打造桌面级学习助手适用平台:HarmonyOS 7.0 (API 26 Beta)一、引言 HarmonyOS 7.0(API 26 Beta)新增了 AgentCard 智能体卡片能力,这是继 HMAF(鸿蒙智能体框架&#x…

2026/8/2 0:23:22 阅读更多 →