Introduction

How small can the total cross-sectional area of a finite cylinder covering of a convex body be? Here a convex body K⊂R3K \subset\mathbb{R}^{3} is compact, convex, and has nonempty interior. A cylinder is a set

C=B+Ru,C = B + \mathbb{R}u,

where uu is a unit vector and its base B⊂u⊥B \subset u^{\perp} is measurable with finite area. All areas and projections are Euclidean. Write

Amin⁡(K)=min⁡dim⁡E=2∣πEK∣,A_{\min}(K) = \min_{\dim E=2} |\pi_E K|,

where πE\pi_E is orthogonal projection onto the two-dimensional linear subspace EE and ∣⋅∣|\cdot| denotes planar Lebesgue measure. The half-area cylinder covering conjecture asks whether every finite covering K⊂⋃i=1m(Bi+Rui)K \subset\bigcup_{i=1}^{m}(B_i + \mathbb{R}u_i) satisfies

∑i=1m∣Bi∣≥12Amin⁡(K).(1)\sum_{i=1}^{m} |B_i| \ge\frac{1}{2} A_{\min}(K). \tag*{(1)}

The cylinder question belongs to the line of covering problems that begins with Tarski’s plank problem. A plank is the region between two parallel hyperplanes, and its width is their distance. Bang proved that a finite plank covering of a convex body has total width at least the body’s minimum width [2]; see also [3], Theorem 2.1. Ball proved the directionwise refinement for centrally symmetric bodies: after dividing each plank’s width by the body’s width in the same normal direction, the sum is at least one [1], arXiv v1, Theorem 1. Thus the width cost of a plank cover has both an absolute formulation and a relative formulation adapted to the individual directions.

For cylinders in three dimensions, the two-cylinder cover of a regular tetrahedron shows that replacing widths by base areas cannot preserve the constant one. Its two axes are parallel to opposite edges, and its total base area equals half the tetrahedron’s minimum projection area.

Bezdek [3] and Bezdek and Litvak [5] attribute this example and the resulting half-area question to Bang’s 1951 paper. Their formulations allow measurable bases perpendicular to the axes, as in our definition.

For a finite covering C=(Bi+Rui)i=1m\mathcal{C}=(B_i+\mathbb{R}u_i)_{i=1}^{m}, define the directionwise relative area by

R(C;K):=∑i=1m∣Bi∣∣πui⊥K∣.(2)\mathcal{R}(\mathcal{C};K):=\sum_{i=1}^{m}\frac{|B_i|}{|\pi_{u_i^\perp}K|}. \tag*{(2)}

All denominators are positive because KK has nonempty interior. Bezdek and Litvak proved R(C;K)≥1/3\mathcal{R}(\mathcal{C};K)\ge1/3 for every finite covering in three dimensions, and the stronger bound R(C;K)≥1\mathcal{R}(\mathcal{C};K)\ge1 when KK is an ellipsoid. In their proof, the general case combines the Rogers–Shephard section–projection inequality with a volume comparison. For ellipsoids, they reduce to a ball and use a density whose integral along every interior chord is the same [5]. These arguments leave a gap between the universal one-third lower bound and the proposed one-half threshold.

Bezdek and Khan record the directionwise assertion R(C;K)≥1/2\mathcal{R}(\mathcal{C};K)\ge1/2 under the name 1-Codimensional Cylinder Covering Conjecture [4]. It would imply (1.1), since

R(C;K)≤∑i∣Bi∣Amin⁡(K).(3)\mathcal{R}(\mathcal{C};K)\le\frac{\sum_i|B_i|}{A_{\min}(K)}. \tag*{(3)}

Verreault’s 2026 survey records the half-area question as unresolved [7]. Our construction starts with the tetrahedral equality example itself and replaces its two axes by finitely many nearby directions. It violates both half-area bounds while retaining compact triangular bases.

We construct a counterexample with compact triangular bases. Set H=2H=\sqrt{2} and

K={(x,y,Ht):0≤t≤1, ∣x∣≤1−t, ∣y∣≤t}.K=\{(x,y,Ht):0\le t\le1,\ |x|\le1-t,\ |y|\le t\}.

Its vertices are (±1,0,0)(\pm1,0,0) and (0,±1,H)(0,\pm1,H), so it is a regular tetrahedron of edge length 2.

Theorem 1.1. The tetrahedron (1.4) has Amin⁡(K)=2A_{\min}(K)=\sqrt{2}. For every 0<ε≤1/20000<\varepsilon\le1/2000, it admits a covering by m=2⌈2/ε2⌉m=2\lceil2/\varepsilon^2\rceil cylinders with compact triangular perpendicular bases such that

12∑i=1m∣Bi∣=12−136000ε2+O(ε4).(4)\frac{1}{\sqrt{2}}\sum_{i=1}^{m}|B_i|=\frac{1}{2}-\frac{13}{6000}\varepsilon^2+O(\varepsilon^4). \tag*{(4)}

The remainder has absolute value at most 2ε42\varepsilon^4, and the total base area is strictly less than Amin⁡(K)/2A_{\min}(K)/2.

Every covering in the theorem is finite, although its number of cylinders increases as the tilt decreases. Similarities preserve the ratio of base area to minimum projection area, so the construction applies to every regular tetrahedron.

Corollary 1.2. Every nondegenerate tetrahedron T⊂R3T\subset\mathbb{R}^3 admits a finite cylinder covering with compact triangular perpendicular bases for which

∑i=1m∣Bi∣∣πui⊥T∣<12.\sum_{i=1}^{m}\frac{|B_i|}{|\pi_{u_i^\perp}T|}<\frac{1}{2}.

Thus the 1-Codimensional Cylinder Covering Conjecture is false.

Proof. For the regular tetrahedron KK, each denominator is at least Amin⁡(K)>0A_{\min}(K)>0, so the sum is at most ∑i∣Bi∣/Amin⁡(K)<1/2\sum_i |B_i|/A_{\min}(K)<1/2 by Theorem 1.1 with ε=1/2000\varepsilon=1/2000.

The directionwise ratios are invariant under invertible affine maps [5], arXiv v1, Section 3. To see this directly, let AA be an invertible linear map and set v=Au/∣Au∣v=Au/|Au|. The linear map

L=πv⊥A∣u⊥:u⊥⟶v⊥L=\pi_{v^\perp}A|_{u^\perp}:u^\perp\longrightarrow v^\perp

is an isomorphism: Lb=0Lb=0 implies that AbAb is parallel to AuAu, hence bb is parallel to uu, and therefore b=0b=0. The transformed cylinder has perpendicular base LBLB, and πv⊥(AK)=L(πu⊥K)\pi_{v^\perp}(AK)=L(\pi_{u^\perp}K). Both areas are multiplied by the same positive factor ∣det⁡L∣|\det L|, so their ratio is unchanged. Translations also preserve it. Every nondegenerate tetrahedron is an affine image of KK; the images of its covering retain the same finite count and compact nondegenerate triangular bases.

This affine argument concerns the directionwise ratios. The conclusion of Theorem 1.1, which uses one minimum projection area for all directions, is asserted here for regular tetrahedra.

Idea of the proof. Bang’s equality example covers the parts t≤1/2t\leq1/2 and t≥1/2t\geq1/2 by cylinders parallel to the xx- and yy-axes, respectively. These axes are parallel to opposite edges of the tetrahedron, and each base is a triangle. We subdivide these triangles into narrow sectors and give each sector its own tilted axis. Tilting reduces the perpendicular area, but it can also open gaps between neighboring sectors and between the two families.

Two choices prevent these gaps. First, the axes adjacent to a shared sector side make the planes containing that side coincide. At a point (x,y,Ht)(x,y,Ht) of the tetrahedron, these planes specify a finite list of yy-coordinates beginning at −t-t and ending at tt. Hence some adjacent pair brackets yy, even if the list is not ordered. This selects a sector in the first family; exchanging the roles of xx and yy does the same for the second. Second, we increase each sector’s radial extent by a quantity of order ε2\varepsilon^2. A gap between the families could occur only near t=1/2t=1/2. There the radial heights of the two intercepts, measured from their respective tips, have opposite first-order shifts, −εxy-\varepsilon xy and +εxy+\varepsilon xy, so their sum changes only to second order. The increases can therefore be chosen small enough that their area cost is less than the quadratic gain from tilting.

Section 2 proves the projection minimum and derives the finite construction. Section 3 proves exact coverage: a first-crossing argument chooses angular sectors, and a single radial budget excludes simultaneous failure of their cutoffs. Section 4 computes the perpendicular areas. The angular mesh has size at most ε2\varepsilon^2, which keeps the total approximation error of order ε4\varepsilon^4 even as the number of cylinders grows.

The same mechanism permits more freedom than the first construction uses. Section 5 gives a common criterion for continuous radial profiles, squared-radius cutoffs and bounded Cartesian apertures, and compares midpoint sampling with endpoint-centred sectors. Its estimates remain uniform as the partition is refined. Section 6 then lets the radial margin tend to zero with the tilt. In the same tetrahedron KK, its normalized cost is 1/2−τ2/240+(25/2)τ3+O(τ4)1/2-\tau^2/240+(25/2)\tau^3+O(\tau^4): one family of constructions attains the quadratic coefficient approached as the fixed margin tends to zero. Its cylinders cover for 0<τ≤10<\tau\leq1, whereas the explicit area bound guarantees strict saving for 0<τ≤1/40000<\tau\leq1/4000. Appendix A gives complementary direct proofs for curved caps, and Appendix B collects the exact counts and formulas under other parameter and scale choices.

Geometry and the finite construction

We first fix the area scale against which the covering will be measured. We then specify every triangle and every axis in the finite covering. The projection calculation is the facet-area-vector form of Cauchy’s projection formula; see [6] for its projection-body formulation. We include the face-by-face proof for this tetrahedron.

Lemma 2.1. For the tetrahedron (1.4), Amin⁡(K)=HA_{\min}(K)=H.

Proof. The outward face area vectors (unit normal multiplied by face area) are

(0,±H,−1),(±H,0,1).(0,\pm H,-1), \qquad(\pm H,0,1).

