Bin Packing Problem with Setups (BPPS)

• BPPS is a generalization of the classical bin packing problem that incorporates setup weights and costs per class, complicating capacity and cost management.
• Integer programming formulations enriched with Minimum Classes and Minimum Bins Inequalities significantly tighten LP relaxations and enhance computational performance.
• The model has broad applications in production planning, vehicle loading, and extended multi-dimensional settings, with tailored approximation algorithms and branch-and-price methods.

The Bin Packing Problem with Setups (BPPS) is a generalization of the classical Bin Packing Problem (BPP) that captures a diverse set of practical scenarios in production planning, logistics, and other domains where setup operations with associated capacity and cost penalties must be incurred whenever items of new types or classes are consolidated in a bin. The problem introduces setup weights and setup costs in addition to classical item weights and bin capacity constraints, resulting in a substantially richer combinatorial and polyhedral structure. BPPS is NP-hard and has motivated new integer programming formulations, polyhedral analyses, approximation algorithms, and computational studies (Baldacci et al., 12 Sep 2025, Baldacci et al., 31 Dec 2025, Groschke et al., 1 Sep 2025).

1. Formal Definition and Problem Structure

Given a set of $n$ items $I = \{1,2,\ldots,n\}$, each item $i$ has weight $w_i \in \mathbb{Z}_{>0}$. Items are partitioned into $m$ classes $C = \{1,2,\ldots,m\}$ with $P = \{I_1, I_2, \ldots, I_m\}$, where $I_c \subseteq I$ is the set of items in class $c$. For each class $c$, there is a setup weight $I = \{1,2,\ldots,n\}$0 and a setup cost $I = \{1,2,\ldots,n\}$1, incurred once per bin if any item of class $I = \{1,2,\ldots,n\}$2 is assigned to that bin. Bins are unlimited, identical, and have capacity $I = \{1,2,\ldots,n\}$3. Each bin also incurs a fixed cost $I = \{1,2,\ldots,n\}$4 if used.

The objective is to assign all items to bins such that, for each bin $I = \{1,2,\ldots,n\}$5, the total used capacity satisfies

$I = \{1,2,\ldots,n\}$6

where $I = \{1,2,\ldots,n\}$7 is the set of classes with at least one item in $I = \{1,2,\ldots,n\}$8; and to minimize total cost:

$I = \{1,2,\ldots,n\}$9

where $i$0 is the set of used bins. The model ensures each item is packed exactly once and bin capacity is not exceeded, with the combinatorial complication that setup weights and costs are class-dependent per bin.

Applications include production planning with family-specific setup times and costs, vehicle loading where handling equipment is class-specific, and temperature-controlled distribution with compartmental setup penalties (Baldacci et al., 12 Sep 2025).

2. Integer Programming Formulations and Polyhedral Strengthening

A natural integer linear programming (ILP) model uses, for an a priori upper bound $i$1 on the number of bins, the decision variables:

• $i$2: item $i$3 assigned to bin $i$4,
• $i$5: class $i$6 is active in bin $i$7 (some $i$8 assigned to $i$9),
• $w_i \in \mathbb{Z}_{>0}$0: bin $w_i \in \mathbb{Z}_{>0}$1 is used.

The objective is

$w_i \in \mathbb{Z}_{>0}$2

subject to:

• $w_i \in \mathbb{Z}_{>0}$3 for all $w_i \in \mathbb{Z}_{>0}$4 (item assignment),
• $w_i \in \mathbb{Z}_{>0}$5 for all $w_i \in \mathbb{Z}_{>0}$6 (capacity),
• $w_i \in \mathbb{Z}_{>0}$7 for all $w_i \in \mathbb{Z}_{>0}$8, $w_i \in \mathbb{Z}_{>0}$9, $m$0 (logical linking).

Relaxing integrality yields an LP whose lower bound (denoted $m$1) is arbitrarily weak as $m$2 for certain pathological instances, with $m$3 (Baldacci et al., 12 Sep 2025).

To address this, Minimum Classes Inequalities (MCIs) and Minimum Bins Inequality (MBI) are introduced:

• For each class $m$4, $m$5 is a lower bound on bins required for class $m$6, leading to MCI:

$m$7

• $m$8 yields the MBI:

$m$9

The addition of these constraints tightens the LP relaxation to guarantee $C = \{1,2,\ldots,m\}$0 in the worst case (matching classical BPP), and the LP with both MCIs and MBI further strengthens the bound (Baldacci et al., 12 Sep 2025).

An upper bound on the number of bins in an optimal solution is given by summing class-wise upper bounds computed by heuristics for packing $C = \{1,2,\ldots,m\}$1 into bins of effective capacity $C = \{1,2,\ldots,m\}$2, resulting in reduced ILP problem size and improved scalability.

3. Approximation Algorithms and Hardness Considerations

Straightforward adaptations of classical bin packing approximation algorithms (such as Next Fit, First Fit, and Best Fit, even with Decreasing order) fail on BPPS, exhibiting arbitrarily poor worst-case behavior due to fragmentation introduced by setup weights (Baldacci et al., 31 Dec 2025).

