1 Introduction

A packing by rectanglesFootnote 1 of a region \(\Omega \subset {\mathbb {R}}^2\) is a finite or countably infinite family of rectangles in \(\Omega \) with disjoint interiors. We say that the packing is perfect if the rectangles cover \(\Omega \) up to null sets. Note that this forces the Lebesgue measure \(m(\Omega )\) of \(\Omega \) to equal the sum \(\sum _{n=1}^\infty m(R_n)\) of the areas of the rectangles.

Meir and Moser [10] posed the question of whether rectangles of dimensions \(n^{-1}\times (n+1)^{-1}\) for \(n\ge 1\) can perfectly pack the unit square \([0,1]^2\), as well as the very similar question of whether squares of sidelength \(n^{-1}\) for \(n\ge 2\) can perfectly pack a rectangle of area \(\sum _{n=2}^\infty n^{-2}=\pi ^2/6-1\). These questions remain open; see for instance [3, pp. 112–113] and [1, Chap. 3] for further discussion. As one measure of partial progress towards these results, Paulhus [11] showedFootnote 2 that one could pack rectangles of dimensions \(n^{-1}\times (n+1)^{-1}\) for \(n\ge 1\) into a square of area \(1+1/(10^9+1)\), and squares of sidelength \(n^{-1}\) for \(n\ge 2\) into a rectangle of area \(\pi ^2/6-1+1/1244918662\). Very recently, it was shown in [14] that the rectangles \(n^{-1}\times (n+1)^{-1}\) for \(1\le n\le 1.35\times 10^{11}\) could be packed into the unit square.

Another direction in which partial progress has been made is to consider whether, for any \(t>1/2\), squares of sidelength \(n^{-t}\) for \(n\ge 1\) can perfectly pack a square or rectangle of area \(\sum _{n=1}^\infty n^{-2t}\) (which is finite when \(t>1/2\)). The goal is then to get t as close as possible to 1, to address the second question of Meir and Moser posed above. Recently an affirmative answer to this question was given in the range \(1/2<t\le 2/3\) by Januszewski and Zielonka [5], building upon previous work in [2, 7, 13], as well as a packing algorithm in the previously mentioned paper [11]. In this note we extend the range of t to almost reach the value \(t=1\) corresponding to the question of Meir and Moser, at the expense of removing the first few squares in the sequence:

Theorem 1.1

Let \(1/2<t<1\) and suppose that \(n_0\) is a natural number that is sufficiently large depending on t. Then squares of sidelength \(n^{-t}\) for \(n\ge n_0\) can perfectly pack a square of area \(\sum _{n=n_0}^\infty n^{-2t}\).

As a corollary, for every \(1/2<t<1\), the squares of sidelength \(n^{-t}\) for \(n\ge 1\) can perfectly pack a finite union of squares; this latter claim was first established in the range \(1/2<t\le 5/6\) in [2], and extended to the range \(1/2<t<2/3\) in [13]. We remark that since the initial release of this preprint, a (somewhat complicated) explicit expression for \(n_0\) in terms of t was given, and the dimensions \(n^{-t}\) also replaced with more general dimensions \(f(n)^{-t}\) for certain classes of function f. Also, a higher dimensional version of this result has since been established in [9].

The strategy of proof is similar to that in the previous works [2, 5, 7, 13], in which one performs a recursive algorithm to pack the first few squares \(n^{-t}\), \(n_0\le n<n_1\) into a square of the indicated area, with the remaining space being described by a union of a family \({\mathcal {R}}_{n_1}\) of rectangles which have a certain controlled size. In previous algorithms, the total perimeter of this family \({\mathcal {R}}_{n_1}\) was comparable to the total perimeter \(\sum _{n=n_0}^{n_1-1}4n^{-t}\) of the squares that one had already packed, and thus (for large \(n_1\)) also comparable to \(n_1^t\) times the total area \(\sum _{n=n_1}^\infty n^{-2t}\) of the remaining rectangles. It is this relationship between the total perimeter and total area of \({\mathcal {R}}_{n_1}\) that prevents t from getting too close to 1, as otherwise one could not eliminate the possibility that all remaining rectangles in \({\mathcal {R}}\) had width less than \(n_1^{-t}\), thus preventing one from continuing the packing. By arranging the squares in near-lattice formations, we are able (for \(n_0\) large enough) to make the total perimeter of \({\mathcal {R}}_{n_1}\) significantly smaller than the perimeter of the squares that one has already packed, and thus significantly smaller than \(n_1^t\) times the total area of \({{\mathcal {R}}}_{n_1}\); this will allow us to take t arbitrarily close to 1. Unfortunately the argument does not seem to extend to the critical case \(t=1\) (or to the supercritical cases \(t>1\)).