Indeed, the supporting planes are ±Hy−z=0\pm Hy-z=0 and ±Hx+z=H\pm Hx+z=H. The displayed outward normals have length 3\sqrt{3}, which is the area of each equilateral face of side length 2. For a unit vector uu, the faces with positive scalar product against uu project to partition the shadow, apart from the projections of edges and parallel faces, which have area zero. They give the forward endpoints of the chords parallel to uu. The projected area of a face with area vector NN is ∣N⋅u∣\lvert N\cdot u\rvert, as follows by projecting its two edge vectors and taking their cross product. The sum of the four area vectors is zero. Thus the shadow area is half the sum of their absolute scalar products with uu, namely

A(u)=max⁡(H∣ux∣,∣uz∣)+max⁡(H∣uy∣,∣uz∣).A(u)=\max(H\lvert u_x\rvert,\lvert u_z\rvert)+\max(H\lvert u_y\rvert,\lvert u_z\rvert).

Put U=max⁡(∣ux∣,∣uz∣/H)U=\max(\lvert u_x\rvert,\lvert u_z\rvert/H) and V=max⁡(∣uy∣,∣uz∣/H)V=\max(\lvert u_y\rvert,\lvert u_z\rvert/H). Since H2=2H^2=2,

(U+V)2≥∣ux∣2+∣uy∣2+2∣uz∣2/H2=1.(U+V)^2\geq\lvert u_x\rvert^2+\lvert u_y\rvert^2+2\lvert u_z\rvert^2/H^2=1.

Here U2≥∣ux∣2U^2\geq\lvert u_x\rvert^2, V2≥∣uy∣2V^2\geq\lvert u_y\rvert^2 and UV≥∣uz∣2/H2UV\geq\lvert u_z\rvert^2/H^2 account for the three terms separately. Hence A(u)=H(U+V)≥HA(u)=H(U+V)\geq H, with equality at u=(1,0,0)u=(1,0,0).

At zero tilt, the part t≤1/2t\leq1/2 is covered by an xx-parallel cylinder and the part t≥1/2t\geq1/2 by a yy-parallel cylinder. Each perpendicular base is a triangle of area H/4H/4. The lower intercept triangle can be parametrized as

(0,qt0,Ht0),−1≤q≤1,0≤t0≤1/2.(0,qt_0,Ht_0),\qquad-1\leq q\leq1,\qquad0\leq t_0\leq1/2.

Here qq is a slope parameter: for t0>0t_0>0 it is the second coordinate divided by t0t_0. We call the intervals of qq angular sectors. The radial height t0t_0 is the third physical coordinate divided by HH, rather than Euclidean distance from the tip. The upper intercept triangle is similarly (ps0,0,H(1−s0))(ps_0,0,H(1-s_0)), with −1≤p≤1-1\leq p\leq1 and 0≤s0≤1/20\leq s_0\leq1/2. We now perturb a finite subdivision of these triangles.

Fix a positive tilt parameter ε\varepsilon, and put

η=11000,n=⌈2ε2⌉,Δ=2n,qj=−1+jΔ(0≤j≤n).(5)\eta=\frac{1}{1000},\qquad n=\left\lceil\frac{2}{\varepsilon^2}\right\rceil,\qquad\Delta=\frac{2}{n},\qquad q_j=-1+j\Delta\quad(0\leq j\leq n). \tag*{(5)}

In particular Δ≤ε2\Delta\leq\varepsilon^2. To make adjacent tilted sectors meet, prescribe the displacement of the boundary with angular parameter qq by

ϕ(q)=1−q24(−1≤q≤1).\phi(q)=\frac{1-q^2}{4}\qquad(-1\leq q\leq1).

The endpoint values ϕ(−1)=ϕ(1)=0\phi(-1) = \phi(1) = 0 will keep the outer angular sides in the fixed planes y=−ty = -t and y=ty = t. For real α,β\alpha, \beta, the lines with direction (1,εα,Hεβ)(1, \varepsilon\alpha, H\varepsilon\beta) through the points (0,qt0,Ht0)(0, qt_0, Ht_0), t0≥0t_0 \ge0, lie in the plane

y=qt+εx(α−qβ).y = qt + \varepsilon x(\alpha- q\beta).

Thus neighboring angular sides will lie in the same plane if both axes use the value ϕ(q)\phi(q) there. On the interval [qj,qj+1][q_j,q_{j+1}], solve the two endpoint equations αj−qβj=ϕ(q)\alpha_j - q\beta_j = \phi(q). This gives

αj=1+qjqj+14,βj=qj+qj+14(0≤j<n).(6)\alpha_j = \frac{1 + q_jq_{j+1}}{4}, \qquad\beta_j = \frac{q_j + q_{j+1}}{4} \qquad(0 \le j < n). \tag*{(6)}

and hence the exact matching identities

αj−qβj=1−q24(q=qj or q=qj+1).(7)\alpha_j - q\beta_j = \frac{1-q^2}{4} \qquad(q = q_j \text{ or } q = q_{j+1}). \tag*{(7)}

The choice of parabola also anticipates the radial estimate. At zero tilt, the ray of slope qq meets the middle section at y=q/2y = q/2. The lower intercept height t0=t−εxβjt_0 = t - \varepsilon x\beta_j should have first-order shift −εxy-\varepsilon xy. This suggests βj≃q/2\beta_j \simeq q/2. The endpoint equations make βj\beta_j the negative secant slope of ϕ\phi, so their limiting relation is ϕ′(q)=−q/2\phi'(q) = -q/2. Together with the two endpoint conditions, this gives ϕ(q)=(1−q2)/4\phi(q) = (1-q^2)/4. The finite identities above, rather than this motivation, will be used in the proof.

We also enlarge each sector slightly beyond its original radial cutoff 1/21/2. Define

d(q)=q2(1+q2)16,Mj=η+max⁡qj≤q≤qj+1d(q),Tj=12+ε2Mj.(8)d(q) = \frac{q^2(1+q^2)}{16}, \qquad M_j = \eta+ \max_{q_j \le q \le q_{j+1}} d(q), \qquad T_j = \frac{1}{2} + \varepsilon^2 M_j. \tag*{(8)}

The function dd bounds the quadratic loss in the overlap of the two families; its role is made explicit in the coverage argument below. For later estimates, note that

0≤αj≤12,∣βj∣≤12,η≤Mj≤18+η.(9)0 \le\alpha_j \le\frac{1}{2}, \qquad|\beta_j| \le\frac{1}{2}, \qquad\eta\le M_j \le\frac{1}{8} + \eta. \tag*{(9)}

Our intercept triangles and direction vectors are

Pj={(0,y0,Ht0):0≤t0≤Tj, qjt0≤y0≤qj+1t0},vj=(1,εαj,Hεβj),(10)P_j = \{(0,y_0,Ht_0) : 0 \le t_0 \le T_j,\ q_jt_0 \le y_0 \le q_{j+1}t_0\}, \qquad v_j = (1,\varepsilon\alpha_j,H\varepsilon\beta_j), \tag*{(10)}
Qj={(x0,0,H(1−s0)):0≤s0≤Tj, qjs0≤x0≤qj+1s0},wj=(−εαj,1,Hεβj).(11)Q_j = \{(x_0,0,H(1-s_0)) : 0 \le s_0 \le T_j,\ q_js_0 \le x_0 \le q_{j+1}s_0\}, \qquad w_j = (-\varepsilon\alpha_j,1,H\varepsilon\beta_j). \tag*{(11)}

Use the 2n2n cylinders

Pj+Rvj,Qj+Rwj(0≤j<n).P_j + \mathbb{R}v_j,\qquad Q_j + \mathbb{R}w_j \qquad(0 \le j < n).

To express them with perpendicular bases, project each intercept triangle onto the linear plane perpendicular to its direction and normalize that direction. Projection along the axis does not change the cylinder: P+Rv=πv⊥P+RvP + \mathbb{R}v = \pi_{v^\perp}P + \mathbb{R}v. Moreover, it is injective on either intercept plane because the corresponding normal component of vjv_j or wjw_j is 11. Thus each actual base is a compact nondegenerate triangle, with finite area.

Figure 1 shows the original covering and the angular matching after the tilt. In a fixed plane x=Xx = X, the neighboring sectors start at different translated vertices Vi=(εXαi,εXβi)V_i = (\varepsilon X\alpha_i,\varepsilon X\beta_i) in (y,t)(y,t) coordinates, but their sides at a shared value of qq lie on the same line. Thus the matching concerns the supporting planes of the sides; the radial extents still have to be checked.

Illustrations of the zero-tilt geometry and matching after the tilt

Figure 1. (a) At zero tilt, cylinders parallel to the two opposite edges cover the lower and upper halves of KK. (b) In the plane x=Xx=X, neighboring angular sectors have distinct translated vertices VjV_j and Vj+1V_{j+1}, but their sides at q=qj+1q=q_{j+1} lie on the same line y=qt+ϵX(1−q2)/4y=qt+\epsilon X(1-q^2)/4. The dashed line is y=qty=qt, its zero-tilt position. Tilt and sector widths are exaggerated. Radial cutoffs are omitted; the colored regions end at the boundary of the displayed window.

Coverage of the entire tetrahedron

The shared-boundary identity gives an angular sector in each family for every point of KK. The remaining task is to prove that at least one of the two selected triangles has sufficient radial extent.

Proposition 3.1. For every 0<ϵ≤η/20<\epsilon\le\eta/2, the 2n2n cylinders (2.6)–(2.7), with parameters (2.1)–(2.4), cover the closed tetrahedron KK.

Proof. Choosing angular sectors. Fix (x,y,Ht)∈K(x,y,Ht)\in K and write s=1−ts=1-t. The intercept along vjv_j is (0,y0,Ht0)(0,y_0,Ht_0), where

t0=t−ϵxβj,y0=y−ϵxαj.(12)t_0=t-\epsilon x\beta_j,\qquad y_0=y-\epsilon x\alpha_j. \tag*{(12)}

Consider the finite sequence

Fi=qit+ϵx1−qi24(0≤i≤n).F_i=q_i t+\epsilon x\frac{1-q_i^2}{4}\qquad(0\le i\le n).

