cs.LGOct 23, 2026

RadonBench: A Benchmark of Classical and Learning-Guided Search Techniques for Problems in Combinatorial Geometry

Authors: Xiaoan Xu, Ethan X. Fang, Jesus A. De Loera, Junwei Lu, Luze Xu

Organizations: Department of Electrical and Computer Engineering, Duke University. · Department of Biostatistics and Bioinformatics, Duke University. · Department of Mathematics, University of California, Davis. · Department of Biostatistics, Harvard T.H. Chan School of Public Health. · 

Abstract

Classical optimization and learning-guided search offer complementary approaches to open problems in discrete mathematics, but their outputs are difficult to compare when search spaces, validity checks, and stopping rules differ. We present RadonBench, a benchmark for constructing maximum-cardinality Radon-free lattice configurations: finite subsets of Zd\mathbb{Z}^d admitting no integral Radon partition. It remains open whether Radon-free configurations with 1111 points in Z3\mathbb{Z}^3 or 2121 points in Z4\mathbb{Z}^4 are possible. If they exist, they would improve the corresponding lower bounds on the integer Radon numbers. We compare exhaustive extension, greedy local search, PatternBoost, AlphaZero, and AlphaEvolve. All outputs use one success predicate: a convex target-sized configuration succeeds exactly when its number of integral Radon points is zero. PatternBoost and AlphaZero construct verified Radon-free configurations of sizes 1010 and 2020, one point short of the respective targets; AlphaZero's best convex target terminals have two integral Radon points in both dimensions. Exhaustive extension over fixed seeds and declared finite search boxes likewise reaches a minimum of two integral Radon points. Greedy local search and AlphaEvolve each produce target-sized near misses with one integral Radon point. No audited run reaches zero. These results establish reproducible near-miss baselines over explicitly scoped search domains.

Figures & tables

Appendix figures & tables7 assets

Supplementary material from the paper’s appendix.

Appendix

Explore similar work

Jun 24, 2026

Geometry-Aware MCTS for Extremal Problems in Combinatorial Geometry

We study certain extremal problems in combinatorial geometry that ask about configurations of points in an n×nn \times n grid that satisfy strict, global geometric constraints. Classical exact solvers suffer from combinatorial explosion for these types of problems, and standard reinforcement learning and transformer-based models struggle with the sparse reward "validity cliff" and quadratic token-consumption limits. To overcome these bottlenecks, we propose a Geometry-Aware Monte Carlo Tree Search (MCTS) framework. Our approach strictly enforces geometric constraints through incremental updates to the feasible action space. For constraints about collections of collinear points, like those that occur in the classic No-Three-in-Line problem (Max-N3IL), this mechanism reduces the constraint checking complexity from O(n3)O(n^3) to O(n2)O(n^2). To improve search efficiency, we exploit geometric symmetries in two ways: canonical pruning during node expansion to reduce the branching factor, and symmetric batch transitions to accelerate the discovery of promising configurations. We perform extensive experiments and establish new best-known computational results on five out of six of the problems that we considered. Notably, for Max-N3IL we find configurations of size roughly 1.8n1.8 n for grids of size 82≤n≤11982 \le n \le 119. For the Smallest Complete Set problem, we find configurations of size roughly 0.95n0.95 n, providing new upper bounds within the tested grids. This work establishes Geometry-Aware MCTS as a highly adaptable framework for discovering novel configurations in combinatorial geometry.
May 1, 2026

New Bounds for Zarankiewicz Numbers via Reinforced LLM Evolutionary Search

The Zarankiewicz number Z(m,n,s,t)\textbf{Z}(m, n, s, t) is the maximum number of edges in a bipartite graph Gm,nG_{m, n} such that there is no complete Ks,tK_{s, t} bipartite subgraph. We determine for the first time the exact values of three Zarankiewicz numbers: Z(11,21,3,3)=116\textbf{Z}(11, 21, 3, 3)=116, Z(11,22,3,3)=121\textbf{Z}(11, 22, 3, 3)=121, and Z(12,22,3,3)=132\textbf{Z}(12, 22, 3, 3)=132. We further establish lower bounds for 41 more Zarankiewicz numbers, including several that are within one edge of the best known upper bound, and we match the established value in four more closed cases. Our results are obtained using OpenEvolve, an open-source evolutionary algorithm based on Large Language Models (LLMs) that iteratively improves algorithms for generating mathematical constructions by optimizing a reward signal which we tailored for this specific problem. These findings provide new extremal graph constructions and demonstrate the potential of LLM-guided evolutionary search to contribute to mathematical research. In addition to presenting the resulting constructions, we report the generation algorithms produced, describe the relevant implementation details, and provide our computational costs. Our costs are remarkably low, at less than $30 for each Zarankiewicz parameter combination, showing that LLM-guided evolutionary search can be an inexpensive, reproducible, and accessible tool for discovering new combinatorial constructions.
Mar 10, 2026

Reinforced Generation of Combinatorial Structures: Ramsey Numbers

We present improved lower bounds for nine classical Ramsey numbers: R(3,13)\mathbf{R}(3, 13) is increased from 6060 to 6161, R(3,18)\mathbf{R}(3, 18) from 9999 to 100100, R(4,13)\mathbf{R}(4, 13) from 138138 to 139139, R(4,14)\mathbf{R}(4, 14) from 147147 to 148148, R(4,15)\mathbf{R}(4, 15) from 158158 to 159159, R(4,16)\mathbf{R}(4, 16) from 170170 to 174174, R(4,18)\mathbf{R}(4, 18) from 205205 to 209209, R(4,19)\mathbf{R}(4, 19) from 213213 to 219219, and R(4,20)\mathbf{R}(4, 20) from 234234 to 237237. These results were achieved using AlphaEvolve, an LLM-based code mutation agent. Beyond these new results, we successfully recovered lower bounds for all Ramsey numbers known to be exact, and matched the best known lower bounds across many other cases. These include bounds for which previous work does not detail the algorithms used. Virtually all known Ramsey lower bounds are derived computationally, with bespoke search algorithms each delivering a handful of results. AlphaEvolve is a single meta-algorithm yielding search algorithms for all of our results.