n = 39 open ★=

13320≤s(39)≤13621444166137320000000000000

The best packing known for 39 squares, side 6.81072208…, Thomas Schadt 2026
6.6506.811
678
6.2457.245
nn+1

Proven

6.650000≤s(39)≤6.810723

  • new result
  • exact

Citation record n-039

lowerwand125 after Tokoharu, Levy et al. 2026, GitHub (confirmed T-110)

upperSchadt 2026, Squares in Squares (confirmed T-101)

Open

  • optimality

Bounds

Best known packing

6.81072208…

6.81072208306864
Found by
Thomas Schadt 2026
Construction
annealing
Minimal polynomial, degree 5
9s5−171s4+999s3−1959s2+1636s+166=0
Source
[Kingbird]
Evidence
E-kingbird-upper-register
Verified upper bound

13621444166137320000000000000

6.81072208306865

The reported value, verified here.

Evidence
E-evand-exact-ceilings-2026-10-05-exact-replay, E-evand-exact-ceilings-2026-10-05-source-replay
Reported lower bound

13320

Proved by
wand125 2026
Kind
counting
Scope
Unrestricted unit-square packing with independent rotations and disjoint interiors.
Note
wand125's square-packing-bounds (6 October 2026) reports s(39)≥133/20 from a density of 600 uniform rectangles of total mass 3899999/100000, on a net the certificate declares: core side 999/1000 and 416 half-angle tangents of step 1/1001, accepted there at every direction by sqverify-proof-net, the source's copy of this repository's sqverify_fast changed to read a declared net (a check2 bundle, with no C++ record). It is above the 1327/200 the record reported (T-068). sqverify-fast, this repository's clean-room measure verifier, decided it here at all 416 directions on 6 October 2026.
Source
[wand125 mixed bounds check2 2026-10-06]
Evidence
E-n039-wand125-mixed-665-report
Verified lower bound

13320

The reported value, verified here.

Evidence
E-n039-wand125-mixed-665-sqverify-fast-replay
Gap

0.16072208…

Verified upper minus verified lower.

Results in the register

Verification

upper: replayed here; lower: replayed here

—

Rigidity

not rigid, numerically checked, numerical multiprecision

Evidence: E-translation-escape-not-rigid

Scope

Square 9 of the retained witness (witness id 10) translates 0.380393 along (1, 0) with the packing still valid, so the configuration admits a non-trivial feasible motion; 13 of its 39 squares do. Every constraint is exactly affine in the slide parameter, so the arithmetic carries no linearization error, but the coordinates are the witness's own finite-precision transcription: this settles the retained configuration, not the true optimum. Rigidity and optimality are independent, and this bears only on the former.

s(39) — open

Exact certificate, 2026-10-06. Evan Daniel’s square-packing published on 5 October 2026 an exact rational certificate of this packing, decided here over ℚ by two exact checkers that share no code with each other or with the source, and by the source’s own two run here: it proves s(39)≤6.81072208306865, the verified upper bound (T-101). Until then the verified upper bound was the trivial grid bound 7.

External intake, 2026-10-06. wand125’s check2 source reports s(39)≥133/20=6.65 (T-110), from a density of 600 rectangles of total mass 3899999/100000<39, on a net the certificate declares: core side 999/1000 and 416 half-angle tangents of step 1/1001. Since 999/1000·(1+1/1001)<1, every unit square contains a core at a net angle strictly in its interior. The source ships no record of its C++ checker for it: it accepts it at every direction with its copy of this repository’s sqverify_fast, changed to read a declared net. It is above the reported 1327/200 (T-068) by 0.015. This repository’s clean-room verifier sqverify-fast decided it here on 6 October 2026 at all 416 directions of its net, and refused two mutants scaled below coverage one: confirmed, re-implemented sharing the producer’s components (the source’s check is a copy of the same crate), so it is also the verified lower bound. wand125’s README says parts of the work were produced with AI assistance under human direction.

External intake, 2026-10-01. wand125’s rectangle-density source reports s(39)≥1327/200=6.635, since superseded (T-110), with total mass 3899/100=38.99<39, accepted by Tokoharu’s unchanged interval checker. The complete 201-direction coverage replay here accepted it again, after this repository’s exact audit checked that the regenerated checker input is the published one and checked the mass and net premises, so it was also verified until 2026-10-06, when the certificate above superseded it in both lanes. wand125’s README says parts of the work were produced with AI assistance under human direction.

External intake, 2026-09-28. wand125’s rectangle-density source reports a direct 663/100=6.63 certificate for this case, whose reported bound the 2026-10-01 intake above raises, with total mass 3899/100=38.99<39, accepted by Tokoharu’s unchanged interval checker. This repository’s exact audit checks that the regenerated checker input is the published one, and checks the mass and net premises; the complete coverage replay has not yet run here, so the verified lower bound is unchanged. wand125’s README says parts of the work were produced with AI assistance under human direction.