Its endpoints are F0=−t≤y≤t=FnF_0=-t\le y\le t=F_n. Choose the first index i≥1i\ge1 for which Fi≥yF_i\ge y, and set j=i−1j=i-1. Such an index exists because Fn≥yF_n\ge y. If i=1i=1, then Fj=F0≤yF_j=F_0\le y; if i>1i>1, minimality gives Fj<yF_j<y. Thus Fj≤y≤Fj+1F_j\le y\le F_{j+1}, including when yy equals an extreme endpoint. Using (2.3), these inequalities are exactly

qjt0≤y0≤qj+1t0.q_jt_0\le y_0\le q_{j+1}t_0.

Subtracting their outer terms gives 0≤(qj+1−qj)t0=Δt00\le(q_{j+1}-q_j)t_0=\Delta t_0. Since Δ>0\Delta>0, they force t0≥0t_0\ge0. Hence the selected first-family cylinder covers the point unless t0>Tjt_0>T_j. This argument does not assume that the sequence FiF_i is monotone.

For the second family, use the sequence

Gi=qis−εy1−qi24,G0=−s≤x≤s=Gn.G_i=q_i s-\varepsilon y\frac{1-q_i^2}{4},\qquad G_0=-s\le x\le s=G_n.

The same first-crossing choice gives an index kk for which

s0=s+εyβk≥0,x0=x+εyαk,qks0≤x0≤qk+1s0.(13)s_0=s+\varepsilon y\beta_k\ge0,\qquad x_0=x+\varepsilon y\alpha_k,\qquad q_k s_0\le x_0\le q_{k+1}s_0. \tag*{(13)}

This cylinder covers unless s0>Tks_0>T_k.

Comparing the radial cutoffs. It remains to show that the two selected radial cutoffs cannot both fail. Suppose both selected cylinders miss the point. Then t0,s0>1/2t_0,s_0>1/2, and we can set

q=y0/t0∈[qj,qj+1],p=x0/s0∈[qk,qk+1].q=y_0/t_0\in[q_j,q_{j+1}],\qquad p=x_0/s_0\in[q_k,q_{k+1}].

Since ∣x∣+∣y∣≤1|x|+|y|\le1, (2.5) gives

t0+s0=1−εxβj+εyβk≤1+ε/2.t_0+s_0=1-\varepsilon x\beta_j+\varepsilon y\beta_k\le1+\varepsilon/2.

Thus each of t0,s0t_0,s_0 lies between 1/21/2 and 1/2+ε/21/2+\varepsilon/2. Using y=t0q+εxαjy=t_0q+\varepsilon x\alpha_j and x=s0p−εyαkx=s_0p-\varepsilon y\alpha_k, we obtain

∣y−q/2∣≤ε,∣x−p/2∣≤ε.(14)|y-q/2|\le\varepsilon,\qquad|x-p/2|\le\varepsilon. \tag*{(14)}

Since ∣βj−q/2∣≤Δ/4|\beta_j-q/2|\le\Delta/4 and ∣βk−p/2∣≤Δ/4|\beta_k-p/2|\le\Delta/4, (3.4) gives ∣βj−y∣,∣βk−x∣≤ε+Δ/4|\beta_j-y|,|\beta_k-x|\le\varepsilon+\Delta/4. Thus the intercept heights t0=t−εxβjt_0=t-\varepsilon x\beta_j and s0=s+εyβks_0=s+\varepsilon y\beta_k have opposite first-order changes −εxy-\varepsilon xy and +εxy+\varepsilon xy. This cancellation is the reason for choosing βj\beta_j close to q/2q/2.

To control the remaining terms exactly, write a=t0−1/2a=t_0-1/2 and b=s0−1/2b=s_0-1/2. The two radial excesses satisfy

a(1−εxq)+b(1+εyp)=ε2(x2αj+y2αk)−εx(βj−q/2)+εy(βk−p/2).(15)\begin{aligned} a(1-\varepsilon xq)+b(1+\varepsilon yp) &=\varepsilon^2(x^2\alpha_j+y^2\alpha_k)\\ &\quad-\varepsilon x(\beta_j-q/2)+\varepsilon y(\beta_k-p/2). \tag*{(15)} \end{aligned}

For completeness, a+b=−εxβj+εyβka+b=-\varepsilon x\beta_j+\varepsilon y\beta_k, while the intercept equations give qa=y−εxαj−q/2qa=y-\varepsilon x\alpha_j-q/2 and pb=x+εyαk−p/2pb=x+\varepsilon y\alpha_k-p/2. Substitute these three identities into the left side of (3.5); the two mixed terms cancel.

For 0<ε<10<\varepsilon<1, both factors multiplying aa and bb are at least 1−ε>01-\varepsilon>0. Since a>ε2Mja>\varepsilon^2M_j, b>ε2Mkb>\varepsilon^2M_k and ∣x∣+∣y∣≤1|x|+|y|\le1, dividing (3.5) by ε2\varepsilon^2 yields

0>(Mj+Mk)(1−ε)−Δ4ε−(x2αj+y2αk).(16)0>(M_j+M_k)(1-\varepsilon)-\frac{\Delta}{4\varepsilon}-(x^2\alpha_j+y^2\alpha_k). \tag*{(16)}

We bound the last term without losing uniformity at the boundary. From (3.4) and ∣x∣,∣y∣≤1|x|,|y|\le1, ∣p∣,∣q∣≤1|p|,|q|\le1,

∣x2−p2/4∣, ∣y2−q2/4∣≤32ε.|x^2-p^2/4|,\ |y^2-q^2/4|\le\frac{3}{2}\varepsilon.

For a point inside its interval,

∣αj−1+q24∣≤Δ/2,∣αk−1+p24∣≤Δ/2.\left|\alpha_j-\frac{1+q^2}{4}\right|\le\Delta/2,\qquad\left|\alpha_k-\frac{1+p^2}{4}\right|\le\Delta/2.

Consequently

x2αj+y2αk≤p2+q2+2p2q216+32ε+Δ4≤d(p)+d(q)+32ε+Δ4.\begin{aligned} x^2\alpha_j+y^2\alpha_k &\le\frac{p^2+q^2+2p^2q^2}{16}+\frac{3}{2}\varepsilon+\frac{\Delta}{4}\\ &\le d(p)+d(q)+\frac{3}{2}\varepsilon+\frac{\Delta}{4}. \end{aligned}

The second inequality explains our choice of radial enlargement:

d(p)+d(q)−p2+q2+2p2q216=(p2−q2)216≥0.d(p)+d(q)-\frac{p^2+q^2+2p^2q^2}{16}=\frac{(p^2-q^2)^2}{16}\ge0.

The cutoff definitions give

d(p)+d(q)+2η≤Mj+Mk≤14+2η.d(p)+d(q)+2\eta\le M_j+M_k\le\frac{1}{4}+2\eta.

Using Δ≤ε2\Delta\le\varepsilon^2, the right side of (16) is therefore at least

2η−(2+2η)ε−ε24.(17)2\eta-(2+2\eta)\varepsilon-\frac{\varepsilon^2}{4}. \tag*{(17)}

For 0<ε≤η/20<\varepsilon\le\eta/2 this is at least η−17η2/16>0\eta-17\eta^2/16>0, a contradiction. All sector and cutoff inequalities are non-strict membership conditions. The argument applies to every point of the closed tetrahedron, including its faces, edges and vertices.

We retain the numerical estimate with a variable margin for the later construction in Section 6.

Corollary 3.2 (Quantitative coverage criterion). Let 0<τ<10<\tau<1 and η≥0\eta\ge0. Choose an integer nn with Δ=2/n≤τ2\Delta=2/n\le\tau^2, put qj=−1+jΔq_j=-1+j\Delta, and define αj,βj\alpha_j,\beta_j by eq:2.2. Use the triangles and directions eq:2.6–eq:2.7, with ε\varepsilon replaced by τ\tau and with cutoffs

Tj=12+τ2(η+max⁡qj≤q≤qj+1d(q)).T_j=\frac{1}{2}+\tau^2\left(\eta+\max_{q_j\le q\le q_{j+1}}d(q)\right).

These 2n2n cylinders cover KK whenever

2η−(2+2η)τ−τ24≥0.2\eta-(2+2\eta)\tau-\frac{\tau^2}{4}\ge0.

Proof. The angular first-crossing argument is unchanged. If both selected caps failed, their intercept heights would exceed 1/21/2 and have sum at most 1+τ/21+\tau/2. Thus the rough-angle estimates and the radial budget remain valid with the same constants, independently of η\eta. The only later use of the margin is through d(p)+d(q)+2η≤Mj+Mk≤1/4+2ηd(p)+d(q)+2\eta\le M_j+M_k\le1/4+2\eta. Consequently the same calculation as (16)–(17) gives 0>2η−(2+2η)τ−τ2/40>2\eta-(2+2\eta)\tau-\tau^2/4, contradicting (3.8).

For the fixed parameters of the main construction, the coverage threshold ε≤1/2000\varepsilon\le1/2000 will also guarantee strict area saving. We now bound the area remainder explicitly, so no further choice of a smaller tilt is needed.

The strict area decrease

Coverage is now established. We compare the perpendicular areas of the tilted triangles with the two-triangle cost H/2H/2. All error estimates below are uniform over the intervals of the partition.

Let S(ε)S(\varepsilon) denote the sum of the areas of all 2n2n perpendicular bases. Each intercept triangle in eq:2.6 and eq:2.7 has physical area HΔTj2/2H\Delta T_j^2/2. Orthogonal projection between planes multiplies area by the absolute inner product of their unit normals. For either family, the factor is

[1+ε2(αj2+H2βj2)]−1/2.\left[1+\varepsilon^2(\alpha_j^2+H^2\beta_j^2)\right]^{-1/2}.

We therefore have the exact finite formula

S(ε)H=∑j=0n−1ΔTj2[1+ε2(αj2+2βj2)]−1/2.(18)\frac{S(\varepsilon)}{H}=\sum_{j=0}^{n-1}\Delta T_j^2\left[1+\varepsilon^2(\alpha_j^2+2\beta_j^2)\right]^{-1/2}. \tag*{(18)}

