REVIEW 3 major objections 5 minor 34 references
A Cycle Walk for Sampling Measures on Spanning Forests for Redistricting
T0 review · 3 major / 5 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read The Cycle Walk, a Metropolis-corrected Markov chain on spanning forests, can sample balanced graph partitions for target measures that penalize compactness directly—no longer relying on spanning-tree counts as a proxy—at parameter values wh
desk verdict The Cycle Walk is a genuinely new sampler with real promise, but the paper as written mis-defines the target measure and the reported 2-Tree acceptance ratio does not match it; the efficiency claims cannot be trusted until that is fixed. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is the added-edge cycle move. In a spanning forest, a 1-Tree step picks a district, adds one non-tree edge inside it, and removes one edge of the resulting unique cycle, changing the tree but not the partition; it is reversible with respect to the tree-count-weighted measure and acts as internal randomisation. A 2-Tree step picks two adjacent districts, adds two boundary edges to create a unique cycle spanning both, and removes two edges subject to population tolerance, yielding a new balanced pair. The reverse probability is computable because forward and backward proposals differ only by a small set of edge choices, which makes Metropolis–Hastings acceptance exact. Th
What would settle it
The decisive check is to compute the maximum total-variation distance, over a wider set of observables such as county splits and district boundary graphs, between independent runs of the Cycle Walk at γ=0.7 with the isoperimetric term; if that distance does not decay to a small value after 10^8 proposals, the empirical convergence claim fails. A secondary check is wall-clock cost per effective sample versus Forest RECOM.
Extended reading notes
Core claim
The discovery the paper tries to establish is that a Metropolised walk on spanning forests—proposals formed by adding edges to make a cycle and removing edges—samples measures νγ(τ) ∝ exp(−γJ_tree(ξ_τ)−J(ξ_τ)) for γ much closer to 1 than existing forest-recombination chains. A 1-Tree step adds a missing edge inside one district and removes an edge from the unique cycle, changing the forest but not the partition, reversibly. A 2-Tree step joins two adjacent districts with two boundary edges, forms a unique cycle, and removes two edges under population tolerance to give a new balanced pair. Since proposal probabilities factor into edge weights and a few choices, reverse probabilities are compu
Load-bearing premise
The load-bearing premise is that convergence can be judged from the maximum total-variation distance between the ranked-order partisan marginals of chains started from different initial plans; if the chain is slow in some feature that those marginals do not see, or if the assumed per-step cost advantage over forest recombination does not hold, the paper's claim of efficient sampling is not established.
Editorial extensions
If this is right
- Redistricting practitioners can target measures with an explicit isoperimetric compactness term, not only the spanning-tree-count proxy, at γ values where earlier tree-based samplers stall.
- The mixture of internal and inter-district moves gives a tunable interpolation between local flip dynamics and full recombination, so the sampler can be adapted to the energy landscape of a specific measure.
- Weighted edges can encode county-preservation or other structural preferences while keeping a known target measure, as demonstrated on South Carolina's 50-district senate maps.
- Because the 1-Tree step preserves the partition, the walk can refresh internal tree structure at low cost, which numerically accelerates convergence of the 2-Tree step at a roughly 9:1 ratio.
- The method's efficiency at γ near 1 suggests that exact-MCMC ensembles for uniform-on-partitions-plus-compactness may be reachable with modest compute, an object that has been effectively inaccessible.
Reading between the lines
- If the paper's convergence picture holds for the reported marginals, the same design should be tested on other observables such as county splits, racial or ethnic composition, and full boundary graphs before relying on it in legal settings; the paper's diagnostic is deliberately marginal-based.
- The 1-Tree step's value suggests a general recipe: add cheap partial resampling of internal tree structure to any global balanced-proposal chain to improve Metropolis acceptance; this could be tested by ablating the internal-to-paired ratio.
- The paper's Balanced Tree Up-Down Walk variant has a uniform stationary measure and may admit rigorous mixing bounds; a proof of rapid mixing there would transfer theoretical support to the Cycle Walk family.
- The paper explicitly forgoes a formal per-step cost analysis; a direct wall-clock comparison of effective sample rates would settle whether the convergence advantage translates into computational savings over forest recombination.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes the Cycle Walk, a Metropolis--Hastings Markov chain on spanning forests of a graph, designed to sample balanced graph partitions for redistricting. The chain mixes 1-tree cycle moves (internal to a district) and 2-tree cycle moves (exchanging population between adjacent districts), with population balance enforced in the proposal mechanism. The stated target family νγ interpolates between the uniform measure on spanning forests and the uniform measure on partitions lifted to forests, and the paper claims that the Cycle Walk mixes efficiently for γ values and compactness score functions that are difficult for existing forest-RECOM methods. Numerical experiments on 4×4 grids, North Carolina, Connecticut, regular lattices, and South Carolina county-weighted graphs are used to support the efficiency claims; the authors also release a Julia implementation.
Significance. If the central claim holds, the Cycle Walk would be a useful contribution: it would allow MCMC sampling of redistricting measures that penalize policy-relevant compactness directly, rather than implicitly through spanning-tree counts, while retaining exact Metropolis--Hastings correction. The paper has several concrete strengths: validation against exact enumeration on a 4×4 grid, comparison with Metropolized Forest RECOM for γ=0 and γ=0.2, release of code and tests, and explicit discussion of the limitations of tree-count-based compactness. However, the target measure and the reported acceptance probability are mutually inconsistent as written, so the numerical results do not currently establish the claimed capability. The issues are specific and correctable, but they affect the central claim.
major comments (3)
- [§6, Eq. (11)] The family νγ is stated inconsistently. With J_tree = -log tree(ξ), Eq. (11) gives νγ(τ) ∝ tree(ξτ)^γ e^{-J(ξτ)}. At γ=1 the partition marginal is proportional to tree(ξ)π(ξ), not π(ξ) as required by (6) and as stated in §10.2.1 ('γ=1 ... uniform on partitions'). To interpolate between ν0 and ν1 as defined in (6), the exponent needs the opposite sign, e^{+γ J_tree - J}, or J_tree must be defined as +log tree. As written, the target measure described in the text is not the measure defined by Eq. (11).
- [§8, Eq. (10) vs Eq. (13)] Equation (10) contains no factor involving treeα(ξτ)/treeα(ξτ'). Substituting the stated target (13) into the Metropolis--Hastings ratio yields an additional factor of the form (treeα(ξτ)/treeα(ξτ'))^{γ} or its inverse, depending on the sign convention used to correct Eq. (11). Since the paper states that 'our rejection probabilities take on the equation reported in (10),' the implemented 2-tree acceptances are inconsistent with the stated νγ for γ≠0. Consequently, the γ-dependence displayed in Figures 5–8 cannot be attributed to sampling the measure νγ as defined, and a reader cannot reproduce the target distribution from the manuscript.
- [§10.2.1, Figs. 5–8] The central efficiency claim rests on convergence diagnostics that are proxies rather than full-distribution checks: maximum total-variation distance between 14 ranked-order marginals from four chains, and single-chain effective sample sizes. The paper also explicitly forgoes a formal cost analysis in §10.2.1. Even after the target-measure and acceptance-ratio issues are corrected, the experiments need to be re-run with the corrected chain before the headline conclusion—efficient sampling for γ near 1 and for non-tree compactness scores—can be considered established. The current numerical evidence does not isolate the stated νγ.
minor comments (5)
- [§4.1, §6, §10.1.2] Typos: 'alos' should be 'also' (§4.1), 'mesure' should be 'measure' (§6), 'intrest' should be 'interest' (§10.1.2).
- [§10.2.4 and Fig. 10] The text says 'We begin by sampling ν0 over a variety of grid sizes,' but the surrounding discussion and Figure 10 caption refer to γ=1 with isoperimetric weights c=0.3,0.5,0.7,0.9. Please reconcile this discrepancy.
- [References] Reference [2] is corrupted ('Micha/suppress l Derezi nski'); please fix the author name and title.
- [§11] The walk is called 'Up-Down Walk' in the first paragraph and 'Down-Up Walk' in the next sentence; please use consistent terminology.
- [Figures 5–6] The labels 'Linked Forest Recom' are used in multiple panels without explicitly distinguishing top/bottom rows; adding '(no Polsby-Popper)' and '(w/ Polsby-Popper)' directly in each panel would improve readability.
Circularity Check
No significant circularity: the Cycle Walk's target measures and MH acceptance ratios are independently specified; self-citations are baselines, not load-bearing. (An internal sign inconsistency in the definition of ν_γ is a correctness issue, not circularity.)
full rationale
The central derivation is not circular. The target measures ν_γ are defined explicitly in Eqs. (6), (11), and (13) in terms of tree counts, tree-weight products, and an independently chosen score J(ξ); they are not fitted to the sampler's output. The Metropolis–Hastings acceptance probabilities in Eqs. (9)–(10) follow from the standard detailed-balance ratio ν(τ')Q(τ',τ)/(ν(τ)Q(τ,τ')), with the proposal kernels defined geometrically by adding edges to form cycles and removing edges to restore a forest. No parameter is calibrated to make a 'prediction' match held-out data. The compactness weight c_γ is tuned in Appendix A to keep mean isoperimetric score roughly constant across γ; this affects only the choice of test measures, not the sampler's correctness or the convergence diagnostic. Self-citations to Forest RECOM [4,6] and multiscale single-node flip [12] are used as experimental baselines, and the sampler is additionally validated against exact enumeration on a 4×4 grid (Sec. 10.1.1), which is independent supporting evidence. No load-bearing argument reduces to a self-citation chain. The one serious concern is an internal algebraic inconsistency: Section 6 defines J_tree = −log tree and then writes ν_γ ∝ e^{−γJ_tree−J}, which equals tree^γ e^{−J}, whereas (6) and the stated γ=1 uniform-on-partitions behavior require ν_γ ∝ tree^{−γ} e^{−J}; Eq. (10) also omits the corresponding tree-count ratio for γ≠1. This would change the measure actually sampled for γ≠0, but it is not circular: the target measure is still independently specified, and the bug does not make the output equivalent to the input by construction.
Assumptions & free parameters
free parameters (3)
- c_gamma (compactness weight) =
linear fit through origin, slope about 0.35 (Appendix A)
- kappa (mixing ratio of 1-tree to 2-tree moves) =
0.99 in the main NC experiments; varied in Section 10.2.2
- county edge upweights for South Carolina =
1, 10, 40
assumptions (4)
- standard math Metropolis-Hastings convergence requires the proposal chain to be irreducible and aperiodic on the relevant state space.
- standard math The 1-Tree Cycle Walk is reversible with respect to the measure m_xi^(alpha) from the cited literature [31, 3].
- domain assumption The population-balance tolerance and score functions define the support of the target measure.
- domain assumption The graphs studied are finite, connected, undirected planar graphs with positive edge weights.
Cite this review
Pith. "Pith review of A Cycle Walk for Sampling Measures on Spanning Forests for Redistricting." pith.science (2026). https://pith.science/paper/VELC3NZO
@misc{pith2026250908629,
author = {Pith},
title = {Pith review of: A Cycle Walk for Sampling Measures on Spanning Forests for Redistricting},
year = {2026},
howpublished = {\url{https://pith.science/paper/VELC3NZO}},
note = {Machine review of arXiv:2509.08629}
}
read the original abstract
We introduce the Cycle Walk, a new Markov chain Monte Carlo method for sampling distributions on balanced graph partitions, motivated by applications in political redistricting. The method operates on spanning forests and combines two types of updates: local "cycle" moves within districts and global moves that exchange population between adjacent districts while preserving balance constraints. This construction enables efficient Metropolis--Hastings correction while allowing proposals at multiple spatial scales. We show that the Cycle Walk naturally interpolates between existing approaches based on local updates and a class of global update methods derived from recombination (RECOM). Through a range of numerical experiments on synthetic graphs and real-world precinct data, we demonstrate that the Cycle Walk exhibits improved empirical convergence diagnostics for distributions that place weaker weight on spanning-tree counts, a regime that is challenging for existing methods. In particular, the algorithm remains effective when incorporating alternative compactness measures that more closely reflect policy-relevant criteria. These results suggest that the Cycle Walk provides a flexible and computationally efficient framework for sampling from a broader class of redistricting distributions than previously accessible with MCMC techniques.
Figures
Figures from the paper (9 more)
Reference graph
Works this paper leans on
-
[1]
Common Cause, No
Rucho v. Common Cause, No. 18-422, 588 U.S. (2019)
2019
-
[2]
Domain sparsifica- tion of discrete distributions using entropic independence.LIPIcs, Volume 215, ITCS 2022, 215:5:1–5:23, 2022
Nima Anari, Micha/suppress l Derezi´ nski, Thuy-Duong Vuong, and Elizabeth Yang. Domain sparsifica- tion of discrete distributions using entropic independence.LIPIcs, Volume 215, ITCS 2022, 215:5:1–5:23, 2022
2022
-
[3]
Log-concave polynomials iv: approximate exchange, tight mixing times, and near-optimal sampling of forests
Nima Anari, Kuikui Liu, Shayan Oveis Gharan, Cynthia Vinzant, and Thuy-Duong Vuong. Log-concave polynomials iv: approximate exchange, tight mixing times, and near-optimal sampling of forests. InProceedings of the 53rd Annual ACM SIGACT Symposium on The- ory of Computing, STOC 2021, page 408–420, New York, NY, USA, 2021. Association for Computing Machinery
2021
-
[4]
Metropolized forest recombination for monte carlo sampling of graph partitions.SIAM Jour- nal on Applied Mathematics, 83(4):1366–1391, 2023
Eric Autry, Daniel Carter, Gregory J Herschlag, Zach Hunter, and Jonathan C Mattingly. Metropolized forest recombination for monte carlo sampling of graph partitions.SIAM Jour- nal on Applied Mathematics, 83(4):1366–1391, 2023
2023
-
[5]
Autry, Daniel Carter, Gregory Herschlag, Zach Hunter, and Jonathan C
Eric A. Autry, Daniel Carter, Gregory Herschlag, Zach Hunter, and Jonathan C. Mattingly. Multi-scale merge-split markov chain monte carlo for redistricting. 2020. 4Seehttps://julialang.org/ 5For GerryChain seehttps://mggg.github.io/GerryChain/index.html 6Seehttps://networkx.org/ A CYCLE W ALK FOR SAMPLING SPANNING FOREST MEASURES 33
2020
-
[6]
Metropolized multiscale forest recombination for redistricting.Multiscale Modeling & Simu- lation, 19(4):1885–1914, 2021
Eric A Autry, Daniel Carter, Gregory J Herschlag, Zach Hunter, and Jonathan C Mattingly. Metropolized multiscale forest recombination for redistricting.Multiscale Modeling & Simu- lation, 19(4):1885–1914, 2021
1914
-
[7]
Models of random spanning trees, 2024
Eric Babson, Moon Duchin, Annina Iseli, Pietro Poggi-Corradini, Dylan Thurston, and Jamie Tucker-Foltz. Models of random spanning trees, 2024
2024
-
[8]
Sachet Bangia, Christy Vaughn Graves, Gregory Herschlag, Han Sung Kang, Justin Luo, Jonathan C. Mattingly, and Robert Ravier. Redistricting: Drawing the line.arXiv, 1704.03360, 2017
arXiv 2017
Show all 34 references
-
[9]
Spanning tree methods for sampling graph partitions, 2022
Sarah Cannon, Moon Duchin, Dana Randall, and Parker Rule. Spanning tree methods for sampling graph partitions, 2022
2022
-
[10]
Sampling balanced forests of grids in polynomial time, 2024
Sarah Cannon, Wesley Pegden, and Jamie Tucker-Foltz. Sampling balanced forests of grids in polynomial time, 2024
2024
-
[11]
On the complexity of sam- pling redistricting plans, 2023
Moses Charikar, Paul Liu, Tianyu Liu, and Thuy-Duong Vuong. On the complexity of sam- pling redistricting plans, 2023
2023
-
[12]
Mattingly
Gabriel Chuang, Gregory Herschlag, and Jonathan C. Mattingly. Multiscale parallel temper- ing for fast sampling on redistricting plans, 2024
2024
-
[13]
Colorado in Context: Congressional Redistricting and Competing Fairness Criteria in Col- orado.Journal of Computational Social Science, 5:180–226, May 2021
Jeanne Clelland, Haley Colgate, Daryl DeFord, Beth Malmskog, and Flavia Sancier-Barbosa. Colorado in Context: Congressional Redistricting and Competing Fairness Criteria in Col- orado.Journal of Computational Social Science, 5:180–226, May 2021
2021
-
[14]
Clelland, Nicholas Bossenbroek, Thomas Heckmaster, Adam Nelson, Peter Rock, and Jade VanAusdall
Jeanne N. Clelland, Nicholas Bossenbroek, Thomas Heckmaster, Adam Nelson, Peter Rock, and Jade VanAusdall. Compactness statistics for spanning tree recombination, 2021
2021
-
[15]
Cormen, Charles E
Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, and Clifford Stein.Introduction to algorithms. MIT Press, Cambridge, MA, third edition, 2009
2009
-
[16]
Redistricting reform in virginia: Districting criteria in context.Virginia Policy Review, 12(2):120–146, 2019
Daryl DeFord and Moon Duchin. Redistricting reform in virginia: Districting criteria in context.Virginia Policy Review, 12(2):120–146, 2019
2019
-
[17]
Recombination: A family of markov chains for redistricting.Harvard Data Science Review, 3(1), 2021
Daryl DeFord, Moon Duchin, and Justin Solomon. Recombination: A family of markov chains for redistricting.Harvard Data Science Review, 3(1), 2021
2021
-
[18]
Outlier analysis for pennsylvania congressional redistricting
Moon Duchin. Outlier analysis for pennsylvania congressional redistricting. https://www.governor.pa.gov/wp-content/uploads/2018/02/md-report.pdf
2018
-
[19]
Comparing Algorithms for Generating Ensembles to Detecting Ger- rymandering
Jonas Eichenlaub. Comparing Algorithms for Generating Ensembles to Detecting Ger- rymandering. https://sites.duke.edu/quantifyinggerrymandering/2023/08/23/comparing- algorithms-for-generating-ensembles-to-detecting-gerrymandering/, 2023
2023
-
[20]
A new automated redis- tricting simulator using Markov chain Monte Carlo.Work
Benjamin Fifield, Michael Higgins, Kosuke Imai, and Alexander Tarr. A new automated redis- tricting simulator using Markov chain Monte Carlo.Work. Pap., Princeton Univ., Princeton, NJ, 2015
2015
-
[21]
Automated redistrict- ing simulation using markov chain monte carlo.Journal of Computational and Graphical Statistics, 29(4):715–728, 2020
Benjamin Fifield, Michael Higgins, Kosuke Imai, and Alexander Tarr. Automated redistrict- ing simulation using markov chain monte carlo.Journal of Computational and Graphical Statistics, 29(4):715–728, 2020
2020
-
[22]
Marzuola, Jonathan C
Greg Herschlag, Jeremy L. Marzuola, Jonathan C. Mattingly, and Andrew Sun. A spectral analysis of tree based markov chains for redistricting problems.Preprint, 2025
2025
-
[23]
Quantifying gerrymandering in north carolina
Gregory Herschlag, Han Sung Kang, Justin Luo, Christy Vaughn Graves, Sachet Bangia, Robert Ravier, and Jonathan C Mattingly. Quantifying gerrymandering in north carolina. Statistics and Public Policy, 7(1):30–38, 2020
2020
-
[24]
Mattingly, Matthias Sachs, and Evan Wyse
Gregory Herschlag, Jonathan C. Mattingly, Matthias Sachs, and Evan Wyse. Non-reversible markov chain monte carlo for sampling of districting maps, 2020
2020
-
[25]
Kruskal, Jr
Joseph B. Kruskal, Jr. On the shortest spanning subtree of a graph and the traveling salesman problem.Proc. Amer. Math. Soc., 7:48–50, 1956
1956
-
[26]
J. C. Mattingly and C. Vaughn. Redistricting and the Will of the People.ArXiv e-prints, October 2014
2014
-
[27]
Mattingly
Jonathan C. Mattingly. Expert report for Common Cause v. Lewis.Common Cause v. Lewis, 2019
2019
-
[28]
Mattingly
Jonathan C. Mattingly. Expert report for Harper v. Hall/Moore.Harper v. Hall/Moore, 2021
2021
-
[29]
Complexity and geometry of sampling con- nected graph partitions, 2019
Elle Najt, Daryl DeFord, and Justin Solomon. Complexity and geometry of sampling con- nected graph partitions, 2019
2019
-
[30]
Empirical sampling of connected graph parti- tions for redistricting.Physical Review E, 104(6):064130, 2021
Elle Najt, Daryl DeFord, and Justin Solomon. Empirical sampling of connected graph parti- tions for redistricting.Physical Review E, 104(6):064130, 2021
2021
-
[31]
Linking and cutting spanning trees.Algorithms, 11(4):53, 2018
Lu´ ıs MS Russo, Andreia Sofia Teixeira, and Alexandre P Francisco. Linking and cutting spanning trees.Algorithms, 11(4):53, 2018. 34 A CYCLE W ALK FOR SAMPLING SPANNING FOREST MEASURES
2018
-
[32]
Swendsen and Jian-Sheng Wang
Robert H. Swendsen and Jian-Sheng Wang. Nonuniversal critical dynamics in monte carlo simulations.Phys. Rev. Lett., 58:86–88, Jan 1987
1987
-
[33]
Generating random spanning trees more quickly than the cover time
David Bruce Wilson. Generating random spanning trees more quickly than the cover time. InProceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing - STOC ’96, pages 296–303, Philadelphia, Pennsylvania, United States, 1996. ACM Press
1996
-
[34]
Mathematically quantifying non-responsiveness of the 2021 georgia congressional districting plan
Zhanzhan Zhao, Cyrus Hettle, Swati Gupta, Jonathan Christopher Mattingly, Dana Randall, and Gregory Joseph Herschlag. Mathematically quantifying non-responsiveness of the 2021 georgia congressional districting plan. InEquity and Access in Algorithms, Mechanisms, and Optimizati...
2021
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.