Abstract 1 Introduction 2 Preliminaries 3 The 13/6-Approximation 4 Conclusion References

A 13/6-Approximation for Strip Packing via the Bottom-Left Algorithm

Stefan Hougardy ORCID Research Institute for Discrete Mathematics, University of Bonn, Germany
Hausdorff Center for Mathematics, University of Bonn, Germany
Bart Zondervan ORCID Faculty of Mathematics and Computer Science, University of Bremen, Germany
Abstract

In the Strip Packing problem, we are given a vertical strip of fixed width and unbounded height, along with a set of axis-parallel rectangles. The task is to place all rectangles within the strip, without overlaps, while minimizing the height of the packing. This problem is known to be NP-hard. The Bottom-Left Algorithm is a simple and widely used heuristic for Strip Packing. Given a fixed order of the rectangles, it places them one by one, always choosing the lowest feasible position in the strip and, in case of ties, the leftmost one. Baker, Coffman, and Rivest proved in 1980 that the Bottom-Left Algorithm has approximation ratio 3 if the rectangles are sorted by decreasing width [1]. For the past 45 years, no alternative ordering has been found that improves this bound. We introduce a new rectangle ordering and show that with this ordering the Bottom-Left Algorithm achieves a 13/6 approximation for the Strip Packing problem.

Keywords and phrases:
Approximation Algorithm, Strip Packing, Bottom-Left Algorithm, Rectangle Packing
Funding:
Stefan Hougardy: funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2047/1 – 390685813.
Copyright and License:
[Uncaptioned image] © Stefan Hougardy and Bart Zondervan; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Packing and covering problems
Related Version:
Full Version: https://arxiv.org/abs/2509.04654
Editors:
Meena Mahajan, Florin Manea, Annabelle McIver, and Nguyễn Kim Thắng

1 Introduction

The Strip Packing problem is a fundamental problem in combinatorial optimization, with applications in cutting-stock manufacturing, VLSI design, and scheduling. In this problem, rectangles must be packed orthogonally into a strip of fixed width, minimizing the total height used.

The Strip Packing problem is a generalization of both Bin Packing and Load Balancing (P||Cmax), which implies that the problem is NP-hard [15] and that there cannot exist a polynomial-time algorithm with absolute approximation ratio 3/2ε, unless P=NP. Additionally, the problem is known to be strongly NP-hard [9].

These hardness results provide strong motivation for the development of efficient approximation algorithms. The earliest of such algorithms for Strip Packing is the Bottom-Left Algorithm (called the BL algorithm for short), introduced in 1980 by Baker, Coffman, and Rivest [1]. This algorithm is conceptually simple yet powerful: it processes the input rectangles in a predetermined order, placing each at the lowest possible position in the strip and breaking ties by choosing the leftmost position. Despite its simplicity, the BL algorithm laid the foundation for decades of research into geometric packing heuristics and many open questions surrounding the algorithm remain. In their pioneering work, Baker et al. [1] proved that the BL algorithm achieves an absolute approximation ratio of 3, and no better bound for this algorithm has been established since. Furthermore, Chazelle [4, 17] established a quadratic runtime implementation of the BL algorithm.

Subsequent research on absolute approximation algorithms for Strip Packing has led to significant advances [5, 20, 19, 21, 11]. In 2014, Harren, Jansen, Prädel, and van Stee [10] presented an algorithm achieving the current best-known absolute approximation ratio of 5/3+ε. For a special case in which all rectangles are skewed (that is, each rectangle has either width less than a δ-fraction of the strip width or height less than a δ-fraction of the optimal packing height), Gálvez, Grandoni, Jabal Ameli, Jansen, Khan, and Rau [8] further designed an almost tight (3/2+ε)-approximation. While these results represent major theoretical progress, the proposed algorithms are intricate and appear to have limited applicability in practical settings.

In contrast, the Bottom-Left Algorithm stands out for its conceptual simplicity and its widespread use in real-world scenarios [13, 18, 6].

As the performance of the BL algorithm strongly depends on the predetermined order of the rectangles, it is tempting to think that among the n! possible orderings of n rectangles, at least one of them produces an optimal packing. Unfortunately, this is not the case: even under the best possible ordering, the absolute approximation ratio of the BL algorithm cannot surpass 4/3ε [12], a fact already foreshadowed in earlier works [1, 3]. Nevertheless, the BL algorithm is a 3-approximation when rectangles are ordered by decreasing width (using an arbitrary order for rectangles with the same width) [1].

Although several orderings of rectangles have been explored over the past 45 years [1, 7, 24], no ordering has been found that improves this 3-approximation. To illustrate this, we present in Table 1 the known approximation ratios achieved by the BL algorithm under different ordering strategies, for both the special case of Square Strip Packing (all rectangles are squares) and for general Strip Packing.

Table 1: Upper bounds on the approximation ratio of the Bottom-Left Algorithm for different orderings of squares and rectangles. Results with a star have a matching lower bound.
Ordering Square Strip Packing Strip Packing
Decreasing width 2 [1] 3 [1]
Decreasing height 2 (as above) unbounded [1]
Increasing width 3 [24] unbounded [1]
Last-Row-Full 2 [24] undefined
Arbitrary 16 [7, 24] unbounded [1]
𝒬𝒲 2 (Section 3.5.1) 1362.167 (Theorem 1)

The main contribution of this paper is to show that, for a carefully chosen ordering, which we call 𝒬𝒲, the BL algorithm attains an approximation ratio of 13/6, improving the best known bound for this classic algorithm.

Theorem 1.

The Bottom-Left Algorithm for Strip Packing has absolute approximation ratio 13/6.

This result substantially narrows the gap between the known lower bound of 4/3ε and the old upper bound of 3. Importantly, it leaves open the intriguing possibility that the BL algorithm could eventually outperform the current state-of-the-art (5/3+ε)-approximation, which serves as a primary motivation for our work.

Related Work.