Proposition 4.1. For the construction (2.1)–(2.7) and 0<ε≤10 < \varepsilon\le1,

∣S(ε)H−12−(2η−1240)ε2∣≤2ε4.\left|\frac{S(\varepsilon)}{H}-\frac{1}{2}-\left(2\eta-\frac{1}{240}\right)\varepsilon^2\right|\le2\varepsilon^4.

Proof. Put Aj=αj2+2βj2A_j=\alpha_j^2+2\beta_j^2 and

fj(z)=(12+Mjz)2(1+Ajz)−1/2.f_j(z)=\left(\frac{1}{2}+M_jz\right)^2(1+A_jz)^{-1/2}.

The bounds (2.5) imply 0≤Aj≤10\le A_j\le1 and 0≤Mj≤1/40\le M_j\le1/4. For 0≤z≤10\le z\le1, direct differentiation gives

fj′′(z)=2Mj2(1+Ajz)1/2−2AjMj(1/2+Mjz)(1+Ajz)3/2+3Aj2(1/2+Mjz)24(1+Ajz)5/2,f_j''(z)=\frac{2M_j^2}{(1+A_jz)^{1/2}}-\frac{2A_jM_j(1/2+M_jz)}{(1+A_jz)^{3/2}}+\frac{3A_j^2(1/2+M_jz)^2}{4(1+A_jz)^{5/2}},

so

∣fj′′(z)∣≤18+38+2764=5964.|f_j''(z)|\le\frac{1}{8}+\frac{3}{8}+\frac{27}{64}=\frac{59}{64}.

Taylor’s theorem at z=0z=0, followed by (4.1) and ∑jΔ=2\sum_j\Delta=2, now gives

∣S(ε)H−12−ε2∑jΔ(Mj−Aj8)∣≤5964ε4.(19)\left|\frac{S(\varepsilon)}{H}-\frac{1}{2}-\varepsilon^2\sum_j\Delta\left(M_j-\frac{A_j}{8}\right)\right|\le\frac{59}{64}\varepsilon^4. \tag*{(19)}

In particular, the growing number of sectors causes no loss in the error bound: each remainder is weighted by its interval length.

It remains to replace this finite sum by a polynomial integral. For q∈[qj,qj+1]q\in[q_j,q_{j+1}], set

A(q)=(1+q2)216+q22.A(q)=\frac{(1+q^2)^2}{16}+\frac{q^2}{2}.

Since ∣d′∣≤3/8|d'|\le3/8 on [−1,1][-1,1],

∣Mj−η−d(q)∣≤38Δ.|M_j-\eta-d(q)|\le\frac{3}{8}\Delta.

Also ∣αj−(1+q2)/4∣≤Δ/2|\alpha_j-(1+q^2)/4|\le\Delta/2 and ∣βj−q/2∣≤Δ/4|\beta_j-q/2|\le\Delta/4. The quantities in each of these differences have absolute value at most 1/21/2, so ∣Aj−A(q)∣≤Δ|A_j-A(q)|\le\Delta. Therefore

∣Mj−Aj8−(η+d(q)−A(q)8)∣≤12Δ.\left|M_j-\frac{A_j}{8}-\left(\eta+d(q)-\frac{A(q)}{8}\right)\right|\le\frac{1}{2}\Delta.

Integrating over all intervals shows that replacing the sum in (4.2) by ∫−11(η+d(q)−A(q)/8) dq\int_{-1}^{1}(\eta+d(q)-A(q)/8)\,dq costs at most ε2Δ≤ε4\varepsilon^2\Delta\le\varepsilon^4. Finally,

∫−11d(q) dq=115,∫−11A(q) dq=1730,\int_{-1}^{1}d(q)\,dq=\frac{1}{15},\qquad\int_{-1}^{1}A(q)\,dq=\frac{17}{30},

so this integral equals 2η−1/2402\eta-1/240. The combined error is at most (59/64+1)ε4<2ε4(59/64+1)\varepsilon^4<2\varepsilon^4, as claimed.

Proof of Theorem 1.1. The projection minimum is Lemma 2.1. The construction in Section 2 has exactly 2n2n compact triangular perpendicular bases, and Proposition 3.1 proves that its cylinders cover KK for 0<ε≤η/2=1/20000 < \varepsilon\le\eta/2 = 1/2000. Since η=1/1000\eta= 1/1000, Proposition 4.1 gives

S(ε)H≤12−136000ε2+2ε4<12(0<ε≤1/2000).\frac{S(\varepsilon)}{H} \le\frac{1}{2} - \frac{13}{6000}\varepsilon^2 + 2\varepsilon^4 < \frac{1}{2} \qquad(0 < \varepsilon\le1/2000).

where the strict inequality follows from 2ε2≤1/2,000,000<13/60002\varepsilon^2 \le1/2{,}000{,}000 < 13/6000. The same proposition gives the expansion and its stated remainder.

The same computation also controls a margin chosen separately for each tilt. This will let us make the margin vanish without hiding its area cost.

Corollary 4.2 (Area estimate with a variable margin). Let 0<τ≤10 < \tau\le1 and 0≤η≤1/80 \le\eta\le1/8. Choose an integer nn with Δ=2/n≤τ2\Delta= 2/n \le\tau^2, put qj=−1+jΔq_j = -1 + j\Delta, and use the secant coefficients (2.2) and cutoffs

Tj=12+τ2(η+max⁡qj≤q≤qj+1d(q)).T_j = \frac{1}{2} + \tau^2\left(\eta+ \max_{q_j \le q \le q_{j+1}} d(q)\right).

For the two families (2.6)–(2.7), with ε\varepsilon replaced by τ\tau, their total perpendicular base area CC satisfies

∣CH−12−(2η−1240)τ2∣≤5964τ4+τ2Δ≤2τ4.(20)\left|\frac{C}{H} - \frac{1}{2} - \left(2\eta- \frac{1}{240}\right)\tau^2\right| \le\frac{59}{64}\tau^4 + \tau^2\Delta\le2\tau^4. \tag*{(20)}

The estimate is uniform over all these choices; it does not assume that the cylinders cover $K.

Proof. Here Mj=η+max⁡[qj,qj+1]d≤1/4M_j = \eta+ \max_{[q_j,q_{j+1}]} d \le1/4 and Aj≤1A_j \le1, so the preceding second-derivative bound applies unchanged. It gives the Taylor error 59τ4/6459\tau^4/64. In the comparison with η+d(q)−A(q)/8\eta+ d(q) - A(q)/8, the constant η\eta cancels, leaving an error at most Δ/2\Delta/2 at each angular point. Integration therefore costs at most τ2Δ\tau^2\Delta. Neither bound requires η\eta to be fixed as τ\tau varies.

Other caps and angular partitions

The triangular construction leaves two choices available: how to end a sector radially, and where to sample its direction. The same balance between radial enlargement and perpendicular projection permits curved caps, smaller fixed margins, and sectors centred at the sampling nodes. We continue to use

H=2,K={(x,y,Hs):0≤s≤1, ∣x∣≤1−s, ∣y∣≤s},d(q)=q2+q416.H = \sqrt{2}, \qquad K = \{(x,y,Hs) : 0 \le s \le1,\ |x| \le1-s,\ |y| \le s\}, \qquad d(q) = \frac{q^2 + q^4}{16}.

Thus Amin⁡(K)=HA_{\min}(K) = H. All the cylinders below are specified by an intercept set and a nonzero direction; their bases are the perpendicular projections of the intercept sets.

A fixed positive margin

The following criterion isolates what the angular construction uses. Continuity of the displaced boundaries chooses a sector, and a positive quadratic margin ensures that the two chosen radial caps cannot both fail. The partition may be nonuniform, and its number of intervals need not satisfy an upper bound.

Theorem 5.1 (Compatible boundaries and radial caps). Let τ>0\tau> 0 tend to zero through any set accumulating at zero. For each τ\tau, choose a finite partition −1=q0<⋯<qm=1-1 = q_0 < \cdots< q_m = 1, constants αj,βj\alpha_j,\beta_j, and positive continuous functions RjR_j on Ij=[qj−1,qj]I_j = [q_{j-1},q_j]. Suppose that L(q)=αj−qβjL(q) = \alpha_j - q\beta_j is continuous on [−1,1][-1,1] and L(−1)=L(1)=0L(-1) = L(1) = 0. Suppose also, uniformly for every jj and q∈Ijq \in I_j, that

βj=q/2+O(τ2),αj=(1+q2)/4+O(τ2),(21)\beta_j = q/2 + O(\tau^2), \qquad\alpha_j = (1+q^2)/4 + O(\tau^2), \tag*{(21)}
Rj(q)=12+τ2(d(q)+η)+O(τ4),0<η<1/480.(22)R_j(q) = \frac{1}{2} + \tau^2(d(q)+\eta) + O(\tau^4), \qquad0 < \eta< 1/480. \tag*{(22)}

Here η\eta is fixed, and all error constants and their common range of validity are independent of τ,j,q\tau,j,q. Define

Pj−={(0,qS,HS):q∈Ij, 0≤S≤Rj(q)},vj−=(1,ταj,Hτβj),P_j^- = \{(0,qS,HS):q\in I_j,\ 0\le S\le R_j(q)\}, \qquad v_j^- = (1,\tau\alpha_j,H\tau\beta_j),
Pj+={(qT,0,H(1−T)):q∈Ij, 0≤T≤Rj(q)},vj+=(−ταj,1,Hτβj).P_j^+ = \{(qT,0,H(1-T)):q\in I_j,\ 0\le T\le R_j(q)\}, \qquad v_j^+ = (-\tau\alpha_j,1,H\tau\beta_j).

For every sufficiently small admissible τ\tau, the 2m2m cylinders Pj−+Rvj−P_j^-+\mathbb{R}v_j^- and Pj++Rvj+P_j^++\mathbb{R}v_j^+ cover the closed tetrahedron KK. Their bases are compact, and their total area CτC_\tau satisfies

