MacLaren-Marsaglia算法故事
30寻求更优随机性：一个算法的故事
5W1H分析
What（是什么）
MacLaren-Marsaglia算法是一种通过组合两个不同周期的线性同余生成器（LCG）来生成高质量伪随机数的算法。它使用一个大小为k的打乱表（shuffle table）来打破LCG的线性相关性，显著改善随机数的统计特性。
Why（为什么）
单个LCG生成的随机数存在明显的缺陷：
- 相邻数值之间存在线性相关性
- 在多维空间中呈现格点结构（lattice structure）
- 低位比特的随机性较差
MacLaren和Marsaglia在1965年提出此算法，目的是在不增加计算复杂度的前提下，通过打乱技术消除这些缺陷，生成统计特性更好的随机数序列。
Who（谁）
- 发明者: Donald E. Knuth在TAOCP第2卷第3.2.2节中详细描述了该算法
- 原始贡献者: G. Marsaglia和T.A. MacLaren在1965年的论文中首次提出
- 使用者: 需要高质量伪随机数的科学计算、模拟、密码学应用开发者
When（何时）
- 算法首次发表：1965年
- 被TAOCP收录：1969年（第2卷第1版）
- 适用场景：当需要比简单LCG更好的随机性，但不需要密码学级安全时
Where（何地）
- 应用场景：蒙特卡洛模拟、数值积分、随机采样、游戏开发、测试数据生成
- 算法地位：作为LCG的改进版本，是许多现代随机数生成器的基础
How（如何）
算法核心思想：
- 维护一个大小为k的表V，用LCG1填充
- 每次生成随机数时：
- 用LCG2生成索引 j ∈ 0, k-1
- 输出Vj作为结果
- 用LCG1生成新值替换Vj
这种"洗牌"机制打破了原始LCG的线性相关性，同时保持了两倍LCG的周期长度。
需求定义
功能需求
-
初始化函数
maclaren_marsaglia_init(k, seed1, seed2)- 参数k：表大小（1到1024）
- 参数seed1：LCG1的种子
- 参数seed2：LCG2的种子
- 返回：成功/失败状态
-
随机数生成函数
maclaren_marsaglia_next()- 返回：32位无符号整数 0, 2\^32-1
-
浮点数生成函数
maclaren_marsaglia_next_double()- 返回：[0.0, 1.0) 范围内的双精度浮点数
-
资源清理函数
maclaren_marsaglia_cleanup()- 释放动态分配的内存
非功能需求
- 性能：生成速度应接近单个LCG
- 内存：使用O(k)的额外空间
- 可移植性：纯C99实现，不依赖平台特性
验收标准
功能验收
- k=64时能正常生成1000个随机数
- 不同种子产生不同序列
- k=1时正确工作（退化为LCG1）
- 资源正确分配和释放
统计特性验收
- 卡方检验显示比单个LCG更均匀
- 自相关测试显示相邻数独立性更好
- 周期长度显著长于单个LCG
- 无明显重复模式
代码质量验收
- 使用C99标准编译无警告
- 包含完整的测试框架
- 代码注释清晰，符合TAOCP风格
- 内存管理无泄漏
参考
- TAOCP Volume 2, Section 3.2.2: "The MacLaren-Marsaglia Method"
- Marsaglia, G. and MacLaren, M.D. (1965). "Uniform random number generators"
- LCG参数来源：Numerical Recipes in C (LCG1), glibc (LCG2)