REVIEW 4 major objections 7 minor 20 references
A learned edge filter can throw away 95–99% of TSP edges and still leave a tour within 1% of optimal.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-14 17:24 UTC pith:4N67NNID
load-bearing objection Solid incremental sparsification pipeline with real empirical gains; threshold opacity and Christofides safety net are the main soft spots, not a collapse of the result. the 4 major comments →
GES-TSP: Graph Edge Sparsification for TSP
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
A GAT trained on Delaunay-prefiltered Euclidean TSP graphs, using four edge features (distance, KNN membership, iterative-MST weight, and local quality ratio) plus a Christofides safety net, produces sparse subgraphs that prune 95% of edges on 100-city MATILDA instances and over 99% on large TSPLIB instances while keeping the optimality gap of the subsequent exact solution under 1%.
What carries the argument
GES: Delaunay coarse graph → four-dimensional edge features → GATv2 edge scoring → threshold filter plus Christofides edges → exact solver on the residual sparse graph.
Load-bearing premise
That a model trained only on a fraction of 100-city MATILDA instances, with a fixed score threshold and Christofides augmentation, will keep a near-optimal tour’s edges on much larger and differently distributed TSPLIB cities.
What would settle it
Take any TSPLIB instance larger than pr2392 for which an optimal tour is known; run the published GES pipeline with the same threshold and check whether the sparse graph still contains a tour whose length is within 1% of that optimum.
If this is right
- Exact MIP solvers can be applied to TSP instances whose complete graphs were previously too large for memory.
- The same learned candidate set can be handed to any heuristic (LKH, Christofides, etc.) to accelerate local search.
- Because the model generalizes from 100-city training data to thousands of cities, a single trained sparsifier can serve as a reusable front-end for many Euclidean TSP solvers.
- Ablation evidence shows that dropping the Delaunay pre-filter or any of the four edge features measurably worsens either pruning rate or optimality gap, so each component is necessary for the claimed performance.
Where Pith is reading between the lines
- The same feature-plus-attention recipe may transfer to other geometric combinatorial problems (vehicle routing, Steiner trees) where short edges dominate high-quality solutions.
- If the score threshold can be made adaptive to graph size or density, the method could avoid the residual infeasibility that currently forces Christofides edges on hard instances.
- Because the model never sees the final tour during inference, it can be viewed as a learned, instance-specific replacement for classic candidate-set heuristics used inside LKH-style solvers.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes GES-TSP, a learning-based edge sparsification pipeline for Euclidean TSP. It first builds a Delaunay triangulation as a coarse candidate graph, constructs four edge features (Euclidean distance, binary KNN membership, a normalized iterative-MST score, and a local quality ratio Q_ij), scores edges with a GATv2 model trained as supervised binary classification against optimal tours, retains edges above a score threshold, and augments the sparse graph with Christofides edges to guarantee a feasible tour. The sparsified instance is then solved with SCIP. On MATILDA (100-node instances), the method reports pruning rates around 94% with optimality gaps typically well below 1%, outperforming SGN and Fitzpatrick baselines on both gap and pruning. On TSPLIB instances up to pr2392, it reports pruning rates above 99% on large instances with gaps remaining under 1%. Ablations address the coarse filter, individual features, and the Christofides augmentation.
Significance. If the empirical claims hold under fully specified and fair evaluation, the work is a useful practical contribution to exact/high-quality TSP solving: it shows that a relatively light GAT, combined with classical geometric structure and a feasibility safety net, can discard the vast majority of edges while preserving near-optimal tours, with apparent transfer from 100-node training instances to larger TSPLIB graphs. Strengths include a clear modular pipeline, comparison against two relevant sparsification baselines, feature and coarse-filter ablations (Table 3, Fig. 8), a feasibility ablation for Christofides (Appendix Table 7), and consistent reporting of both pruning rate and optimality gap. The contribution is primarily empirical and engineering-oriented rather than theoretical; its value for the community depends on reproducibility of the sparsification rule and on clarifying how much of the reported quality comes from learned scores versus Delaunay structure and Christofides edges.
major comments (4)
- §4.4 (Sparse Graph) never states the numerical score threshold used to retain edges, nor whether it is fixed across instances, chosen by validation, or adjusted per graph size. This threshold is load-bearing for every pruning-rate and optimality-gap number in Tables 1–2, Figs. 4–5, and Appendix Tables 4–5. Without the value (or a selection procedure), the central results are not reproducible and the comparison of pruning rates against SGN (top-k) and Fitzpatrick is not fully interpretable.
- §4.4 and Appendix Table 7: without Christofides augmentation, the pipeline already yields infeasible graphs on 8/127 CLKhard and 8/127 LKCChard instances—the hard distributions used in training. The TSPLIB “strong generalization” claim (pruning >99%, gap <1% after training only on ~100-node MATILDA subsets) therefore needs a quantitative breakdown of how many retained tour edges come from GAT scores versus forced Christofides edges, and of the gap when Christofides edges are excluded or replaced by a weaker connectivity fix. Otherwise it is unclear whether the headline generalization is driven by the learned sparsifier or by the 3/2-approximation safety net.
- §5.1 comparative study: the paper does not fully specify how the SGN and Fitzpatrick baselines were configured for a fair head-to-head (e.g., k for SGN, decision threshold or feature set for Fitzpatrick, whether all three methods feed the same SCIP formulation and time limits, and whether GNN inference time is included in speedups in Fig. 6). Because pruning rate and gap trade off, unequal operating points could inflate GES’s advantage. Please report the exact baseline settings and, ideally, a gap-vs-pruning Pareto comparison.
- §5 training setup and free parameters: the model is trained on one-third of CLKhard and LKCChard (100 nodes) with a fixed positive class weight of 30, GATv2 hyperparameters, and an unspecified KNN neighborhood size k in the edge features. The generalization premise to TSPLIB instances up to ~2400 nodes is central to the abstract and §6. A sensitivity analysis of threshold, k, and class weight—and at least one experiment training or validating on a different size/distribution—would substantially strengthen the claim that the GAT, not only Delaunay+Christofides, transfers.
minor comments (7)
- Abstract and title block: missing space after period (“benchmark.In some”); affiliation typo “Lanzhou Universtiy” (twice).
- §3.2 Eq. (1)–(2): notation mixes h_v and h^{(l)}_v; edge feature r_vu is introduced without defining how the four-dimensional e_ij of §4.2 is mapped into the attention input.
- §4.2 Eq. (4): Q_ij uses min_k d_ik; clarify whether the min is over neighbors in the Delaunay graph or over all other nodes, and state the KNN parameter k used in Eq. (3).
- §5: “we consider instances ranging from ch130 to pr2392, since optimal solution files are not available…” is a run-on sentence; also clarify that MATILDA categories each have 190 instances of size 100.
- Figures 4–6 and 9: axis labels, units, and whether times include model inference vs. SCIP only should be stated in captions for readability.
- Related work: end-to-end neural TSP methods are cited but the positioning relative to other learning-to-prune / candidate-set methods could briefly note differences in supervision and solver coupling.
- Appendix Table 6 vs. main Fig. 8: feature-ablation gaps on TSPLIB are useful; ensure the main text points readers to the full table for exact numbers.
Circularity Check
No circularity: supervised edge classification trained on optimal tours and evaluated via independent external metrics (pruning rate, optimality gap) on held-out and out-of-distribution instances.
full rationale
The paper's pipeline is standard supervised learning for sparsification: Delaunay coarse graph + hand-crafted geometric features (distance, KNN, Q, normalized IMST) feed a GAT that is trained by binary cross-entropy against known optimal-tour edges on a MATILDA training split; a score threshold plus Christofides edges then produce a sparse graph that is solved by SCIP. All reported quantities (optimality gap to true optima, pruning rate relative to the complete graph, speed-up) are computed on held-out MATILDA categories and on external TSPLIB instances never seen during training. None of the six circularity patterns applies: there is no self-definitional equation, no fitted parameter re-labeled as a prediction of a dependent quantity, no load-bearing self-citation, no uniqueness theorem imported from the authors, no ansatz smuggled via self-citation, and no renaming of a known empirical pattern. Geometric facts (Delaunay, Christofides) are classical background, not circular inputs. The method is therefore self-contained against external benchmarks; any concerns about threshold choice or the Christofides safety net are generalization/correctness issues, not circularity.
Axiom & Free-Parameter Ledger
free parameters (5)
- edge score threshold for sparsification
- BCE positive class weight
- GAT architecture and training hyperparameters
- KNN neighborhood size k
- IMST iteration depth R=ceil(log2 n) and normalization
axioms (5)
- domain assumption Delaunay triangulation retains many edges that appear in high-quality Euclidean TSP tours and is a sufficient O(n) coarse candidate set.
- standard math Christofides tour edges guarantee the sparsified graph contains at least one Hamiltonian cycle.
- domain assumption Optimal-tour edge labels on training instances are a valid supervision signal for predicting useful candidate edges on unseen instances.
- domain assumption Euclidean distances and the four hand-crafted edge features (d, KNN, IMST, Q) are informative enough for GAT message passing to rank near-optimal edges.
- ad hoc to paper SCIP exact solutions on the sparsified graph are comparable across methods for reporting optimality gaps.
invented entities (2)
-
GES sparsification pipeline (Delaunay → GAT edge scores → threshold + Christofides)
no independent evidence
-
Simplified local Q-value edge feature Q_ij = (1+d_ij)/(1+min_k d_ik)
no independent evidence
read the original abstract
Solving large-scale instances of the Traveling Salesman Problem (TSP) exactly is computationally expensive. Researchers often employ graph sparsification methods to improve computational efficiency. Traditional sparsification methods typically rely on fixed heuristics and fail to fully exploit instance-specific structural information. In this paper, we propose Graph Edge Sparsification (GES), a learning-based sparsification approach for Euclidean TSP. By incorporating geometric structural information and combinatorial optimization technology, our proposed method adaptively generates a sparsification graph for different instances, significantly reducing the graph size and accelerating the solving process. Experimental results demonstrate that our sparsification method can prune up to 95% of edges on the MATILDA dataset, while keeping the solution gap within 1% of the optimal value. Moreover, our approach exhibits strong generalization capability on the TSPLIB benchmark.In some large-scale instances, the pruning rate exceeds 99%, while the optimality gap remains below 1%.
Figures
Reference graph
Works this paper leans on
-
[1]
The traveling salesman problem
Michael Jünger, Gerhard Reinelt, and Giovanni Rinaldi. The traveling salesman problem. Handbooks in operations research and management science, 7:225–330, 1995
1995
-
[2]
Spectral sparsification of graphs.SIAM Journal on Computing, 40(4):981–1025, 2011
Daniel A Spielman and Shang-Hua Teng. Spectral sparsification of graphs.SIAM Journal on Computing, 40(4):981–1025, 2011
2011
-
[3]
Springer, 2008
Mark De Berg, Otfried Cheong, Marc Van Kreveld, and Mark Overmars.Computational geometry: algorithms and applications. Springer, 2008
2008
-
[4]
Xiangping Xu, Jun Li, and MengChu Zhou. Delaunay-triangulation-based variable neighbor- hood search to solve large-scale general colored traveling salesman problems.IEEE Transactions on Intelligent Transportation Systems, 22(3):1583–1593, 2020
2020
-
[5]
The transformer network for the traveling salesman problem.arXiv preprint arXiv:2103.03012, 2021
Xavier Bresson and Thomas Laurent. The transformer network for the traveling salesman problem.arXiv preprint arXiv:2103.03012, 2021
Pith/arXiv arXiv 2021
-
[6]
Chaitanya K Joshi, Thomas Laurent, and Xavier Bresson. An efficient graph convolutional network technique for the travelling salesman problem.arXiv preprint arXiv:1906.01227, 2019
Pith/arXiv arXiv 1906
-
[7]
Pointer networks.Advances in neural information processing systems, 28, 2015
Oriol Vinyals, Meire Fortunato, and Navdeep Jaitly. Pointer networks.Advances in neural information processing systems, 28, 2015
2015
-
[8]
Irwan Bello, Hieu Pham, Quoc V Le, Mohammad Norouzi, and Samy Bengio. Neural combina- torial optimization with reinforcement learning.arXiv preprint arXiv:1611.09940, 2016
Pith/arXiv arXiv 2016
-
[9]
Pomo: Policy optimization with multiple optima for reinforcement learning.Advances in neural information processing systems, 33:21188–21198, 2020
Yeong-Dae Kwon, Jinho Choo, Byoungjip Kim, Iljoo Yoon, Youngjune Gwon, and Seungjai Min. Pomo: Policy optimization with multiple optima for reinforcement learning.Advances in neural information processing systems, 33:21188–21198, 2020
2020
-
[10]
Learning to sparsify travelling salesman problem instances
James Fitzpatrick, Deepak Ajwani, and Paula Carroll. Learning to sparsify travelling salesman problem instances. InInternational Conference on Integration of Constraint Programming, Artificial Intelligence, and Operations Research, pages 410–426. Springer, 2021
2021
-
[11]
Liang Xin, Wen Song, Zhiguang Cao, and Jie Zhang. Neurolkh: Combining deep learning model with lin-kernighan-helsgaun heuristic for solving the traveling salesman problem.Advances in Neural Information Processing Systems, 34:7472–7483, 2021
2021
-
[12]
Combhelper: a neural approach to reduce search space for graph combinatorial problems
Hao Tian, Sourav Medya, and Wei Ye. Combhelper: a neural approach to reduce search space for graph combinatorial problems. InProceedings of the AAAI Conference on Artificial Intelligence, volume 38, pages 20812–20820, 2024
2024
-
[13]
A comprehensive survey on graph reduction: Sparsification, coarsening, and condensation
Mohammad Hashemi, Shengbo Gong, Juntong Ni, Wenqi Fan, B Aditya Prakash, and Wei Jin. A comprehensive survey on graph reduction: Sparsification, coarsening, and condensation. 2024
2024
-
[14]
The graph neural network model.IEEE transactions on neural networks, 20(1):61–80, 2008
Franco Scarselli, Marco Gori, Ah Chung Tsoi, Markus Hagenbuchner, and Gabriele Monfardini. The graph neural network model.IEEE transactions on neural networks, 20(1):61–80, 2008
2008
-
[15]
Open graph benchmark: Datasets for machine learning on graphs
Weihua Hu, Matthias Fey, Marinka Zitnik, Yuxiao Dong, Hongyu Ren, Bowen Liu, Michele Catasta, and Jure Leskovec. Open graph benchmark: Datasets for machine learning on graphs. Advances in neural information processing systems, 33:22118–22133, 2020
2020
-
[16]
Graph attention networks.stat, 1050:4, 2018
Petar Velickovic, Guillem Cucurull, Arantxa Casanova, Adriana Romero, Pietro Lio, and Yoshua Bengio. Graph attention networks.stat, 1050:4, 2018
2018
-
[17]
Worst-case analysis of a new heuristic for the travelling salesman problem
Nicos Christofides. Worst-case analysis of a new heuristic for the travelling salesman problem. InOperations Research Forum, volume 3, page 20. Springer, 2022
2022
-
[18]
Scip: solving constraint integer programs.Mathematical Programming Computation, 1(1):1–41, 2009
Tobias Achterberg. Scip: solving constraint integer programs.Mathematical Programming Computation, 1(1):1–41, 2009
2009
-
[19]
Solving the equality generalized traveling salesman problem using the lin– kernighan–helsgaun algorithm.Mathematical Programming Computation, 7(3):269–287, 2015
Keld Helsgaun. Solving the equality generalized traveling salesman problem using the lin– kernighan–helsgaun algorithm.Mathematical Programming Computation, 7(3):269–287, 2015. 10
2015
-
[20]
Yuan Sun, Andreas Ernst, Xiaodong Li, and Jake Weiner. Generalization of machine learn- ing for problem reduction: a case study on travelling salesman problems.arXiv preprint arXiv:2005.05847, 2020. 11 A Technical Appendices and Supplementary Material This appendix supplements the line charts in the main text by reporting the corresponding detailed numeri...
Pith/arXiv arXiv 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.