Skip to content

圆与球的致密填充问题(Circle & Sphere Packing)

这组问题涉及几何学中经典的圆堆积(Circle Packing)与球堆积(Sphere Packing)理论。在数学上,这些问题多属于具有极高挑战性的优化问题。以下是详细的研究与解答。


1. 相同尺寸圆的填充数量数列 $a_n$

问题: 半径为 $R$的圆,以半径分别为$R, R/2, R/3 \dots R/n$的圆最大数量填充,求数列$a_n$ 的通项公式。

结论与分析

在数学上,该数列没有精确的、封闭形式的通用通项公式。因为小圆在大圆内的致密排布取决于极其复杂的几何边界重组(属于 NP-hard 问题),目前学术界仅对特定的小 $n$ 值有确切的几何证明。

我们可以给出该数列的前几项精确值以及当 $n \to \infty$ 时的渐近通项公式

前几项精确值

  • $a_1 = 1$:小圆半径为 $R$,即大圆自身,只能放下 1 个。
  • $a_2 = 2$:小圆半径为 $R/2$,直径为 $R$。在大圆的直径上恰好可以并排放置 2 个,剩余空间无法塞下第 3 个。
  • $a_3 = 7$:小圆半径为 $R/3$。最佳排布为蜂窝状的“中心一个,周围一圈环绕六个”,最大数量为 $1 + 6 = 7$。
  • $a_4 = 11$:小圆半径为 $R/4$。经严密数学证明,此时大圆内部会产生松动空隙,但空间不足以放下第 12 个,因此最大值为 11。

渐近通项公式

当 $n \to \infty$ 时,大圆边界对内部排布的影响趋于 0。此时,小圆在大圆内部的排布完美趋近于平面最密堆积——正六边形(蜂窝状)排布,其填充密度(面积占比)极限为: $$\eta = \frac{\pi}{2\sqrt{3}} \approx 90.69\%$$

根据面积占比关系: $$a_n \cdot \frac{\pi R^2}{n^2} \approx \frac{\pi}{2\sqrt{3}} \cdot \pi R^2$$

因此,该数列的渐近通项公式为: $$a_n \sim \frac{\pi}{2\sqrt{3}} n^2 \approx 0.9069 n^2 \quad (n \to \infty)$$


2. 单一尺寸小圆填充的最大比例

问题: 半径为 $R$的圆,用任意$R' < R/2$ 的单一固定尺寸圆最大数量填充,所填充的最大比例为?

结论与分析

为了使填充比例达到最大化,我们应当让挑选的单一小圆半径 $R'$无限趋近于 0(即令$n \to \infty$)。

  • 如果选择较大的 $R'$(例如接近 $R/2$),由于边界处会留下巨大的、无法利用的死角,填充比例会非常低(如 $n=2$时仅为$50\%$)。
  • 当 $R' \to 0$ 时,大圆内部的排布完全逼近平面等圆最密堆积的几何极限。

最大比例

$$\text{Maximum Ratio} = \frac{\pi}{2\sqrt{3}} = \frac{\pi}{\sqrt{12}} \approx 90.69\%$$


3. 多尺寸任意小圆填充的比例

问题: 半径为 $R$ 的圆,用半径任意小的多尺寸圆填充,所填充的比例为?

结论与分析

注意本题与第 2 问的区别:第 2 问限定使用“某一个固定尺寸”的圆;而本题允许使用无数种大小不一、半径任意小的圆进行“嵌套套娃”填充。

在几何学中,这对应著名的分形几何结构——阿波罗尼奥斯垫片(Apollonian Gasket)

  1. 先在大圆内部放入几个尽可能大的、互切的圆;
  2. 在这些圆与大圆、或圆与圆之间留下的弧形孔隙中,塞入能容纳的次大圆;
  3. 在次大圆留下的更小孔隙中,继续塞入更小的圆,如此进行无限次级联填充。

测度论结论

测度论已经严格证明,通过这种无限的多尺度圆填充,未被覆盖的残差空隙的勒贝格测度(即面积)在极限状态下趋于 0。

因此,所填充的极限比例为 $1$(即 $100\%$)


4. 三维空间(球体)的拓展

问题: 将上述问题中的圆换成球体,又如何?

当维度从 2D(圆)拓展到 3D(球体)时,相关结论修正如下:

针对第 1 问:等球堆积数列 $b_n$

  • $b_1 = 1$
  • $b_2 = 2$:在一条直径线上刚好放下 2 个。
  • $b_3 = 13$:半径为 $R/3$的小球。这在三维中等价于接吻数问题(Kissing Number)。在 3D 空间中,一个球的周围最多可以同时外切 12 个等大球。因此,最佳排布是 1 个中心球外围环绕 12 个球,总共$1 + 12 = 13$个,它们恰好能被半径为$3$ 的大球包裹。
  • 渐近通项公式($n \to \infty$):根据三维等球最密堆积(开普勒猜想,如面心立方 FCC 堆积),其空间最高体积分数为 $\eta_{3D} = \frac{\pi}{3\sqrt{2}} \approx 74.05\%$。结合体积公式,其渐近通项为: $$b_n \sim \frac{\pi}{3\sqrt{2}} n^3 \approx 0.7405 n^3$$

针对第 2 问:单尺寸小球填充最大比例

当挑选的单一小球半径 $R' \to 0$ 时,空间利用率逼近三维等球最密堆积极限: $$\text{最大比例} = \frac{\pi}{3\sqrt{2}} = \frac{\pi}{\sqrt{18}} \approx 74.05\%$$

针对第 3 问:多尺寸任意小球填充比例

在三维空间中,允许使用无数种不同尺寸、任意小的球体去级联填充大球内部留下的三维孔隙(三维阿波罗尼奥斯球体堆积),剩余孔隙的总体积同样可以被无限逼近于 0。

  • 填充极限比例同样可以达到 $1$(即 $100\%$)

相关笔记

基于 Obsidian 整理 · 由 VitePress 构建