优化滑动窗口中位数:使用延迟删除和双堆方法解决Python TLE问题

优化滑动窗口中位数:使用延迟删除和双堆方法解决Python TLE问题

本文深入探讨了滑动窗口中位数问题,并针对传统双堆方法中因低效移除操作导致的超时(TLE)问题,提出了一种基于延迟删除策略的优化方案。通过将元素与索引绑定并利用自定义堆实现,该方案避免了昂贵的O(K)移除操作,将时间复杂度从O(NK)有效降低至O(N log K),从而在大规模数据集上实现了高性能。

1. 问题概述:滑动窗口中位数

滑动窗口中位数问题要求在一个整数数组 nums 上,维护一个大小为 k 的滑动窗口。当窗口从左向右移动时,每次移动一个位置,都需要计算并返回当前窗口内的中位数。中位数是排序后位于中间的那个数;如果窗口大小 k 是偶数,则中位数是中间两个数的平均值。

例如,对于 nums = [1,3,-1,-3,5,3,6,7] 和 k = 3:

窗口 [1,3,-1],中位数是 1窗口 [3,-1,-3],中位数是 -1…依此类推

2. 传统双堆方法的性能瓶颈

解决滑动窗口中位数问题的一个常见且高效的思路是使用两个堆(一个最大堆 small 存储较小的一半元素,一个最小堆 large 存储较大的一半元素)。通过维护这两个堆,可以O(1)或O(log K)地获取中位数,并在O(log K)时间内添加新元素。

然而,当需要从窗口中“移除”一个旧元素时,问题就出现了。原始代码中的 popNum 方法采用了以下逻辑:

def popNum(self, num):    if num > (self.small[0] * -1): # 假设small[0]是最大堆的堆顶,其真实值是负数        self.large.remove(num)        heapq.heapify(self.large)    else:        self.small.remove(num * -1)        heapq.heapify(self.small)    self.balance()

这里的瓶颈在于 list.remove(num) 和 heapq.heapify(heap)。

立即学习“Python免费学习笔记(深入)”;

list.remove(num):在Python列表中查找并移除一个元素的时间复杂度是O(N),其中N是列表的当前大小(在这里是K)。heapq.heapify(heap):在移除元素后,堆的结构被破坏,需要调用 heapify 来重新构建堆,其时间复杂度也是O(N)(或O(K))。

因此,每次移除操作的复杂度高达O(K)。在一个包含 N 个元素的数组上进行 N-K+1 次窗口滑动,每次滑动都涉及一次移除和一次添加,导致总时间复杂度飙升至 O(NK)。当 N=100000 且 K=50000 时,`NK的量级将达到5 * 10^9`,这远远超出了通常的计算限制,从而导致“时间限制超出”(TLE)错误。

3. 优化的双堆方法:延迟删除策略

为了解决移除操作的效率问题,我们可以采用“延迟删除”(Lazy Deletion)策略。这种方法的核心思想是:当一个元素离开窗口时,我们不立即从堆中物理删除它,而是对其进行“标记”。当我们需要从堆中获取元素(如peek或pop)时,如果遇到已标记为“删除”的元素,就跳过它,直到找到一个有效元素。

实现延迟删除需要解决几个关键问题:

唯一标识元素: 窗口中可能存在重复的数值。为了区分它们,我们需要将每个元素与其在原始数组中的索引绑定,形成 (值, 索引) 对。例如,[1, 5, 1] 中有两个 1,通过 (1, 0) 和 (1, 2) 可以区分。标记删除: 如何高效地标记一个元素为“已删除”?最简单的方法是维护一个窗口的起始索引 lowindex。任何 (值, 索引) 对中,如果 索引 自定义堆封装: Python的 heapq 模块提供了堆的基本操作,但不支持延迟删除。我们需要封装 heapq,创建自定义的堆类来处理 (值, 索引) 对和延迟删除逻辑。

4. 核心实现细节

我们将构建两个自定义的堆类:MinWindowHeap(最小堆)和 MaxWindowHeap(最大堆),以及一个 Solution 类来协调它们。

4.1 自定义堆结构:MinWindowHeap 和 MaxWindowHeap

这两个类将封装 heapq 的功能,并添加延迟删除的逻辑。

MinWindowHeap (最小堆):

