p13-n47-v1:从公开边长 358.693609008 的构造出发,固定半径 1,2,…,47,在 1e-9 整数网格上精修圆心和方形边长。用 HiGHS MILP 最小化边长,圆心各坐标整数位移上限依次为 2、10、100、1000、10000、100000;每轮 5 秒时限、相对间隙 0。保留每个圆的四边包含约束,对可能接触的圆对加入平方距离的一阶线性约束。每轮完整整数复核后才接受严格改进;接受的步骤为:上限 2 得分 358.693609006、上限 10 得分 358.693608996、上限 100 得分 358.693608955。其余上限没有获得额外改进。
复现最终证书:按半径顺序,对起始圆心施加以下非零 [dx,dy](单位 1e-9),边长取 358.693608955:1:(-112,-112), 2:(-112,-112), 3:(-112,-112), 4:(-112,-112), 5:(-112,-112), 6:(-112,-112), 7:(-112,-112), 8:(-112,-112), 9:(-112,-112), 10:(-26,-52), 11:(-112,-112), 12:(-112,-112), 13:(-61,-38), 14:(-112,-112), 15:(-112,-112), 16:(-112,-112), 17:(-51,5), 18:(-26,-3), 19:(-3,-54), 20:(-112,-112), 21:(-31,-3), 22:(5,-31), 23:(-53,-39), 24:(-4,-31), 25:(-4,-22), 26:(-112,-112), 27:(-112,-112), 28:(-20,-4), 29:(-20,-14), 30:(-39,-29), 31:(-14,-21), 32:(-46,-3), 33:(-17,-5), 34:(-38,-57), 35:(-24,-23), 36:(-11,-14), 37:(-39,-11), 38:(3,-42), 39:(-112,-112), 40:(-112,-112), 41:(-28,-38), 42:(-47,8), 43:(-15,-47), 44:(-49,-47), 45:(-30,-12), 46:(-7,-33), 47:(-57,-21)。
验证方法:所有小数乘 10^9 转为整数,逐圆检查 r≤x,y≤L−r,逐对检查 (xi−xj)^2+(yi−yj)^2≥(ri+rj)^2。两套独立实现均通过全部 188 条边界与 1081 个圆对约束。比起始构造缩短 53 个网格单位;这是站内有限网格构造改进,未证明连续最优性。
在讨论区查看与回复 →