Kenyon and Rémila [16] proved that there exists an asymptotic fully polynomial-time approximation scheme for Strip Packing. Their result has been sharpened in [14, 22]. Furthermore, Harren and van Stee [11] provided a polynomial-time approximation scheme for Strip Packing where rotations of rectangles by 90 degrees are permitted. This also implies a PTAS for Square Strip Packing (all rectangles are squares).

Beyond classical Strip Packing, the Bottom-Left Algorithm has also been studied in other contexts. In Multiple Strip Packing, where there are several strips of different widths, Zhuk [23] constructs a 10-approximation that combines an online algorithm for assigning the rectangles to strips, with the BL algorithm using decreasing-width order to place the rectangles within the strip. Furthermore, Fekete, Kamphans, and Schweer [7] considered the online variant of Square Strip Packing with Tetris and gravity constraints, proving a competitive ratio of 3.5. The methods developed in this paper may also be useful in these contexts, offering the potential for improved bounds or deeper understanding of the BL algorithm in its various extensions.

Additionally, the Bottom-Left Algorithm fits within the framework of fixed priority algorithms as pioneered by Borodin, Nielsen and Rackoff [2], since it irrevocably places rectangles in a predetermined order without knowledge of future rectangles. Hence granting the algorithm the power to choose the ordering of the rectangles from an online instance yields constant competitiveness.

2 Preliminaries

A Strip Packing instance =(,W) consists of a vertical strip of fixed width W and infinite height, together with a set of n closed axis-aligned rectangles. Each rectangle r has a given width wr and height hr. The maximum height of a rectangle in is denoted by hmax() or simply by hmax if the set is clear from the context. We assume that each rectangle r of the instance fits into the strip, i.e., wrW.

A packing of the rectangles into the strip is defined by specifying the lower left coordinate (xr,yr) of each rectangle r which we call the position of the rectangle r. The packing is feasible if all rectangles lie within the strip and no two rectangles overlap within their interior, i.e., for each rr the packing satisfies the two conditions:

xr0,xr+wrW,yr 0,
(xr,xr+wr)×(yr,yr+hr)(xr,xr+wr)×(yr,yr+hr) =.

Define the height of a feasible packing as max{yr+hr:r}. The goal of Strip Packing is to compute a feasible packing of minimum height for the given instance . Denote this value by hOPT() or simply by hOPT. We do not allow rotation of rectangles in this paper.

2.1 The Bottom-Left Algorithm

The Bottom-Left Algorithm is a simple and widely used heuristic for Strip Packing. The algorithm was introduced by Baker et al. [1]. Given a Strip Packing instance (,W) together with an ordering of the rectangles r1,,rn, the BL algorithm places each rectangle in the specified order at the lowest available position in the strip, choosing the leftmost such position in case of a tie. More formally, the BL algorithm places the rectangle r1 at position (0,0), which is a feasible packing of the first rectangle. Next, assuming that the BL algorithm has obtained a feasible packing of the first i1 rectangles into the strip, it chooses a position (xi,yi) that results in a feasible packing for the first i rectangles, such that (yi,xi) is lexicographically minimum among all possible choices for the position (xi,yi). Observe that this BL packing heavily depends on the ordering of the rectangles, which is part of the input rather than being computed by the algorithm. To be precise, for a Strip Packing instance =(,W) and a sorting algorithm 𝒜 of the rectangles in , denote by 𝒜 the ordered instance where the rectangles are sorted according to 𝒜. We call the packing obtained by the BL algorithm a BL packing and denote it by BL(𝒜) for an ordered instance 𝒜. The height of a BL packing on an instance 𝒜 is denoted by hBL(𝒜). If the ordering is clear from the context, then we might drop the additional notation.

An important property of the BL algorithm is that for every rectangle ri not touching the left strip boundary, there is a rectangle rj whose right face touches the left face of ri and that is before ri in the ordering, i.e., j<i (as otherwise ri could be placed further to the left). We call the rectangle rj a left supporter of ri. A left supporter prevents the rectangle ri from being placed more to the left. Analogously, every rectangle that does not touch the bottom strip boundary has a bottom supporter.

The (absolute) approximation ratio achieved by the BL algorithm on an instance ordered by 𝒜 is defined as the ratio hBL(𝒜)/hOPT().

2.2 The Horizontal Strip Partition

For our analysis of the BL algorithm we construct a partition of the strip based on the BL packing of the rectangles. This partition might be of interest independent of our results, and has already been used by Baker et al. [1], although they did not explicitly mention it. The main idea is to partition the strip into horizontal regions, where a region is the space under the bottom face of a rectangle and above the bottom faces of all rectangles placed before by the BL algorithm. The crucial property of this partition is that a region does not contain the bottom faces of any rectangle placed prior to placing the first rectangle above the region, and hence each unoccupied gap in this region has width strictly less than the width of this first rectangle that is placed above it, as otherwise the BL algorithm could have placed the rectangle at a lower position. We make this formal in Lemma 2.

Formally, let ={r1,,rn} such that the BL algorithm places ri before rj if i<j. We define the horizontal strip partition of a BL packing as the partition of the strip [0,W]×[0,hBL(,W)] into the regions