__init__(self, conv=lambda x: x): 初始化内部列表 self.heap。conv 是一个转换函数,用于处理元素入堆时的形式(例如,MaxWindowHeap 会用它来取反值)。self.lowindex 记录当前窗口的起始索引。peek(self): 这是延迟删除的关键。它会循环检查堆顶元素:如果堆顶元素的索引小于 self.lowindex,说明该元素已过期,将其弹出并继续检查下一个堆顶元素,直到找到一个有效的(未过期的)元素或者堆为空。push(self, item): 将 item((值, 索引) 对)通过 conv 函数转换后推入堆。pop(self): 先调用 peek() 确保获取到的是有效堆顶,然后从 self.heap 中实际弹出该元素。

MaxWindowHeap (最大堆):

继承自 MinWindowHeap。__init__(self): 调用 super().__init__(negate),其中 negate 函数用于将 (值, 索引) 对的 值 取反,从而将 heapq 的最小堆行为模拟成最大堆。

4.2 Solution 类中的核心逻辑

Solution 类将协调 small (MaxWindowHeap) 和 large (MinWindowHeap) 两个堆,实现滑动窗口中位数的功能。

rebalance(self, add):

维护一个 self.balance 计数器,表示 large 堆比 small 堆多出的元素数量(或负数表示 small 堆多出)。当 abs(self.balance) 大于1时,说明两个堆不平衡,需要将一个堆的顶部元素移动到另一个堆,直到平衡。此方法在每次插入或逻辑删除元素后被调用。

insert(self, item):

item 是 (值, 索引) 对。根据 item 的值与 large 堆的堆顶元素比较,决定将其插入 small 堆还是 large 堆。插入后调用 rebalance 保持堆的平衡。

remove(self, item):

item 是 (值, 索引) 对,表示要移除的旧元素。关键优化: 它不实际从堆中移除元素,而是通过更新 self.large.lowindex 和 self.small.lowindex 为 item[1] + 1 来标记所有索引小于等于 item[1] 的元素为“过期”。然后调用 rebalance 来调整平衡计数。

getMedian(self):

根据 self.balance 的值和两个堆的堆顶元素(通过 peek() 获取有效元素)来计算当前窗口的中位数。如果 balance 为0,取两个堆顶的平均值;否则取元素较多那个堆的堆顶。

medianSlidingWindow(self, nums, k):

主函数,用于解决问题。首先将 nums 转换为 (值, 索引) 对的列表 items。初始化 small 和 large 堆,并用前 k 个元素填充第一个窗口。计算第一个窗口的中位数并添加到结果列表。然后通过循环滑动窗口:zip(items, items[k:]) 用于同时获取离开窗口的 olditem 和进入窗口的 item。对 olditem 调用 remove(),对 item 调用 insert()。每次滑动后计算当前窗口的中位数并添加到结果。

5. 完整解决方案代码

