[论文解读] Integer Programming for Classifying Orthogonal Arrays

Dursun A. Bulutoglu, Kenneth J. Ryan|arXiv (Cornell University)|Jan 9, 2015

Formal Methods in Verification参考文献 29被引用 9

一句话总结

本文提出了两种新颖的整数规划（IP）公式，用于分类此前未分类的非同构正交数组 OA(N,k,2,4)，通过同构剪枝和一种新的正交设计等价性（OD-等价性）来减轻计算负担。该方法成功实现了对 OA(160,k,2,4)（k=9,10）和 OA(176,k,2,4)（k=5–10）的分类，这些情况此前因对称性和组合爆炸而难以处理。

ABSTRACT

Classifying orthogonal arrays is a well known important class of problems that asks for finding all non-isomorphic, non-negative integer solutions to a class of systems of constraints. Solved instances are scarce. We develop two new methods based on finding all non-isomorphic solutions of two novel integer linear programming formulations for classifying all non-isomorphic OA(N,k,s,t) given a set of all non-isomorphic OA(N,k-1,s,t). We also establish the concept of orthogonal design equivalence of OA(N,k,2,t) to reduce the number of integer linear programs (ILPs) whose all non-isomorphic solutions need to be enumerated by our methods. For each ILP, we determine the largest group of permutations that can be exploited with the branch-and-bound (B&B) with isomorphism pruning algorithm of Margot [Discrete Optim~4 (2007), 40-62] without losing isomorphism classes of OA(N,k,2,t). Our contributions brought the classifications of all non-isomorphic OA(160,k,2,4) for k=9,10 and OA(176,k,2,4) for k=5,6,7,8,9,10 within computational reach. These are the smallest s=2, t=4 cases for which classification results are not available in the literature.

研究动机与目标

• 解决正交数组（OAs）在 s=2、t=4 且 N=160,176 时分类结果稀缺的问题。
• 通过正交设计等价性（OD-等价性）利用对称性，降低分类OAs的计算复杂度。
• 开发高效的整数规划（IP）公式，通过带同构剪枝的分支定界法实现非同构OAs的枚举。
• 通过系统性地向已知低指标OAs添加列，扩展现有OAs分类方法。

提出的方法

• 将OAs分类问题建模为两种新颖的整数线性规划（ILP）模型，用于在现有 OA(N,k-1,2,4) 的基础上添加新列。
• 引入正交设计等价性（OD-等价性），将符号翻转和置换对称性下的等价2-符号OAs进行分组，从而减少需求解的ILP数量。
• 应用Margot提出的带同构剪枝的分支定界算法，枚举每个ILP的所有非同构解，同时保持同构类的完整性。
• 采用混合策略，结合完整枚举（Full）、ILP扩展（ILP-Ext.）以及一种混合方法，在计算效率与解的完备性之间取得平衡。
• 通过识别与同构剪枝兼容的最大置换群，实现对每个ILP公式的对称性约减。
• 使用已知的分类结果验证方法，并与MCS算法及其他先进方法进行性能比较。

实验结果

研究问题

• RQ1 新颖的ILP公式能否高效分类此前文献中未分类的非同构OA(160,k,2,4)和OA(176,k,2,4)（k=5–10）？
• RQ2 正交设计等价性（OD-等价性）在OAs分类中，能在多大程度上减少需求解的ILP数量？
• RQ3 所提出的基于ILP的方法与MCS算法及其他现有方法相比，性能如何？
• RQ4 同构剪枝与对称感知ILP公式相结合，能否扩展至此前因λ>1且t=4而难以处理的OAs实例？
• RQ5 当λ=1且k>8时，OA(N,k,2,4)在N=160和N=176时的分类在计算上是否可行？

主要发现

• 成功完成OA(160,9,2,4)和OA(160,10,2,4)的分类，分别得到15,083个和0个非同构类，采用ILP-Ext.方法。
• 对于OA(176,k,2,4)，该方法成功分类了k=5至10的所有非同构数组，其中k=9时有26个非同构类，k=10时为0个。
• OD-等价性的使用通过将对称设计分组，显著减少了需求解的ILP数量，大幅提升计算效率。
• ILP-Ext.方法在约378,241分钟内实现了OA(176,9,2,4)的实时分类，证明了其在高k情况下的可扩展性。
• 在多个案例中，该方法优于MCS算法，尤其在N和k较大时：OA(160,k,2,4)的ILP方法耗时711,298分钟，而MCS方法耗时未知（不可用）。
• 混合方法在中等k值时表现最优，例如OA(160,8,2,4)耗时1,184分钟，得到99,618个同构类。