Hi={[0,W]×[max{yr1,,yri1},yri)if 1in,[0,W]×[max{yr1,,yrn},hBL]if i=n+1.

Here [x,x)= and max=0. Thus, Hi is the space below the bottom face of ri and above the bottom faces of r1,,ri1 and Hn+1 is all the space above the highest bottom face. See Figure 1 for an example.

2.3 Covering Proper Horizontal Lines

A key component in our analysis of the approximation ratio of the BL algorithm is to bound the total area occupied by rectangles in the BL packing in terms of WhOPT. Baker et al. [1] show that when the BL algorithm places the rectangles in order of decreasing width, then each region Hi is at least half occupied by rectangles for all 1in. Moreover, as the height of Hn+1 is at most hOPT, this implies that the BL algorithm is a 3-approximation. Using similar arguments, we will derive a slightly stronger result in Lemmas 2 and 3.

For this, we want to analyze the fraction of a horizontal line in the strip that is occupied by rectangles. We call a horizontal line proper if it does not intersect the top or bottom face of any rectangle. We do not need to consider lines that are not proper, because their total area is zero. A proper horizontal line alternates between parts that are occupied by rectangles and parts that are unoccupied. An unoccupied gap is a maximal connected unoccupied part on the line.

Lemma 2.

Let i{1,,n} and be a proper horizontal line in the region Hi. Then each unoccupied gap in has width strictly less than wri.

Proof.

By definition of Hi, all rectangles placed prior to ri have their bottom face below Hi. Hence, if an unoccupied gap in has size at least wri when ri is placed, then this contradicts the BL algorithm, because the gap is a lower position where we can place ri.

Lemma 3.

Let i{1,,n}. Suppose there exists a proper horizontal line in the region Hi that intersects k rectangles just before ri is placed. If each of these k rectangles has width at least wri, then any proper horizontal line in Hi is at least k/(2k+1) occupied. In particular, it is at least 1/2 occupied if an endpoint of is occupied by a rectangle prior to placing ri.

Proof.

We consider the iteration of the BL algorithm in which rectangle ri is placed. By Lemma 2, in this iteration, each occupied part of is at least as wide as each unoccupied gap, because the rectangles intersecting have width at least wri by assumption. If the line contains k1 occupied parts in this iteration then it can contain at most k+1 unoccupied parts. Thus, at least a k/(2k+1) fraction of the line is occupied. If the leftmost or rightmost part of the line is occupied by a rectangle, then there are at most as many unoccupied parts as there are occupied parts. Hence in this case at least a k/(2k)=1/2 fraction of the line is occupied. In later iterations of the BL algorithm more rectangles may be placed that intersect the line. However, this will only increase the fraction of the line that is occupied by rectangles.

It remains to show that the statement holds for any horizontal line in the region Hi. The k rectangles that intersect have their bottom face below Hi by definition of the horizontal strip partition. And the k rectangles have their top faces above Hi, as else ri can be placed lower on top of such a rectangle, since the k rectangles each have width at least wri. Hence, any proper horizontal line in Hi intersects the k rectangles that intersect , and thus each such line satisfies the statements as desired.

Observe that when rectangles are ordered by decreasing width, then every region satisfies the condition in Lemma 3. Additionally, in this case, Baker et al. [1] show that the leftmost part of a proper horizontal line is always occupied by a rectangle.

3 The 13/6-Approximation

This section is devoted to proving that the Bottom-Left Algorithm is a 13/6-approximation when placing the rectangles in the so-called 𝒬𝒲-ordering. We start in Section 3.1 by constructing this novel ordering that is based on the 𝒬𝒲-partition of the rectangles, and establish several key properties of the partition. After that, in Section 3.2, we use the horizontal strip partition to derive lower bounds on the area occupied by rectangles in the BL packing. Most regions in this partition are at least half occupied by rectangles, though a few require special attention, being in total a little less than half occupied. Section 3.3 develops tools to distinguish different types of horizontal lines, which leads to the construction of a quadratic program to bound the maximum unoccupied area in the special regions. This analysis culminates in the proof of Theorem 1 in Section 3.4, followed by final remarks on our approach in Section 3.5.

3.1 The Algorithm

Baker, Coffman, and Rivest [1] showed that the Bottom-Left Algorithm is a 3-approximation for Strip Packing when the rectangles are placed in order of decreasing width. Moreover, they show that the analysis is tight by constructing a lower bound instance. Their example consists of a so-called checkerboard followed by a very tall rectangle with small width that is placed on top of the checkerboard construction. The checkerboard is essentially half occupied by rectangles and half unoccupied, which alone would result in a lower bound of 2 on the approximation ratio. However, the final tall rectangle enhances their lower bound to 3. Therefore, we aim to place as many tall rectangles as possible first on the strip bottom. Indeed this improves the approximation guarantee of the Bottom-Left Algorithm.

With this in mind, we partition the rectangles into three sets, and use this partition to design an ordering in which the BL algorithm packs the rectangles. First of all, let be a maximal set of the tallest possible rectangles that fit next to each other on the bottom of the strip. To be more precise, go over the rectangles in order of decreasing height and add a rectangle to if the sum of widths inside remains less or equal to the strip width W. In case two rectangles have the same height, we do not care in which order they are considered. Second of all, let 𝒲 be the rectangles in with width greater than half of the strip width. Last of all, define 𝒬=(𝒲) to be the set of remaining rectangles. We call the constructed partition of rectangles the 𝒬𝒲-partition. Pseudocode for computing the partition is given in Algorithm 1.

Algorithm 1 𝒬𝒲-partition of rectangles.

Given: Strip Packing instance (,W).
Result: Partition of =𝒬𝒲.

An important property of the 𝒬𝒲-partition is that rectangles in 𝒬𝒲 are not too tall. For a number h, we denote by h all rectangles in the set of height at least h, i.e., h={fhfh}.

Lemma 4.

Let ,𝒬,𝒲 be the sets of the 𝒬𝒲-partition of a Strip Packing instance (,W) according to Algorithm 1. For every r𝒬𝒲, we have hr12hOPT.

Proof.

For r𝒬𝒲 we have wr+fhrwf>W by the definition of . Hence there are at least two rectangles in hr{r} that cannot be next to each other in the optimal packing. As each such rectangle has height at least hr, it follows that 2hr is a lower bound on the height hOPT of an optimum packing.

For a Strip Packing instance , we order the rectangles based on the 𝒬𝒲-partition as follows: first sort the rectangles from by decreasing height, next order the rectangles from 𝒬 by decreasing width, and finally take any ordering of 𝒲. We break ties arbitrarily. We denote the BL packing obtained from this 𝒬𝒲-ordering by BL(𝒬𝒲). Figure 1 illustrates an example.

Figure 1: The packing BL(𝒬𝒲) together with the horizontal strip partition. Blue rectangles are in , brown in 𝒬, and green in 𝒲.

3.2 Setup of Analysis

We dedicate this part to the analysis of the approximation factor of the Bottom-Left Algorithm following the 𝒬𝒲-ordering. For the regions of the horizontal strip partition of this packing, we derive lower bounds on the fraction that is occupied by rectangles. These bounds imply an approximation ratio of 13/6 as stated in Theorem 1.

We index the rectangles from the instance (,W) by r1,,rn such that the BL algorithm places ri before rj if i<j. Furthermore, we denote a=|| and b=|𝒬|, which implies that ={r1,,ra}, 𝒬={ra+1,,ra+b} and 𝒲={ra+b+1,,rn} (see Figure 1).

Let H1,,Hn+1 be the horizontal strip partition as defined in Section 2.2. Observe that the regions H1,,Ha are empty, because all rectangles from are placed on the bottom of the strip. Furthermore, in case that 𝒬= or the highest top face of a rectangle from 𝒬 is placed below height hmax, then we even have a 2-approximation.

Lemma 5.

If 𝒬= or max{yr+hrr𝒬}hmax, then the Bottom-Left Algorithm of 𝒬𝒲 has approximation ratio 2.

Proof.

All rectangles from are placed on the bottom of the strip and have height at most hmax. If there are no rectangles in 𝒬, or if the highest top face in 𝒬 is below hmax, then in the worst case all rectangles from 𝒲 are stacked on top of each other at height hmax. Thus it holds that hBLhmax+r𝒲hr. Now both hmax and r𝒲hr are lower bounds on the height of an optimal packing, which implies hBL2hOPT.

Thus, from now on, we may assume that 𝒬 and that there exists some rectangle r𝒬 for which it holds that yr+hr>hmax. The next result explains why the rectangles from 𝒲 are placed after rather than before 𝒬.

Lemma 6.

Every proper horizontal line in Ha+1 is at least half occupied.

Proof.

The rectangles in are placed next to each other at the bottom of the strip and are ordered from left to right by decreasing height. The first rectangle ra+1 from 𝒬 is placed according to the BL algorithm on top of this first row at height yra+1. Let f be the rightmost rectangle in whose top face is at height yra+1. Such a rectangle always exists as ra+1 must have a bottom supporter (cf. Section 2.1). As ra+1 cannot be placed at a lower position, it follows that the gap between the right side of f and the right strip boundary is strictly less than wra+1, which is less than 12W by definition of 𝒬. Therefore, each horizontal line in Ha+1 is at least half occupied.

We continue by analyzing the regions corresponding to rectangles from 𝒲.

Lemma 7.

For a+b+2in+1, every proper horizontal line in the region Hi is at least half occupied.

Proof.

Each region Hi contains at least one rectangle from 𝒲. Thus, a proper horizontal line in Hi intersects a rectangle with width at least W/2.

Next, we consider the region Ha+b+1 between the highest bottom face of a rectangle from 𝒬 and the bottom face of the first rectangle that is placed from 𝒲. Define the top rectangle of 𝒬 as the rectangle rT in 𝒬 whose top face is placed highest by the BL algorithm; and in case of ties we let rT be one with highest bottom face (in the example shown in Figure 1 the rectangle r11 is the rectangle rT). Then any proper horizontal line in Ha+b+1 intersects rT. The rectangle rT might be the only rectangle that such a line intersects, and wrT can be small, hence an arbitrarily small fraction of the line might be occupied by rectangles. Thus, the region Ha+b+1 can be very sparsely occupied. However, for a superset of Ha+b+1 the following lemma shows that we get again an at least half occupied space. Note that we have Ha+b+1[0,W]×[yrT,yrT+hrT].

Lemma 8.

The space H=[0,W]×([0,hrT][yrT,yrT+hrT]) is at least half occupied.

Proof.

It holds that yrThrT by construction of . Moreover, the total width of the rectangles in {rT}hrT exceeds the width of the strip. Thus the BL algorithm occupies inside H an area of at least WhrT by rectangles. As H has area 2WhrT it is at least half occupied.

It remains to study the regions Hi for i{a+2,,a+b}. For this, we define the left rectangle of 𝒬 to be the first rectangle rL from 𝒬 that touches the left strip boundary (in the example shown in Figure 1 the rectangle r8 is the rectangle rL). If no rectangle in 𝒬 touches the left strip boundary, then we set rL:=rT. In Section 3.3 we will develop a technique to prove that the union Ha+2HL is almost half occupied by rectangles. But first, we consider the regions HL+1,,Ha+b. Since the rectangles in 𝒬 are ordered by decreasing width we can apply Lemma 3.

Lemma 9.

For L+1ia+b, either HiH or every proper horizontal line in the region Hi is at least half occupied.

Proof.

If rL=rT, then HiH. We thus may assume that rLrT. A proper horizontal line in Hi is above the bottom face of rL. Therefore, prior to placing ri, the line only intersects rectangles from 𝒬 and these rectangles have width at least wri. Additionally, we will show that the leftmost part of is occupied by a rectangle. Then Lemma 3 implies that is at least half occupied.

This argument is based on [1]. Let r be the highest rectangle placed before ri that touches the left strip boundary. Such r exists and is part of 𝒬, because rL𝒬 touches the left strip boundary and no rectangle from is above rL. Since r is placed before ri, it follows that the bottom face of r must be placed below the proper line by definition of the horizontal strip partition. If r intersects , then the first part of is occupied by a rectangle. Otherwise, we show that there is no rectangle placed before ri that is above r and whose left face is strictly to the left of the right face of r. Namely, suppose there exists such a rectangle, then let r′′ be the leftmost such rectangle. With the same argument as before, as is a proper horizontal line in Hi, the bottom face of r′′ is also below . Furthermore, r′′ must be supported on the left, but as r′′ is the leftmost rectangle, such a supporter is the left strip boundary, a contradiction to the definition of r. This implies that just before ri is placed, the space above r is unoccupied, and wrwri by the decreasing width ordering of 𝒬, hence ri could have been placed lower on top of r by the BL algorithm. This contradiction implies that the leftmost part of the line must be occupied by a rectangle.

3.3 The Space 𝑯𝒂+𝟐𝑯𝑳

A proper horizontal line in the regions Ha+2,,HL might be less than half occupied prior to placing the first rectangle above the line, because, contrary to before, the leftmost part of the line is either occupied by rectangles from or is unoccupied. However, if the line intersects k rectangles from 𝒬 prior to placing the first rectangle above the line, then, according to Lemma 3, the line is at least a k/(2k+1)1/3-fraction occupied. The central idea of this section is to establish a bound on the number of lines intersecting a fixed number of rectangles, which will subsequently imply that the union of the regions Ha+2HL is occupied by at least the amount 12hW112hOPTW where h=yrLyra+1 is the height of the union of these regions.

We begin with two lemmas that bound the height of any rectangle in 𝒬, and consequently, the height of the space Ha+2HL. For this, define rB to be the leftmost bottom supporter of ra+1 (cf. Section 2.1), that is, of the first rectangle that is placed from 𝒬.

Lemma 10.

For each rectangle r𝒬, we have hrhrB.

Proof.

Since ra+1 is the first rectangle that is placed from 𝒬, it follows that the bottom supporter rB is a rectangle from . It follows that hra+1hrB, as otherwise ra+1 would have been part of . As 𝒬 is ordered by decreasing width, it holds for any rectangle r𝒬 that wrwra+1. Therefore, it must hold that hrhrB, as else r would belong to .

This yields an upper bound on the height at which the left rectangle rL can be placed.

Lemma 11.

It holds that yrLhmax+hrB.

Proof.

The statement is certainly true if rL=ra+1. Thus, we may assume that there is at least one rectangle in 𝒬 that is placed before rL; let r𝒬 be the leftmost such rectangle whose left face touches the right face of a rectangle from . If there are multiple such rectangles, then let r be the top one. Observe that by Lemma 10 the top face of r is at most at height hmax+hrB. There are two options. If there exists a rectangle in 𝒬 that is placed to the left of r, then the first such rectangle must be rL because the only available left supporter left of r is the left strip boundary by definition of r. Otherwise, the space left of r is not occupied by rectangles from 𝒬. Now, the first rectangle that is placed above r (if it exists) must be placed on top of r and must touch the left strip boundary; because the width of such a rectangle is less or equal to the width of r. In conclusion, either there is a first rectangle in 𝒬 touching the left strip boundary, which then is placed at height at most hmax+hrB; or there is no such rectangle, and then the bottom face of rL=rT is placed below the top face of r, implying that yrL<hmax+hrB.

Next, we study the properties of horizontal lines that intersect a fixed number of rectangles just before placing a rectangle above the line. To this end, define the type of a proper horizontal line as the set T() of rectangles from 𝒬 that intersect prior to placing the first rectangle above . If the cardinality of T() equals k, then we say that the line has order k. Moreover, denote by LM() the leftmost rectangle from T().

Lemma 12.

Let be a proper horizontal line in Ha+2HL. Then LM() has a rectangle from as left supporter.

Proof.

Suppose not, then there is a rectangle from 𝒬 that is the left supporter of LM(). Let r be such a left supporter with highest yr and let r be the first rectangle that is placed above . Then it holds that wrwr. Prior to placing r, no rectangle intersects to the left of LM(), and no rectangle is above . As is a proper line and the space to the left of LM() is larger than wr, it follows that r can be placed with its bottom face below to the left of LM(), contradicting that r is above .

For two horizontal lines and , we denote < when is below .

Lemma 13.

Let < be proper horizontal lines. Let r𝒬 be a rectangle such that rT() and let r𝒬 be a rectangle such that both rT() and rT(). Then wrwr.

Proof.

Suppose that wr<wr, then the BL algorithm places r before r. As rT() is placed before the first rectangle is placed above , it follows that r must intersect . However, this contradicts rT().

Let k be the set of all proper horizontal lines of order k inside of the space Ha+2HL. Define k={r𝒬k:r=LM()}, that is, the set of all rectangles in 𝒬 that are the leftmost rectangle of some line of order k. Furthermore, we denote by LMk the leftmost rectangle from k, and if there are multiple such rectangles, then let LMk be the highest among them.

Lemma 14.

Let k. Then either LM()=LMk, or is at least half occupied.

Proof.

Suppose that LM()LMk. Lemma 12 states that each rectangle in k has a rectangle from as left supporter. As the rectangles in are placed from left to right by decreasing height, it follows that LMk also has the highest bottom face among all the rectangles in k. Let k be a line intersecting LMk.

Let T()={ri1,,rik} and T()={ri1,,rik} ordered from left to right. There can be at most k+1 gaps on line respectively just before the first rectangle is placed above the line, namely, at most one to the left of ri1 (resp. ri1), at most one between each consecutive pair of rectangles from the type, and at most one to the right of rik (resp. rik). Denote these gaps from left to right (with possibly length 0) by g0,,gk respectively g0,,gk. We have that xri1xri1g0, because xri1xri1 and both have a rectangle from as left supporter by Lemma 12. This implies that

j=0kgj+j=1kwrijj=1k(gj+wrij).

Now, Lemma 13 implies that j=1kwrijj=1kwrij. Thus it follows that

j=0kgjj=1kgj.

Moreover, we know that each gap on has length less than the width of any rectangle adjacent to the gap just after the first rectangle is placed above . As those rectangles have width less than or equal to the rectangles intersecting , it follows that

j=0kgjj=1kgjj=1kwrijj=1kwrij.

Thus the line is at least half occupied, which proves the first statement. Figure 2 illustrates this for k=2.

Figure 2: Example showing that the lower horizontal line of order 2 is at least half occupied. Blue rectangles are in and brown ones are in 𝒬.

Corollary 15.

Every k that intersects LMk is at least k/(2k+1) occupied. All other lines of order k are at least 1/2 occupied.

Proof.

Consequences of Lemma 3 respectively Lemma 14.

Lemma 16.

Let k, k and suppose that <. Then kk. Moreover, if k<k, then the downward projection of the gap between the rightmost rectangle of that intersects and LM() is disjoint from the downward projection of the gap between the rightmost rectangle of that intersects and LM().

Proof.

Suppose that the first statement is false, then there exist k and k with < and k>k, such that for every ′′k′′ with <′′< it holds that T(′′)=T() or T(′′)=T(). This means that all rectangles in T()T() are placed on top of rectangles of T()T(). Consider the following assignment: go over the rectangles in rT()T() from left to right, and assign a rectangle rT()T() to r if r is the first rectangle in this ordering of T()T() that is above r. By definition, a rectangle from T()T() cannot be matched with two rectangles, because then it cannot have a left supporter as its width is less than the width of the rectangles from T()T() by Lemma 13. Thus each rectangle is matched with at most one other rectangle, but then as k>k, there must be a rectangle in T()T() that has no rectangle from T()T() on top of it, and therefore, the first rectangle that is placed above can now be placed lower on top of this free rectangle. Contradiction.

Next suppose that kk. Let r be the rightmost rectangle from intersecting . Then xLM()xr+wr, because is placed in order of decreasing height, is above and both have a left supporter in by Lemma 12. This immediately implies that the downward projections are disjoint, since the first gap on ends at xLM(), and the first gap on begins at xr+wr.

To measure the distance between two horizontal lines, we denote by d(,) the difference in y-coordinates.

Theorem 17.

The union of the regions Ha+2HL consists of at least 12hW112hOPTW occupied space with h=yrLyra+1.

Proof.

By Corollary 15 the space Ha+2HL is maximally unoccupied when we have for all k that |k|1, because only lines of order k intersecting the leftmost rectangle of k can be (k+1)/(2k+1) unoccupied, while other lines of order k are only at most half unoccupied.

Define αk=sup{d(,),k s.t. LMkT()T()}, that is, the largest distance between two lines of order k that both intersect the leftmost rectangle of k. If no line of order k exists, then define αk=0. From Lemma 14 it follows that all lines of order k that are less than half occupied are considered in this supremum. It holds that αkhLMk12hOPT by Lemma 4. Furthermore, the height of the space Ha+2HL is h by definition, thus k=1αk=h. Notice that for an instance consisting of n rectangles, it holds that αk=0 for any k>n, hence k=1αk equals the finite sum k=1nαk. In conclusion, we have the linear constraints αk12hOPT for every k and k=1αk=h.

Next, consider k that is less than half occupied. This line might intersect some rectangles r1,,rz from and it intersects k rectangles ri1,,rik from 𝒬 prior to placing the first rectangle above (see Figure 3). We define x=xrz+wrz if rz exists, and otherwise we define x=0. By Lemma 2 at this moment, each gap between rectangles has width strictly less than the width of the smallest rectangle from 𝒬 that intersects the line . Hence, there exists an x-coordinate x[x,xri1] such that the line is exactly half occupied to the right of x. Define the length of the remaining gap on the line by β=xx.

Figure 3: Example of a horizontal line of order 2. To the right of x the line is exactly half occupied prior to placing the first rectangle above , the remaining gap between x and x is defined as β.

Furthermore, for all lines of order k that are less than W/2 occupied, consider the largest such β-gap

βk=max{βk is less than W/2 occupied}.

Again define βk=0 if k is empty. It holds that the gap βk is upper bounded by 12k+1W, because by Lemma 3 a line of order k is at most k2k+1W occupied and each rectangle on the line is larger than any unoccupied gap on the line. Moreover, it holds that k=1βkW, because for different values of k the downward projections of the gaps defined by βk are disjoint by Lemma 16.

Next we describe the relation between the amount of unoccupied space and the values α and β. A line of order k has exactly β+12(Wx) unoccupied space (see Figure 3). By definition it holds that ββk. Furthermore, we have xj=kβj, because the rightmost point of the βj+1-gap is left of the leftmost point of the βj-gap for all j by Lemma 16. Therefore, has at most βk+12(Wj=kβj) unoccupied space. At most αk of the lines of order k are more than half unoccupied by Lemma 14. Hence the maximum amount of unoccupied space in Ha+2HL is a solution to the following quadratic program

max k=1(αkβk+12αk(Wj=kβj))subject to,
0αk12hOPTfor all k,
k=1αk=h,
0βk12k+1Wfor all k,
k=1βkW.

The objective of this quadratic program can be rewritten as

max12 k=1αkγkwhere γk=βk+Wj=k+1βj.

Let k1 and k2 be different indices such that γk1γk2γk for all k{k1}. Then setting αk1=min{12hOPT,h}, αk2=max{0,h12hOPT} and αk=0 for k{k1,k2} maximizes the objective (provided that γk1,γk20), because by Lemma 11 we have that h=yrLyra+1hmaxhOPT, as rB is the leftmost bottom supporter of ra+1. If h12hOPT, then αk2=0 and the maximum objective value is 12hγk1=23hW when k1=1, β1=13W and βk=0 for all k2. And observe that 23hW=12hW+16hW12hW+112hOPTW as h12hOPT. Otherwise, if h>12hOPT, then the objective is 14hOPTγk1+12(h12hOPT)γk2. There exists a maximum solution with βk=0 for k{k1,k2}, because βk>0 only increases the objective when αk>0, which is only the case for k{k1,k2}. Thus if k1<k2, then, as hhOPT, the objective becomes

14hOPT(βk1+Wβk2)+12(h12hOPT)(βk2+W) 14hOPTβk1+12hW

and in case k1>k2, the objective is

14hOPT(βk1+W)+12(h12hOPT)(βk2+Wβk1)14hOPTβk2+12hW.

This is maximum for min{k1,k2}=1 and βmin{k1,k2}=13W. This also gives the desired result 112hOPTW+12hW.

3.4 The Proof of the 𝟏𝟑/𝟔-Approximation Ratio

We now combine the above results to establish Theorem 1.

Theorem 1. [Restated, see original statement.]

The Bottom-Left Algorithm for Strip Packing has absolute approximation ratio 13/6.

Proof.

Consider the packing BL(𝒬𝒲). If 𝒬= or yr+hrhmax for all r𝒬, then Lemma 5 implies that the BL algorithm is a 2-approximation. Hence assume that 𝒬 and max{yr+hrr𝒬}>hmax. The space H defined in Lemma 8 is contained in Ha+1HLHa+b+1 by Lemma 10. Thus, Lemma 6, 7, 8 and 9 imply that the space Ha+1HLHa+b+1 is at least half occupied by rectangles. Furthermore, for h=yrLyra+1, Theorem 17 states that the space Ha+2HL is at least 12hW112hOPTW occupied. Thus the BL algorithm covers at least an area of

12(hBLh)W+12hW112hOPTW=12hBLW112hOPTW.

As the total area covered by rectangles divided by W is a lower bound for hOPT we get that

hOPT12hBL112hOPT.

This implies the desired approximation guarantee

hBL136hOPT.

3.5 Final Remarks on the 𝓕𝓠𝓦-Ordering

We conclude this section with remarks on the special case involving only squares, followed by lower bounds on the performance of our approach. Moreover, we briefly discuss an alternative 𝒬𝒲-ordering of rectangles, where is ordered by increasing height, rather than decreasing height.

3.5.1 The Square Case

In the special case where all rectangles are squares, the BL algorithm on 𝒬𝒲 has an approximation ratio of 2. The reason for this is that packing 𝒬 using the 𝒬𝒲-ordering is identical to packing those squares in decreasing order of size, because placing the squares in 𝒬 by decreasing width coincides with placing them by decreasing height, thus the squares that are packed at the bottom of the strip are exactly the squares from , and 𝒬 is placed by decreasing size on top of . Baker et al. [1] showed that a BL packing of squares in decreasing size has at least half of its area occupied by squares. Adding further squares of width at least W/2 on top preserves this property. Consequently, hBL(𝒬𝒲)2hOPT() for squares.

3.5.2 Lower Bounds

The checkerboard instance of [1] gives a lower bound of 2 on the approximation ratio of the BL algorithm when squares are packed in order of decreasing size. The same lower bound holds for the 𝒬𝒲-ordering as the squares of size 2iε constitute the set .

There remains a small gap between the lower bound of 2 and our upper bound of 13/6 for the approximation ratio of the BL algorithm with 𝒬𝒲-ordering. Closing this gap cannot be achieved merely by analyzing the horizontal strip partition and bounding the area occupied by rectangles in each region, as the following example demonstrates.

Figure 4: The BL packing of an instance following the 𝒬𝒲-ordering on a strip of width W=3w with rectangles ={(w,h+1),(w+1,h)}, 𝒬={(w+1,h),(w+1,1)(w+1,h)} and 𝒲=.

Consider the BL packing in 𝒬𝒲-order for an instance with strip width 3w consisting of one rectangle of size (w,h+1), three rectangles of size (w+1,h), and one of size (w+1,1) (see Figure 4). As h,w, the combined space H=H3H6 is half occupied, while H4 is only one-third occupied. The total height of H3H6 approaches hOPT, and the height of H4 approaches 12hOPT. Using the area bound this implies that

12(hBLhr3)+13hr3=12hBL16hr3hOPT.

This example shows that narrowing the gap in the approximation ratio of the Bottom-Left Algorithm with 𝒬𝒲-ordering requires new techniques that analyze the union of the regions Ha+2HL jointly with the other regions. In particular, we have

hBL2hOPT+13hr3(2+16)hOPT.

3.5.3 Ordering 𝓕 by Increasing Height

An alternative natural ordering arises from the 𝒬𝒲-partition by arranging in order of increasing height (rather than decreasing), followed by sorting 𝒬 and 𝒲 as before. The advantage of this ordering is that every horizontal line in Ha+2,,Ha+b has its leftmost part occupied by a rectangle. However, the drawback is that there may be a gap between the rightmost rectangle from and the right strip boundary, which complicates bounding the occupied area inside [0,W]×[0,hmax].

Fortunately, using Lemma 3, the regions Ha+2,,Ha+b that are above hmax are at least half occupied. Also, Lemma 7 and 8 still holds, hence the regions Ha+b+2,,Hn+1 and the space H are at least half occupied. Now under the most pessimistic assumption, the space [0,W]×[0,hmax] is completely unoccupied. This results in an upper bound of 3 on the approximation ratio for the BL algorithm under this ordering of the rectangles, because with an area argument it holds that 12(hBLhmax)hOPT.

Whenever there is no gap between and the right strip boundary, we obtain a 2-approximation when ordering by increasing height, because, following the reasoning in Lemma 6, the width of the unoccupied gap left of is at most the width of a rectangle from 𝒬, which is bounded by W/2.

4 Conclusion

In this paper, we presented a new ordering of rectangles under which the Bottom-Left Algorithm achieves an approximation ratio of 13/6, improving the previously best-known bound of 3 for the BL algorithm given by [1]. A key ingredient in our analysis is the detailed study of horizontal lines in the packing and the fraction of each line that is covered by rectangles. For this, we developed a technique based on formulating and solving a quadratic program to determine the number of lines of a given order. This method may also prove useful for refining the analysis of other packing algorithms.

Determining the exact approximation ratio of the BL algorithm remains an open and intriguing problem. There are instances for which no ordering of rectangles can achieve a ratio better than 4/3ε for the BL algorithm [12], even when all rectangles are squares. Closing the gap between 4/3ε and 13/6 is interesting, because of the possibility for the BL algorithm to surpass the current best-known (5/3+ε)-approximation for Strip Packing [10].

References

  • [1] Brenda S. Baker, Edward G. Coffman Jr., and Ronald L. Rivest. Orthogonal packings in two dimensions. SIAM J. Comput., 9(4):846–855, 1980. doi:10.1137/0209064.
  • [2] Allan Borodin, Morten N. Nielsen, and Charles Rackoff. (Incremental) priority algorithms. Algorithmica, 37(4):295–326, 2003. doi:10.1007/S00453-003-1036-3.
  • [3] Donna J. Brown. An improved BL lower bound. Information Processing Letters, 11(1):37–39, 1980. doi:10.1016/0020-0190(80)90031-9.
  • [4] Bernard Chazelle. The bottom-left bin-packing heuristic: An efficient implementation. IEEE Transactions on Computers, C-32(8):697–707, 1983. doi:10.1109/TC.1983.1676307.
  • [5] E. G. Coffman, Jr., M. R. Garey, D. S. Johnson, and R. E. Tarjan. Performance bounds for level-oriented two-dimensional packing algorithms. SIAM Journal on Computing, 9(4):808–826, 1980. doi:10.1137/0209062.
  • [6] Jie Fang, Xubing Chen, Yunqing Rao, Yili Peng, and kuan Yan. A new approach for bin packing problem using knowledge reuse and improved heuristic. Scientific Reports, 14(1):31708, 2024. doi:10.1038/s41598-024-81749-5.
  • [7] Sándor P. Fekete, Tom Kamphans, and Nils Schweer. Online square packing with gravity. Algorithmica, 68(4):1019–1044, 2014. doi:10.1007/s00453-012-9713-8.
  • [8] Waldo Gálvez, Fabrizio Grandoni, Afrouz Jabal Ameli, Klaus Jansen, Arindam Khan, and Malin Rau. A tight (3/2+ε)-approximation for skewed strip packing. Algorithmica, 85(10):3088–3109, 2023. doi:10.1007/s00453-023-01130-2.
  • [9] M. R. Garey and D. S. Johnson. “Strong” NP-completeness results: motivation, examples, and implications. Journal of the Association for Computing Machinery, 25(3):499–508, 1978. doi:10.1145/322077.322090.
  • [10] Rolf Harren, Klaus Jansen, Lars Prädel, and Rob van Stee. A (5/3+ε)-approximation for strip packing. Computational Geometry, 47(2):248–267, 2014. doi:10.1016/j.comgeo.2013.08.008.
  • [11] Rolf Harren and Rob van Stee. Improved absolute approximation ratios for two-dimensional packing problems. In Irit Dinur, Klaus Jansen, Joseph Naor, and José Rolim, editors, Approximation, Randomization, and Combinatorial Optimization. Proceedings of the 12th International Workshop, APPROX 2009 and 13th International Workshop, RANDOM 2009, LNCS 5687, pages 177–189, 2009. doi:10.1007/978-3-642-03685-9_14.
  • [12] Stefan Hougardy and Bart Zondervan. The bottom-left algorithm for the strip packing problem. In Adele Anna Rescigno and Ugo Vaccaro, editors, Combinatorial Algorithms. 35th International Workshop, IWOCA 2024, LNCS 14764, pages 433–445, 2024. doi:10.1007/978-3-031-63021-7_33.
  • [13] Stefan Jakobs. On genetic algorithms for the packing of polygons. European Journal of Operational Research, 88(1):165–181, 1996. doi:10.1016/0377-2217(94)00166-9.
  • [14] Klaus Jansen and Roberto Solis-Oba. Rectangle packing with one-dimensional resource augmentation. Discrete Optimization, 6:310–323, 2009. doi:10.1016/j.disopt.2009.04.001.
  • [15] Richard M. Karp. Reducibility among combinatorial problems. In Raymond E. Miller and James W. Thatcher, editors, Complexity of Computer Computations (Proc. Sympos., IBM Thomas J. Watson Res. Center, Yorktown Heights, N.Y., 1972), pages 85–103. Plenum Press, New York-London, 1972. doi:10.1007/978-1-4684-2001-2_9.
  • [16] Claire Kenyon and Eric Rémila. A near-optimal solution to a two-dimensional cutting stock problem. Mathematics of Operations Research, 25(4):645–656, 2000. doi:10.1287/moor.25.4.645.12118.
  • [17] Stefan Michel. Revisiting Chazelle’s implementation of the Bottom-Left heuristic: A corrected and rigorous analysis, 2025. doi:10.48550/arXiv.2511.07008.
  • [18] T. Rolich, D. Domović, and M. Golub. Bottom-left and sequence pair for solving packing problems. In 39th International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO), pages 1318–1323, 2016. doi:10.1109/MIPRO.2016.7522344.
  • [19] Ingo Schiermeyer. Reverse-fit: A 2-optimal algorithm for packing rectangles. In Second Annual European Symposium on Algorithms (ESA ’94), LNCS 855, pages 290–299, 1994. doi:10.1007/BFb0049416.
  • [20] Daniel D. K. D. B. Sleator. A 2.5 times optimal algorithm for packing in two dimensions. Information Processing Letters, 10(1):37–40, 1980. doi:10.1016/0020-0190(80)90121-0.
  • [21] A Steinberg. A strip-packing algorithm with absolute performance bound 2. SIAM Journal on Computing, 26(2):401–409, 1997. doi:10.1137/S0097539793255801.
  • [22] Maxim Sviridenko. A note on the Kenyon–Remila strip-packing algorithm. Information Processing Letters, 112:10–12, 2012. doi:10.1016/j.ipl.2011.10.003.
  • [23] S. N. Zhuk. Approximate algorithms to pack rectangles into several strips. Discrete Mathematics and Applications, 16(1):73–85, 2006. doi:10.1515/156939206776241264.
  • [24] Bart Zondervan. The Bottom-Left algorithm for the strip packing problem. Master’s thesis, Research Institute for Discrete Mathematics, University of Bonn, 2023.