pith. sign in

arxiv: 2604.12070 · v1 · submitted 2026-04-13 · 🧮 math.OC

Accelerating Column Generation in Highly Degenerate Integer Programming Problems with Template Pricing

Pith reviewed 2026-05-10 14:49 UTC · model grok-4.3

classification 🧮 math.OC
keywords column generationtemplate pricinginteger programmingdegeneracygeneralized assignment problempricing subproblemLagrangian relaxationdual smoothing
0
0 comments X

The pith

Template pricing accelerates column generation by selecting columns similar to a template vector rather than those with the absolute lowest reduced cost.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces Template pricing as a new strategy for solving pricing subproblems within column generation for integer programs. Instead of always picking the column with the best reduced cost, the method maximizes similarity to a chosen template vector among columns that still have acceptable reduced costs. This coordination across subproblems aims to reduce the slow progress caused by degeneracy in the restricted master problem. Computational tests on Generalized Assignment Problem instances show large speedups in degenerate cases, along with stronger bounds and better feasible solutions on a broad set of benchmarks.

Core claim

Template pricing replaces standard Dantzig pricing by solving a modified subproblem that maximizes column similarity to a template vector subject to a reduced-cost threshold. The subproblem is solved either exactly or via a Lagrangian relaxation heuristic. On benchmark GAP instances the approach yields convergence speedups exceeding 1000 times in highly degenerate cases and over 100 times when paired with dual smoothing, while achieving optimal column-generation bounds on every one of 1735 ISA test instances and producing stronger bounds in 43 percent and better integer solutions in 9 percent of them.

What carries the argument

Template pricing: the pricing subproblem that maximizes similarity of generated columns to a fixed template vector while restricting attention to columns whose reduced costs satisfy a suitability threshold, solved exactly or by Lagrangian relaxation.

If this is right

  • Several benchmark GAP instances reach optimality more than 1000 times faster than with Dantzig pricing.
  • The same instances reach optimality more than 100 times faster when template pricing is combined with adaptive dual smoothing.
  • Column-generation bounds become optimal on every one of the 1735 ISA instances tested.
  • Stronger bounds are obtained in 43 percent of the ISA instances and improved integer feasible solutions appear in 9 percent.
  • The method produces both faster convergence and higher-quality primal solutions without changing the master problem formulation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The template-selection step could be automated by learning representative vectors from previously solved subproblems.
  • Coordinated pricing may lessen dual oscillation enough to reduce reliance on separate stabilization techniques.
  • The Lagrangian heuristic for the template subproblem could transfer to other constrained similarity tasks inside decomposition algorithms.
  • Similar speedups might appear in other degenerate problems such as vehicle routing or cutting-stock formulations that share the same pricing structure.

Load-bearing premise

A suitable template vector exists or can be generated so that maximizing similarity among acceptable-reduced-cost columns reliably speeds overall convergence without introducing new bottlenecks or lowering bound quality.

What would settle it

Running template pricing against Dantzig pricing on the same degenerate GAP and ISA instance sets and observing no reduction in iteration count, runtime, or bound quality would show the claimed acceleration does not hold.

Figures

Figures reproduced from arXiv: 2604.12070 by Luke Marshall, Prachi Shah, Santanu S. Dey.

Figure 1
Figure 1. Figure 1: The optimal 𝛼 weight and delta change per iteration for a single machine in LT pricing for the instance c10400. 0 20 40 60 80 100 120 140 160 180 200 220 240 260 280 300 10−3 100 103 Feasible RMP # iterations log(Weight) avg min/max fit [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Average of the optimal 𝛼 weights across all machines per iteration in LT pricing for instance c10400. Min/Max values of 𝛼 across machine is shown by shaded region. The fit plotted above corresponds to: 0.014𝑒0.0234𝑥 [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Convergence of c10400 to optimality with various initialization methods (seed = 0). See Table [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Chosen age threshold policy vs. degeneracy and algorithm. Policy curves are determined by minimizing [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Yagiura: solved vs solve time [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: ISA: Solved vs solve time. Instances where LR does not reach optimal bound are treated as timed-out. [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Geometric view of Pessoa’s directional smoothing and limited [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Column management heatmap for Dantzig pricing on instance c10400. Solve time is shown as a function [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Age-threshold sweep for Dantzig pricing on instance a05200. Solve time is shown as a function of the [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
read the original abstract