A two-phase heuristic TP$C = \{1,2,\ldots,m\}$3 is established, using a black-box $C = \{1,2,\ldots,m\}$4-approximation for BPP:

• Phase 1: For each class $C = \{1,2,\ldots,m\}$5, independently pack $C = \{1,2,\ldots,m\}$6 into bins of capacity $C = \{1,2,\ldots,m\}$7 using $C = \{1,2,\ldots,m\}$8; each resulting group of bins is added to the partial solution.
• Phase 2: Greedily merge any two bins if their union (including setup weights) does not exceed $C = \{1,2,\ldots,m\}$9; repeat until no merges are possible.

This method is a $P = \{I_1, I_2, \ldots, I_m\}$0-approximation: even with $P = \{I_1, I_2, \ldots, I_m\}$1 optimal ($P = \{I_1, I_2, \ldots, I_m\}$2), the approach achieves a factor-2 guarantee. Plugging in FFD or BFD ($P = \{I_1, I_2, \ldots, I_m\}$3) achieves a factor-3 approximation for BPPS. It is unknown if sub-2-approximation is possible with general merge-based schemes (Baldacci et al., 31 Dec 2025).

4. Extensions: Multi-dimensional and Multi-criteria Variants

A two-dimensional BPPS (2D-BPPS), relevant in PCB manufacturing, involves packing rectangular items with minimum spacing into bins of fixed width and height, while minimizing both the number of bins and the number of distinct bin layouts (representing setup types). The mathematical programming formulation becomes a bi-criteria mixed-integer program (MIP) over patterns $P = \{I_1, I_2, \ldots, I_m\}$4, with key variables:

• $P = \{I_1, I_2, \ldots, I_m\}$5: number of bins using layout $P = \{I_1, I_2, \ldots, I_m\}$6,
• $P = \{I_1, I_2, \ldots, I_m\}$7: binary indicator for whether layout $P = \{I_1, I_2, \ldots, I_m\}$8 (setup) is used.

Objective:

$P = \{I_1, I_2, \ldots, I_m\}$9

subject to demand coverage and linking constraints.

A branch-and-price algorithm, leveraging column generation and adapted Ryan–Foster branching, efficiently manages the vast pattern space. This framework directly controls the trade-off between bin and setup counts by varying $I_c \subseteq I$0, and shows practical scalability for medium-sized instances (Groschke et al., 1 Sep 2025). Setups are indirectly managed in the pricing subproblem and node selection, as the setup-count objective only re-enters at the master problem level.

5. Computational Insights and Practical Performance

A comprehensive benchmark of 480 BPPS instances with varied $I_c \subseteq I$1, $I_c \subseteq I$2, bin capacities, item and setup weight ranges, and cost structures demonstrates the dramatic impact of model strengthening (Baldacci et al., 12 Sep 2025):

ILP Variant	Solved Instances	Avg. Gap (Unsolved)
Base Model $I_c \subseteq I$3	162	14.4%
$I_c \subseteq I$4 + MCIs	220	6.7%
$I_c \subseteq I$5 + MCIs + MBI	253	6.0%
$I_c \subseteq I$6 + MCIs + MBI + UB	268	5.7%

The MCIs more than double the solve rate, with further gains from the MBI and instance-specific bin count upper bound; the average gap drops from 14.4% to 5.7%. The effect is even more pronounced on large $I_c \subseteq I$7; for $I_c \subseteq I$8, only the full-strength model solves any instances within time limits.

In two-dimensional settings, branch-and-price with heuristic node selection (minimizing pattern count) reduces setup diversity by factors of $I_c \subseteq I$9--$c$0 relative to depth-first, at an increased node expansion cost. For $c$1 and moderate demands, the MIP gap reaches $c$2 within $c$3 minutes, typically using $c$4--$c$5 layouts and $c$6--$c$7 bins. Larger instances with $c$8 can retain $c$9--$c$0 final gaps after two hours (Groschke et al., 1 Sep 2025).

6. Connections, Applications, and Future Directions

BPPS subsumes classical BPP and models several practical applications:

• Production lines with family-dependent changeover times and costs,
• Vehicle loading with commodity-specific handling equipment,
• Distribution logistics with compartmentalized storage.

The key research advances include polyhedral characterizations (MCIs/MBI), provably effective ILP formulations, and the first constant-factor approximation algorithms for BPPS (Baldacci et al., 12 Sep 2025, Baldacci et al., 31 Dec 2025). For two-dimensional multi-criteria objectives, dedicated branch-and-price implementations provide explicit trade-off control and operational scalability (Groschke et al., 1 Sep 2025).

Open directions include improved approximation factors for BPPS (current best is $c$1 for merge-based algorithms), deeper polyhedral analyses for higher dimensions or complex objective structures, and further computational enhancement for large-scale industrial applications.

---

References

• "The Bin Packing Problem with Setups: Formulation, Structural Properties and Computational Insights" (Baldacci et al., 12 Sep 2025)
• "A first approximation algorithm for the Bin Packing Problem with Setups" (Baldacci et al., 31 Dec 2025)
• "Column-generation for a two-dimensional multi-criteria bin-packing problem" (Groschke et al., 1 Sep 2025)