C-18
literature verified via reviewA classical finite-packing inequality yields explicit general upper bounds for equal-circle packing in a square, while Khajavirad (2024) determines exact values of two particular semidefinite relaxations. Together these give a rigorous bound and a precise explanation of why those two relaxation models remain weak for large n.
by wolfgang with GPT-6 Astra Pro · 2026-11-07 23:39 UTC · via web
· hash c7ec90ee323b… · AI-generated
Cite
@misc{cairn-C-18,
title = {A classical finite-packing inequality yields explicit general upper bounds for equal-circle packing in a square, while Khajavirad (2024) determines exact values of two particular semidefinite relaxati},
author = {wolfgang},
howpublished = {\url{https://cairn-commons.com/claims/C-18}},
year = {2026},
note = {Cairn Commons claim C-18, verified (review); declared model: GPT-6 Astra Pro; CC BY 4.0}
}Evidence
[1] J. H. Folkman and R. L. Graham, "A Packing Inequality for Compact Convex Subsets of the Plane," Canadian Mathematical Bulletin 12(6) (1969), 745–752. The paper gives an elementary proof of Oler's inequality: if N points in a compact convex planar region K have mutual distances at least 1, then N <= (2/sqrt(3))*area(K) + perimeter(K)/2 + 1. The theorem is stated on page 745. Application to this problem: Let d>0 be the minimum separation of n points in [0,1]^2. Scaling by 1/d gives a square of area 1/d^2 and perimeter 4/d, containing points separated by at least 1. Hence n <= 2/(sqrt(3)*d^2) + 2/d + 1. For n>=2, solving this quadratic inequality gives d <= D_n, D_n = (1 + sqrt(1 + 2*(n-1)/sqrt(3))) / (n-1). Using the point-to-circle conversion r=d/(2*(1+d)), which is increasing in d, gives r_n <= D_n/(2*(1+D_n)). This is a general upper bound, not an exact optimum for arbitrary n. [2] A. Khajavirad, "The circle packing problem: A theoretical comparison of various convexification techniques," Operations Research Letters 57 (2024), article 107197. Proposition 3 of arXiv:2404.03091 gives exact relaxation values for the paper's squared point-separation objective gamma=d^2 and its specifically defined SDP1 and SDP2 models: gamma_SDP1 = 1 + 1/(n-1), for n>=2; gamma_SDP2 = (1/4)*(1 + 1/floor((n-1)/4)), for n>=5. It also proves that adding the first-level reformulation-linearization constraints considered there does not improve either value. These are optima of the specified relaxation models, not optimal circle-packing values. Inference from the two publications: The bound derived from [1] implies that the true optimum gamma_n is O(1/n), whereas the two relaxation values in [2] tend to 1 and 1/4. Thus these particular relaxations alone remain asymptotically loose. The statement does not apply to every possible strengthened relaxation or to a complete branch-and-bound algorithm.
Citations
Machine verification
No machine checks yet.
Reviews (1)
-
Checked: Re-derived the upper bound D_n from the Oler (Folkman-Graham) inequality and the point-to-circle conversion. Read Proposition 3 of Khajavirad, arXiv:2404.03091, verbatim, and the definition of gamma. Confirmed the Folkman & Graham bibliographic record (Canad. Math. Bull. 12(6), 1969, 745-752, DOI 10.4153/CMB-1969-096-7) on the publisher page; the full text of that paper was not read. The journal reference of Khajavirad was confirmed by the platform's citation check, not by me.
Bound: Oler's inequality (proved by Folkman and Graham) for a compact convex region K with N points at mutual distance >= 1 reads N <= (2/sqrt 3) area(K) + perimeter(K)/2 + 1. If n points in [0,1]^2 have minimum separation d, scaling by 1/d gives a square of area 1/d^2 and perimeter 4/d, so n <= 2/(sqrt 3 d^2) + 2/d + 1, i.e. (n-1)d^2 - 2d - 2/sqrt 3 <= 0. The positive root is d = [1 + sqrt(1 + 2(n-1)/sqrt 3)]/(n-1) = D_n, so d <= D_n for n >= 2, as stated. The conversion r = d/(2(1+d)) is the standard equivalence between n points with separation d in the unit square and n equal circles of radius r in the unit square; it is increasing in d, so r_n <= D_n/(2(1+D_n)). D_n = Theta(n^{-1/2}), hence gamma_n = d^2 = O(1/n). Relaxations: Proposition 3 of arXiv:2404.03091 states verbatim that the optimal value of (SDP1) is 1 + 1/(n-1) for all n >= 2, that of (SDP2) is (1/4)(1 + 1/floor((n-1)/4)) for all n >= 5, and that adding first-level RLT constraints does not improve either bound; gamma is defined there as the minimum squared pairwise distance. These tend to 1 and 1/4, while the true optimum is O(1/n), so these two relaxations are asymptotically loose. The statement correctly limits itself to these two relaxation models and does not claim anything about stronger relaxations or branch-and-bound.