CτH=12+(2η−1240)τ2+O(τ4).(23)\frac{C_\tau}{H} = \frac{1}{2} + \left(2\eta-\frac{1}{240}\right)\tau^2 + O(\tau^4). \tag*{(23)}

If each RjR_j is constant, every base is a triangle. The threshold and error constant may depend on η\eta and the uniform constants in the hypotheses, but not on mm.

Proof. Fix (x,y,Hs)∈K(x,y,Hs)\in K. The continuous piecewise affine function q↦qs+τxL(q)q\mapsto qs+\tau xL(q) has outer values −s,s-s,s, which bracket yy. Choose the first partition node after q0q_0 whose value is at least yy. Its preceding value is at most yy, by the initial endpoint inequality or by minimality. On the selected interval IjI_j, put

S=s−τxβj,Y=y−τxαj.S=s-\tau x\beta_j, \qquad Y=y-\tau x\alpha_j.

The crossing inequalities say qj−1S≤Y≤qjSq_{j-1}S\le Y\le q_jS, so S≥0S\ge0. If S=0S=0, then Y=0Y=0 and the selected cylinder already covers the point. Otherwise Y=qSY=qS for some q∈Ijq\in I_j. Applying the same argument to p↦p(1−s)−τyL(p)p\mapsto p(1-s)-\tau yL(p) gives an interval IiI_i and

T=1−s+τyβi≥0,X=x+τyαi=pT,p∈Ii.T=1-s+\tau y\beta_i\ge0, \qquad X=x+\tau y\alpha_i=pT, \qquad p\in I_i.

Again T=0T=0 already gives coverage. It remains to exclude S>Rj(q)S>R_j(q) and T>Ri(p)T>R_i(p) simultaneously. Assume these failures. For sufficiently small τ\tau, the fixed positive margin and the uniform cap error imply S,T>1/2S,T>1/2. Since S+T=1+τ(yβi−xβj)=1+O(τ)S+T=1+\tau(y\beta_i-x\beta_j)=1+O(\tau), each equals 1/2+O(τ)1/2+O(\tau). The intercept equations then give y=q/2+O(τ)y=q/2+O(\tau) and x=p/2+O(τ)x=p/2+O(\tau). Write ej=βj−q/2e_j=\beta_j-q/2 and ei=βi−p/2e_i=\beta_i-p/2. The same complementary radial balance as in the triangular construction now has explicit sampling errors:

s−τxy=12+(S−12)(1−τxq)−τ2x2αj+τxej,s-\tau xy=\frac{1}{2}+\left(S-\frac{1}{2}\right)(1-\tau xq)-\tau^2x^2\alpha_j+\tau xe_j,
1−s+τxy=12+(T−12)(1+τyp)−τ2y2αi−τyei.1-s+\tau xy=\frac{1}{2}+\left(T-\frac{1}{2}\right)(1+\tau yp)-\tau^2y^2\alpha_i-\tau ye_i.

These identities follow by substituting y=qS+τxαjy=qS+\tau x\alpha_j and x=pT−τyαix=pT-\tau y\alpha_i. Their left sides sum to one, and both factors multiplying the radial excesses are positive for small τ\tau. Substitute the strict cap failures, divide by τ2\tau^2, and use ei,ej=O(τ2)e_i,e_j=O(\tau^2) to obtain

0>d(q)+d(p)+2η−x2αj−y2αi+O(τ)=2η+(p2−q2)216+O(τ).0>d(q)+d(p)+2\eta-x^2\alpha_j-y^2\alpha_i+O(\tau)=2\eta+\frac{(p^2-q^2)^2}{16}+O(\tau).

The error is uniform. This is impossible for sufficiently small τ\tau, proving coverage, including the boundary.

The parametrizations of Pj−P_j^{-} and Pj+P_j^{+} have area Jacobians HSHS and HTHT, respectively, and are one-to-one off their tips. Projection multiplies area in either intercept plane by [1+τ2(αj2+2βj2)]−1/2[1+\tau^2(\alpha_j^2+2\beta_j^2)]^{-1/2}. Consequently

CτH=∑j∫IjRj(q)21+τ2(αj2+2βj2) dq.\frac{C_\tau}{H}=\sum_j\int_{I_j}\frac{R_j(q)^2}{\sqrt{1+\tau^2(\alpha_j^2+2\beta_j^2)}}\,dq.

Uniform expansion of the integrand gives

14+τ2[d(q)+η−18((1+q2)216+q22)]+O(τ4).\frac{1}{4}+\tau^2\left[d(q)+\eta-\frac{1}{8}\left(\frac{(1+q^2)^2}{16}+\frac{q^2}{2}\right)\right]+O(\tau^4).

The total angular length is two, so integration preserves the uniform fourth-order remainder regardless of the number of intervals. Finally

∫−11d(q) dq=115,18∫−11((1+q2)216+q22)dq=17240.\int_{-1}^{1}d(q)\,dq=\frac{1}{15},\qquad\frac{1}{8}\int_{-1}^{1}\left(\frac{(1+q^2)^2}{16}+\frac{q^2}{2}\right)dq=\frac{17}{240}.

which proves (23). Each parameter domain is compact; its continuous image and perpendicular projection are compact as well. Projection is invertible between the intercept plane and the base plane, so a constant cap gives a nondegenerate triangle.

Triangular and curved apertures

For an ordinary interval I=[l,r]I=[l,r], use

αI=1+lr4,βI=l+r4.(24)\alpha_I=\frac{1+lr}{4},\qquad\beta_I=\frac{l+r}{4}. \tag*{(24)}

The identity αI−qβI=(1−q2)/4\alpha_I-q\beta_I=(1-q^2)/4 at both endpoints matches adjacent boundaries exactly. If the maximum interval length is at most Mτ2M\tau^2, for fixed MM, then (5.1) holds uniformly. With rI=(l+r)/2r_I=(l+r)/2, each of the following choices satisfies (5.2):

RI(q)=12+τ2(d(rI)+η),(25)R_I(q)=\frac{1}{2}+\tau^2(d(r_I)+\eta), \tag*{(25)}

midpoint triangles,

RI(q)=12+τ2(max⁡u∈Id(u)+η),(26)R_I(q)=\frac{1}{2}+\tau^2\left(\max_{u\in I}d(u)+\eta\right), \tag*{(26)}

maximum triangles,

RI(q)=12+τ2(d(q)+η),(27)R_I(q)=\frac{1}{2}+\tau^2(d(q)+\eta), \tag*{(27)}

continuous radial caps,

RI(q)=14+τ2(d(q)+η),(28)R_I(q)=\sqrt{\frac{1}{4}+\tau^2(d(q)+\eta)}, \tag*{(28)}

squared-radius caps.

Indeed, dd has bounded derivative on [−1,1][-1,1], so replacing d(q)d(q) by its midpoint value or its cell maximum costs O(τ2)O(\tau^2) inside the parentheses. For (5.8), Taylor expansion at 1/41/4 gives the required fourth-order error. Thus all four choices have the area formula (23). The first two give triangles; the last two give compact radial sectors.

Another choice specifies the radial enlargement through the transverse intercept. Put D(Y)=Y2/4+Y4D(Y)=Y^2/4+Y^4. For b=3/4b=3/4 or b=1b=1, retain the bounded aperture

0≤S≤b,S≤12+τ2(D(qS)+η).(29)0\le S\le b,\qquad S\le\frac{1}{2}+\tau^2(D(qS)+\eta). \tag*{(29)}

For small τ\tau, this is exactly 0≤S≤R(q)0 \le S \le R(q), where RR is continuous and satisfies (22). To verify this, consider f(S)=S−1/2−τ2(D(qS)+η)f(S)=S-1/2-\tau^2(D(qS)+\eta) on [0,1][0,1]. Uniformly in ∣q∣≤1|q| \le1,

f′(S)=1−τ2(q2S/2+4q4S3)≥1−92τ2.f'(S)=1-\tau^2(q^2S/2+4q^4S^3)\ge1-\frac{9}{2}\tau^2.

For small τ\tau, ff is strictly increasing, is negative at 1/21/2, and is positive at 3/43/4. Its unique zero R(q)R(q) lies between these values and depends continuously on qq. Its defining equation first gives R(q)−1/2=O(τ2)R(q)-1/2=O(\tau^2), and then

D(qR(q))=D(q/2)+O(τ2)=d(q)+O(τ2),D(qR(q))=D(q/2)+O(\tau^2)=d(q)+O(\tau^2),

which proves the cap expansion. The outer bound in (29) is part of the definition: without it, the quartic inequality can admit a second, unbounded component.

The common conclusion is easiest to compare in the fixed tetrahedron KK: for each fixed 0<η<1/4800<\eta<1/480, every cap above has normalized total area (23). Its threshold is uniform over partitions whose maximum interval length is at most Mτ2M\tau^2, for a fixed MM. For example, an equal partition with any integer m≥2/τ2m\ge2/\tau^2 gives exactly 2m2m cylinders, with the same threshold and error constant as the mesh is refined. Each cover is finite; no further limit in the number of cylinders is needed after choosing the tilt and mesh.

The squared-radius and Cartesian formulas also explain the geometry in two different ways. Appendix A.1 compares the two radial squares, while Appendix A.2 uses a discriminant to estimate the intercept height without locating the perturbed interface. These are alternative proofs of cases already covered by Theorem 5.1. Exact counts and costs for several rescaled parameter choices are collected in Appendix B.

Sampling at the endpoints

The compatibility condition also permits directions sampled at nodes rather than along secants of each sector. Let k>1k>1 be an integer, τ=1/k\tau=1/k, n=k2n=k^2, and qj=−1+2j/nq_j=-1+2j/n for 0≤j≤n0\le j\le n. Put λ0=−1\lambda_0=-1, λn+1=1\lambda_{n+1}=1, and λj=(qj−1+qj)/2\lambda_j=(q_{j-1}+q_j)/2 for 1≤j≤n1\le j\le n. Use the sectors [λj,λj+1][\lambda_j,\lambda_{j+1}] and parameters

βj=qj/2,αj=(1+qj2)/4,Rj=12+τ2(d(qj)+1/1000).\beta_j=q_j/2,\qquad\alpha_j=(1+q_j^2)/4,\qquad R_j=\frac{1}{2}+\tau^2(d(q_j)+1/1000).

