如何高效随机选取多个不重叠的等长子序列

陌丽君_3179

陌丽君_3179

2026-10-28

687人浏览

原创

如何高效随机选取多个不重叠的等长子序列

本文介绍一种时间与空间复杂度均优、且能保证均匀随机性的算法:通过“虚拟分隔点”映射法,从给定序列中无偏地抽取 m 个长度为 n 的互不重叠连续子序列。

本文介绍一种时间与空间复杂度均优、且能保证均匀随机性的算法:通过“虚拟分隔点”映射法,从给定序列中无偏地抽取 m 个长度为 n 的互不重叠连续子序列。

在序列分析、数据采样、强化学习经验回放或测试用例生成等场景中,常需从原始序列中随机、均匀、无偏地选出若干个固定长度、彼此不重叠(即索引区间无交集)的连续子序列。朴素的贪心逐次采样(如先随机选一个子序列,再在剩余片段中递归选取)会因路径依赖导致概率失衡——某些合法解可能永远无法生成(例如首段过早占据中间位置,导致后续无法容纳足够子序列),因而不满足均匀性要求。

更优解法的核心思想是:将“选择 M 个不重叠长度为 N 的子序列”转化为“在压缩后的索引空间中均匀随机选取 M 个起始位置”,再通过偏移校正还原为真实起始索引。

? 原理简析

设原序列长度为 L,需选 M 个长度为 N 的子序列。每个子序列占据 N 个连续位置,M 个共占 M × N 位;剩余 L − M × N 个“空闲位”需被分配到 M + 1 个间隙中(子序列前、之间、之后),每个间隙可为空(≥0)。这等价于求 M + 1 个非负整数之和为 L − M × N 的组合数。

但直接枚举间隙大小效率低。巧妙的替代方案是:

  • 构造一个长度为 L − M × (N − 1) 的“逻辑索引空间”(即预留出每个子序列后 N−1 位作为其“专属延伸区”);
  • 在该空间中均匀随机选取 M 个不重复起始位置(即 random.sample(range(L - M*(N-1)), M));
  • 将选中的位置排序后,第 i 个位置 d_i 对应真实起始索引为 d_i + i × (N − 1) —— 这一偏移恰好为前 i 个子序列“腾出”了它们各自的 N−1 个后续位置,从而天然保证子序列互不重叠且覆盖完整。

该映射是一一对应(bijection),因此均匀采样逻辑位置 → 均匀采样所有合法子序列组合。

✅ 实现代码

import random

def sample_nonoverlapping_subsequences(seq, M, N):
    """
    从序列 seq 中均匀随机选取 M 个长度为 N 的不重叠连续子序列。

    Args:
        seq: 输入序列(list/tuple)
        M: 子序列数量
        N: 每个子序列长度

    Returns:
        list[list]: M 个子序列组成的列表,每个为 seq 的切片

    Raises:
        ValueError: 若 M*N > len(seq),无解
    """
    L = len(seq)
    if M * N > L:
        raise ValueError(f"Impossible: need {M}*{N}={M*N} positions, but sequence has only {L}")

    # 逻辑空间大小:预留 M*(N-1) 位用于子序列内部延伸
    logical_size = L - M * (N - 1)

    # 随机选取 M 个逻辑起始点(无序,自动去重)
    dividers = sorted(random.sample(range(logical_size), M))

    # 映射为真实起始索引:第 i 个子序列向后偏移 i*(N-1)
    starts = [d + i * (N - 1) for i, d in enumerate(dividers)]

    # 提取子序列
    return [seq[start : start + N] for start in starts]

# 示例使用
seq = list(range(1, 17))  # [1, 2, ..., 16]
samples = sample_nonoverlapping_subsequences(seq, M=3, N=3)
print(samples)
# 可能输出:[[2, 3, 4], [7, 8, 9], [12, 13, 14]]

⚠️ 注意事项与验证要点

  • 时间复杂度:O(M log M)(主要开销在 sorted(),random.sample 为 O(M));空间复杂度:O(M + N×M)(仅存储结果与临时索引);远优于暴力枚举(指数级)或回溯(最坏多项式但常数大)。
  • 均匀性保障:因 random.sample 在逻辑空间中严格均匀,且映射为双射,故所有 C(L − M×(N−1), M) 种合法解等概率出现。
  • 边界安全:函数内置长度校验,避免无效输入;range(logical_size) 自动处理 logical_size == 0(即刚好填满)的情形。
  • 可复现性:如需调试或测试,可在调用前设置 random.seed(42)。

? 总结

该算法摒弃了易导致偏差的贪心策略,转而利用组合数学中的“星与条”(stars and bars)思想进行空间压缩与坐标映射,以极简代码实现高效、正确、可证明均匀的随机采样。它适用于任意支持索引切片的序列类型(list, array, str 等),是处理不重叠子序列随机选取任务的推荐标准解法。

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1831

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4884

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1789

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

29117

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3487

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3487

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1223

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

616

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2703

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
热门推荐
/
最新课程
phpStudy极速入门视频教程
phpStudy极速入门视频教程

共6课时 | 54.7万人学习

独孤九贱(4)_PHP视频教程
独孤九贱(4)_PHP视频教程

共89课时 | 134.1万人学习