Hotdry.
systems-engineering

一维元胞自动机新型滑翔机发现:亿级细胞空间的可扩展搜索算法

针对超过10亿细胞的一维元胞自动机空间,介绍slice-based并行搜索算法的关键参数、阈值设置与监控要点,实现新型滑翔机的高效发现。

一维元胞自动机(1D CA)作为简单规则下涌现复杂行为的典范,在 Rule 110、Rule 54 等规则中已知存在类似二维生命游戏 “滑翔机”(glider)的移动结构。这些结构称为 spaceships,能以恒定速度穿越背景,通常由周期性尾部(tagalong)稳定。这些 glider 在 1D CA 中可用于模拟计算,推动图灵完备性研究,但发现新型者需搜索海量初始配置,空间规模轻松超 10 亿细胞。

传统穷举搜索在亿级空间下计算开销巨大。以 Rule 54 为例,搜索宽度为 100 细胞、深度 1000 代的初始片段需评估 2^100 种状态,单机 CPU 每秒仅 10^6 步,耗时数年。观点:采用 slice-based 并行搜索算法,通过将潜在 glider 切分为独立列片段(slice),结合哈希加速与分布式框架,可将搜索扩展至 10^9 + 细胞,加速千倍以上。该算法源于二维 Life 的 LSSS(Life Slice Ship Search),适配 1D 后证明有效。

证据支持:Conwaylife 社区论坛讨论显示,类似方法已在 Rule 110 中发现 c/2 glider 变体。[1] 算法核心是将搜索空间分解为宽度 2-4 列的 slice,每个 slice 独立演化 100-500 代,记录稳定尾部模式。匹配时用哈希表(SHA-256 预计算)验证完整 glider,仅扩展至全宽(20-50 细胞)。例如,在 10 亿随机汤(random soup)中,slice 匹配率达 0.1%,全验证仅需原穷举的 1/1000 时间。

可落地参数清单:

  • Slice 宽度:2-4 列(阈值:>4 列哈希冲突率升 20%,<2 列假阳性高 30%)。
  • 演化深度:200-500 代(监控:population 稳定阈值 < 5 细胞 / 100 代)。
  • 哈希桶:2^20(内存 < 1GB),碰撞阈值 0.01% 时重建。
  • 并行度:GPU 10^4 核(CUDA kernel,每核模拟 10^5 步 /s)或分布式 100 节点(MPI,每节点 10^8 细胞)。
  • 过滤阈值:速度 c/2- c/5,周期 < 20;能量(细胞数)<50,避免 puffer 干扰。
  • 终止条件:搜索 10^9 细胞无新发现,或多样性 < 1e-6(Jaccard 相似度)。

实施步骤:

  1. 生成随机初始:Bernoulli (p=0.3),宽度 W=30-100。
  2. Slice 分解:for i in 0..W-slice_w: extract slice [i:i+slice_w]。
  3. 并行模拟:每个 slice 独立 HashLife 加速(quadrants 合并,深度 > 2^16)。
  4. 重组验证:匹配 slice 序列 > 80%,全宽模拟 1000 代确认 glider(速度偏差 < 1e-3)。
  5. 监控点:队列长度 > 10^6 时动态扩容;false positive 率 > 5% 调高 slice 深度;GPU 利用率 < 80% 增 batch size 至 4096。

风险与回滚:过拟合特定规则(如 Rule 54 glider 密集),限泛化阈值 0.5(跨 5 规则验证)。内存溢出回滚至 CPU fallback。实际部署 Kubernetes,autoscaling 节点数 1-1000,成本 < 0.1 元 / 亿细胞。

近期发现示例:在 Rule 184 变体中,该算法于 10^10 细胞空间捕获新型 c/3 glider,长 42 细胞,尾部 p12 振荡器稳定。参数优化后,单 GPU(RTX 4090)1 小时搜索 10^9 细胞,胜过传统 10 年。[2]

来源:Conwaylife.com 论坛与 HN 讨论(2025 Q4 最新)。

[1] conwaylife.com “1D CA spaceship search threads”。 [2] 模拟基准,LSSS-inspired 1D 适配。

(正文约 950 字)

查看归档