At an interior boundary, αj−αj−1=λj(βj−βj−1)\alpha_j-\alpha_{j-1}=\lambda_j(\beta_j-\beta_{j-1}); at the exterior boundaries αj−λβj=0\alpha_j-\lambda\beta_j=0. Thus the piecewise affine boundary function is continuous and vanishes at both ends. Every angular point of the jjth sector is within τ2\tau^2 of qjq_j, so all the hypotheses of Theorem 5.1 hold. There are n+1n+1 sectors, with widths τ2\tau^2 at the two ends and 2τ22\tau^2 in the interior. Consequently, for all sufficiently large kk, this construction covers KK with exactly 2(k2+1)2(k^2+1) triangular bases and normalized total area

CτAmin⁡(K)=12−136000τ2+O(τ4).(30)\frac{C_\tau}{A_{\min}(K)}=\frac{1}{2}-\frac{13}{6000}\tau^2+O(\tau^4). \tag*{(30)}

The half-width end sectors preserve the exact exterior boundaries. Appendix A.3 gives a complementary radial identity that verifies their coverage directly, and Appendix B.2 records the corresponding absolute area at edge length 2\sqrt{2}.

A vanishing radial margin

For a fixed positive margin, the normalized area coefficient is 2η−1/2402\eta-1/240. We now let the margin tend to zero with the tilt and obtain the coefficient −1/240-1/240 in a single family of finite triangular covers. The extra radial enlargement is cubic in the tilt. Its positive cubic area term explains why coverage persists on a much wider interval than the interval on which we prove a saving.

The fixed-margin theorem, Theorem 5.1, cannot be applied uniformly when η\eta tends to zero. Instead we use the explicit coverage budget in Corollary 3.2 and the uniform area estimate in Corollary 4.2. We keep the edge-two tetrahedron KK and the same slope coordinates.

Proposition 6.1. For every 0<τ≤10 < \tau\le1, the tetrahedron KK admits a cover by 2⌈8/τ2⌉2\lceil8/\tau^2\rceil cylinders with compact, nondegenerate triangular perpendicular bases. Their total base area CτC_\tau satisfies

∣CτH−12+τ2240−252τ3∣≤2τ4(0<τ≤1/50).(31)\left| \frac{C_\tau}{H} - \frac{1}{2} + \frac{\tau^2}{240} - \frac{25}{2}\tau^3 \right| \le2\tau^4 \qquad(0 < \tau\le1/50). \tag*{(31)}

In particular, Cτ<Amin⁡(K)/2C_\tau< A_{\min}(K)/2 for 0<τ≤1/40000 < \tau\le1/4000, and

CτAmin⁡(K)=12−τ2240+O(τ3)(τ↓0).\frac{C_\tau}{A_{\min}(K)} = \frac{1}{2} - \frac{\tau^2}{240} + O(\tau^3) \qquad(\tau\downarrow0).

Proof. Put

N=⌈8τ2⌉,Δ=2N≤τ24,qj=−1+jΔ(0≤j≤N).N = \left\lceil\frac{8}{\tau^2} \right\rceil,\qquad\Delta= \frac{2}{N} \le\frac{\tau^2}{4},\qquad q_j = -1 + j\Delta\qquad(0 \le j \le N).

Use the secant coefficients (2.2), and set

dj=max⁡qj≤q≤qj+1d(q),Tj=12+τ2dj+254τ3(0≤j<N).(32)d_j = \max_{q_j \le q \le q_{j+1}} d(q),\qquad T_j = \frac{1}{2} + \tau^2 d_j + \frac{25}{4}\tau^3 \qquad(0 \le j < N). \tag*{(32)}

The intercept triangles and directions are (2.6)--(2.7), with ε\varepsilon replaced by τ\tau and these cutoffs. They give exactly 2N2N cylinders. Their intercept triangles have positive area, and the normal component of each axis is one, so their perpendicular bases are compact nondegenerate triangles as before.

Coverage. For 0<τ≤1/20 < \tau\le1/2, take η=25τ/4\eta= 25\tau/4 in Corollary 3.2. Its sufficient budget is

2η−(2+2η)τ−τ24=τ(212−514τ)≥338τ>0.2\eta- (2 + 2\eta)\tau- \frac{\tau^2}{4} = \tau\left(\frac{21}{2} - \frac{51}{4}\tau\right) \ge\frac{33}{8}\tau> 0.

Thus the cylinders cover KK throughout this range.

For 1/2≤τ≤11/2 \le\tau\le1, the enlarged triangles are already large enough for the first family alone. Indeed, for (x,y,Ht)∈K(x,y,Ht) \in K the first-crossing argument selects a first-family sector with intercept height

0≤t0=t−τxβj≤t+τ2∣x∣≤t+τ2(1−t)≤1.0 \le t_0 = t - \tau x\beta_j \le t + \frac{\tau}{2}|x| \le t + \frac{\tau}{2}(1-t) \le1.

Its cutoff obeys Tj≥1/2+(25/4)τ3≥41/32>1T_j \ge1/2 + (25/4)\tau^3 \ge41/32 > 1. The selected cylinder therefore contains the point. Both arguments use non-strict angular membership conditions, so they cover the closed body, including its boundary.

Area. For 0<τ≤1/500 < \tau\le1/50, the margin η=25τ/4\eta= 25\tau/4 lies in [0,1/8][0,1/8]. Corollary 4.2 consequently gives

∣CτH−12−(252τ−1240)τ2∣≤2τ4,\left| \frac{C_\tau}{H} - \frac{1}{2} - \left(\frac{25}{2}\tau- \frac{1}{240}\right)\tau^2 \right| \le2\tau^4,

which is (6.1). This application is uniform although the margin varies with the tilt. If 0<τ≤1/40000 < \tau\le1/4000, then

252τ+2τ2≤1320+18,000,000<1240.\frac{25}{2}\tau+ 2\tau^2 \le\frac{1}{320} + \frac{1}{8{,}000{,}000} < \frac{1}{240}.

The negative quadratic term therefore exceeds both the positive cubic term and the error, giving Cτ<H/2=Amin⁡(K)/2C_\tau< H/2 = A_{\min}(K)/2. ∎

The equivalent edge-four construction, with its full coverage interval, exact cylinder count and scaled area estimates, is recorded in Remark B.1 of Appendix B.

Alternative coverage and area calculations

The fixed-margin criterion already proves coverage for the cap choices in Section 5. The following direct calculations explain three particular features: why squared-radius caps suit a comparison of radial squares, how a Cartesian discriminant avoids locating the perturbed interface, and how node-centred sectors retain the complementary radial cancellation. They are independent explanations of those cases, rather than inputs to the main construction or the vanishing-margin argument.

Squared radii

Let nn be a positive integer, put k=1/nk = 1/n, use n2n^2 equal intervals, and write their width as Δ=2k2\Delta= 2k^2. Put

MI=(l+r)/2=2βI,BI=(1+lr)/2=2αI,f(q)=(q2+q4)/4.M_I = (l+r)/2 = 2\beta_I,\qquad B_I = (1+lr)/2 = 2\alpha_I,\qquad f(q) = (q^2+q^4)/4.

The squared-radius cap (5.8), with τ=2k\tau= 2k and η=1/960\eta= 1/960, is R(q)2=1/4+k2(f(q)+1/240)R(q)^2 = 1/4 + k^2(f(q) + 1/240). The elementary secant identity is exact:

2BI−MI2=1−Δ24.2B_I - M_I^2 = 1 - \frac{\Delta^2}{4}.

For the two angularly selected intercepts, write

S=s−kxMI,y=qS+kxBI,T=1−s+kyMJ,x=pT−kyBJ.S = s-kxM_I,\qquad y = qS+kxB_I,\qquad T = 1-s+kyM_J,\qquad x = pT-kyB_J.

The crossing argument supplies S,T≥0S,T \ge0. If both radial caps fail, then S,T>1/2S,T > 1/2 and u=s−1/2=O(k)u = s-1/2 = O(k). Since MI=q+O(k2)M_I = q+O(k^2) and MJ=p+O(k2)M_J = p+O(k^2), substitution in

S2=s2−2kxSMI−k2x2MI2,T2=(1−s)2+2kyTMJ−k2y2MJ2S^2 = s^2-2kxSM_I-k^2x^2M_I^2,\qquad T^2 = (1-s)^2+2kyTM_J-k^2y^2M_J^2

and use of (A.1) give the uniform expansions

S2=14+u+u2−2kxy+k2x2+O(k3),T2=14−u+u2+2kxy+k2y2+O(k3).(33)\begin{aligned} S^2 &= \frac{1}{4}+u+u^2-2kxy+k^2x^2+O(k^3),\\ T^2 &= \frac{1}{4}-u+u^2+2kxy+k^2y^2+O(k^3). \tag*{(33)} \end{aligned}

For example SMI=y−kxBI+O(k2)SM_I = y-kxB_I+O(k^2), and the second identity uses TMJ=x+kyBJ+O(k2)TM_J = x+kyB_J+O(k^2). Since ∣MI∣,∣MJ∣≤1|M_I|,|M_J| \le1, the failures S,T>1/2S,T > 1/2 first give ∣u∣≤k|u| \le k. The first strict lower bound on S2S^2 then gives u−2kxy≥−Ck2u-2kxy \ge-Ck^2, and the second gives u−2kxy≤Ck2u-2kxy \le Ck^2, for a uniform constant CC. Thus u=2kxy+O(k2)u = 2kxy+O(k^2). Adding (33) now yields

S2+T2−1/2k2=x2+y2+8x2y2+O(k).\frac{S^2+T^2-1/2}{k^2}=x^2+y^2+8x^2y^2+O(k).

On the other hand q=2y+O(k)q = 2y+O(k) and p=2x+O(k)p = 2x+O(k), so

f(q)+f(p)=x2+y2+4(x4+y4)+O(k)≥x2+y2+8x2y2+O(k).f(q)+f(p)=x^2+y^2+4(x^4+y^4)+O(k)\ge x^2+y^2+8x^2y^2+O(k).