import heapq# 辅助函数:将 (值, 索引) 对的值取反,用于模拟最大堆def negate(item):    return -item[0], item[1]class MinWindowHeap(object):    def __init__(self, conv=lambda x: x):        self.heap = []        self.conv = conv  # 转换函数 (例如,用于MaxHeap取反值)        self.lowindex = 0  # 窗口下限索引,用于识别已删除项    def peek(self):  # 返回 (值, 索引) 或 None (如果堆为空或仅包含已删除项)        while self.heap:            # 转换堆顶元素,例如 MaxWindowHeap 会将值取反            item = self.conv(self.heap[0])             if item[1] >= self.lowindex: # 如果索引在当前窗口内,则有效                return item            # 元素已过期(索引小于lowindex),从堆中弹出            heapq.heappop(self.heap)        return None # 堆中没有有效元素    def push(self, item):        # 将 (值, 索引) 对通过转换函数推入堆        heapq.heappush(self.heap, self.conv(item))    def pop(self):        item = self.peek()  # 获取有效堆顶,同时清除所有过期的堆顶        if item:            heapq.heappop(self.heap) # 实际弹出有效堆顶        return item # 返回被弹出的有效元素class MaxWindowHeap(MinWindowHeap):    def __init__(self):        # Python 3 中 super() 可以不带参数        super(MaxWindowHeap, self).__init__(negate) # 使用negate函数将最小堆模拟为最大堆class Solution(object):    def rebalance(self, add):        """        调整两个堆的平衡。        add > 0 表示 large 堆增加了元素,或 small 堆移除了元素。        add < 0 表示 small 堆增加了元素,或 large 堆移除了元素。        """        self.balance += add        if abs(self.balance)  1: # large 堆元素过多,需要移动一个到 small 堆            self.small.push(self.large.pop())        elif self.balance = pivot[0])        heap = self.large if islarge else self.small        heap.push(item)        # 更新平衡计数并重新平衡        self.rebalance(1 if islarge else -1)    def remove(self, item):        """        逻辑删除一个元素(通过更新 lowindex 标记过期)并调整平衡。        item 是 (值, 索引) 对,表示离开窗口的旧元素。        """        # 获取 large 堆的堆顶作为枢轴        pivot = self.large.peek()        # 判断旧元素是属于 large 堆还是 small 堆        islarge = (pivot is not None) and (item[0] >= pivot[0])        # 关键步骤:更新两个堆的 lowindex,标记所有索引小于等于 item[1] 的元素为过期        # 这意味着窗口向右移动了,item[1] 及其之前的所有元素都可能已过期        self.large.lowindex = self.small.lowindex = item[1] + 1        # 更新平衡计数并重新平衡        self.rebalance(-1 if islarge else 1)    def getMedian(self):        """        获取当前窗口的中位数。        """        if self.balance == 0: # 两个堆大小相同,中位数是两个堆顶的平均值            return (self.large.peek()[0] + self.small.peek()[0]) * 0.5        # 否则,元素较多的那个堆的堆顶就是中位数        return self.large.peek()[0] if self.balance > 0 else self.small.peek()[0]    def medianSlidingWindow(self, nums, k):        """        滑动窗口中位数的主函数。        :type nums: List[int]        :type k: int        :rtype: List[float]        """        self.small = MaxWindowHeap() # 存储较小一半元素的堆 (最大堆)        self.large = MinWindowHeap() # 存储较大一半元素的堆 (最小堆)        self.balance = 0             # 记录 large 堆与 small 堆的元素数量差        # 将原始数组转换为 (值, 索引) 对,以便唯一标识元素        items = [(val, i) for i, val in enumerate(nums)]        # 填充第一个窗口        for item in items[:k]:            self.insert(item)        result = [self.getMedian()] # 记录第一个窗口的中位数        # 滑动窗口        # olditem 是离开窗口的元素,item 是进入窗口的元素        for olditem, item in zip(items, items[k:]):            self.remove(olditem) # 逻辑删除旧元素            self.insert(item)    # 插入新元素            result.append(self.getMedian()) # 计算并记录当前窗口的中位数        return result

6. 注意事项与总结

时间复杂度优化: 通过延迟删除策略,remove 操作不再需要遍历列表或重建堆。更新 lowindex 是 O(1) 操作。peek 和 pop 操作在遇到过期元素时,每次弹出都是 O(log K),但每个元素只会因过期而被弹出一次。因此,平均来看,insert、remove 和 getMedian 的操作都保持在 O(log K)。整个算法的时间复杂度从 O(N*K) 优化到了 O(N log K),在大规模数据集上性能显著提升。空间复杂度: 两个堆最多存储 K 个 (值, 索引) 对,所以空间复杂度为 O(K)。处理重复值: 将值与索引绑定 (值, 索引) 是处理数组中重复值的关键,确保每个元素都能被唯一标识和跟踪。Python heapq 模块: heapq 默认是最小堆。为了实现最大堆,我们通过存储元素的负值来实现。平衡维护: rebalance 方法确保 small 和 large 两个堆的大小差异不超过1,这是高效获取中位数的基础。

这种优化的双堆方法提供了一个健壮且高效的解决方案,适用于处理大规模数据下的滑动窗口中位数问题。它巧妙地利用了堆的特性和延迟删除机制,避免了传统方法中的性能瓶颈。

以上就是优化滑动窗口中位数:使用延迟删除和双堆方法解决Python TLE问题的详细内容,更多请关注创想鸟其它相关文章!

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。
如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件至 chuangxiangniao@163.com 举报,一经查实,本站将立刻删除。
发布者:程序猿,转转请注明出处:https://www.chuangxiangniao.com/p/1374986.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
在 CPU 上运行任何量化的 GGUF 模型进行本地推理的教程
上一篇 2025年12月14日 14:37:08
Python实现弗洛伊德三角形:从基础到优化
下一篇 2025年12月14日 14:37:31

