[论文解读] The area generating function for simplex-duplex polyominoes

Svjetlan Feretić|ArXiv.org|Oct 26, 2009

Advanced Combinatorial Mathematics参考文献 5被引用 3

一句话总结

本文引入了单峰双组件多格骨牌——一种列双组件多格骨牌，其相邻列中不会同时出现两个连通分量——并使用改进的Temperley方法推导出其面积生成函数。该方法通过逐次添加一个或两个列来避免双组件列带来的复杂性，并求解涉及生成函数一阶与二阶导数的函数方程，最终得到n格单峰双组件多格骨牌的渐近增长率为约0.119443 × 3.522020^n。

ABSTRACT

Back in the early days of polyomino enumeration, a model called column-convex polyominoes was introduced and its area generating function was found. That generating function is rational: the numerator has degree four and the denominator has degree three. Let a column-duplex polyomino be a polyomino whose columns can have either one or two connected components. A simplex-duplex polyomino is a column-duplex polyomino in which there is no occurrence of two adjacent columns each having two connected components. Simplex-duplex polyominoes are not easy to deal with, but their area generating function can still be found. To find this generating function, we use an upgraded version of the Temperley method. Though that technique is widely used in these times, our application presents two interesting features. Firstly, we add one or two columns at a time, thus bypassing those simplex-duplex polyominoes which end with a two-component column. (It is somewhat more usual to add just one column at a time.) Secondly, we obtain a functional equation that involves both the first and the second derivatives of the sought-for generating function. (Such equations usually involve the first derivative only. In some cases, no derivative is involved at all.) Right because of this latter interesting feature, the Temperley method produces a very complicated formula for the generating function. Anyway, from that formula it is easy to compute Taylor polynomials. Thus we get plenty of evidence that the number of n-celled simplex-duplex polyominoes behaves asymptotically as 0.119443*3.522020^n. For comparison, the number of n-celled column-convex polyominoes behaves asymptotically as 0.180916*3.205569^n.

研究动机与目标

• 开发一种可解模型，以在原本难以精确计数的列凸多格骨牌与列双组件多格骨牌之间实现插值。
• 定义并分析单峰双组件多格骨牌，即列双组件多格骨牌的一个子类，其特征为任意相邻两列不会同时具有两个连通分量。
• 推导出双变量生成函数G(q,w)，用于按面积（q）和双组件列数量（w）计数单峰双组件多格骨牌。
• 计算面积生成函数G(q,1)，并确定其渐近行为，以获得n格单峰双组件多格骨牌的总数。

提出的方法

• 采用升级版Temperley方法，通过逐次添加一个或两个列来避免双组件列带来的复杂性。
• 构建一个涉及生成函数一阶与二阶导数的函数方程，此特征非标准，显著增加了复杂度。
• 求解该函数方程，得到以单组件列结尾的单峰双组件多格骨牌的生成函数。
• 从解中重构完整的生成函数G(q,1)，从而可计算泰勒级数系数。
• 通过高阶泰勒多项式（最高达320次）进行数值计算，以分析渐近行为。
• 通过考察相邻系数之比并以增长常数进行归一化，执行渐近分析。

实验结果

研究问题

• RQ1 n格单峰双组件多格骨牌的数量的渐近增长速率是多少？与列凸多格骨牌相比有何差异？
• RQ2 尽管模型复杂且缺乏简单递推结构，是否仍可推导出单峰双组件多格骨牌的面积生成函数？
• RQ3 在相邻列中双组件列受限制的条件下，该限制如何影响生成函数的结构与奇点性质？
• RQ4 在函数方程中使用二阶导数在多大程度上影响了生成函数的可解性与复杂度？
• RQ5 单峰双组件多格骨牌的最终生成函数是否比列凸多格骨牌的更复杂？这对未来推广意味着什么？

主要发现

• 基于G(q,1)的320次泰勒多项式分析，n格单峰双组件多格骨牌的数量渐近行为为0.119442870405 × 3.522019812882^n。
• 面积生成函数G(q,1)的主奇点为位于约0.283927988236处的一阶极点。
• 当w=0时，列凸多格骨牌的面积生成函数作为双变量生成函数的特例被恢复，证实与已知结果一致。
• 推导中使用的函数方程涉及生成函数的一阶与二阶导数，此特征极为罕见且非平凡，显著增加了公式的复杂度。
• 增长常数3.522020远高于列凸多格骨牌的增长常数（3.205569），表明该类多格骨牌集合更为密集。
• 尽管生成函数复杂，但通过高精度数值计算泰勒系数，其渐近行为得到高度可信的验证。