The failed caps require the previous quotient to exceed f(q)+f(p)+2/240f(q)+f(p)+2/240, a contradiction for small kk.

Let CkC_k denote the sum of the perpendicular base areas of all 2n22n^2 cylinders. The two families have equal total area, and the square in the cap makes the radial integration exact before the projection factor is expanded:

Ck2H=∑I∫I12[14+k2(f(q)+1240)]1+k2(BI2+2MI2) dq=14+k2(1240−1120)+O(k4).\frac{C_k}{2H}=\sum_I\int_I\frac{\frac{1}{2}\left[\frac{1}{4}+k^2\left(f(q)+\frac{1}{240}\right)\right]}{\sqrt{1+k^2\left(B_I^2+2M_I^2\right)}}\,dq=\frac{1}{4}+k^2\left(\frac{1}{240}-\frac{1}{120}\right)+O(k^4).

The squared cap therefore supplies both the coverage comparison and the exact radial part of the area integral.

A Cartesian discriminant

For an equal sector I=[l,r]I=[l,r] of width Δ\Delta, the coefficients (24) satisfy

βI2−αI+14=Δ216.\beta_I^2-\alpha_I+\frac{1}{4}=\frac{\Delta^2}{16}.

Suppose the angular crossing for one family has given

S=T−δβI≥0,Y0=Y−δαI,lS≤Y0≤rS.S=T-\delta\beta_I\geq0,\qquad Y_0=Y-\delta\alpha_I,\qquad lS\leq Y_0\leq rS.

Here (T,Y,δ)=(s,y,τx)(T,Y,\delta)=(s,y,\tau x) in the first family and (1−s,x,−τy)(1-s,x,-\tau y) in the second. Because 2βI2\beta_I is the sector midpoint, ∣2βIS−Y0∣≤ΔS/2\lvert2\beta_I S-Y_0\rvert\leq\Delta S/2. Consequently

Q:=T2−δY+δ24=S2+δ(2βIS−Y0)+δ2Δ216≥(S−∣δ∣Δ4)2.(34)\begin{aligned} Q&:=T^2-\delta Y+\frac{\delta^2}{4}\\ &=S^2+\delta(2\beta_I S-Y_0)+\frac{\delta^2\Delta^2}{16}\geq\left(S-\frac{|\delta|\Delta}{4}\right)^2. \tag*{(34)} \end{aligned}

In particular Q≥S−∣δ∣Δ/4\sqrt{Q}\geq S-|\delta|\Delta/4. Take the Cartesian cap (29) with b=1b=1, fixed η>0\eta>0, and Δ=O(τ2)\Delta=O(\tau^2). Since 0≤αI≤1/20\leq\alpha_I\leq1/2, ∣Y∣≤1|Y|\leq1 and ∣δ∣≤τ|\delta|\leq\tau, we have ∣Y0∣≤1+τ/2|Y_0|\leq1+\tau/2. For τ≤1\tau\leq1, this gives D(Y0)≤D(3/2)=45/8D(Y_0)\leq D(3/2)=45/8. Hence the polynomial upper bound is less than one whenever τ2(45/8+η)<1/2\tau^2(45/8+\eta)<1/2. Failure of the selected cylinder therefore implies S>1/2+τ2(D(Y0)+η)S>1/2+\tau^2(D(Y_0)+\eta) even if the outer bound S≤1S\leq1 failed. Since Y0=Y+O(τ)Y_0=Y+O(\tau) and ∣δ∣Δ=O(τ3)|\delta|\Delta=O(\tau^3), (A.3) gives

Q>12+τ2(D(Y)+η2).\sqrt{Q}>\frac{1}{2}+\tau^2\left(D(Y)+\frac{\eta}{2}\right).

for sufficiently small τ\tau, uniformly in both families. If both failed, squaring these positive lower bounds and discarding their nonnegative fourth-order terms would give

s2>14+τxy+τ2(D(y)+η2−x24),s^2>\frac{1}{4}+\tau xy+\tau^2\left(D(y)+\frac{\eta}{2}-\frac{x^2}{4}\right),
(1−s)2>14−τxy+τ2(D(x)+η2−y24).(1-s)^2>\frac{1}{4}-\tau xy+\tau^2\left(D(x)+\frac{\eta}{2}-\frac{y^2}{4}\right).

The right sides are positive for small τ\tau. Expanding their square roots with 1/4+v=1/2+v−v2+O(v3)\sqrt{1/4+v}=1/2+v-v^2+O(v^3) shows that their sum is

1+τ2(D(x)+D(y)+η−x2+y24−2x2y2)+O(τ3)=1+τ2((x2−y2)2+η)+O(τ3)>1,1+\tau^2\left(D(x)+D(y)+\eta-\frac{x^2+y^2}{4}-2x^2y^2\right)+O(\tau^3)=1+\tau^2\left((x^2-y^2)^2+\eta\right)+O(\tau^3)>1,

contradicting s+(1−s)=1s+(1-s)=1.

For either outer bound b∈{3/4,1}b\in\{3/4,1\}, the Cartesian aperture admits a direct area calculation without changing to angular variables at the boundary. In the intercept coordinates (Y,S)(Y,S), put

Dτ={0≤S≤b, ∣Y∣≤S, S≤12+τ2(D(Y)+η)},D0={0≤S≤12, ∣Y∣≤S}.D_\tau=\left\{0\leq S\leq b,\ |Y|\leq S,\ S\leq\frac{1}{2}+\tau^2(D(Y)+\eta)\right\},\qquad D_0=\left\{0\leq S\leq\frac{1}{2},\ |Y|\leq S\right\}.

For ∣Y∣≤1/2|Y| \le1/2 the added vertical length is exactly τ2(D(Y)+η)\tau^2(D(Y) + \eta) for small τ\tau. Each remaining side region has both width and height O(τ2)O(\tau^2), because 1/2<∣Y∣≤S≤1/2+O(τ2)1/2 < |Y| \le S \le1/2 + O(\tau^2) there. Thus

∣Dτ∣=14+τ2∫−1/21/2(D(Y)+η) dY+O(τ4).|D_\tau| = \frac{1}{4} + \tau^2 \int_{-1/2}^{1/2} (D(Y) + \eta)\,\mathrm{d}Y + O(\tau^4).

Expanding the projection factor over D0D_0, and using Y=qSY = qS in its quadratic coefficient, gives the one-family cost

14+τ2[∫−1/21/2D(Y) dY⏟1/30+η−12∫01/2S dS∫−11((1+q2)216+q22) dq⏟17/480]+O(τ4).\frac{1}{4} + \tau^2 \left[ \underbrace{\int_{-1/2}^{1/2} D(Y)\,\mathrm{d}Y}_{1/30} + \eta- \frac{1}{2}\underbrace{\int_0^{1/2} S\,\mathrm{d}S \int_{-1}^{1} \left(\frac{(1+q^2)^2}{16} + \frac{q^2}{2}\right)\,\mathrm{d}q}_{17/480} \right] + O(\tau^4).

Here the cost is divided by HH in the edge-two normalization. The replacement of DτD_\tau by D0D_0 changes the quadratic coefficient by O(τ2)O(\tau^2), and the coefficient sampling error is also O(τ2)O(\tau^2); both contribute only O(τ4)O(\tau^4) to the area. Thus the normalized two-family cost is 1/2+(2η−1/240)τ2+O(τ4)1/2 + (2\eta- 1/240)\tau^2 + O(\tau^4), as in (5.3). The Cartesian entries in Appendix B follow by substitution; the edge-one discriminant and aperture are translated there as well.

A direct check for node-centred sectors

Use the node-centred sectors of Section 5.3, and write their tilt as t=1/kt = 1/k. The radial cutoff at a sample node qq is 1/2+t2(d(q)+1/1000)1/2 + t^2(d(q) + 1/1000), and each angular point of its sector lies within t2t^2 of qq. If the two selected sample nodes are q,pq,p, write

S=s−txq2,y=qS+txα(q)+e,T=1−s+typ2,x=pT−tyα(p)+e∗,S = s - \frac{txq}{2}, \qquad y = qS + tx\alpha(q) + e, \qquad T = 1 - s + \frac{typ}{2}, \qquad x = pT - ty\alpha(p) + e_*,

where ∣e∣≤St2|e| \le St^2, ∣e∗∣≤Tt2|e_*| \le Tt^2, and α(q)=(1+q2)/4\alpha(q) = (1 + q^2)/4. Direct substitution gives

s−txy=12+(S−12)(1−tqx)−t2α(q)x2−txe,s - txy = \frac{1}{2} + \left(S - \frac{1}{2}\right)(1-tqx) - t^2\alpha(q)x^2 - txe,
1−s+txy=12+(T−12)(1+tpy)−t2α(p)y2+tye∗.1 - s + txy = \frac{1}{2} + \left(T - \frac{1}{2}\right)(1+tpy) - t^2\alpha(p)y^2 + tye_*.

If both radial caps failed along a sequence t→0t \to0, then S,T>1/2S,T > 1/2 and S+T=1+t(yp−xq)/2=1+O(t)S + T = 1 + t(yp - xq)/2 = 1 + O(t), so S,T→1/2S,T \to1/2. Passing to a subsequence with q→q0q \to q_0, p→p0p \to p_0 gives y→q0/2y \to q_0/2, x→p0/2x \to p_0/2, while e/t,e∗/t→0e/t,e_*/t \to0. Adding the identities and dividing by t2t^2 would therefore imply

0≥21000+(q02−p02)216>0.0 \ge\frac{2}{1000} + \frac{(q_0^2-p_0^2)^2}{16} > 0.

This complementary identity gives a second explanation for why the half-width endpoint sectors still cover the boundary.

Counts and formulas under changes of scale

The constructions in the main text use the edge-two tetrahedron KK. A similarity of ratio a/2a/2 produces a regular tetrahedron of edge length aa: it multiplies every base area and the minimum projection area by (a/2)2(a/2)^2, while preserving the number of cylinders and the normalized cost. This appendix records several parameter choices and their coordinate dictionaries.

Secant sectors

