n = 41 open ★=

27140≤s(41)≤692669309446881100000000000000

The best packing known for 41 squares, side 6.92669309…, Thomas Schadt 2025
6.7756.927
678
6.4037.403
nn+1

Proven

6.775000≤s(41)≤6.926694

  • new result
  • exact

Citation record n-041

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

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

Open

  • optimality

Bounds

Best known packing

6.92669309…

6.92669309446880
Found by
Thomas Schadt 2025
Construction
annealing
Source
[Kingbird]
Evidence
E-kingbird-upper-register
Minimal polynomial, degree 42

144s42−33248s41+3740531s40−273229120s39+14568177368s38−604345329616s37+20303247278518s36−567715628580628s35+13476130642163772s34−275622556171657148s33+4913118839607229315s32−77021965442580593792s31+1069597207525632250760s30−13234280063158548374864s29+146588403144234109714492s28−1459056537531761947694412s27+13090219490685085049164304s26−106115059640167069135194108s25+778709778173545540562913112s24−5180160724110239826615572336s23+31267085211757278545052994144s22−171331300125735569491450805184s21+852412555299622931388971786184s20−3849639590878015114188275848896s19+15771592794879254264477226832440s18−58556476540137831140983424890112s17+196742347065286547712609208667628s16−597075318553361988026330293830592s15+1632807555219691500831155576662224s14−4011703707846363075271797958908992s13+8823126218415607049547609565313808s12−17292499393880618294971765830702496s11+30033784585675389426408059928238624s10−45904080196423917967770013765165584s9+61200148145342575539111090507985440s8−70370773645277985938951580858638528s7+68756165329665893470887878785349152s6−55949600498940958320310297578555360s5+36867848125763978702951438802849616s4−18873332426700882570902855047275200s3+7023595398126017089078028623797120s2−1682258751137636203725622554061120s+192930128676231207430057837613968=0

Verified upper bound

692669309446881100000000000000

6.92669309446881

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

27140

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(41)≥271/40 from a density of 376 uniform rectangles of total mass 4099999/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 169/25 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-n041-wand125-mixed-6775-report
Verified lower bound

27140

The reported value, verified here.

Evidence
E-n041-wand125-mixed-6775-sqverify-fast-replay
Gap

0.15169309…

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 6 of the retained witness (witness id 7) translates 0.085821 along (0.707107, -0.707107) with the packing still valid, so the configuration admits a non-trivial feasible motion; 4 of its 41 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(41) — 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(41)≤6.92669309446881, 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(41)≥271/40=6.775 (T-111), from a density of 376 rectangles of total mass 4099999/100000<41, 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 169/25 (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(41)≥169/25=6.76, since superseded (T-111), with total mass 4099/100=40.99<41, 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 1351/200=6.755 certificate for this case, whose reported bound the 2026-10-01 intake above raises, with total mass 4099/100=40.99<41, 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 1349/200=6.745 certificate for this case, whose reported bound the 2026-09-28 intake above raises, with total mass 4099/100=40.99<41, 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 source reports s(40)≥13/2=6.5 and identifies the monotone consequence at n=41. The complete exact n40 replay verifies the source premise, so monotonicity verifies the n41 bound.

Open. The best known packing gives s(41)≤6.92669310, and the verified lower bound is s(41)≥271/40=6.775, from wand125’s check2 certificate on a declared net (T-111, 2026-10-06, V3/C3), which superseded its rectangle-density certificate of 1 October at 169/25=6.76 (T-074), leaving a gap of about 0.1517.

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 41 squares, each a rational centre and a rational t=tan(θ/2), in a square of side 6.9266930944688022317…, 2.2×10−15 above the side the Kingbird catalogue prints, 6.92669309446880. Square for square, its pose lies within 1.5×10−4 of the binary64 pose the atlas pictures for this count. 2 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(41)≤6.92669309446881, 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 2025, via simulated annealing. Its side length is algebraic of degree 42 over ℚ.

The lower bound

The verified lower bound is 271/40=6.775, from wand125’s check2 certificate of 6 October on a declared net (T-111), decided here by sqverify-fast at every direction of its net. Before it, 169/25=6.76, from its own n41 rectangle-density certificate of 1 October, replayed in full here (T-074); its certificate of 28 September, 1351/200=6.755 (T-070), held the field earlier on 2 October. Until then it was 13/2: the complete n40 exact replay verifies that bound, and deleting one square from any hypothetical packing proves it for n41 too. A second interval implementation or generalized native importer would add confirmation and reuse, but it is not a missing premise of either 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-n041-wand125-mixed-6775-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)