We propose a new pricing strategy for column generation (CG), referred to as Template pricing. This method is motivated by the desire to coordinate solutions of different pricing subproblems in order to accelerate the convergence of the CG process and simultaneously obtain good quality integer feasible solutions. Instead of finding a column with the optimal reduced cost, Template pricing tries to maximize the similarity of columns with a given template vector, while restricting the search to columns with suitable reduced cost. We present an exact and heuristic method (based on Lagrangian relaxation) to efficiently solve the Template pricing problem. We conduct extensive computational experiments on benchmark instances of the Generalized Assignment Problem (GAP). Our results demonstrate that Template pricing can significantly accelerate the CG algorithm, especially in the presence of significant degeneracy, where several benchmark GAP instances solved over 1000x faster than Dantzig pricing, and over 100x with adaptive dual-smoothing. Template pricing allows us to achieve CG optimal bounds on all 1735 ISA instances, finding stronger bounds in 43% and improved integer solutions in 9% of these instances than previously released.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript introduces Template pricing for column generation (CG) in integer programs, particularly highly degenerate cases such as the Generalized Assignment Problem (GAP). Rather than selecting the column with minimum reduced cost (Dantzig pricing), the approach maximizes similarity to a supplied template vector among columns whose reduced cost lies below an acceptable threshold. Exact and Lagrangian-relaxation heuristic solvers are presented for the resulting pricing subproblem. Extensive experiments on GAP and ISA benchmark sets report speedups exceeding 1000× versus Dantzig pricing (and >100× with adaptive dual smoothing) on selected instances, attainment of CG-optimal bounds on all 1735 ISA instances, stronger bounds in 43% of cases, and improved integer solutions in 9%.

Significance. If the template-generation procedure can be made fully automatic, instance-independent, and cost-accounted, the method would constitute a practical advance for CG on degenerate master problems, simultaneously accelerating convergence and improving primal solution quality. The provision of both exact and Lagrangian heuristic pricing solvers, together with large-scale empirical validation on standard GAP and ISA sets, supplies concrete evidence that coordinated column selection can mitigate degeneracy effects.

major comments (3)
  1. [Abstract and §3] Abstract and §3 (Template pricing formulation): the central claims of 1000× speedups and guaranteed CG-optimal bounds on all 1735 ISA instances rest on the assumption that a suitable template vector can be supplied or generated without dominating runtime, biasing master duals, or requiring instance-specific tuning. The manuscript treats the template as 'given' and does not detail a general, zero-overhead procedure that satisfies conditions (a)–(c) in the skeptic note; this is load-bearing for the portability of the reported accelerations.
  2. [§5] §5 (Computational experiments): the reported timings and bound-quality improvements must explicitly subtract the cost of template construction and selection; if templates are chosen post-hoc per instance class or via an auxiliary solve whose time is omitted, the 1000× and 100× speedup figures and the 'all 1735 instances' optimality claim cannot be taken at face value.
  3. [§4.2] §4.2 (Lagrangian heuristic): the heuristic pricing solver is claimed to preserve the CG-optimal bound; however, no formal argument or empirical verification is supplied showing that the similarity-maximization step, even when restricted to acceptable reduced-cost columns, cannot produce a different set of columns that alters the master-problem dual trajectory or final bound quality relative to exact Dantzig pricing.
minor comments (2)
  1. [§3] Notation for the similarity measure and the 'acceptable reduced-cost threshold' should be introduced with explicit symbols and related to the standard reduced-cost definition used in Dantzig pricing.
  2. [§5] Table captions and axis labels in the experimental figures should state whether reported times include template generation and whether the same template is reused across multiple CG iterations.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment in turn below and specify the changes to be made in the revised version.

