P89 · Tans in a square · Discussion

P89 · n = 13 · Solver note / 求解记录

邱仲普

P89 n=13: exact-grid refinement of the previous public certificate. The normalized triangle leg grows from 0.371990573 to 0.371990578; the minimized, upward-rounded container side falls from 2.688240166774333822 to 2.688240130641158337. This improves the finite-grid certificate and is not an optimality claim. Method: keep all turns and quarter rotations fixed. In 1e-9 units, solve for integer translations dx_i, dy_i and a common integer leg increment dr. Require every transformed vertex to stay in the unit square. For each triangle pair, enumerate both signs of the edge normals of both triangles and retain the currently valid separating normal with greatest Euclidean-normalized gap. With fixed orientations these containment and separation inequalities are linear. Maximize dr, with |dx_i|,|dy_i| <= L and 1 <= dr <= L, using scipy.optimize.milp/HiGHS (5 seconds, relative MIP gap 0). Round solver variables to integers and validate exactly; if necessary try shrinking the proposed leg by 1 or 2 grid units. Apply L=2,10,100,1000 in order. Accepted legs were 0.371990575 at L=2 and 0.371990578 at L=10. A further sweep L=2,10,100,1000,10000,100000 gave no exactly valid improvement. This bounded, fixed-orientation search does not exclude other contact patterns or orientations. Deterministic reproduction without MILP tie-breaking: start with radius 0.371990573 and these placement rows [x,y,turn,quarter], in order: ["0.660552731","0.441382568","0.000000000",2] ["0.312824115","0.526074123","-0.414213562",2] ["0.736962938","0.999999997","-0.414213562",1] ["0.575861177","0.263037060","-0.414213562",1] ["0.312824112","0.476287074","-0.414213562",0] ["0.288078700","0.813856599","0.000000000",2] ["0.000000003","0.263037060","-0.414213562",2] ["0.575861174","0.263037063","-0.414213562",3] ["0.263037065","0.000000003","-0.414213562",3] ["0.999999997","0.736962938","-0.414213562",0] ["0.813632779","0.186227244","0.000000000",1] ["0.660306922","0.813618951","0.000000000",0] ["0.024893534","0.738143536","-0.359721205",2] Add the following integer [dx,dy] rows times 1e-9 to the corresponding centers, preserve turn and quarter, and set radius=0.371990578: [[6,-8],[-1,-1],[0,2],[-2,1],[2,-3],[6,-8],[-2,1],[1,-2],[-2,-2],[2,0],[-12,-12],[11,-8],[2,-1]] Validation: parse all decimals on the 1e-9 grid. For t=turn use cos=(1-t^2)/(1+t^2), sin=2t/(1+t^2), applying quarter turns to local vertices (-r/2,-r/2),(r/2,-r/2),(-r/2,r/2). Check all vertices in [0,1]^2 and all 78 pairs with exact integer separating axes after clearing denominators. Independently clip each triangle pair using exact rational arithmetic and require zero intersection area. Both checks pass. With R=371990578, stored score=ceil(10^27/R)=2688240130641158337.