External intake, 2026-09-27. wand125’s rectangle-density source reports a direct 331/50=6.62 certificate for this case, whose reported bound the 2026-09-28 intake above raises, with total mass 3899/100=38.99<39, accepted by Tokoharu’s unchanged interval checker. This repository’s exact audit checks that the regenerated checker input is the published one, and checks the mass and net premises; the complete coverage replay has not yet run here, so the verified lower bound is unchanged.

External intake, 2026-09-22. wand125’s retained point-certificate source reports s(39)≥13/2=6.5, which was the verified lower bound until 2026-10-02. The complete exact replay verifies the retained certificate.

Open. The best known packing gives s(39)≤6.81072209, and the verified lower bound is s(39)≥133/20=6.65, from wand125’s check2 certificate on a declared net (T-110, 2026-10-06, V3/C3), which superseded its rectangle-density certificate of 1 October at 1327/200=6.635 (T-074), leaving a gap of about 0.1607.

The exact certificate

Evan Daniel’s square-packing published on 5 October 2026 exact rational certificates of the packings this register lists as best known, solved from its own witnesses, and offered those that do not lower a printed side as a replay of the existing bounds, asking for nothing to be registered from them. The certificate for this count holds the same 39 squares, each a rational centre and a rational t=tan(θ/2), in a square of side 6.8107220830686488518…, 8.9×10−15 above the side the Kingbird catalogue prints, 6.81072208306864. Square for square, its pose lies within 1.5×10−4 of the binary64 pose the atlas pictures for this count. 16 squares move by more than 1×10−8, all of them squares the source lists as carrying no force; every other square’s centre rounds to the witness’s. Evan Daniel’s solver moves the binary64 pose to a nearby exact KKT point of the problem of minimizing the side under non-overlap, computed at 80 digits, and rounds it outward to rationals.

This repository decides the certificate exactly. Converted without rounding, every pair and every wall is decided over ℚ twice, by sqpack’s exact separating-axis test and by an independent checker that shares no code with it, and the source’s own two checkers, run here as retained, accept it as well (receipts). That proves s(39)≤6.81072208306865, the verified upper bound: the certificate’s side rounded up at the fourteen decimals the catalogue prints, as the record writes any certificate of a printed side. It says nothing about optimality. On jlevy/squares#375 its author wrote that the solver, its checkers and the batch “were written with Claude (Anthropic) as a coding and research agent, directed and reviewed by me.”

The verified upper bound and the printed side now agree to one unit of the printed side’s last place, the precision at which the record compares them. The printed side itself is not certified here: the certificate’s side lies above it.

The packing

Found by Thomas Schadt in 2026, via simulated annealing. Its side length is algebraic of degree 5 over ℚ. Its Galois group over ℚ is S5, computed here from the recorded minimal polynomial with sympy.galois_group, and S5 is not solvable. The side length is therefore not expressible in radicals, so the empty exact_form is a permanent property of this case rather than a transcription gap.

The lower bound

The earlier record from [Friedman DS7] gives the lower-bound expression 103/37+22+113/37 for s(39) (approximately 6.350603018684). Friedman’s DS7 survey, Theorem 9, k=6, reports this bound at n=37; reference [8] is Green’s private communication (2000). The source proof has not been recovered. The unavoidable-set argument DS7’s Figure 34 illustrates does not prove it: at k=6 that point pattern leaves a unit square empty (review). Inherited at n=39 by monotonicity. That report remains as historical evidence; the verified point certificate gave the stronger bound 13/2, and the replayed rectangle-density certificate now gives 663/100. The source audit compares the exact theorem expressions separately from opaque table decimals.

The exact replay verifies all five covering-certificate conditions over the complete rational direction net. A second interval implementation or generalized native importer would add confirmation and reuse, but it is not a missing premise of this fixed bound.

Corrected 2 October 2026: Nagamochi’s Lemma 1 is false (Karakuş 2026), so his closed form below is now a reported bound whose published proof is incomplete (review). This register had recorded that proof as verified, its own error, logged as defect D-516. Nagamochi’s earlier general closed form remains part of the evidence history and applies to every N≥4:

s(N)≥min{⌈N⌉,N−2⌊N⌋+1+1}

Verification Code

The programs behind this case’s verified bounds, by their evidence. The code column says how the code that ran stands to the code its producer used. VERIFIERS.md says what each program is and whose it is.

bound evidence run code programs
verified lower E-n039-wand125-mixed-665-sqverify-fast-replay replayed here shared components V-sqverify-fast (first-party)
verified upper E-evand-exact-ceilings-2026-10-05-exact-replay replayed here independent V-sqpack-verify, V-check-rational-witness-independent (first-party); V-evand-exact-certificates (first-party, premises)
verified upper E-evand-exact-ceilings-2026-10-05-source-replay replayed here producer’s code V-evand-verify-cert-py, V-evand-verify-cert2-py (external); V-evand-exact-certificates (first-party, premises)