We remark that the same argument (with minor notational changes) would also allow one to pack rectangles of dimensions \(n^{-t}\times (n+1)^{-t}\) for \(n\ge n_0\) perfectly into a square of area \(\sum _{n=n_0}^\infty n^{-t}(n+1)^{-t}\); we leave the details of this modification to the interested reader. The quantity \(n_0\) could be calculated explicitly as a function of t, but we have not attempted to optimize this quantity. In principle, one could combine the arguments here with some initial packing of the first \(n_0\) squares, located for instance by computer search, in order to be able to replace \(n_0\) by 1 for certain values of t that are sufficiently far from 1, but we will not attempt to do so here (among other things, it would require \(n_0\) to be reduced to a magnitude suitable for computer assistance to be viable).

2 Initial Reductions

Throughout this paper we fix the parameter \(1/2<t<1\), and then introduce the exponent

$$\begin{aligned} \delta \,\,{:}{=}\,\, 1-t \end{aligned}$$

Note that because we are in the regime \(1/2<t<1\), we have \(0<\delta <1\) and

$$\begin{aligned} t+\delta t<1. \end{aligned}$$
(2.1)

In fact, these are the only two properties of \(\delta \) that we will need in the sequel. We will use this exponent \(\delta \) to define a certain technical modification of the concept of the total perimeter of a family of rectangles.

We adopt the asymptotic notation \(X=O(Y)\), \(X\ll Y\), or \(Y\gg X\) to denote the estimate \(|X|\le C_tY\) for some constant \(C_t\) that is allowed to depend only on t (or equivalently, on \(\delta \)); in particular, these constants will be independent of the parameters M or \(N_0\) that we shall shortly introduce. We write \(X \asymp Y\) for \(X\ll Y\ll X\). Next, we select two large parameters:

  • We pick a natural number M which is sufficiently large depending on \(\delta ,t\). (One can for instance take \(M\,\,{:}{=}\,\,\lfloor \delta ^{-C/\delta }\rfloor \) for a suitably large absolute constant.Footnote 3) Roughly speaking, we will pack our squares in groups of cardinality \(\asymp M^2\) at a time, arranged into approximate lattices with \(\asymp M\) squares in each row and column.

  • Finally, we pick a number \(N_0\) that is sufficiently large depending on \(M,\delta ,t\). (For instance, one can check that \(N_0\,{:}{=}\,M^{10/\delta }\) would work in the arguments below, though this choice is far from best possible.) This will be our lower bound for the parameter \(n_0\) in Theorem 1.1; in particular, \(n_0\) will be far larger than M or \(M^2\).

Given a rectangle R, we define the width w(R) to be the smaller of the two sidelengths, and the height h(R) to be the larger of the two sidelengths (with \(w(R)=h(R)\) when R is a square), thus the area m(R) is equal to . Given a finite family \({\mathcal {R}}\) of rectangles with disjoint interiors, we can thus define the total area

and unweighted total perimeter

For technical reasons we will often work instead with the weighted total perimeter

One should think of this weighted total perimeter as a slight modification of the unweighted total perimeter, in which narrower rectangles are given slightly less weight than wider rectangles. This modification is convenient for technical induction purposes; our algorithms will at one point replace a wide rectangle with several narrower rectangles, with a favorable control on the weighted total perimeter of the latter, despite having unfavorable control on the unweighted total perimeter.

In previous literature, proofs of results such as Theorem 1.1 were given by detailing a specific recursive algorithm for generating the desired packing, and then verifying that the algorithm produced a packing with all the required properties. Here we will arrange the argument slightly differentlyFootnote 4 by using induction instead of recursion, and more precisely by using a downward induction to establish the following more technical proposition, that allows us to perfectly pack any family of rectangles that has well controlled weighted total perimeter (and also obeys some other minor conditions), and which easily implies Theorem 1.1:

Proposition 2.1

(perfectly packing some families of rectangles)    Let \(n_{\max }\ge n_0\ge N_0\) and suppose that \({\mathcal {R}}\) is a finite family of rectangles with disjoint interiors, with total area

$$\begin{aligned} {{\,\textrm{area}\,}}({\mathcal {R}})=\sum _{n=n_0}^\infty \frac{1}{n^{2t}}, \end{aligned}$$
(2.2)

and obeying the weighted total perimeter bound

$$\begin{aligned} {{\,\textrm{perim}\,}}_\delta ({\mathcal {R}})\le M^{-1+\delta /2}\sum _{n=1}^{n_0-1}\frac{1}{n^{t+\delta t}} \end{aligned}$$
(2.3)

and the crude height bound

$$\begin{aligned} \sup _{R\in {\mathcal {R}}}h(R)\le 1. \end{aligned}$$
(2.4)

Then one can pack \(\bigcup _{R \in {\mathcal {R}}}R\) by squares of sidelength \(n^{-t}\) for \(n_0\le n<n_{\max }\).

Indeed, if \(n_0\ge N_0\) and we take \({\mathcal {R}}\) to consist solely of a square S of area \(\sum _{n=n_0}^\infty n^{-2t}\), then S has sidelength \(O(n_0^{1/2-t})\) (here we use the hypothesis \(t>1/2\)), and hence

$$\begin{aligned} {{\,\textrm{perim}\,}}_\delta ({\mathcal {R}})\ll n_0^{(1/2-t)(1+\delta )}. \end{aligned}$$

On the other hand, from (2.1) we have

$$\begin{aligned} \sum _{n=1}^{n_0-1}\frac{1}{n^{t+\delta t}}\gg n_0^{1-t-\delta t}=n_0^{(1-\delta )/2}n_0^{(1/2-t)(1+\delta )}. \end{aligned}$$

Since \(n_0\ge N_0\) and \(N_0\) is sufficiently large depending on \(M,\delta ,t\), we conclude that the condition (2.3) holds. Also it is clear that S has height at most 1. Applying Proposition 2.1, we conclude that we can pack S by the squares of sidelength \(n^{-t}\) for \(n_0\le n<n_{\max }\) for any \(n_{\max }\). Sending \(n_{\max }\rightarrow \infty \) and using a standard compactness argument (see e.g., [8]) we can then pack S by squares of sidelength \(n^{-t}\) for \(n\ge n_0\), which is then a perfect packing by comparison of areas. Theorem 1.1 follows.

The key step in establishing Proposition 2.1 will be to prove the following assertion.

Proposition 2.2

(efficiently packing a small rectangle of bounded eccentricity) Let \(n_0\ge N_0\), and suppose that R is a rectangle whose dimensions w(R), h(R) obey the inequalities

$$\begin{aligned} M n_0^{-t}\le w(R)\le h(R)\le 3Mn_0^{-t}. \end{aligned}$$
(2.5)