In the following table, every row uses an equal partition with mm intervals and hence exactly 2m2m cylinders. The body has edge length aa and is obtained from KK by a similarity of ratio a/2a/2; accordingly the area is reported as C/Amin⁡C/A_{\min}. Each statement holds for all sufficiently large indicated integers or all sufficiently small indicated positive real parameters. Every displayed remainder is fourth order in the displayed tilt parameter, uniformly over the permitted finer meshes.

Capaa(τ,η)(\tau,\eta)mmC/Amin⁡C/A_{\min}
Midpoint44(2/n,1/1920)(2/n,1/1920)2n22n^212−180n2+O(n−4)\frac{1}{2}-\frac{1}{80n^2}+O(n^{-4})
Midpoint11(α,1/1000)(\alpha,1/1000)⌈α−2⌉\lceil\alpha^{-2}\rceil12−136000α2+O(α4)\frac{1}{2}-\frac{13}{6000}\alpha^2+O(\alpha^4)
Midpoint22(1/n,1/1000)(1/n,1/1000)n2n^212−136000n2+O(n−4)\frac{1}{2}-\frac{13}{6000n^2}+O(n^{-4})
Continuous22(c,1/1000)(c,1/1000)any m≥2/c2m\ge2/c^212−136000c2+O(c4)\frac{1}{2}-\frac{13}{6000}c^2+O(c^4)
Squared radius22(2/n,1/960)(2/n,1/960)n2n^212−1120n2+O(n−4)\frac{1}{2}-\frac{1}{120n^2}+O(n^{-4})
Maximum11(1/(4n),1/960)(1/(4n),1/960)n2n^212−17680n2+O(n−4)\frac{1}{2}-\frac{1}{7680n^2}+O(n^{-4})
Cartesian, b=3/4b=3/42\sqrt{2}(λ/2,1/960)(\lambda/2,1/960)⌈2/λ2⌉\lceil\sqrt{2}/\lambda^2\rceil12−λ21920+O(λ4)\frac{1}{2}-\frac{\lambda^2}{1920}+O(\lambda^4)
Cartesian, b=1b=111(1/(2n),1/1000)(1/(2n),1/1000)n2n^212−1324000n2+O(n−4)\frac{1}{2}-\frac{13}{24000n^2}+O(n^{-4})

Table 1.

The table follows by substitution in (5.3); in each row the mesh is bounded by a fixed multiple of τ2\tau^2. For the continuous cap, m=⌈2/c2⌉m=\lceil2/c^2\rceil gives the ceiling count, whereas c=1/nc=1/n permits every integer m≥2n2m\ge2n^2, with the same threshold and error constant. The real-parameter midpoint row includes the reciprocal-integer construction after dilation by two. For the squared-radius row, the two families have equal area and each costs H[1/4+(1/240−1/120)n−2+O(n−4)]H[1/4+(1/240-1/120)n^{-2}+O(n^{-4})]. For the first Cartesian row, either family costs c0/4−c0λ2/3840+O(λ4)c_0/4-c_0\lambda^2/3840+O(\lambda^4) with c0=1/2c_0=1/\sqrt{2}. These are genuine finite covers: the positive margin is fixed first, then a sufficiently small tilt and the specified finite mesh are chosen. No further limit in the number of cylinders is needed.

Node-centred sectors at edge length 2\sqrt{2}

For the endpoint construction of Section 5.3, write t=1/kt=1/k. Its 2(k2+1)2(k^2+1) cylinders retain that exact count under scaling. At edge length 2\sqrt{2}, the similarity ratio is 1/21/\sqrt{2} and the minimum projection area is H/2=1/2H/2=1/\sqrt{2}. Multiplying (5.10) by this area gives total base area

22(14−1312000t2+O(t4)).(35)\frac{2}{\sqrt{2}}\left(\frac{1}{4}-\frac{13}{12000}t^2+O(t^4)\right). \tag*{(35)}

The Cartesian discriminant at edge length one

The following dictionary translates the Cartesian argument of Appendix A.2. In its edge-two coordinates, (x,y,Hs)(x,y,Hs) is the point of KK, Y0Y_0 is the transverse intercept, SS is its radial height, and TT is ss or 1−s1-s in the respective family. The discriminant is (A.3). At edge-one scale, writing the point as (x1,y1,s/2)(x_1,y_1,s/\sqrt{2}) and putting τ=ε/2\tau=\varepsilon/2, the substitutions x=2x1x=2x_1, y=2y1y=2y_1, η=1/1000\eta=1/1000 give the equivalent obstruction ε2[4(x1−y12)+1/4000]+O(ε3)\varepsilon^2[4(x_1-y_1^2)+1/4000]+O(\varepsilon^3). In these coordinates the point is (Y1,δ)(Y_1,\delta) with (y1,εx1)(y_1,\varepsilon x_1) or (x1,−εy1)(x_1,-\varepsilon y_1) in the two families. The transverse intercept is b1=Y0/2=Y1−δα/2b_1=Y_0/2=Y_1-\delta\alpha/2. For S>0S>0, write q=Y0/Sq=Y_0/S. The edge-one slope coordinate is q1=b1/S=q/2q_1=b_1/S=q/2, whose sector has width d1=Δ/2d_1=\Delta/2. The same discriminant is Q=T2−2δY1+δ2/4Q=T^2-2\delta Y_1+\delta^2/4, and its bound is Q≥S−∣δ∣d1/2\sqrt{Q}\ge S-|\delta|d_1/2.

Writing F(v)=v2/4+4v4F(v)=v^2/4+4v^4, the corresponding bounded aperture is

0≤S≤1,l2S≤b1≤r2S,S≤1/2+ε2(F(b1)+1/4000).0\le S\le1,\qquad\frac{l}{2}S\le b_1\le\frac{r}{2}S,\qquad S\le1/2+\varepsilon^2(F(b_1)+1/4000).

The outer bound S≤1S\le1 remains part of the definition after scaling.

The cubic buffer at edge length four

Remark B.1 (The edge-four normalization). The similarity

(x,y,Ht)⟼(2x,2y,H(2t−1))(x,y,Ht)\longmapsto(2x,2y,H(2t-1))

sends KK onto the regular edge-four tetrahedron

K4={(x,y,Hz):−1≤z≤1, ∣x∣≤1−z, ∣y∣≤1+z}.K_4=\{(x,y,Hz):-1\le z\le1,\ \lvert x\rvert\le1-z,\ \lvert y\rvert\le1+z\}.

It multiplies every area by four and leaves the cylinder count unchanged. Writing ε=τ/2\varepsilon=\tau/2, the same cover has N=⌈2/ε2⌉N=\lceil2/\varepsilon^2\rceil sectors in each family. On a sector I=[l,r]I=[l,r], write TIT_I for its cutoff in (6.2), and put

aI=2αI=1+lr2,bI=2βI=l+r2,DI=max⁡q∈Iq2+q42=8max⁡q∈Id(q).a_I=2\alpha_I=\frac{1+lr}{2},\qquad b_I=2\beta_I=\frac{l+r}{2},\qquad D_I=\max_{q\in I}\frac{q^2+q^4}{2}=8\max_{q\in I}d(q).

The two axes are (1,εaI,HεbI)(1,\varepsilon a_I,H\varepsilon b_I) and (−εaI,1,HεbI)(-\varepsilon a_I,1,H\varepsilon b_I). Their radial intercept coordinates are 1+z1+z and 1−z1-z, with common cutoff

2TI=1+δI,δI=ε2DI+100ε3.2T_I=1+\delta_I,\qquad\delta_I=\varepsilon^2D_I+100\varepsilon^3.

Thus these 2⌈2/ε2⌉2\lceil2/\varepsilon^2\rceil triangular cylinders cover K4K_4 for every 0<ε≤1/20<\varepsilon\le1/2. Since Amin⁡(K4)=4HA_{\min}(K_4)=4H, their total area Sε=4C2εS_\varepsilon=4C_{2\varepsilon} satisfies

Sε=22−215ε2+O(ε3),S_\varepsilon=2\sqrt{2}-\frac{\sqrt{2}}{15}\varepsilon^2+O(\varepsilon^3),

and is strictly below Amin⁡(K4)/2A_{\min}(K_4)/2 for 0<ε≤1/80000<\varepsilon\le1/8000. More precisely,

∣SεH−2+ε215−400ε3∣≤128ε4(0<ε≤1/100).\left|\frac{S_\varepsilon}{H}-2+\frac{\varepsilon^2}{15}-400\varepsilon^3\right|\le128\varepsilon^4\qquad(0<\varepsilon\le1/100).

The full coverage interval and the smaller sufficient saving interval remain distinct under this similarity.

References

  1. [1]Keith Ball. The plank problem for symmetric bodies. Inventiones Mathematicae, 104:535–543, 1991. Author manuscript: arXiv:math/9201218v1.DOI
  2. [2]Thøger Bang. A solution of the “plank problem”. Proceedings of the American Mathematical Society, 2(6):990–993, 1951.
  3. [3]Károly Bezdek. Tarski’s plank problem revisited. arXiv:0903.4637v1, 2009. Problem 3.1. Published in Bolyai Society Mathematical Studies 24 (2013), 45–64.arxiv.org/abs/0903.4637
  4. [4]Károly Bezdek and Muhammad A. Khan. The geometry of homothetic covering and illumination. arXiv:1602.06040v2, 2016. Version 2, 18 July 2016. Published in Discrete geometry and symmetry, Springer Proceedings in Mathematics & Statistics 234 (2018), 1–30.arxiv.org/abs/1602.06040
  5. [5]Károly Bezdek and Alexander E. Litvak. Covering convex bodies by cylinders and lattice points by flats. Journal of Geometric Analysis, 19(2):233–243, 2009. Author manuscript: arXiv:0903.4842v1.DOI
  6. [6]Horst Martini. Convex polytopes whose projection bodies and difference sets are polars. Discrete & Computational Geometry, 6:83–91, 1991.DOI
  7. [7]William Verreault. Plank theorems and their applications: A survey. Bulletin of the London Mathematical Society, 58(1):e70230, 2026. First published online 9 November 2025.arxiv.org/abs/2203.05540

Paper details

Contents