read point-by-point responses
  1. Referee: [Abstract and §3] Abstract and §3 (Template pricing formulation): the central claims of 1000× speedups and guaranteed CG-optimal bounds on all 1735 ISA instances rest on the assumption that a suitable template vector can be supplied or generated without dominating runtime, biasing master duals, or requiring instance-specific tuning. The manuscript treats the template as 'given' and does not detail a general, zero-overhead procedure that satisfies conditions (a)–(c) in the skeptic note; this is load-bearing for the portability of the reported accelerations.

    Authors: We agree that the manuscript would benefit from an explicit description of how a suitable template can be obtained. The current version treats the template vector as given in order to concentrate on the template pricing formulation itself. We will revise the abstract and Section 3 to incorporate a general, automatic template generation procedure that uses information from the master problem duals. This procedure will be instance-independent, require no per-instance tuning, and add negligible computational cost, thereby supporting the reported speedups and bound improvements without bias. revision: yes

  2. Referee: [§5] §5 (Computational experiments): the reported timings and bound-quality improvements must explicitly subtract the cost of template construction and selection; if templates are chosen post-hoc per instance class or via an auxiliary solve whose time is omitted, the 1000× and 100× speedup figures and the 'all 1735 instances' optimality claim cannot be taken at face value.

    Authors: We acknowledge that the reported timings should account for all computational costs, including template construction. In the revised manuscript, we will update Section 5 to include the time spent on template generation and selection in all timing results. We will also provide details on the template selection process used for the ISA instances to substantiate the claims regarding CG-optimal bounds on all 1735 instances. revision: yes

  3. Referee: [§4.2] §4.2 (Lagrangian heuristic): the heuristic pricing solver is claimed to preserve the CG-optimal bound; however, no formal argument or empirical verification is supplied showing that the similarity-maximization step, even when restricted to acceptable reduced-cost columns, cannot produce a different set of columns that alters the master-problem dual trajectory or final bound quality relative to exact Dantzig pricing.

    Authors: We thank the referee for this important remark. The Lagrangian heuristic is restricted to columns with acceptable reduced costs, so any column it returns is a valid improving column for the master problem. Consequently, the column generation algorithm reaches the same CG-optimal bound as with exact Dantzig pricing, although the sequence of columns and the dual trajectory may differ. We will add both a formal argument establishing the invariance of the final bound and empirical verification on a representative subset of instances to Section 4.2. revision: yes

Circularity Check

0 steps flagged

No circularity: algorithmic proposal with direct empirical validation against baselines

full rationale

The paper introduces Template pricing as a new column-generation pricing rule that maximizes similarity to a supplied template vector subject to an acceptable reduced-cost threshold, then supplies both exact and Lagrangian-relaxation heuristics for solving the resulting subproblem. All reported speed-ups (1000× on some GAP instances) and bound-quality improvements are obtained by running the full CG procedure on 1735 ISA instances and comparing wall-clock times and final bounds directly against Dantzig pricing and adaptive dual-smoothing; no equation, parameter, or uniqueness claim is fitted to the same data and then re-used as a “prediction.” The template vector is treated as an exogenous input whose generation cost is not subtracted from the reported timings, but this is an explicit modeling assumption rather than a self-referential loop. Consequently the derivation chain remains open and externally falsifiable.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard column generation theory without introducing new mathematical axioms or postulated entities. Any free parameters (e.g., weights in the similarity objective or Lagrangian multipliers) are not detailed in the abstract.

axioms (1)
  • standard math Linear programming duality and reduced-cost optimality conditions hold for the restricted master problem.
    Foundational to all column generation methods and implicitly used throughout.

pith-pipeline@v0.9.0 · 5489 in / 1235 out tokens · 56606 ms · 2026-05-10T14:49:10.510491+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