Then one can find \(n'_0\ge n_0\) with \(n'_0-n_0\asymp M^2\) and a perfect packing of R by the squares of sidelength \(n^{-t}\) for \(n_0\le n<n'_0\), together with an additional finite family \(\mathcal R\) of rectangles with disjoint interiors and widths \(O(n_0^{-t})\), obeying the unweighted total perimeter bound

$$\begin{aligned} {{\,\textrm{perim}\,}}({\mathcal {R}})\ll Mn_0^{-t}. \end{aligned}$$
(2.6)

The point here is that the unweighted total perimeter of the rectangles \({\mathcal {R}}\) is only \(O(Mn_0^{-t})\), as compared against the unweighted total perimeter of the squares of sidelength \(n^{-t}\) for \(n_0\le n<n'_0\) which is comparable to \(M^2n_0^{-t}\). This gain of \(O(M^{-1})\) is superior to the factor of \(M^{-1+\delta /2}\) which appears in (2.3), which in turn is superior to the factor \(M^{-1+\delta }\) which is what would be needed to ensure the condition (2.5) is satisfied for certain rectangles \(R_i\) that we will construct shortly.

We prove Proposition 2.2 in the next section. Assuming it for now, we conclude the proof of Proposition 2.1 and hence Theorem 1.1. We fix \(n_{\max }\) and perform a downward induction on \(n_0\); that is to say, we assume inductively that Proposition 2.1 holds for any larger choice of \(n_0\) and any family \({\mathcal {R}}\) of rectangles obeying the various hypotheses of that proposition. Proposition 2.1 is vacuously true for \(n_0=n_{\max }\) (in this case there are no rectangles to pack), so suppose that \(n_0<n_{\max }\) and that the claim has already been proven for larger values of \(n_0\). Let \({\mathcal {R}}\) obey the hypotheses of the proposition. From (2.3) and (2.1) we have

On the other hand, from (2.2) we have

From the pigeonhole principle, we conclude that there exists \(R\in {\mathcal {R}}\) with

which simplifies (using \((1-\delta /2)/(1-\delta )>1+\delta /2\)) to

$$\begin{aligned} w(R)\gg M^{1+\delta /2}n_0^{-t}. \end{aligned}$$

Since M is assumed to be sufficiently large depending on \(\delta \) (and t), this implies

$$\begin{aligned} h(R)\ge w(R)\ge 2Mn_0^{-t}. \end{aligned}$$
(2.7)

We can then partition R into a rectangle \(R_0\) of dimensions \((w(R)-Mn_0^{-t})\times h(R)\) and a rectangle \(R_*\) of dimensions \(Mn_0^{-t}\times h(R)\). By cutting off squares of sidelength \(Mn_0^{-t}\) from \(R_*\) until the height of the remaining rectangle dips below \(2Mn_0^{-t}\), we see from (2.7) that one can partition \(R_*\) into rectangles \(R_1,\dots ,R_m\) of dimensions \(Mn_0^{-t}\times h(R_i)\) with

$$\begin{aligned} Mn_0^{-t}\le h(R_i)<2Mn_0^{-t} \end{aligned}$$

for \(i=1,\dots ,m\), and

$$\begin{aligned} \sum _{i=1}^mh(R_i)=h(R). \end{aligned}$$

From (2.4) we conclude in particular the crude upper bound

$$\begin{aligned} m\le n_0^t \end{aligned}$$
(2.8)

and we have the perfect packing

$$\begin{aligned} R=R_0\cup R_*=R_0\cup R_1\cup \dots \cup R_m. \end{aligned}$$
(2.9)

Applying Proposition 2.2m times, we can then find natural numbers

$$\begin{aligned} n_0=n'_0<n'_1<\ldots <n'_m \end{aligned}$$

with

$$\begin{aligned} n'_{i+1}-n'_i\asymp M^2 \end{aligned}$$
(2.10)

for all \(0\le i\le m-1\), which by (2.8) (and the hypothesis that \(n_0\ge N_0\) is large depending on \(M,\delta ,t\)) implies in particular that

(2.11)

(say) for all \(0\le i\le m\), and a perfect packing of each \(R_i\), \(i=1,\dots ,m\), by squares of sidelength \(n^{-t}\) for \(n'_{i-1}\le n<n'_i\), together with an additional family \({\mathcal {R}}_i\) of rectangles of disjoint interiors, widths \(O(n_0^{-t})\), and with

$$\begin{aligned} {{\,\textrm{perim}\,}}({\mathcal {R}}_i)\ll Mn_0^{-t}. \end{aligned}$$
(2.12)

If we then define the new family of rectangles

$$\begin{aligned} {\mathcal {R}}'\,{:}{=}\,(\mathcal R\setminus \{R\})\cup \{R_0\}\cup \bigcup _{i=1}^m{\mathcal {R}}_i \end{aligned}$$

then we see that the rectangles in \({{\mathcal {R}}}'\) have disjoint interiors, and \(\bigcup _{R'\in {\mathcal {R}}}R'\) is perfectly packed by squares of sidelength \(n^{-t}\) for \(n_0\le n<n'_m\), together with the rectangles in \({\mathcal {R}}'\). If \(n'_m\ge n_{\max }\) then we are now done, so assume that \(n'_m<n_{\max }\). We compute (using \(w(R_0)\le w(R)\), \(h(R_0)=h(R)\), (2.12), (2.11), (2.10), and (2.3) in turn, and using the size hypotheses on M and \(n_0\))

that is to say, \({\mathcal {R}}'\) obeys the condition (2.3) (with \(n_0\) replaced by \(n'_m\)). Also, the total area of \({\mathcal {R}}'\) can be computed to be

$$\begin{aligned} {{\,\textrm{area}\,}}({\mathcal {R}}')={{\,\textrm{area}\,}}(\mathcal R)-\sum _{n=n_0}^{n'_m-1}\frac{1}{n^{2t}}=\sum _{n=n'_m}^\infty \frac{1}{n^{2t}} \end{aligned}$$

and from (2.4) we easily see that all rectangles in \(\mathcal R'\) have height at most 1. Thus by induction hypothesis (with \({\mathcal {R}}\) replaced by \({\mathcal {R}}'\)), we can pack \(\bigcup _{R'\in {\mathcal {R}}'}R'\) by squares of sidelength \(n^{-t}\) for \(n'_m\le n<n_{\max }\). This gives the desired packing of \({\mathcal {R}}\) by squares of sidelength \(n^{-t}\) for \(n_0\le n<n_{\max }\), closing the induction.

It remains to establish Proposition 2.2. This is the purpose of the next section.

Remark 2.3

The above analysis can be converted into the following algorithm for constructing the perfect packing in Theorem 1.1:

  1. (i)

    Select a sufficiently large natural number M, initialize \(n_0\) to be the quantity in Theorem 1.1, and let \({\mathcal {R}}\) consist of a single square S of area \(\sum _{n=n_0}^\infty n^{-2t}\).

  2. (ii)

    Let R be a rectangle in \({{\mathcal {R}}}\) of maximal width w(R), and perform the subdivision (2.9) of R into rectangles \(R_0,R_1,\dots ,R_m\) as indicated above. (This assumes that \(w(R)\ge 2Mn_0^{-t}\); if this is not the case, terminate with error.)

  3. (iii)

    For each \(i=1,\dots ,m\) in turn, apply Proposition 2.2 to \(R_i\) to subdivide that rectangle into squares of sidelength \(n^{-t}\) for \(n_0\le n<n'_0\), together with an additional family of rectangles \({\mathcal {R}}_i\); then replace \(n_0\) with \(n'_0\) and continue iterating in i.

  4. (iv)

    Replace the rectangle R in \({\mathcal {R}}\) by \(R_0\) together with the rectangles in \({\mathcal {R}}_1\cup \dots \cup \mathcal R_m\), then return to step (ii).

The above analysis then ensures (for \(n_0\) large enough) that this algorithm never terminates and produces a perfect packing of the original square S.

3 Efficiently Packing a Small Rectangle of Bounded Eccentricity

We now prove Proposition 2.2. Without loss of generality we may take R to be the rectangle

$$\begin{aligned} R=[0,w(R)]\times [0,h(R)]. \end{aligned}$$

From (2.5) we may find natural numbers

$$\begin{aligned} M\le M_1\le M_2<3M \end{aligned}$$

such that

We will now take \(n'_0\,{:}{=}\,n_0+M_1M_2\), then clearly \(n'_0-n_0\asymp M^2\). We index the set \(\{n:n_0\le n<n_0+M_1M_2\}\) “lexicographically” as \(\{n_{i,j}:0\le i<M_1;\,0\le j<M_2\}\), where

$$\begin{aligned} n_{i,j}\,{:}{=}\,n_0+jM_1+i. \end{aligned}$$

Our task is then to perfectly pack R by \(M_1M_2\) squares \(S_{i,j}\) of sidelength \(n_{i,j}^{-t}\) for \(0\le i<M_1\) and \(0\le j<M_2\), together with some additional finite family \({\mathcal {R}}\) of rectangles with disjoint interiors and heights \(O(n_0^{-t})\) obeying (2.6).

To motivate the construction, suppose temporarily that the squares \(S_{i,j}\) were required to have sidelength \(n_0^{-t}\) instead of \(n_{i,j}^{-t}\). Then we could simply use the lattice packing

(3.1)

for \(0\le i<M_1\), \(0\le j<M_2\), as these squares perfectly pack the rectangle

$$\begin{aligned} \left[ 0,M_1n_0^{-t}\right] \times \left[ 0,M_2n_0^{-t}\right] \end{aligned}$$

and the remaining portion of the original rectangle R can then be perfectly packed by the two rectangles

$$\begin{aligned} \left[ 0,M_1n_0^{-t}\right] \times \left[ M_2n_0^{-t},h(R)\right] \quad \text {and}\quad \left[ M_1n_0^{-t},w(R)\right] \times \left[ 0,h(R)\right] \end{aligned}$$

which have widths \(O(n_0^{-t})\) and heights \(O(Mn_0^{-t})\) (and thus perimeters \(O(Mn_0^{-t})\)), giving the claim.

In our actual problem, the squares \(S_{i,j}\) are slightly smaller, being required to have sidelength \(n_{i,j}^{-t}\) instead of \(n_0^{-t}\). If one attempts to position the bottom left corners of the \(S_{i,j}\) in the same location \((in_0^{-t},jn_0^{-t})\) as in the lattice packing (3.1), thus

$$\begin{aligned} S_{i,j}\,{:}{=}\,\bigl [in_0^{-t},in_0^{-t}+n_{i,j}^{-t}\bigr ]\times \bigl [jn_0^{-t},jn_0^{-t}+n_{i,j}^{-t}\bigr ] \end{aligned}$$

then this would still form a packing of the rectangle R, but there would now be a large number of gaps between the squares, necessitating \({\mathcal {R}}\) to consist of something like \(\asymp M^2\) rectangles of perimeter \(\asymp n_0^{-t}\) each, which would not give the desired bound (2.6). However, it is possible to close most of these gaps by sliding the squares \(S_{i,j}\) closer together, thus reducing the perimeter of \({{\mathcal {R}}}\) substantially. More precisely, our actual construction of the \(S_{i,j}\) will take the form

$$\begin{aligned} S_{i,j}\,{:}{=}\,\big [x_{i,j},x_{i,j}+n_{i,j}^{-t}\big ]\times \big [y_{i,j},y_{i,j}+n_{i,j}^{-t}\big ] \end{aligned}$$

where

$$\begin{aligned} x_{i,j}\,{:}{=}\,w(R)-\sum _{i'=i}^{M_1-1}n_{i',j}^{-t}\quad \text {and}\quad y_{i,j}\,{:}{=}\,\sum _{j'=0}^{j-1}n_{i,j'}^{-t}; \end{aligned}$$

see Fig. 1. Here we adopt the usual convention that an empty sum such as \(\sum _{i=a}^{a-1}x_i\) vanishes. Note from the mean value theorem, the triangle inequality, and the hypothesis \(n_0\ge N_0\) that

(3.2)

and

(3.3)

Thus, up to errors of \(O(M^3n_0^{-t}/N_0)\), the points \((x_{i,j},y_{i,j})\) are arranged in a lattice of spacing \(n_0^{-t}\). Note that for any \(0\le i<M_1\) and \(0\le j<M_2\) we have

$$\begin{aligned} 0\le w(R)-M_1n_0^{-t}\le x_{i,j}\le x_{i,j}+n_{i,j}^{-t}\le w(R) \end{aligned}$$

and

$$\begin{aligned} 0\le y_{i,j}\le y_{i,j}+n_{i,j}^{-t}\le M_2n_0^{-t}\le h(R) \end{aligned}$$

and so all the squares \(S_{i,j}\) are contained in R. Next, for any \(0\le i,i'<M_1\) and \(0\le j,j'<M_2\) with \((i,j)\ne (i',j')\), we argue that the squares \(S_{i,j},S_{i',j'}\) have disjoint interiors as follows.

  • If \(j'<j\) and \(i'\ge i\), then \(y_{i',j'}+n_{i',j'}^{-t}\le y_{i,j}\), and hence the interior of \(S_{i',j'}\) lies below the interior of \(S_{i,j}\), giving disjointness. By symmetry, one also has disjointness if \(j<j'\) and \(i\ge i'\).

  • If \(i'<i\) and \(j'\le j\), then \(x_{i',j'}+n_{i',j'}^{-t}\le x_{i,j}\), and hence the interior of \(S_{i',j'}\) lies to the left of the interior of \(S_{i,j}\), giving disjointness. By symmetry, one also has disjointness if \(i<i'\) and \(j\le j'\). This covers all the possible cases for \(i,j,i',j'\).

If \(0\le i<M_1-1\) and \(0\le j<M_2-1\), then (using (3.2), (3.3) as necessary, as well as the size hypotheses on \(n_0\)) we have the relations

$$\begin{aligned} x_{i+1,j}&=x_{i,j}+n_{i,j}^{-t},\\ y_{i,j}&<y_{i+1,j}+n_{i+1,j}^{-t}< y_{i,j} + n_{i,j}^{-t},\\ x_{i,j}&<x_{i,j+1}<x_{i,j}+n_{i,j}^{-t} ,\\ y_{i,j+1}&=y_{i,j}+n_{i,j}^{-t},\\ x_{i+1,j+1}&=x_{i,j+1}+n_{i,j+1}^{-t} ,\\ y_{i+1,j+1}&<y_{i,j+1}<y_{i+1,j+1} + n_{i+1,j+1}^{-t} ,\\ x_{i+1,j}&<x_{i+1,j+1}<x_{i+1,j}+n_{i+1,j}^{-t} ,\\ y_{i+1,j+1}&=y_{i+1,j}+n_{i+1,j}^{-t} \end{aligned}$$

(see Fig. 1). As a consequence, the squares \(S_{i,j},S_{i+1,j},S_{i,j+1},S_{i+1,j+1}\) surround the rectangle

$$\begin{aligned} {[}x_{i+1,j},x_{i+1,j+1}]\times [y_{i+1,j+1},y_{i,j+1}] \end{aligned}$$
(3.4)

which by (3.2), (3.3) has width and height , and hence perimeter also.

Fig. 1
figure 1

A rectangle R (with \(M_1=3\) and \(M_2=4\)), which is perfectly packed by \(M_1M_2=12\) squares \(S_{i,j}\) with \(0\le i<3\) and \(0\le j<4\) (the square \(S_{i,j}\) depicted is for \((i,j)=(1,1)\)), together with \((M_1-1)(M_2-1)=6\) small rectangles of the form (3.4) between the squares \(S_{i,j}\), \(M_2=4\) rectangles of the form (3.5) on the left side of R, \(M_1=3\) rectangles of the form (3.6) on the upper side of R, and one rectangle (3.7) on the upper left of R. This becomes a reasonably efficient packing of the rectangle R by squares once M (and hence \(M_1,M_2\)) gets large, and \(n_0\) is extremely large compared to M

From Fig. 1 we now see that the rectangle R can be packed by the squares \(S_{i,j}\) for \(0\le i<M_1\), \(0\le j<M_2\) together with the rectangles (3.4) for \(0\le i<M_1-1\), \(0\le j<M_2-1\), as well as the additional rectangles

$$\begin{aligned} {[}0,x_{0,j}]\times \bigl [y_{0,j},y_{0,j}+n_{0,j}^{-t}\bigr ] \end{aligned}$$
(3.5)

for \(0\le j<M_2\), the rectangles

$$\begin{aligned} \bigl [x_{i,M_2-1},x_{i,M_2-1}+n_{i,M_2-1}^{-t}\bigr ]\times \bigl [y_{i,M_2-1}+n_{i,M_2-1}^{-t},h(R)\bigr ] \end{aligned}$$
(3.6)

for \(0\le i<M_1\), and the rectangle

$$\begin{aligned} {[}0,x_{0,M_2-1}]\times \bigl [y_{0,M_2-1}+n_{0,M_2-1}^{-t},h(R)\bigr ]. \end{aligned}$$
(3.7)

All of these rectangles have width and height \(O(n_0^{-t})\), thanks to (3.2), (3.3), and hence perimeter \(O(n_0^{-t})\) also. Collecting these rectangles into a family \({{\mathcal {R}}}'\), we see that

which gives (2.6) since \(N_0\) is large compared with M. The claim follows.