Ray projection for optimizing polytopes with prohibitively many constraints in set-covering column generation

• 成果类型: Article
• 署名作者: Porumbel, Daniel
• 署名单位: heSam Universite; Conservatoire National Arts & Metiers (CNAM)
• 刊物名称: MATHEMATICAL PROGRAMMING
• ISSN/ISSBN: 0025-5610
• DOI: 10.1007/s10107-014-0840-7
• 发表日期: 2016
• 页码: 147-197
• 关键词: cutting-stock problem; arc routing problem; integer programs; algorithms; Knapsack; packing; plane; price

摘要

A recurrent task in mathematical programming consists of optimizing polytopes with prohibitively many constraints, e.g., the primal polytope in cutting-planes methods or the dual polytope in column generation (CG). We present a method to optimize such polytopes using the following intersection sub-problem: given ray , find the maximum such that is feasible. We interpret as a direction of growth from the origin ; is the intersection point between and the polytope boundary. To ease the exposition, the first part of the paper is devoted to the general context of a (primal or dual) polytope with rational data such that: (i) the intersection sub-problem can be effectively solved in practice for input rays and (ii) is feasible. By iterating over such rays , our method produces a sequence of feasible solutions that we prove to be finitely convergent. From Section 3 on, we focus on dual polytopes in CG formulations. Typically, if the CG (separation) sub-problem can be solved by dynamic programming (DP), so can be the intersection sub-problem. The proposed integer ray method (IRM) only uses integer rays, and so, the intersection sub-problem can be solved using a DP scheme based on states indexed by integer ray values. We show that under such conditions, the intersection sub-problem can be even easier than the CG sub-problem, especially when no other integer data is available to index states in DP, i.e., if the CG sub-problem input consists of fractional (or large-range) values. As such, the IRM can tackle scaled instances (with large-range weights) of capacitated problems that seem prohibitively hard for classical CG. This is confirmed by numerical experiments on various capacitated Set-Covering problems: Capacitated Arc-Routing, Cutting-Stock and other three versions of Elastic Cutting-Stock (i.e., a problem class that includes Variable Size Bin Packing).