相关推荐

  • Matplotlib 地图中多类型图例的创建与优化

    Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化Matplotlib 地图中多类型图例的创建与优化

    本教程旨在解决matplotlib地图可视化中,如何在一个图例中同时展示颜色块(如区域分类)和自定义标记(如特定兴趣点)的问题。文章详细介绍了当传统`patch`对象无法正确显示标记时,如何利用`matplotlib.lines.line2d`创建标记图例句柄,并将其与颜色块图例句柄合并,从而生成一…

    2026年5月10日 用户投稿
    300
  • 利用海象运算符简化条件赋值:Python教程与最佳实践

    本文旨在探讨Python中海象运算符(:=)在条件赋值场景下的应用。通过对比传统if/else语句与海象运算符,以及条件表达式,分析海象运算符在简化代码、提高可读性方面的优势与局限性。并通过具体示例,展示如何在列表推导式等场景下合理使用海象运算符,同时强调其潜在的复杂性及替代方案,帮助开发者更好地掌…

    2026年5月10日
    100
  • 怎么在PHP代码中实现图片上传功能_PHP图片上传功能实现与安全处理教程

    首先创建含enctype的HTML表单,再用PHP接收文件,检查目录、移动临时文件,验证类型与大小,生成唯一文件名,并调整php.ini限制以确保上传成功。 如果您尝试在PHP项目中添加图片上传功能,但服务器无法正确接收或保存文件,则可能是由于表单配置、文件处理逻辑或安全限制的问题。以下是实现该功能…

    2026年5月10日
    300
  • 比特币新手教程 比特币交易平台有哪些

    比特币是一种去中心化的数字货币,基于区块链技术实现点对点交易,具有匿名性、有限发行和不可篡改等特点;新手可通过交易所购买,P2P交易获得比特币,常用平台包括Binance、OKX和Huobi;交易流程包括注册账户、实名认证、绑定支付方式、充值法币并下单购买,可选择市价单或限价单;比特币存储方式有交易…

    2026年5月10日
    000
  • RichHandler与Rich Progress集成:解决显示冲突的教程

    在使用rich库的`richhandler`进行日志输出并同时使用`progress`组件时,可能会遇到显示错乱或溢出问题。这通常是由于为`richhandler`和`progress`分别创建了独立的`console`实例导致的。解决方案是确保日志处理器和进度条组件共享同一个`console`实例…

    2026年5月10日
    000
  • 使用 Jupyter Notebook 进行探索性数据分析

    Jupyter Notebook通过单元格实现代码与Markdown结合,支持数据导入(pandas)、清洗(fillna)、探索(matplotlib/seaborn可视化)、统计分析(describe/corr)和特征工程,便于记录与分享分析过程。 Jupyter Notebook 是进行探索性…

    2026年5月10日
    000
  • 深入理解 Express.js 中 next() 参数的作用与中间件机制

    本文深入探讨 express.js 中间件函数中的 `next()` 参数。它负责将控制权传递给请求-响应周期中的下一个中间件或路由处理程序。文章将详细解释 `next()` 的工作原理、中间件的注册与执行顺序,以及不正确使用 `next()` 可能导致请求挂起的风险,并通过代码示例和实际应用场景,…

    2026年5月10日
    000
  • Python命令怎样使用profile分析脚本性能 Python命令性能分析的基础教程

    使用Python的cProfile模块分析脚本性能最直接的方式是通过命令行执行python -m cProfile your_script.py,它会输出每个函数的调用次数、总耗时、累积耗时等关键指标,帮助定位性能瓶颈;为进一步分析,可将结果保存为文件python -m cProfile -o ou…

    2026年5月10日
    000
  • Python递归函数追踪与性能考量:以序列打印为例

    本文深入探讨了Python中一种递归打印序列元素的方法,并着重演示了如何通过引入缩进参数来有效追踪递归函数的执行流程和参数变化。通过实际代码示例,文章揭示了递归调用可能带来的潜在性能开销,特别是对调用栈空间的需求,以及Python默认递归深度限制可能导致的错误,为读者提供了理解和优化递归算法的实用见…

    2026年5月10日
    000
  • python中zip函数详解 python多序列压缩zip函数应用场景

    zip函数的应用场景包括:1) 同时遍历多个序列,2) 合并多个列表的数据,3) 数据分析和科学计算中的元素运算,4) 处理csv文件,5) 性能优化。zip函数是一个强大的工具,能够简化代码并提高处理多个序列时的效率。 在Python中,zip函数是一个非常有用的工具,它能够将多个可迭代对象打包成…

    2026年5月10日
    000
  • c++如何实现UDP通信_c++基于UDP的网络通信示例

    UDP通信基于套接字实现,适用于实时性要求高的场景。1. 流程包括创建套接字、绑定地址(接收方)、发送(sendto)与接收(recvfrom)数据、关闭套接字;2. 服务端监听指定端口,接收客户端消息并回传;3. 客户端发送消息至服务端并接收响应;4. 跨平台需处理Winsock初始化与库链接,编…

    2026年5月10日
    100
  • Python中怎样使用pymongo?

    在python中使用pymongo可以轻松地与mongodb数据库进行交互。1)安装pymongo:pip install pymongo。2)连接到mongodb:from pymongo import mongoclient; client = mongoclient(‘mongod…

    2026年5月10日
    000
  • Golang空接口如何应用在项目中

    空接口可用于接收任意类型值,常见于日志函数、通用数据结构、JSON动态解析及配置驱动逻辑,提升代码灵活性,但需配合类型断言确保安全,避免滥用以降低维护成本。 空接口 interface{} 在 Go 语言中是一个非常灵活的类型,它可以存储任何类型的值。虽然它牺牲了一部分类型安全,但在实际项目中合理使…

    2026年5月10日
    100
  • JavaScript计算器开发:解决数值显示与初始化问题

    本教程深入探讨了使用JavaScript构建计算器时常见的数值显示异常问题,特别是由于类属性未初始化导致的`Cannot read properties of undefined`错误。我们将详细分析问题根源,并通过在构造函数中调用初始化方法来解决该问题,同时优化显示逻辑,确保计算器功能稳定且界面显…

    2026年5月10日
    000
  • Python 函数参数类型:如何使用可变参数和动态参数?

    python 中的参数类型:关键词参数、可变参数和动态参数 在 python 中,函数的参数可以分为以下几种类型: 关键词参数(kw)**:这些参数具有名称,并且在调用函数时明确指定。可变参数(*args):这些参数没有名称,允许函数接受任意数量的位置参数。它们将被收集到一个元组中。动态参数(kwa…

    2026年5月10日
    000
  • Circle为何在凌晨向Solana新增铸造5亿枚USDC?USDC增发原因与对SOL生态影响深度解析

    近日,链上数据显示,Circle 在凌晨向 Solana 链新增铸造了 5亿枚USDC。此次大规模增发引起市场关注,投资者需要了解背后的原因以及对 Solana 生态的潜在影响。 USDC增发原因分析 增发 USDC 的主要原因可能包括: 满足市场需求:近期 Solana 上交易活动活跃,USDC …

    2026年5月10日
    000
  • 使用 Ajax 和 FormData 实现文件上传及文本数据提交的完整教程

    本文旨在解决在使用 Ajax 和 FormData 进行文件上传时,遇到的 $_POST 和 $_FILES 为空的问题。通过详细的代码示例和解释,我们将展示如何正确地构建 FormData 对象,并通过 Ajax 将文件和文本数据发送到服务器端,同时避免常见的错误配置,确保数据能够成功地被 PHP…

    2026年5月10日
    000
  • pycharm解析器怎么添加 解析器添加详细流程

    在pycharm中添加解析器的步骤包括:1) 打开pycharm并进入设置,2) 选择project interpreter,3) 点击齿轮图标并选择add,4) 选择解析器类型并配置路径,5) 点击ok完成添加。添加解析器后,选择合适的类型和版本,配置环境变量,并利用解析器的功能提高开发效率。 在…

    2026年5月10日
    100
  • python中numpy的用法

    NumPy是Python中用于科学计算的强大库,它提供了以下功能:多维数组处理矩阵运算快速傅里叶变换(FFT)线性代数随机数生成 NumPy在Python中的强大功能 NumPy是Python中用于科学计算的一个强大且灵活的库。它提供了用于处理多维数组和矩阵的一组高效工具,是数据分析和机器学习项目的…

    2026年5月10日
    100
  • 深入理解MQTT多级通配符#的用法限制与Paho-MQTT订阅实践

    本文旨在解析mqtt多级通配符`#`在订阅主题时的严格使用规则,尤其是在paho-mqtt库中遇到的`valueerror: ‘invalid subscription filter.’`问题。我们将详细阐述mqtt规范中关于`#`必须作为主题过滤器最后一个字符的规定,并通过…

    2026年5月10日
    000

发表回复

登录后才能评论
关注微信