[论文解读] More Kolakoski Sequences

Bernd Sing|arXiv (Cornell University)|Sep 21, 2010

Fractal and DNA sequence analysis参考文献 1被引用 5

一句话总结

本文将经典的Kolakoski序列推广至任意双字母字母表{r,s}，重点关注其中一字母为奇数、另一字母为偶数的情形——此类情形无法约化为已知的替换系统。通过广义替换规则分析字母频率、回文、平方子串及子串复杂度，并推导出C∞-词数量的界，表明复杂度随指数δ = ln(r+s)/ln((r+s)/2)呈多项式增长，且{1,4}与{2,3}的数值数据充分支持该结论。

ABSTRACT

Our goal in this article is to review the known properties of the mysterious Kolakoski sequence and at the same time look at generalizations of it over arbitrary two letter alphabets. Our primary focus will here be the case where one of the letters is odd while the other is even, since in the other cases the sequences in question can be rewritten as (well-known) primitive substitution sequences. We will look at word and letter frequencies, squares, palindromes and complexity.

研究动机与目标

• 将经典的Kolakoski序列扩展至任意双字母字母表{r,s}（r ≠ s），特别关注奇偶性混合的情形。
• 分析广义Kolakoski序列的结构性质，包括字母频率、回文、平方子串及子串复杂度。
• 为不同字母表（特别是{1,4}与{2,3}）建立C∞-词（无限Kolakoski序列的子串）数量的界。
• 探讨Kolakoski序列的子串集合是否与C∞-词集合一致，这是该领域的一个关键猜想。

提出的方法

• 通过交替的替换规则σ₀与σ₁构造广义Kolakoski序列，其中σ₀作用于偶数位置，σ₁作用于奇数位置，定义为σ₀(r) = r^r，σ₀(s) = r^s，σ₁(r) = s^r，σ₁(s) = s^s。
• 对于奇偶性混合的字母表{r,s}，序列无法约化为原生替换系统，因此需直接分析游程长度与子串结构。
• 通过原生词及其对称性质研究回文子串，特别关注以奇数字母为中心的奇数长度回文。
• 利用替换动力学与矩阵谱性质推导的对数增长指数，对复杂度γ(n)（长度为n的C∞-词数量）进行界定。
• 在双对数坐标系中绘制{1,4}与{2,3}的数值数据，以估计增长指数δ = ln(r+s)/ln((r+s)/2)。
• 基于Dekking的结果并推广至任意{r,s}，推导出复杂度的理论界，得到指数为α与β的上下多项式界。

实验结果

研究问题

• RQ1 广义Kolakoski序列在字母表{r,s}下，长度为n的子串（C∞-词）数量的增长率如何？
• RQ2 当一字母为奇数、一字母为偶数时，广义Kolakoski序列中的字母频率行为如何？
• RQ3 广义Kolakoski序列中，回文子串能否系统性地生成与分类？
• RQ4 广义Kolakoski序列在何种程度上可通过替换规则表示？何时会约化为已知的原生系统？
• RQ5 对于{1,4}与{2,3}等字母表，猜想γ(n) ≈ n^δ（δ = ln(r+s)/ln((r+s)/2)）是否得到数值与理论证据的支持？

主要发现

• 对于经典字母表{1,2}，复杂度γ(n)满足C₁·n².⁷⁰⁸⁷ < γ(n) < C₂·n².⁷¹⁰²，其中C₁、C₂为正常数，表明复杂度呈紧密的多项式增长。
• 所猜想的复杂度指数δ = ln(r+s)/ln((r+s)/2)得到数值数据的充分支持，δ ≈ 2.7095（{1,2}），与经验界一致。
• 对于{1,4}与{2,3}，lnγ(n)与ln n的数值图显示增长趋势与所猜想的δ一致，图1中浅灰与深灰曲线与预测的幂律趋势吻合。
• 在{r,s}中一奇一偶的广义Kolakoski序列中，字母频率不相等，可通过推导出的替换规则计算，但分布不均匀。
• 回文子串源于奇数长度的基本回文词，且中心为奇数字母；对每个长度，此类基本词的数量最多为两倍。
• 当两字母均为奇数时，{r,s}的替换规则产生一个三次Pisot-Vijaraghavan数，当2(r+s) ≥ (r−s)²时成立，表明存在深层代数结构，尽管此情形并非本文分析重点。