不等圆packing问题.pdfVIP

不等圆packing 问题 问题: 对于具有 NP 难度的不等圆 packing 问题，讨论其难解性和可计算性，并为此问题 设计一个高效的求解算法。 1 不等圆 packing 问题难解性和可计算性 我们是这样描述不等圆 packing 问题：对于n 个圆（ ）是否可以两两不重叠地 放进到指定的一个大圆 R 中，且不和R 重叠，如果可以应该如何放置。 首先我们定义“状态”是 n 个圆在大圆 R 中分布的一种情况，每个状态一一对应着n 个圆 和大圆 R 的圆心位置。当找到一种状态满足n 个圆和大圆 R 互不重叠且所有 n 个圆都在大圆 R 内，则找到这个packing 问题的一个解。因为大圆内空间的连续性，n 个不等圆中每个圆圆 心在大圆 R 中的位置取值理论上是无限的，加上n 个圆的排列，要找到满足互不重叠的状态是 困难的。下面给出几个求解不等圆packing 问题的一些解法： 最大穴度算法 ：[3] 中形式化定义了人工经验方法为最大穴度算法。在不考虑n 个圆的各种 排列情况下首先将第一个圆放进与大圆 R 相切的一个位置，接着把第二个圆放到与前两个圆 相切的位置之一，然后把第三个圆也放置到与前三个圆中的两个圆相切的位置之一，就这样将 下一个圆放到前 i 个圆中任意两个圆相切的位置，如果这个位置导致了新放入的圆与其他圆有 重叠则将新圆放置在满足与其他两个圆相切的另一位置，如果发现所有位置都不满足不重叠的 条件，则回溯到之前的一个圆，并重新选取之前圆的位置。按照这样的规则满足条件则继续添 加新圆，如果不满足则回溯到上一个圆放置的位置，直到所有的圆都放入到了大圆 R 中，则找 到问题的解，或者遍历所有情况仍旧找不到解，则认为该问题没有解。n 个圆的排列复杂度为 ，将新圆放置在i 个圆和大圆中有 个位置，所以总计计算时间为 ∏ ， 即 ∏ 。 网格算法 ：在实际编程计算中，网格算法将大圆R 用间隔为 step 正交的直线进行划分，横 向和纵向直线的交点为格点，粗略估计这样的格点有 个。这样就可以通过遍历所有 n 个圆在全部格点的位置找到满足 packing 问题的解，而对于每个圆都有 个位置选择，那么 算法时间为 ，即 拟人拟物算法 ：这是一个求解不等圆packing 问题效率较高的方法，具体在 [1] 有详细介 绍。拟物算法将大圆 R 看做封闭刚性不可改变空的圆形容器，将每个小圆看做光滑弹性的小 1 球，开始将这 n 个小球强行放到容器中，我们定义小球之间因为弹性挤压产生的能量为弹性势 能（ ），接着小球由于相互间以及和容器间的作用力不断运动，每个小 球都向受挤压的方向运动，直至所有小球受力到达平衡或者不再受力。拟物方法就是模拟这样 一个过程。拟人方法是为了解决拟物方法中可能出现的一种“死锁”的情况，即 n 个小球在局 部最小解的情况中反复移动无法跳出的情况（这时的弹性势能 U 不再变化且不为 0 ），“人为 地”从容器中取出最拥挤的小球重新放入容器内进行拟物计算，这样就解决了拟物算法可能遇 到的僵持平衡的状态。整个过程直到所有的小球都不受到弹性力，弹性势能为零（ ）则认 为找到了解。 2 高效的不等圆 packing 问题算法实可能会遇到的问题 在 [1, 4] 中给出了拟人拟物