70 extracted references · 70 canonical work pages

  1. [1]

    Mathematical Programming Computation 1, 1–41 (2009)

    Achterberg, T.: Scip: solving constraint integer programs. Mathematical Programming Computation 1, 1–41 (2009)

  2. [2]

    Les Cahiers du GERAD ISSN 711, 2440 (2004)

    Amor, H.B., Desrosiers, J., Frangioni, A.: Stabilization in column generation. Les Cahiers du GERAD ISSN 711, 2440 (2004)

  3. [3]

    Discrete Applied Mathematics 157(6), 1167–1184 (2009)

    Amor, H.M.B., Desrosiers, J., Frangioni, A.: On the choice of explicit stabilizing terms in column generation. Discrete Applied Mathematics 157(6), 1167–1184 (2009)

  4. [4]

    Mathematical Programming 115(2), 351–385 (2008)

    Baldacci, R., Christofides, N., Mingozzi, A.: An exact algorithm for the vehicle routing problem based on the set partitioning formulation with additional cuts. Mathematical Programming 115(2), 351–385 (2008)

  5. [5]

    Operations research 46(3), 316–329 (1998)

    Barnhart, C., Johnson, E.L., Nemhauser, G.L., Savelsbergh, M.W., Vance, P.H.: Branch-and-price: Column generation for solving huge integer programs. Operations research 46(3), 316–329 (1998)

  6. [6]

    https://people.brunel.ac.uk/~mastjjb/jeb/orlib/gapinfo.html, accessed: 2024-12-05

    Beasely, J.E.: OR Library. https://people.brunel.ac.uk/~mastjjb/jeb/orlib/gapinfo.html, accessed: 2024-12-05

  7. [7]

    Operations Research 54(3), 454–463 (2006)

    Ben Amor, H., Desrosiers, J., Valério de Carvalho, J.M.: Dual-optimal inequalities for stabilized column generation. Operations Research 54(3), 454–463 (2006)

  8. [8]

    European Journal of Operational Research 223(2), 360–371 (2012)

    Benchimol, P., Desaulniers, G., Desrosiers, J.: Stabilized dynamic constraint aggregation for solving set partitioning problems. European Journal of Operational Research 223(2), 360–371 (2012)

  9. [9]

    Bixby, R.E., Gregory, J.W., Lustig, I., Marsten, R.E., Shanno, D.F.: Very large-scale linear programming: A case study in combining interior point and simplex methods. Oper. Res. 40, 885–897 (1992), https://api.semanticscholar.org/CorpusID: 31277854

  10. [10]

    Operations research 40(5), 885–897 (1992)

    Bixby, R.E., Gregory, J.W., Lustig, I.J., Marsten, R.E., Shanno, D.F.: Very large-scale linear programming: A case study in combining interior point and simplex methods. Operations research 40(5), 885–897 (1992)

  11. [11]

    In: International Conference on Integer Programming and Combinatorial Optimization

    Black, A.E.: Exponential lower bounds for many pivot rules for the simplex method. In: International Conference on Integer Programming and Combinatorial Optimization. pp. 86–99. Springer (2025)

  12. [12]

    European Journal of Operational Research 262(3), 835–850 (Nov 2017)

    Bouarab, H., El Hallaoui, I., Metrane, A., Soumis, F.: Dynamic constraint and variable aggregation in column generation. European Journal of Operational Research 262(3), 835–850 (Nov 2017). https://doi.org/10.1016/j.ejor.2017.04.049

  13. [13]

    Mathematical programming 113(2), 299–344 (2008)

    Briant, O., Lemaréchal, C., Meurdesoif, P., Michel, S., Perrot, N., Vanderbeck, F.: Comparison of bundle and classical column generation. Mathematical programming 113(2), 299–344 (2008)

  14. [14]

    Annals of Operations Research 86(0), 629–659 (1999)

    Valério de Carvalho, J.: Exact solution of bin-packing problems using column generation and branch-and-bound. Annals of Operations Research 86(0), 629–659 (1999)

  15. [15]

    European journal of operational research 60(3), 260–272 (1992)

    Cattrysse, D.G., Van Wassenhove, L.N.: A survey of algorithms for the generalized assignment problem. European journal of operational research 60(3), 260–272 (1992)

  16. [16]

    European Journal of Operational Research 72(1), 167–174 (Jan 1994)

    Cattrysse, D., Salomon, M., Van Wassenhove, L.N.: A set partitioning heuristic for the generalized assignment problem. European Journal of Operational Research 72(1), 167–174 (Jan 1994). https://doi.org/10.1016/0377-2217(94)90338-7

  17. [17]

    Computers & Operations Research 24(1), 17–23 (1997)

    Chu, P.C., Beasley, J.E.: A genetic algorithm for the generalised assignment problem. Computers & Operations Research 24(1), 17–23 (1997)

  18. [18]

    Computers & Operations Research 24(1), 17–23 (Jan 1997)

    Chu, P., Beasley, J.: A genetic algorithm for the generalised assignment problem. Computers & Operations Research 24(1), 17–23 (Jan 1997). https://doi.org/10.1016/S0305-0548(96)00032-9

  19. [19]

    INFORMS Journal on Computing 34(2), 1141–1156 (2022)

    Costa, L., Contardo, C., Desaulniers, G., Yarkony, J.: Stabilized column generation via the dynamic separation of aggregated rows. INFORMS Journal on Computing 34(2), 1141–1156 (2022)

  20. [20]

    In: Proceedings of the 50th Annual ACM SIGACT Symposium on Theory of Computing

    Dadush, D., Huiberts, S.: A friendly smoothed analysis of the simplex method. In: Proceedings of the 50th Annual ACM SIGACT Symposium on Theory of Computing. pp. 390–403 (2018)

  21. [21]

    Operations research 8(1), 101–111 (1960)

    Dantzig, G.B., Wolfe, P.: Decomposition principle for linear programs. Operations research 8(1), 101–111 (1960)

  22. [22]

    European Journal of Operational Research 97(2), 245–259 (1997)

    Desaulniers, G., Desrosiers, J., Dumas, Y., Marc, S., Rioux, B., Solomon, M.M., Soumis, F.: Crew pairing at air france. European Journal of Operational Research 97(2), 245–259 (1997)

  23. [23]

    Desaulniers, G., Desrosiers, J., Solomon, M.M.: Column generation, vol. 5. Springer Science & Business Media (2006)

  24. [24]

    SPRINGER INTERNATIONAL PU, S.l

    Desrosiers, J., Lübbecke, M., Desaulniers, G., Gauthier, J.B.: BRANCH-AND-PRICE. SPRINGER INTERNATIONAL PU, S.l. (2025), oCLC: 1520195899

  25. [25]

    Discrete Mathematics 194(1-3), 229–237 (1999)

    Du Merle, O., Villeneuve, D., Desrosiers, J., Hansen, P.: Stabilized column generation. Discrete Mathematics 194(1-3), 229–237 (1999)

  26. [26]

    Mathematical Programming 123, 345–370 (2010)

    Elhallaoui, I., Metrane, A., Soumis, F., Desaulniers, G.: Multi-phase dynamic constraint aggregation for set partitioning type problems. Mathematical Programming 123, 345–370 (2010)

  27. [27]

    Operations Research 53(4), 632–645 (2005)

    Elhallaoui, I., Villeneuve, D., Soumis, F., Desaulniers, G.: Dynamic aggregation of set-partitioning constraints in column generation. Operations Research 53(4), 632–645 (2005)

  28. [28]

    Management science 27(1), 1–18 (1981)

    Fisher, M.L.: The lagrangian relaxation method for solving integer programming problems. Management science 27(1), 1–18 (1981)

  29. [29]

    Mathematical programming 57(1), 341– 374 (1992)

    Forrest, J.J., Goldfarb, D.: Steepest-edge simplex algorithms for linear programming. Mathematical programming 57(1), 341– 374 (1992)

  30. [30]

    Trans- portation Science 50(1), 23–34 (2016)

    Fukasawa, R., He, Q., Song, Y.: A branch-cut-and-price algorithm for the energy minimization vehicle routing problem. Trans- portation Science 50(1), 23–34 (2016)

  31. [31]

    In: Di Gaspero, L., Festa, P., Nakib, A., Pavone, M

    Geibinger, T., Kletzander, L., Musliu, N.: Instance space analysis for the generalized assignment problem. In: Di Gaspero, L., Festa, P., Nakib, A., Pavone, M. (eds.) Metaheuristics. pp. 421–435. Springer International Publishing, Cham (2023)

  32. [32]

    In: Approaches to integer programming, pp

    Geoffrion, A.M.: Lagrangean relaxation for integer programming. In: Approaches to integer programming, pp. 82–114. Springer (2009)

  33. [33]

    Optimization Letters 3(1), 123–136 (Jan 2009)

    Ghoniem, A., Sherali, H.D.: Complementary column generation and bounding approaches for set partitioning formulations. Optimization Letters 3(1), 123–136 (Jan 2009). https://doi.org/10.1007/s11590-008-0097-2

  34. [34]

    Mathematical Programming 12, 361–371 (1977)

    Goldfarb, D., Reid, J.K.: A practicable steepest-edge simplex algorithm. Mathematical Programming 12, 361–371 (1977)

  35. [35]

    INFORMS Journal on Computing 28(1), 175–194 (2016) 18 Marshall, Shah, and Dey

    Gschwind, T., Irnich, S.: Dual inequalities for stabilized column generation revisited. INFORMS Journal on Computing 28(1), 175–194 (2016) 18 Marshall, Shah, and Dey

  36. [36]

    Trans- portation science 35(3), 286–303 (2001)

    Haase, K., Desaulniers, G., Desrosiers, J.: Simultaneous vehicle and crew scheduling in urban mass transit systems. Trans- portation science 35(3), 286–303 (2001)

  37. [37]

    Mathematical Programming 5, 1–28 (1973), https://api

    Harris, P.M.J.: Pivot selection methods of the devex lp code. Mathematical Programming 5, 1–28 (1973), https://api. semanticscholar.org/CorpusID:45065141

  38. [38]

    Mathematical Pro- gramming Computation 4(4), 363–381 (2012)

    Held, S., Cook, W., Sewell, E.C.: Maximum-weight stable sets and safe lower bounds for graph coloring. Mathematical Pro- gramming Computation 4(4), 363–381 (2012)

  39. [39]

    Technical report, Optimization Online (November 2025), https://optimization-online.org/2025/11/ the-scip-optimization-suite-10-0/

    Hojny, C., Besançon, M., Bestuzheva, K., Borst, S., Chmiela, A., Dionísio, J., Eifler, L., Ghannam, M., Gleixner, A., Göß, A., Hoen, A., van der Hulst, R., Kamp, D., Koch, T., Kofler, K., Lentz, J., Maher, S.J., Mexi, G., Mühmer, E., Pfetsch, M.E., Pokutta, S., Serrano, F., Shinano, Y., Turner, M., Vigerske, S., Walter, M., Weninger, D., Xu, L.: The SCIP ...

  40. [40]

    Huangfu and J

    Huangfu, Q., Hall, J.A.J.: Parallelizing the dual revised simplex method. Mathematical Programming Computation 10(1), 119–142 (Mar 2018). https://doi.org/10.1007/s12532-017-0130-5

  41. [41]

    In: International Conference on Integer Programming and Combinatorial Optimization

    Kukharenko, K., Sanità, L.: On the number of degenerate simplex pivots. In: International Conference on Integer Programming and Combinatorial Optimization. pp. 252–264. Springer (2024)

  42. [42]

    Transportation Research Part C: Emerging Technologies 129, 103217 (2021)

    Li, Y., Clarke, J.P., Dey, S.S.: Using submodularity within column generation to solve the flight-to-gate assignment problem. Transportation Research Part C: Emerging Technologies 129, 103217 (2021)

  43. [43]

    Liu, D.C., Nocedal, J.: On the limited memory bfgs method for large scale optimization. Math. Program. 45(1–3), 503–528 (Aug 1989)

  44. [44]

    In: Wiley Encyclopedia of Operations Research and Manage- ment Science

    Lübbecke, M.E.: Column Generation. In: Wiley Encyclopedia of Operations Research and Manage- ment Science. John Wiley & Sons, Ltd (2011). https://doi.org/10.1002/9780470400531.eorms0158, _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/9780470400531.eorms0158

  45. [45]

    Lübbecke and Jacques Desrosiers

    Lübbecke, M.E., Desrosiers, J.: Selected Topics in Column Generation. Operations Research 53(6), 1007–1023 (Dec 2005). https://doi.org/10.1287/opre.1050.0234

  46. [46]

    Marshall, L.: Template pricing (2026), https://github.com/mathgeekcoder/template-pricing

  47. [47]

    informs Journal on Computing 8(4), 344–354 (1996)

    Mehrotra, A., Trick, M.A.: A column generation approach for graph coloring. informs Journal on Computing 8(4), 344–354 (1996)

  48. [48]

    University of Melbourne, Department of Mathematics and Statistics (2000)

    Neame, P.: Nonsmooth Dual Methods in Integer Programming. University of Melbourne, Department of Mathematics and Statistics (2000)

  49. [49]

    INFOR: Information Systems and Operational Research 45(3), 123–141 (2007)

    Öncan, T.: A survey of the generalized assignment problem and its applications. INFOR: Information Systems and Operational Research 45(3), 123–141 (2007)

  50. [50]

    Mathematical Programming Computation 9(1), 61–100 (2017)

    Pecin, D., Pessoa, A., Poggi, M., Uchoa, E.: Improved branch-cut-and-price for capacitated vehicle routing. Mathematical Programming Computation 9(1), 61–100 (2017)

  51. [51]

    INFORMS Journal on Computing 30(2), 339–360 (May 2018)

    Pessoa, A., Sadykov, R., Uchoa, E., Vanderbeck, F.: Automation and Combination of Linear-Programming Based Stabilization Techniques in Column Generation. INFORMS Journal on Computing 30(2), 339–360 (May 2018). https://doi.org/10.1287/ijoc.2017.0784

  52. [52]

    In: Experimental Algorithms: 12th International Symposium, SEA 2013, Rome, Italy, June 5-7, 2013

    Pessoa, A., Sadykov, R., Uchoa, E., Vanderbeck, F.: In-out separation and column generation stabilization by dual price smoothing. In: Experimental Algorithms: 12th International Symposium, SEA 2013, Rome, Italy, June 5-7, 2013. Proceedings

  53. [53]

    pp. 354–365. Springer (2013)

  54. [54]

    INFORMS Journal on Computing 30(2), 339–360 (2018)

    Pessoa, A., Sadykov, R., Uchoa, E., Vanderbeck, F.: Automation and combination of linear-programming based stabilization techniques in column generation. INFORMS Journal on Computing 30(2), 339–360 (2018)

  55. [55]

    Electron

    Pigatti, A., De Aragao, M.P., Uchoa, E.: Stabilized branch-and-cut-and-price for the generalized assignment problem. Electron. Notes Discret. Math. 19, 389–395 (2005)

  56. [56]

    Operations Research Letters 35(5), 660–668 (2007)

    Rousseau, L.M., Gendreau, M., Feillet, D.: Interior point stabilization for column generation. Operations Research Letters 35(5), 660–668 (2007)

  57. [57]

    Operations research 45(6), 831–841 (1997)

    Savelsbergh, M.: A branch-and-price algorithm for the generalized assignment problem. Operations research 45(6), 831–841 (1997)

  58. [58]

    Communications of the ACM 52(10), 76–84 (2009)

    Spielman, D.A., Teng, S.H.: Smoothed analysis: an attempt to explain the behavior of algorithms in practice. Communications of the ACM 52(10), 76–84 (2009)

  59. [59]

    Uchoa, E., Pessoa, A., Moreno, L.: Optimizing with Column Generation: Advanced branch-cut-and-price algorithms (Part I). Tech. Rep. L-2024-3, Cadernos do LOGIS-UFF, Universidade Federal Fluminense, Engenharia de Produção (August 2024)

  60. [60]

    Springer (2010)

    Vanderbeck, F., Wolsey, L.A.: Reformulation and decomposition of integer programs. Springer (2010)

  61. [61]

    Vanderbeck, F.: Decomposition and column generation for integer programs. Ph.D. thesis, UCL-Université Catholique de Louvain (1994)

  62. [62]

    Vanderbeck, F.: Decomposition and column generation for integer programs. Ph.D. thesis, Universite’ Catholique do Lovain (1994)

  63. [63]

    In: Desaulniers, G., Desrosiers, J., Solomon, M.M

    Vanderbeck, F.: Implementing Mixed Integer Column Generation. In: Desaulniers, G., Desrosiers, J., Solomon, M.M. (eds.) Column Generation, pp. 331–358. Springer US, Boston, MA (2005)

  64. [64]

    Opera- tions Research Letters 34(3), 296–306 (May 2006)

    Vanderbeck, F., Savelsbergh, M.W.P.: A generic view of Dantzig–Wolfe decomposition in mixed integer programming. Opera- tions Research Letters 34(3), 296–306 (May 2006)

  65. [65]

    International Transactions in Op- erational Research 4(2), 151–162 (1997)

    Wentges, P.: Weighted dantzig-wolfe decomposition for linear mixed-integer programming. International Transactions in Op- erational Research 4(2), 151–162 (1997)

  66. [66]

    http://www.al.cm.is.nagoya-u.ac.jp/~yagiura/gap/, ac- cessed: 2024-12-05

    Yagiura, M.: GAP (generalized assignment problem) instances. http://www.al.cm.is.nagoya-u.ac.jp/~yagiura/gap/, ac- cessed: 2024-12-05

  67. [67]

    Yagiura, M., Ibaraki, T., Glover, F.: A path relinking approach with ejection chains for the generalized assignment problem. European journal of operational research 169(2), 548–569 (2006) 7 Statements and Declarations 7.1 Funding The authors declare that no funds, grants, or other support were received during the preparation of this manuscript. Template ...

  68. [68]

    If non-negative reduced costs, then we have 𝑆𝑖 (𝜋,𝜇)= ∅ and we stop

    Verify 𝑆𝑖 (𝜋,𝜇)≠ ∅, that is, solve the Dantzig pricing problem ( 4). If non-negative reduced costs, then we have 𝑆𝑖 (𝜋,𝜇)= ∅ and we stop

  69. [69]

    𝑥 ∈ 𝑆 𝑖 (𝜋,𝜇)≡ (3)

    Else, optimize template objective over the good columns, that is, solve ( 7): OPT ∶= max 𝑑 (𝑥, 𝑦𝑖) s.t. 𝑥 ∈ 𝑆 𝑖 (𝜋,𝜇)≡ (3)

  70. [70]

    𝑑(𝑥, 𝑦𝑖) ≥OPT 𝑥 ∈ 𝑃 𝑖 (11) Due to numerical issues, ( 11) can fail at times to find a column with negative reduced costs, even when one exists

    Finally, break ties among the optimal solutions of ( 7), that is, solve ( 6): min ∑ 𝑗∈𝐽 (𝑐𝑖𝑗 − 𝜋𝑗) 𝑥𝑗 s.t. 𝑑(𝑥, 𝑦𝑖) ≥OPT 𝑥 ∈ 𝑃 𝑖 (11) Due to numerical issues, ( 11) can fail at times to find a column with negative reduced costs, even when one exists. In these cases, we reduce the value of OPT by 1 and re-solve ( 11) until a suitable column is found. C Col...