Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

Graph Theory

345 missions · 144 completed

Missions

Open201Completed144All345
Linear OptimizationOperations ResearchOptimization·Captain: mikedeng1

Minimum Cost Flows in Graphs with Unit Capacities: When the Second Phase of a Cost Scale Starts, at Most 3λ Excess Is Left for λ = √m and 36λ for λ = n^(2/3)Research Paper

Motivation

Minimum-cost flow is a standard model for routing a quantity through a network while paying a cost on each arc. Unit capacities arise when each input arc can carry at most one unit; they include flow formulations of several matching and assignment problems. Goldberg, Kaplan, Hed and Tarjan analyze a cost-scaling algorithm for this setting. Their second-phase argument depends on a bound on the excess left in a pseudoflow when every excess vertex has undergone enough potential increases. This mission formalizes that bound, Lemma 5 of their STACS 2015 paper. The published paper gives an O(min⁡{m1/2,n2/3} mlog⁡(nC))O(\min\{m^{1/2},n^{2/3}\}\,m\log(nC))O(min{m1/2,n2/3}mlog(nC)) running-time bound for its unit-capacity algorithm; the running-time statement itself is outside this mission.

Setting

A network has a finite vertex set VVV and an arc set EEE that contains both orientations of each input arc. Write n=∣V∣n=|V|n=∣V∣. Every input arc has capacity one; its reverse has capacity zero and opposite cost. Write mmm for the number of input arcs, so ∣E∣=2m|E|=2m∣E∣=2m. The paper assumes no parallel or antiparallel input arcs, allowing an arc to be identified by its ordered endpoints. Costs are antisymmetric between orientations.

A pseudoflow fff assigns an antisymmetric quantity to each arc and respects its capacity. Its excess at vvv is the sum ef(v)e_f(v)ef​(v) of incoming flow. Positive excess means that vvv has flow to discharge; negative excess is a deficit. A circulation is a pseudoflow with zero excess at every vertex. An arc is residual for fff when its residual capacity uf(v,w)=u(v,w)−f(v,w)u_f(v,w)=u(v,w)-f(v,w)uf​(v,w)=u(v,w)−f(v,w) is positive. For vertex prices ppp, its reduced cost is cp(v,w)=c(v,w)−p(v)+p(w)c_p(v,w)=c(v,w)-p(v)+p(w)cp​(v,w)=c(v,w)−p(v)+p(w). A pseudoflow is ε\varepsilonε-optimal at ppp when every residual arc has reduced cost at least −ε-\varepsilon−ε. These are the paper's conventions on p. 408.

At the start of a cost scale, let f′f'f′ be a circulation that is 2ε2\varepsilon2ε-optimal at prices p′p'p′. During the scale, fff is a pseudoflow that is ε\varepsilonε-optimal at prices ppp. The potential increase in scale units is d(v)=(p(v)−p′(v))/εd(v)=(p(v)-p'(v))/\varepsilond(v)=(p(v)−p′(v))/ε. Put λ=min⁡{m,n2/3}\lambda=\min\{\sqrt m,n^{2/3}\}λ=min{m​,n2/3}. The second phase begins when every excess vertex has d(v)≥λd(v)\ge\lambdad(v)≥λ; the paper also states that every deficit then has d(v)=0d(v)=0d(v)=0 and that the levels d(v)d(v)d(v) are integral. The remaining excess is X(f)=∑v:ef(v)>0ef(v)X(f)=\sum_{v:e_f(v)>0}e_f(v)X(f)=∑v:ef​(v)>0​ef​(v).

The comparison graph has arcs E+={(v,w)∈E:f′(v,w)>f(v,w)}E^+=\{(v,w)\in E:f'(v,w)>f(v,w)\}E+={(v,w)∈E:f′(v,w)>f(v,w)} and vertices VVV. An excess–deficit cut Y⊆E+Y\subseteq E^+Y⊆E+ meets every path in this graph from an excess vertex to a deficit vertex. Thus a cut is measured among comparison arcs, not among arbitrary arcs of the original network. Lemmas 3 and 4 of the paper establish the level and cut bounds on which the excess estimate rests.

Formalization targets

The goal is Lemma 5 on p. 410, stated with both constants obtained in its proof. Under the cost-scale conditions above, and the divisibility conventions in footnotes 5–6, its two cases are

λ=m⟹X(f)≤3λ,\lambda=\sqrt m\quad\Longrightarrow\quad X(f)\le 3\lambda,λ=m​⟹X(f)≤3λ, λ=n2/3⟹X(f)≤36λ.\lambda=n^{2/3}\quad\Longrightarrow\quad X(f)\le36\lambda.λ=n2/3⟹X(f)≤36λ.

The second case also assumes that six divides the integer λ\lambdaλ; the first uses divisibility by three. The goal keeps the two conclusions separate so that the stronger 3λ3\lambda3λ bound is retained when m\sqrt mm​ determines the scale. Its milestone list follows the paper: Lemma 3 bounds d(v)−d(w)d(v)-d(w)d(v)−d(w) by three on an arc of E+E^+E+; Lemma 4 bounds total excess by a cut; and the proof of Lemma 5 supplies the cut bands and the two size estimates. The bands use the exact three-level and six-level ranges printed on p. 410.

Significance

The result connects a condition on price increases to a numerical bound on the amount of flow still needing discharge. Combined with the paper's time per unit of remaining excess, this is the bound needed for its second-phase running-time estimate. The two regimes show how the same cost-scale state can be controlled by either the number of input arcs or the number of vertices. Without the excess bound, the paper's argument gives no comparable limit on the remaining work at the phase transition.

The paper proves Lemma 5, but this particular unit-capacity excess argument has no matching theorem in the platform catalog checked for this mission. The network, circulation and pseudoflow interfaces are already published as reusable definitions from earlier cost-scaling work. A completed machine-checked development would add the cut argument and band counts to those interfaces. Those facts can also support other analyses that compare an entry circulation with an intermediate pseudoflow; this mission does not claim to formalize the algorithm's entire execution or its running-time model.

Difficulty

The main obstacle is linking local price inequalities to a global bound on excess. A reduced-cost inequality controls a single residual arc, while excess may move through a path and is measured over many vertices. The comparison graph G+G^+G+ is essential: on a residual arc outside E+E^+E+, the level difference used by Lemma 3 is not available. Counting all residual arcs as if they were comparison arcs would therefore not justify the cut bound. There is a separate counting issue: the represented network has both directions of each input arc, so using ∣E∣|E|∣E∣ for mmm changes the constant in the square-root case.

Formalization scope

Lean represents the network by a finite symmetric arc set with real capacities and antisymmetric real costs. A local unit-capacity predicate marks each reverse pair as capacities (1,0)(1,0)(1,0) or (0,1)(0,1)(0,1); mmm counts the capacity-one orientations. Flows are real valued as on p. 408. Integer costs on that page are generalized to real costs because Lemmas 3–5 use only antisymmetry and reduced-cost inequalities. No result here depends on cost integrality.

The formalization takes the cost-scale state as hypotheses: ε>0\varepsilon>0ε>0, entry circulation and 2ε2\varepsilon2ε-optimality, current ε\varepsilonε-optimality, integral levels, level at least λ\lambdaλ at excess vertices, and zero level at deficits. It takes λ\lambdaλ as a natural number equal to min⁡{m,n2/3}\min\{\sqrt m,n^{2/3}\}min{m​,n2/3} and records the paper's divisibility simplifications. The indexed pigeonhole milestones also require λ>0\lambda>0λ>0 so their index ranges are nonempty; the goal includes the zero-arc boundary without requiring an index. Cuts use reachability after deleting a subset of E+E^+E+, equivalent to meeting every excess-to-deficit path. The proof of Lemma 5 prints “cut in GfG_fGf​”; the formalized claim uses G+G^+G+, the graph for which Lemma 3 gives the required level control and Lemma 4 states its bound.

The definition layer, path-cut lemma and finite counting arguments are the reusable infrastructure. Contributions that establish these exact statements, their necessary elementary flow identities, or a machine-checked goal proof are in scope. Counting mmm with reverse arcs, allowing cuts outside E+E^+E+, using arbitrary real levels in place of integral ones, or replacing both cases by the weaker 36λ36\lambda36λ bound does not express this target.

Selected references

  • A. V. Goldberg, H. Kaplan, S. Hed and R. E. Tarjan, Minimum Cost Flows in Graphs with Unit Capacities, 32nd International Symposium on Theoretical Aspects of Computer Science, LIPIcs 30, pp. 406–419, 2015. doi:10.4230/LIPIcs.STACS.2015.406.
9 thms1 active userReviewed
CombinatoricsGeometry & Topology·Captain: mikedeng1

Planar Graphs Have Bounded Queue-Number 2: Every Graph of Euler Genus g Has Queue-Number at Most 4g + 49Research Paper

Motivation

A queue layout of a graph places its vertices on a line and splits its edges into classes, called queues, so that inside one queue no edge is nested inside another. The minimum number of queues is the queue-number. Queue layouts were introduced by Heath, Leighton and Rosenberg (1992) and Heath and Rosenberg (1992) as the dual of stack layouts (book embeddings): reading the vertex order from left to right, the edges of one queue are opened and closed in first-in first-out order. Queue layouts are closely tied to track layouts and, through them, to three-dimensional grid drawings of graphs (§9 of the paper).

Whether planar graphs have bounded queue-number was asked by Heath, Leighton and Rosenberg in 1992 and stayed open for 27 years. The best general bounds before 2019 grew with the number of vertices: O(n)O(\sqrt n)O(n​) by Heath, Leighton and Rosenberg (1992), O(log⁡2n)O(\log^2 n)O(log2n) by Di Battista, Frati and Pach (2013), and O(log⁡n)O(\log n)O(logn) by Dujmović (2015); for Euler genus ggg, Dujmović, Morin and Wood (2017) had O(g+log⁡n)O(g + \log n)O(g+logn). Dujmović, Joret, Micek, Morin, Ueckerdt and Wood (arXiv:1904.04791, J. ACM 2020) proved that every planar graph has queue-number at most 494949 (Theorem 1), and extended the result to every graph embeddable in a fixed surface (Theorem 2) and to every proper minor-closed class (Theorem 3).

This mission is Theorem 2: graphs of Euler genus ggg have queue-number at most 4g+494g + 494g+49. It is the second mission of a two-mission series; the first poses Theorem 1.

Setting

All graphs are finite and simple. A linear ordering ⪯\preceq⪯ of the vertex set V(G)V(G)V(G) is fixed. Two edges vwvwvw and xyxyxy with four distinct ends nest if v≺x≺y≺wv \prec x \prec y \prec wv≺x≺y≺w (after renaming the ends). A queue is a set of edges no two of which nest. A kkk-queue layout of GGG is a linear ordering together with a partition of E(G)E(G)E(G) into kkk queues, and qn⁡(G)≤k\operatorname{qn}(G) \le kqn(G)≤k means that a kkk-queue layout exists.

A layering of GGG is an ordered partition (V0,V1,… )(V_0, V_1, \dots)(V0​,V1​,…) of V(G)V(G)V(G) such that every edge joins two vertices of the same or of consecutive layers. A BFS layering takes one root rrr in each connected component and puts vvv in layer dist⁡(r,v)\operatorname{dist}(r, v)dist(r,v). A BFS spanning tree of a connected graph GGG rooted at rrr is a spanning tree TTT with dist⁡T(r,v)=dist⁡G(r,v)\operatorname{dist}_T(r, v) = \operatorname{dist}_G(r, v)distT​(r,v)=distG​(r,v) for every vvv. A path in TTT is vertical if the distance to rrr grows by exactly one at every step.

The Euler genus of the orientable surface with hhh handles is 2h2h2h, and that of the non-orientable surface with ccc cross-caps is ccc; the Euler genus of a graph is the least Euler genus of a surface it embeds in (footnote 2 of the paper). Planar graphs are exactly the graphs of Euler genus 000. Combinatorially, an embedding of a graph is described by an embedding scheme: a cyclic order of the edges at each vertex (a rotation system) together with a sign on every edge (twisted or not). Tracing faces through the scheme gives a face count fff, and with nnn vertices, mmm edges and ccc components the scheme has Euler genus 2c−n+m−f2c - n + m - f2c−n+m−f.

Formalization targets

Goal: Theorem 2 with the paper's bound

Euler genus of G≤g ⟹ qn⁡(G)≤4g+49.\text{Euler genus of } G \le g \ \Longrightarrow\ \operatorname{qn}(G) \le 4g + 49 .Euler genus of G≤g ⟹ qn(G)≤4g+49.

The paper states Theorem 2 as qn⁡(G)=O(g)\operatorname{qn}(G) = O(g)qn(G)=O(g) and proves the explicit bound 4g+494g + 494g+49 in §5. The goal poses that explicit bound, for every finite graph, connected or not.

Milestones

  1. Lemma 18 (p. 18). Every BFS layering (V0,V1,… )(V_0, V_1, \dots)(V0​,V1​,…) of a planar graph GGG admits a 494949-queue layout whose vertex ordering lists V0V_0V0​, then V1V_1V1​, then V2V_2V2​, and so on.
  2. Lemma 21 (p. 19). If GGG is connected with Euler genus at most ggg and TTT is a BFS spanning tree with layering (Vi)(V_i)(Vi​), there is a connected (or empty) subgraph ZZZ with at most 2g2g2g vertices in each layer such that G−V(Z)G - V(Z)G−V(Z) is planar and sits inside a connected planar graph G+G^+G+ whose own BFS layering (Wi)(W_i)(Wi​) satisfies Wi∩(V(G)∖V(Z))=Vi∖V(Z)W_i \cap (V(G)\setminus V(Z)) = V_i \setminus V(Z)Wi​∩(V(G)∖V(Z))=Vi​∖V(Z); vertical paths of the BFS tree of G+G^+G+ restrict to vertical paths of TTT.
  3. The 4g4g4g-queue step (p. 23). If GGG has a layering with at most 2g2g2g vertices of ZZZ per layer and G−ZG - ZG−Z has a 494949-queue layout ordered layer by layer, then qn⁡(G)≤4g+49\operatorname{qn}(G) \le 4g + 49qn(G)≤4g+49.

Significance

The result. Theorem 2 settles boundedness of the queue-number for every class of graphs embeddable in a fixed surface, with a bound linear in the Euler genus. Through the connection between queue-number and track-number, the paper derives from it (§9, Theorems 54 and 58) that a graph of Euler genus ggg with nnn vertices has track-number gO(g4/7)g^{O(g^{4/7})}gO(g4/7) and a three-dimensional grid drawing of volume gO(g4/7)ng^{O(g^{4/7})} ngO(g4/7)n, linear in nnn for every fixed surface. It is also a step on the paper's route to Theorem 3 (proper minor-closed classes), since graphs on surfaces are the basic pieces of the Robertson–Seymour graph minor structure theorem. Lemma 21 is of independent use: it reduces a graph on a surface to a planar graph while controlling BFS layers.

Formalizing it. The result is proved; as far as is known it has no machine-checked proof. A formalization needs queue layouts, layerings, BFS layerings and vertical paths, and, independently of queues, a usable combinatorial theory of graphs on surfaces, which Lean's Mathlib does not have. The Euler genus definition of this mission (embedding schemes and face tracing) and Lemma 21 are reusable well beyond queue layouts.

Difficulty

The obvious reduction, deleting few vertices to make the graph planar, is not enough: a planar subgraph alone does not control how the deleted vertices interact with a queue layout. The deleted set must be thin in every BFS layer, and the planar remainder must carry a BFS layering compatible with the original one, which is what Lemma 21 provides. Its proof cuts the surface along a subgraph made of tree paths and must show that the cut graph is planar by an Euler-characteristic count over face-tracing data, and that distances from a new root agree with the old BFS distances. In Lean, the main obstacle is the surface side: Euler's formula for embedding schemes, the cutting operation on rotation systems, and the passage between the combinatorial notion of Euler genus and the drawing-based notion of planarity. The planar input, Lemma 18, is itself the whole of mission 1's argument (Theorem 15 with Lemmas 5 and 8).

Formalization scope

  • Graphs are SimpleGraph V on a finite type V : Type with decidable equality and adjacency.
  • "qn⁡(G)≤k\operatorname{qn}(G) \le kqn(G)≤k" is HasQueueLayout G k: an injective map ord : V → ℕ (the linear ordering) and an assignment of the edges of GGG (G.edgeSet, not all unordered pairs) to Fin k such that no two edges in one queue nest; nesting uses strict inequalities, so edges sharing an end never nest. Queues may be empty.
  • A layering is a layer function L : V → ℕ; a BFS layering takes one root per component. Vertical paths are nonempty lists indexed from 000 with distance d+id + id+i, d≥0d \ge 0d≥0.
  • Euler genus is combinatorial: EulerGenusLE G g asks for an embedding scheme (rotation permutation of the darts, cyclic at every vertex, and a symmetric edge signature) whose Euler genus 2c−n+m−f2c - n + m - f2c−n+m−f, computed in Z\mathbb ZZ, is at most ggg. The face count is half the number of face-tracing orbits plus the number of isolated vertices. Equivalence with footnote 2's topological definition is the Heffter–Edmonds–Ringel rotation principle and its non-orientable extension, with additivity over components; it is not part of the mission.
  • Planarity is the published RobertsonSeymour1986.GM5.IsPlanar: a crossing-free drawing in R2\mathbb R^2R2 by simple arcs. "Euler genus 000 iff planar" is a theorem connecting the two definitions, not built into either.
  • "Euler genus ggg" is read as "at most ggg" throughout; every bound is monotone in ggg. In Lemma 21, "ZZZ is connected" reads "empty or connected" (the paper takes Z=∅Z = \emptysetZ=∅ when g=0g = 0g=0), and the restriction of a vertical path may be empty.
  • Trivializing formalizations are ruled out: queues indexed on all pairs Sym2 V would deny edgeless graphs a 000-queue layout; non-strict nesting would forbid adjacent edges in one queue; a drawing that tolerates crossings would make every graph planar; a face count without halving, or without isolated vertices, would make the Euler genus too small or too large (the sanity checks confirm K3K_3K3​ has a scheme of Euler genus 000, a twisted K3K_3K3​ one of genus 111, and two isolated vertices genus 000).
  • Welcome contributions: Euler's formula for embedding schemes, face-tracing lemmas, the planar case (Lemma 18, via mission 1's results), the cutting construction of Lemma 21, and the purely combinatorial 4g4g4g-queue step.

Selected references

  • V. Dujmović, G. Joret, P. Micek, P. Morin, T. Ueckerdt, D. R. Wood, Planar graphs have bounded queue-number, J. ACM 67(4), 2020; arXiv:1904.04791v5. https://arxiv.org/abs/1904.04791 , https://doi.org/10.1145/3385731
  • L. S. Heath, F. T. Leighton, A. L. Rosenberg, Comparing queues and stacks as mechanisms for laying out graphs, SIAM J. Discrete Math. 5(3), 1992. https://doi.org/10.1137/0405031
  • L. S. Heath, A. L. Rosenberg, Laying out graphs using queues, SIAM J. Comput. 21(5), 1992. https://doi.org/10.1137/0221055
  • V. Dujmović, Graph layouts via layered separators, J. Combin. Theory Ser. B 110, 2015. https://doi.org/10.1016/j.jctb.2014.07.005
  • V. Dujmović, P. Morin, D. R. Wood, Layered separators in minor-closed graph classes with applications, J. Combin. Theory Ser. B 127, 2017. https://doi.org/10.1016/j.jctb.2017.05.006
  • B. Mohar, C. Thomassen, Graphs on Surfaces, Johns Hopkins University Press, 2001.
9 thms1 active userReviewed
ProbabilityRandom Matrix Theory·Captain: mikedeng1

Spectral Statistics of Erdős–Rényi Graphs I: Local Semicircle Law: Local Semicircle Law for the Sparse Noncentered Matrix A = H + f|e⟩⟨e| Down to Spectral Scale η ≈ 1/NResearch Paper

Motivation

The adjacency matrix of an Erdős–Rényi random graph contains a large deterministic mean as well as random fluctuations. The mean affects the largest eigenvalue, while the centered fluctuations determine much of the remaining spectrum. A local spectral law must describe both pieces at a resolution fine enough to distinguish eigenvalues within a narrow energy window. Erdős, Knowles, Yau and Yin established such estimates for a broader class of sparse symmetric matrices and their noncentered rank-one deformations in Spectral statistics of Erdős–Rényi graphs I. This mission formalizes the paper's Theorem 2.9, with the centered-matrix law and eleven intermediate results of the paper as milestones.

The paper treats the graph matrix through a model that also allows other sparse distributions. Its parameters qNq_NqN​ and fNf_NfN​ separate sparsity from the size of the deterministic mean. That separation matters: the normalized resolvent trace can obey a local semicircle law even when the deformation produces a distinguished top eigenvalue. The formal target therefore concerns the deformed matrix itself, rather than replacing it with the centered matrix.

Setting

For each size NNN, let HN=(hij)H_N=(h_{ij})HN​=(hij​) be a real symmetric random matrix. Entries indexed by the upper triangle are independent, centered, and have variance 1/N1/N1/N. Their higher moments obey

E∣hij∣p≤CpNqNp−2,3≤p≤(log⁡N)A0log⁡log⁡N.\mathbf E|h_{ij}|^p\le\frac{C^p}{Nq_N^{p-2}},\qquad 3\le p\le(\log N)^{A_0\log\log N}.E∣hij​∣p≤NqNp−2​Cp​,3≤p≤(logN)A0​loglogN.

The sparsity parameter qNq_NqN​ lies between (log⁡N)3ξN(\log N)^{3\xi_N}(logN)3ξN​ and CNC\sqrt NCN​, and the logarithmic parameter satisfies 1+a0≤ξN≤A0log⁡log⁡N1+a_0\le\xi_N\le A_0\log\log N1+a0​≤ξN​≤A0​loglogN, with a0>0a_0>0a0​>0 and A0≥10A_0\ge10A0​≥10. Define eN=N−1/2(1,…,1)e_N=N^{-1/2}(1,\ldots,1)eN​=N−1/2(1,…,1) and the noncentered matrix AN=HN+fN∣eN⟩⟨eN∣A_N=H_N+f_N|e_N\rangle\langle e_N|AN​=HN​+fN​∣eN​⟩⟨eN​∣, where 0≤fN≤NC0\le f_N\le N^C0≤fN​≤NC. The projection ∣eN⟩⟨eN∣|e_N\rangle\langle e_N|∣eN​⟩⟨eN​∣ has every entry equal to 1/N1/N1/N. These are Definitions 2.1 and 2.2 of the source paper, pp. 5–6.

For a complex number z=E+iηz=E+i\etaz=E+iη with η>0\eta>0η>0, the resolvents are G(z)=(HN−zI)−1G(z)=(H_N-zI)^{-1}G(z)=(HN​−zI)−1 and G~(z)=(AN−zI)−1\widetilde G(z)=(A_N-zI)^{-1}G(z)=(AN​−zI)−1. Their normalized traces are m(z)=N−1Tr⁡G(z)m(z)=N^{-1}\operatorname{Tr}G(z)m(z)=N−1TrG(z) and m~(z)=N−1Tr⁡G~(z)\widetilde m(z)=N^{-1}\operatorname{Tr}\widetilde G(z)m(z)=N−1TrG(z). The comparison object is the semicircle density ρsc(x)=(2π)−1[4−x2]+\rho_{\mathrm{sc}}(x)=(2\pi)^{-1}\sqrt{[4-x^2]_+}ρsc​(x)=(2π)−1[4−x2]+​​ and its Stieltjes transform msc(z)=∫Rρsc(x)/(x−z) dxm_{\mathrm{sc}}(z)=\int_{\mathbb R}\rho_{\mathrm{sc}}(x)/(x-z)\,dxmsc​(z)=∫R​ρsc​(x)/(x−z)dx. Put κE=∣∣E∣−2∣\kappa_E=\bigl||E|-2\bigr|κE​=​∣E∣−2​ and D={E+iη:∣E∣≤Σ, 0<η≤3}D=\{E+i\eta:|E|\le\Sigma,\ 0<\eta\le3\}D={E+iη:∣E∣≤Σ, 0<η≤3} for a fixed Σ≥3\Sigma\ge3Σ≥3. These definitions appear in Section 2, pp. 7–9.

An NNN-dependent event has (ξ,ν)(\xi,\nu)(ξ,ν)-high probability when its failure probability is at most exp⁡[−ν(log⁡N)ξN]\exp[-\nu(\log N)^{\xi_N}]exp[−ν(logN)ξN​] for all sufficiently large NNN. The spectral estimates below refer to single events on which the claimed bound holds at every z∈Dz\in Dz∈D simultaneously.

Formalization targets

Main goal: the local law for A

Theorem 2.9 assumes ξN∼(A0/2)log⁡log⁡N\xi_N\sim(A_0/2)\log\log NξN​∼(A0​/2)loglogN and qN≥(log⁡N)C1ξNq_N\ge(\log N)^{C_1\xi_N}qN​≥(logN)C1​ξN​. Its first conclusion is that, with (ξ,ν)(\xi,\nu)(ξ,ν)-high probability, every z=E+iη∈Dz=E+i\eta\in Dz=E+iη∈D satisfies

∣m~(z)−msc(z)∣≤(log⁡N)C2ξN(min⁡{1qN2κE+η,1qN}+1Nη).|\widetilde m(z)-m_{\mathrm{sc}}(z)|\le(\log N)^{C_2\xi_N}\left(\min\left\{\frac1{q_N^2\sqrt{\kappa_E+\eta}},\frac1{q_N}\right\}+\frac1{N\eta}\right).∣m(z)−msc​(z)∣≤(logN)C2​ξN​(min{qN2​κE​+η​1​,qN​1​}+Nη1​).

Under the additional bound fN≤C0Nf_N\le C_0\sqrt NfN​≤C0​N​, its second conclusion controls every entry at every z∈Dz\in Dz∈D:

∣G~ij(z)−δijmsc(z)∣≤(log⁡N)C2ξN(1qN+Im⁡msc(z)Nη+1Nη).|\widetilde G_{ij}(z)-\delta_{ij}m_{\mathrm{sc}}(z)|\le(\log N)^{C_2\xi_N}\left(\frac1{q_N}+\sqrt{\frac{\operatorname{Im}m_{\mathrm{sc}}(z)}{N\eta}}+\frac1{N\eta}\right).∣Gij​(z)−δij​msc​(z)∣≤(logN)C2​ξN​(qN​1​+NηImmsc​(z)​​+Nη1​).

The constants C1,C2C_1,C_2C1​,C2​ are universal and positive. The probability exponent ν>0\nu>0ν>0 may depend on the fixed parameters named in the theorem. Both conclusions, including the extra condition on fNf_NfN​ in the second, are retained from Theorem 2.9, pp. 9–10.

Milestones

The milestones follow the paper's own sequence of numbered results.

  1. Properties of mscm_{\mathrm{sc}}msc​ (Lemma 3.2): ∣msc∣∼1|m_{\mathrm{sc}}|\sim1∣msc​∣∼1, ∣1−msc2∣∼κE+η|1-m_{\mathrm{sc}}^2|\sim\sqrt{\kappa_E+\eta}∣1−msc2​∣∼κE​+η​, and the size of Im⁡msc\operatorname{Im}m_{\mathrm{sc}}Immsc​ inside and outside the spectrum.
  2. Resolvent identities for minors (Lemma 3.4) and the self-consistent equation Gii=1/(−z−msc−([v]−Υi))G_{ii}=1/(-z-m_{\mathrm{sc}}-([v]-\Upsilon_i))Gii​=1/(−z−msc​−([v]−Υi​)) (Lemma 3.10); both are deterministic.
  3. Large deviation bounds for linear and quadratic forms in independent variables with the sparse moment decay (Lemma 3.8 (ii)).
  4. The weak local semicircle law for HHH (Theorem 3.1) on DL={∣E∣≤Σ, (log⁡N)LN−1≤η≤3}D_L=\{|E|\le\Sigma,\ (\log N)^LN^{-1}\le\eta\le3\}DL​={∣E∣≤Σ, (logN)LN−1≤η≤3}, L≥8ξL\ge8\xiL≥8ξ.
  5. Fluctuation averaging (Lemma 4.1 (4.2)): the average [Z]=N−1∑iZi[Z]=N^{-1}\sum_iZ_i[Z]=N−1∑i​Zi​ of the fluctuation terms is much smaller than each ZiZ_iZi​.
  6. The a priori norm bound ∥H∥≤2+(log⁡N)ξq−1/2\|H\|\le2+(\log N)^\xi q^{-1/2}∥H∥≤2+(logN)ξq−1/2 (Lemma 4.3) and the strong local law for HHH (Theorem 2.8, the analogue of the goal for f=0f=0f=0).
  7. Interlacing of the eigenvalues of HHH and AAA (Lemma 6.1) and the near-orthogonality ∑α≠N∣⟨vα,e⟩∣2=O(f−2)\sum_{\alpha\ne N}|\langle v_\alpha,e\rangle|^2=O(f^{-2})∑α=N​∣⟨vα​,e⟩∣2=O(f−2) of eee to the non-top eigenvectors of AAA (Corollary 6.7).
  8. The trace comparison ∣Λ~−Λ∣≤π/(Nη)|\widetilde\Lambda-\Lambda|\le\pi/(N\eta)∣Λ−Λ∣≤π/(Nη) (Lemma 7.1), which yields the first conclusion of the goal from Theorem 2.8.
  9. Entrywise control of G~\widetilde GG on the event Ω~(η)\widetilde\Omega(\eta)Ω(η) that Λ~d+Λ~o≤(log⁡N)−ξ\widetilde\Lambda_d+\widetilde\Lambda_o\le(\log N)^{-\xi}Λd​+Λo​≤(logN)−ξ on the line Im⁡z=η\operatorname{Im}z=\etaImz=η: the off-diagonal bound Λ~o≤CΦ\widetilde\Lambda_o\le C\PhiΛo​≤CΦ (Proposition 7.6 (7.14)) and the diagonal bounds of Lemma 7.7, with Φ=(log⁡N)ξ/q+(log⁡N)2ξ(Im⁡m~/(Nη)+1/(Nη))\Phi=(\log N)^\xi/q+(\log N)^{2\xi}(\sqrt{\operatorname{Im}\widetilde m/(N\eta)}+1/(N\eta))Φ=(logN)ξ/q+(logN)2ξ(Imm/(Nη)​+1/(Nη)).

Significance

The trace estimate gives quantitative control over spectral mass at local scales, including near the edges E=±2E=\pm2E=±2. The entrywise estimate says more: each resolvent entry is close to the diagonal semicircle approximation with an error that records both sparsity and finite-size effects. The paper's later results derive density-of-states and eigenvector statements from these local controls. Those consequences are outside this mission's goal.

The intended formal development supplies reusable definitions of sparse symmetric ensembles, high-probability events, semicircle transforms, and complex resolvents, followed by exact statements of the paper's estimates. Mathlib contains no local semicircle law, for Wigner or for sparse matrices, and no Stieltjes transform of the semicircle distribution; the local law, the large deviation bounds and the interlacing lemma are each reusable beyond this paper.

Difficulty

The deterministic rank-one comparison controls the normalized trace but does not by itself control each entry of G~\widetilde GG. The entrywise estimate must account for the deformation and for fluctuations caused by sparse entries. Near the spectral edges, the stability scale of the semicircle transform changes with κE+η\kappa_E+\etaκE​+η; a bound that is uniform in the bulk does not automatically retain the paper's edge behavior. The theorem also asks for one event covering a continuum of zzz values, so a collection of separate fixed-zzz probability bounds has the wrong quantifier order. These are substantive parts of Sections 3–7, beyond algebraic manipulation of resolvents.

Formalization scope

Lean represents all HNH_NHN​ on one probability space and takes their entries to be real, symmetric, measurable random variables, independent in the upper triangle. Conditions involving log⁡log⁡N\log\log NloglogN and the moments hold for all sufficiently large NNN; this avoids impossible small-NNN parameter bounds and matches the asymptotic theorem. The moment order ppp is real, as in the paper. Each finite expectation has an explicit integrability condition. A single positive constant CCC is used for the moment, sparsity, and deformation bounds: replacing several fixed constants by their maximum preserves all three inequalities.

Matrix indices use Fin N rather than the paper's 1,…,N1,\ldots,N1,…,N. The resolvent is a complex matrix formed by coercing the real matrix before subtracting zIzIzI. Every result that uses it restricts η>0\eta>0η>0 eventually or explicitly; its inverse at a singular input is never used to make a claim. The semicircle transform is defined by the paper's integral, and 4−x2\sqrt{4-x^2}4−x2​ uses Lean's nonnegative square root, implementing the positive part. The high-probability predicate evaluates the measure of arbitrary sets, so it can state the intersection over all zzz without a separate measurability assumption on that intersection. Positive qNq_NqN​ and positive square-root arguments follow eventually from the model and z∈Dz\in Dz∈D.

The deterministic Lemmas 3.2, 3.4, 3.10, 6.1 and 7.1 are stated for every real symmetric matrix and every zzz with η>0\eta>0η>0 (Lemma 3.2 on DDD); the paper's DLD_LDL​ is contained there. Minors keep the original index names, and Gkl(T)G^{(\mathbb T)}_{kl}Gkl(T)​ is only used for k,l∉Tk,l\notin\mathbb Tk,l∈/T. ZiZ_iZi​ is the explicit expression of (3.15). Ordered eigenvalues are the sorted roots of the characteristic polynomial. Statements the paper makes "for fixed zzz" (Lemmas 4.1, 7.7, Proposition 7.6) are uniform in zzz: one threshold N0N_0N0​ serves every zzz. The unnamed constants CCC, ν\nuν, KKK are chosen after the fixed parameters and before the random matrices; the universal C1,C2C_1,C_2C1​,C2​ come first.

Trivializing encodings are ruled out: every resolvent is taken at η>0\eta>0η>0, so the junk inverse of a singular matrix never enters; the model is satisfiable (independent Rademacher signs divided by N\sqrt NN​ with qN=Nq_N=\sqrt NqN​=N​ and ξN=(A0/2)log⁡log⁡N\xi_N=(A_0/2)\log\log NξN​=(A0​/2)loglogN satisfy it for large NNN); each event of the goal is one intersection over all z∈Dz\in Dz∈D, not a family of pointwise bounds; the constants appear in the paper's order; and the goal is about A=H+f∣e⟩⟨e∣A=H+f|e\rangle\langle e|A=H+f∣e⟩⟨e∣ with fff allowed to be large, so it does not reduce to Theorem 2.8.

Part (iii) of Lemma 3.8 is not a milestone: the printed (3.22) omits the diagonal variance term (N−2∑i∣Bii∣2)1/2(N^{-2}\sum_i|B_{ii}|^2)^{1/2}(N−2∑i​∣Bii​∣2)1/2, without which it fails for Rademacher variables with q≍Nq\asymp\sqrt Nq≍N​. Contributions of proofs of any milestone, and of the missing groundwork (Stieltjes transforms of measures, Schur complements of resolvents, moment bounds for sums of independent variables), are welcome.

Selected references

  • L. Erdős, A. Knowles, H.-T. Yau and J. Yin, Spectral statistics of Erdős–Rényi graphs I: Local semicircle law, Annals of Probability 41(3B), 2013, pp. 2279–2375. arXiv:1103.1919v5; DOI:10.1214/11-AOP734.
16 thms1 active userReviewed
Combinatorics·Captain: mikedeng1

Planar Graphs Have Bounded Queue-Number 1: Every Planar Graph Has a 49-Queue LayoutResearch Paper

Motivation

A queue layout represents a graph by placing its vertices in a line and assigning edges to queues. The order and assignment must prevent nested edges from sharing a queue. This gives a graph parameter that measures how much queue capacity is needed for a particular graph. Heath, Leighton, and Rosenberg asked whether one fixed number of queues suffices for every planar graph. The question remained open while bounds depending on the graph's number of vertices, degree, or treewidth were known. Dujmović, Joret, Micek, Morin, Ueckerdt, and Wood resolved it by proving a uniform bound of 49. Their paper also develops layered partitions, a structural tool used beyond queue layouts. These claims and the historical account appear in the source paper, pp. 3–5 and 16.

Setting

A finite simple graph GGG has vertices V(G)V(G)V(G) and unordered edges E(G)E(G)E(G). Choose a linear order of V(G)V(G)V(G). Two disjoint edges nest when their endpoints appear as a<b<c<da<b<c<da<b<c<d, with one edge joining aaa to ddd and the other joining bbb to ccc. A queue is a set of edges with no nested pair. A kkk-queue layout is a vertex order together with an assignment of every edge to one of kkk queues; qn⁡(G)≤k\operatorname{qn}(G)\le kqn(G)≤k means such a layout exists. Edges sharing an endpoint cannot nest. Empty queues are allowed, so an edgeless graph can use zero queues. These are the paper's definitions in §1, pp. 3–4.

A layering assigns vertices to layers V0,V1,…V_0,V_1,\ldotsV0​,V1​,… such that endpoints of an edge lie in the same or adjacent layers. A BFS layering starts from one selected root in each connected component and places each vertex in the layer equal to its graph distance from that root. A vertex partition divides V(G)V(G)V(G) into disjoint nonempty parts. Its quotient graph G/PG/\mathcal PG/P has those parts as vertices; distinct parts are adjacent if an edge of GGG joins them. An HHH-partition may index the parts by vertices of another graph HHH: the quotient's edges must be represented by edges of HHH, while an indexed part may be empty. A partition has layered width at most ℓ\ellℓ with respect to a fixed layering if every part contains at most ℓ\ellℓ vertices of every layer. The exact definitions are in §2, pp. 7–9.

A graph is planar when it can be drawn in the plane with vertices at distinct points and edges as simple arcs that meet only at shared endpoints. A tree decomposition covers every vertex and edge of a graph by bags indexed by a tree, with the bags containing each vertex forming a connected subtree. Treewidth at most three means that a tree decomposition exists with every bag of size at most four. The paper uses these standard notions in §2 and §4, pp. 7, 12, and 16.

Formalization targets

Queue bound for planar graphs

The goal is the explicit upper bound the paper derives for its Theorem 1:

∀G finite and planar,qn⁡(G)≤49.\forall G\text{ finite and planar},\qquad \operatorname{qn}(G)\le49.∀G finite and planar,qn(G)≤49.

The printed theorem says that planar queue-number is bounded, and the authors state 49 immediately after it. A bound that depends on ∣V(G)∣|V(G)|∣V(G)∣ would not express the result. The goal uses the paper's own constant, as given in Theorem 1 and the following sentence, p. 5.

Layered partition and transfer bounds

The structural target is Theorem 15, p. 16: for every BFS layering of a planar graph, there is a partition of layered width at most three whose quotient remains planar and has treewidth at most three. The queue transfer target is Lemma 8, p. 10:

qn⁡(G)≤3ℓqn⁡(H)+⌊3ℓ2⌋\operatorname{qn}(G)\le 3\ell\operatorname{qn}(H)+ \left\lfloor\frac{3\ell}{2}\right\rfloorqn(G)≤3ℓqn(H)+⌊23ℓ​⌋

when GGG has an HHH-partition of layered width ℓ\ellℓ. Its layout orders vertices layer by layer. The milestone list also includes the rainbow characterization, the blowup lemma, the complete-graph queue number, and the five-queue bound for planar graphs of treewidth at most three. These are the stated supporting results along the route to 49 in §2–§4.1, pp. 7 and 10–16.

Significance

The goal establishes a constant queue bound for all planar graphs, independent of their size, degree, and treewidth. Theorem 15 supplies a reusable structural statement: even if the original graph has large treewidth, its quotient by small layer intersections has treewidth at most three while remaining planar. The paper goes on to treat graphs of bounded Euler genus and proper minor-closed classes, and describes connections with graph products and other graph parameters in its abstract and introduction, pp. 1 and 5.

The result is proved in the cited paper. This mission records Lean statements of the goal and six supporting milestones; their proofs remain to be formalized. A completed development would connect finite graph orders, nested-edge combinatorics, layered partitions, planar drawings, and tree decompositions in one machine-checked chain. The queue and partition definitions can also support later work on graph classes beyond the planar case.

Difficulty

A vertex order chosen without regard to a graph's structure can contain large rainbows, forcing many queues in that order. Planarity alone does not specify an order or an edge assignment, and a direct bound from ordinary treewidth does not cover planar graphs with arbitrarily large treewidth. The central structural requirement is to find a partition whose parts remain small in each BFS layer while its quotient has low treewidth. The queue transfer lemma must retain a layer-by-layer vertex order and account for edges within parts and between parts. The paper, pp. 10–11 and 16, identifies these as distinct constraints in its stated results.

Formalization scope

Vertices form a finite type and edges are those of a Mathlib simple graph. A vertex order is an injective map into N\mathbb NN; strict inequalities on four endpoints encode nesting. Queue assignments are functions on the graph's edge set, so the empty graph's zero-queue case is present. The statements use an existence predicate for qn⁡(G)≤k\operatorname{qn}(G)\le kqn(G)≤k, without taking a numerical infimum. Layerings are functions to N\mathbb NN, with BFS roots chosen separately for disconnected components. Theorem 15 uses a partition into nonempty parts, and the quotient has exactly those parts as vertices. Treewidth at most three uses a published tree-decomposition definition with bags of size at most four.

The locally checked drawing definition has the same fields as the published plane-drawing definition; the published module was unavailable from the local index in this run. Planarity means existence of a crossing-free drawing, not a drawing that permits crossings. The formulation does not use a universal edge carrier that would rule out zero queues, a non-strict nesting test, or a default numerical treewidth. The near-triangulation, Sperner, and tripod statements from the paper's deeper proof are outside this proposal's milestone scope; the partition theorem remains a full milestone. Solvers may contribute the omitted geometric and tree-decomposition infrastructure when proving it.

Selected references

  • V. Dujmović, G. Joret, P. Micek, P. Morin, T. Ueckerdt, and D. R. Wood, Planar graphs have bounded queue-number, Journal of the ACM 67(4), 2020. arXiv:1904.04791v5; DOI:10.1145/3385731.
9 thms1 active userReviewed
CombinatoricsLinear algebra·Captain: mikedeng1

A Harary-Sachs Theorem for Hypergraphs: The Codegree-d Coefficient of the Characteristic Polynomial of a k-Graph as a Weighted Sum over Its Veblen InfragraphsResearch Paper

Motivation

The characteristic polynomial of an ordinary graph records spectral information in coefficients that can also be described by finite graph configurations. The classical Harary–Sachs theorem is an instance of this connection. For a uniform hypergraph, the adjacency object is a tensor rather than a matrix, and its characteristic polynomial is defined through a multivariate resultant. Clark and Cooper's hypergraph theorem gives a corresponding coefficient formula. It identifies the configurations that replace the elementary subgraphs of the graph case and assigns them weights determined by directed Eulerian structures. The result makes a tensor spectral invariant accessible through finite combinatorial objects.

The article appeared as Gregory J. Clark and Joshua N. Cooper, A Harary-Sachs theorem for hypergraphs, Journal of Combinatorial Theory, Series B (2021), DOI 10.1016/j.jctb.2021.01.002. This mission follows the pinned arXiv version 2, whose page and theorem numbers are used below.

Setting

A simple kkk-graph H=([n],E)\mathcal H=([n],E)H=([n],E) has a fixed vertex set of nnn labeled vertices and a set EEE of distinct edges, each containing exactly kkk vertices. Isolated vertices remain part of [n][n][n]. The normalized adjacency hypermatrix AH\mathcal A_{\mathcal H}AH​ gives an ordered kkk-tuple the value 1/(k−1)!1/(k-1)!1/(k−1)! when its underlying set is an edge, and zero otherwise. This normalization makes every ordering of the other k−1k-1k−1 vertices contribute the intended total weight.

The characteristic polynomial ϕ(H)\phi(\mathcal H)ϕ(H) is the multivariate resultant of the coordinates of (λI−AH)xk−1(\lambda\mathcal I-\mathcal A_{\mathcal H})x^{k-1}(λI−AH​)xk−1. The resultant is characterized by three conditions: it vanishes exactly when homogeneous coordinate polynomials have a nonzero common complex root, it evaluates to one on the coordinate monomials, and it is irreducible over the complex coefficient field. Its degree in λ\lambdaλ is t=n(k−1)n−1t=n(k-1)^{n-1}t=n(k−1)n−1. The codegree-ddd coefficient ϕd(H)\phi_d(\mathcal H)ϕd​(H) is the coefficient of λt−d\lambda^{t-d}λt−d; this degree is fixed from the source theorem and is not read from an implementation of ϕ\phiϕ.

A multi-hypergraph SSS is a multiset of kkk-element edges. It is Veblen when the degree of every vertex is divisible by kkk. Its flattening discards repeated edges, so SSS is an infragraph of H\mathcal HH when every edge of its flattening belongs to EEE. Components use only vertices covered by edges; a labeled object may still sit in a larger ambient vertex set with unused vertices.

A rooting of SSS is a sequence of its edges, each assigned a root vertex in that edge, with roots in nondecreasing order. An edge rooted at uuu contributes a directed arc u→vu\to vu→v to every other vertex vvv of the edge. Summing these arcs with multiplicity produces a directed multigraph DRD_RDR​. An Euler rooting has balanced in- and out-degrees and connected positive-degree vertices. Its arborescence number τ(DR)\tau(D_R)τ(DR​) counts spanning trees of these active vertices directed toward a specified root, with parallel arcs distinguished. For connected SSS, its associated coefficient is

CS=∑R: DR Eulerianτ(DR)∏v∈V(DR)deg⁡−(v).C_S=\sum_{R:\,D_R\text{ Eulerian}}\frac{\tau(D_R)}{\prod_{v\in V(D_R)}\deg^-(v)}.CS​=R:DR​ Eulerian∑​∏v∈V(DR​)​deg−(v)τ(DR​)​.

For disconnected SSS, CSC_SCS​ is the product of the coefficients of its connected components. These conventions follow Definitions 6–10 in the source.

Formalization targets

Differential trace and Euler rootings

The supporting targets identify when a monomial differential operator contributes to a matrix trace, and show that an Euler rooting determines a connected Veblen multi-hypergraph. The operator identity is counted by Euler tours with distinguishable parallel arcs. The resultant-to-trace identity quoted in the paper is also a milestone:

ϕd(H)=Pd(−Tr⁡1(H)/1,…,−Tr⁡d(H)/d).\phi_d(\mathcal H)=P_d\left(-\operatorname{Tr}_1(\mathcal H)/1,\ldots,-\operatorname{Tr}_d(\mathcal H)/d\right).ϕd​(H)=Pd​(−Tr1​(H)/1,…,−Trd​(H)/d).

Here PdP_dPd​ is the paper's sum over positive compositions of ddd, and Tr⁡j\operatorname{Tr}_jTrj​ is its differential trace from Equation (1).

Harary–Sachs coefficient formula

The goal uses the labeled form of Theorem 14. Let Tm(d,H)\mathscr T_m(d,\mathcal H)Tm​(d,H) be the set of ordered mmm-tuples (S1,…,Sm)(S_1,\ldots,S_m)(S1​,…,Sm​) of connected labeled Veblen infragraphs whose edge counts sum to ddd. Then

ϕd(H)=∑m=1d(−(k−1)n)mm!∑(S1,…,Sm)∈Tm(d,H)∏i=1mCSi.\phi_d(\mathcal H)=\sum_{m=1}^{d}\frac{\bigl(-(k-1)^n\bigr)^m}{m!}\sum_{(S_1,\ldots,S_m)\in\mathscr T_m(d,\mathcal H)}\prod_{i=1}^{m}C_{S_i}.ϕd​(H)=m=1∑d​m!(−(k−1)n)m​(S1​,…,Sm​)∈Tm​(d,H)∑​i=1∏m​CSi​​.

The paper prints the equivalent sum over isomorphism classes, using automorphism groups and the copy weight (#H⊆H)(\#H\subseteq\mathcal H)(#H⊆H). The labeled expression is the result before that grouping. The isomorphism-class statement is outside this mission's Lean target and remains a separate formalization task.

Significance

The theorem expresses a coefficient of a resultant, defined algebraically over complex common roots, as a finite rational weighted sum over hypergraph configurations. It extends the graph case while preserving the exact normalization, sign, component count, and coefficient index. The formula can be used to study low codegree coefficients by enumerating finite Veblen infragraphs and their Euler rootings. Clark and Cooper also use associated coefficients to compare different hypergraph structures in the paper.

The formal development separates the resultant, the differential trace, and the combinatorial counts into reusable definitions. It would provide machine-checked statements of the paper's main equality and several links between those definitions. The Lean theorems in this draft are open targets with sorry placeholders; compiling their statements does not constitute a machine-checked proof of the paper's claims.

Difficulty

The ordinary matrix determinant has a direct expansion by permutations. A tensor characteristic polynomial is a multivariate resultant, so that expansion does not directly reveal which hypergraphs contribute to a coefficient. The trace side uses repeated partial derivatives, and repeated hyperedges can create several distinct rootings with the same directed arc multiplicities. Counting only a multiset of root assignments loses those orderings and changes the coefficient. The Eulerian condition, arborescence multiplicities, and factorial factors must agree simultaneously for the algebraic and combinatorial expressions to match.

Formalization scope

The Lean representation uses Fin n vertices, finite sets of simple edges, multisets of repeated edges, rational adjacency weights, and complex polynomials only for the resultant criterion. It assumes k≥2k\ge2k≥2, n≥1n\ge1n≥1, and 1≤d≤n(k−1)n−11\le d\le n(k-1)^{n-1}1≤d≤n(k−1)n−1. The source fixes k≥2k\ge2k≥2, and the latter two bounds expose the domain in which the printed coefficient expansion and its proof use ddd. An isolated vertex of H\mathcal HH counts in nnn and in (k−1)n(k-1)^n(k−1)n, while isolated vertices do not pad a multi-hypergraph component.

The formalization takes a polynomial satisfying the resultant's common-root, normalization, and irreducibility properties as an input. Separate sanity proofs give witnesses in the one-variable degree-zero and degree-one cases; degree one covers the first goal-domain degree when k=2k=2k=2. The characteristic polynomial is evaluated from this polynomial, never defined by the trace identity. Rootings are sequences, arborescences are directed toward the least positive-degree vertex, and the associated coefficient is calculated on covered vertices. The goal uses the labeled sum, so no automorphism is computed on an ambient padded vertex set. Contributions formalizing the omitted isomorphism-class equivalence, the BEST circuit count, and proofs of the stated targets are welcome.

Selected references

  • Gregory J. Clark and Joshua N. Cooper, A Harary-Sachs theorem for hypergraphs, Journal of Combinatorial Theory, Series B, 2021. arXiv:1812.00468v2; DOI 10.1016/j.jctb.2021.01.002.
8 thms1 active userReviewed
CombinatoricsGroup Theory·Captain: mikedeng1

Edge-Transitive Bi-Cayley Graphs 3: Every Connected Semisymmetric Bi-Dihedrant Has Valency at Least 6, and Semisymmetric Bi-Dihedrants of Valency 2k Exist for Each Odd k ≥ 3Research Paper

Motivation

A graph may have enough symmetry for every edge to look the same while its vertices still fall into two distinct classes. Such a graph is semisymmetric when it has constant valency, is edge-transitive, and is not vertex-transitive. Conder, Zhou, Feng and Zhang study this phenomenon among bi-Cayley graphs, graphs equipped with a group acting freely on two vertex orbits. Their dihedral case supplies examples related to questions of Marušič and Potočnik about semisymmetric tetracirculants. It also asks how small the valency of a semisymmetric graph with this specific group action can be. The paper states that six is the minimum and that examples occur at every valency 2k2k2k with odd k≥3k\ge3k≥3 (Conder et al., Theorem 1.5 and its introduction, p. 3).

The authors develop the bi-Cayley setting and the relevant automorphism actions before reaching the dihedral result. The lower bound appears as Theorem 6.1, while Example 6.2 defines a graph family and Proposition 6.4 gives a semisymmetry criterion for its odd-kkk members (Conder et al., pp. 11–13). This mission targets those claims within a single formal graph model. The final sentence of Theorem 1.5, concerning an edge-regular family, depends on the separate Theorem 6.5 and is outside this goal.

Setting

Let HHH be a finite group. Choose subsets R,L,S⊆HR,L,S\subseteq HR,L,S⊆H with R−1=RR^{-1}=RR−1=R, L−1=LL^{-1}=LL−1=L, 1∉R∪L1\notin R\cup L1∈/R∪L, and ∣R∣=∣L∣|R|=|L|∣R∣=∣L∣. The graph BiCay⁡(H,R,L,S)\operatorname{BiCay}(H,R,L,S)BiCay(H,R,L,S) has two labelled copies of each element of HHH, written h0h_0h0​ and h1h_1h1​. For each h∈Hh\in Hh∈H, it joins h0h_0h0​ to (xh)0(xh)_0(xh)0​ when x∈Rx\in Rx∈R, joins h1h_1h1​ to (yh)1(yh)_1(yh)1​ when y∈Ly\in Ly∈L, and joins h0h_0h0​ to (zh)1(zh)_1(zh)1​ when z∈Sz\in Sz∈S. These are unoriented edges. Right multiplication by HHH preserves the graph and acts regularly on each of its two vertex copies (Conder et al., pp. 1 and 4).

A bi-dihedrant is a bi-Cayley graph whose group is dihedral. Write Dn=⟨a,b∣an=b2=(ab)2=1⟩D_n=\langle a,b\mid a^n=b^2=(ab)^2=1\rangleDn​=⟨a,b∣an=b2=(ab)2=1⟩ for the dihedral group of degree n≥3n\ge3n≥3. Graph valency is the number of neighbours of a vertex. Connectedness is required throughout the paper. For a semisymmetric graph the common valency exists by definition, but the formal statements quantify over every vertex so that the numerical claims refer to the actual graph degree (Conder et al., definitions on p. 4 and Theorem 6.1 on p. 11).

The constructed family Γ(n,λ,2k)\Gamma(n,\lambda,2k)Γ(n,λ,2k) uses n≥5n\ge5n≥5, k≥2k\ge2k≥2, and a residue λ\lambdaλ of multiplicative order 2k2k2k modulo nnn whose even-power sum vanishes: ∑j=0k−1λ2j=0\sum_{j=0}^{k-1}\lambda^{2j}=0∑j=0k−1​λ2j=0. For 0≤i<k0\le i<k0≤i<k, put ci=∑j=0iλ2jc_i=\sum_{j=0}^{i}\lambda^{2j}ci​=∑j=0i​λ2j and di=λcid_i=\lambda c_idi​=λci​. Set S(n,λ,2k)={aci:0≤i<k}∪{badi:0≤i<k}S(n,\lambda,2k)=\{a^{c_i}:0\le i<k\}\cup\{ba^{d_i}:0\le i<k\}S(n,λ,2k)={aci​:0≤i<k}∪{badi​:0≤i<k}. Then Γ(n,λ,2k)=BiCay⁡(Dn,∅,∅,S(n,λ,2k))\Gamma(n,\lambda,2k)=\operatorname{BiCay}(D_n,\varnothing,\varnothing,S(n,\lambda,2k))Γ(n,λ,2k)=BiCay(Dn​,∅,∅,S(n,λ,2k)) (Conder et al., Example 6.2, pp. 12–13).

Right multiplication R(g):hi↦(hg)iR(g):h_i\mapsto(hg)_iR(g):hi​↦(hg)i​ gives a subgroup R(H)≤Aut⁡(Γ)R(H)\le\operatorname{Aut}(\Gamma)R(H)≤Aut(Γ), and X=NAut⁡(Γ)(R(H))X=N_{\operatorname{Aut}(\Gamma)}(R(H))X=NAut(Γ)​(R(H)) is its normaliser. For α∈Aut⁡(H)\alpha\in\operatorname{Aut}(H)α∈Aut(H) and g,x,y∈Hg,x,y\in Hg,x,y∈H the paper's permutations are σα,g:h0↦(hα)0, h1↦(ghα)1\sigma_{\alpha,g}:h_0\mapsto(h^\alpha)_0,\ h_1\mapsto(gh^\alpha)_1σα,g​:h0​↦(hα)0​, h1​↦(ghα)1​ and δα,x,y:h0↦(xhα)1, h1↦(yhα)0\delta_{\alpha,x,y}:h_0\mapsto(xh^\alpha)_1,\ h_1\mapsto(yh^\alpha)_0δα,x,y​:h0​↦(xhα)1​, h1​↦(yhα)0​; the sets F\mathrm FF and I\mathrm II of (2.2) collect those that preserve the edge rules. A graph is normal edge-transitive when XXX, not merely Aut⁡(Γ)\operatorname{Aut}(\Gamma)Aut(Γ), is transitive on edges (Conder et al., pp. 4–5).

Formalization targets

The milestones, in attack order, are Proposition 2.1(b) (one may assume 1∈S1\in S1∈S), Proposition 2.2 (the normaliser X=R(H)(F∪I)X=R(H)(\mathrm F\cup\mathrm I)X=R(H)(F∪I)), the opening step of the proof of Theorem 6.1 (R=L=∅R=L=\varnothingR=L=∅ for a semisymmetric graph), Theorem 6.1, the claims of Example 6.2, and Proposition 6.4.

Minimum valency

For every connected semisymmetric bi-dihedrant over DnD_nDn​, n≥3n\ge3n≥3, and every vertex vvv,

deg⁡Γ(v)≥6.\deg_{\Gamma}(v)\ge6.degΓ​(v)≥6.

This is Theorem 6.1 and the first claim of Theorem 1.5. The milestone for the proof's opening step records that semisymmetry forces R=L=∅R=L=\varnothingR=L=∅; the theorem itself still ranges over all valid R,L,SR,L,SR,L,S (Conder et al., Theorem 6.1, p. 11).

Examples at every odd index

For each odd integer k≥3k\ge3k≥3, there is a connected semisymmetric bi-dihedrant Γ\GammaΓ such that

deg⁡Γ(v)=2kfor every vertex v.\deg_{\Gamma}(v)=2k\qquad\text{for every vertex }v.degΓ​(v)=2kfor every vertex v.

Example 6.2 is the first milestone toward this claim: under its hypotheses Γ(n,λ,2k)\Gamma(n,\lambda,2k)Γ(n,λ,2k) is connected, ∣S∣=2k|S|=2k∣S∣=2k, every vertex has degree 2k2k2k, and σα,b\sigma_{\alpha,b}σα,b​ for α:(a,b)↦(aλ,ba)\alpha:(a,b)\mapsto(a^\lambda,ba)α:(a,b)↦(aλ,ba) lies in XXX and permutes the neighbours of 101_010​ cyclically, so the graph is normal edge-transitive (Conder et al., p. 13). Proposition 6.4 is the second: under the defining conditions for Γ(n,λ,2k)\Gamma(n,\lambda,2k)Γ(n,λ,2k), odd kkk and λk=−1\lambda^k=-1λk=−1 modulo nnn imply semisymmetry (Conder et al., Proposition 6.4, p. 13). An auxiliary arithmetic theorem states the parameter existence needed to apply it for every odd k≥3k\ge3k≥3.

Significance

The lower bound excludes semisymmetric bi-dihedrants of valency at most five. The existence assertion says that six is attained and places examples at an infinite sequence of larger valencies. Together these claims distinguish the restrictions imposed by a dihedral bi-Cayley action from edge transitivity alone. The construction also gives a concrete family on which questions about stronger automorphism properties can be posed (Conder et al., §6, pp. 11–14).

The result is proved in the cited paper; this mission seeks a machine-checked version of its first two Theorem 1.5 claims. The definition of a bi-Cayley graph, the precise graph symmetry predicates, and the arithmetic parameter family are reusable beyond this particular lower bound. The statements here are open proof targets in Lean, with their mathematical content fixed before proof development. The edge-regular claim from the theorem's final sentence requires additional source results and remains a separate formalization target.

Difficulty

Edge transitivity alone does not make a graph vertex-transitive: a semisymmetric graph is precisely the regular case where that inference fails. The lower bound therefore cannot be obtained by treating every vertex as lying in one automorphism orbit. In the dihedral case, the proof must account for both the rotations and reflections in the connection set, and for the constraints that connectedness and graph automorphisms impose on low-valency possibilities (Conder et al., proof of Theorem 6.1, pp. 11–12).

The existence statement has a distinct arithmetic requirement. Proposition 6.4 is conditional on n,λ,kn,\lambda,kn,λ,k satisfying Example 6.2 and λk=−1\lambda^k=-1λk=−1. The paper asserts examples for every odd k≥3k\ge3k≥3 without spelling out parameters for all such kkk. The auxiliary arithmetic statement records that remaining obligation explicitly; a proof of the conditional proposition alone would not establish the quantified existence claim.

Formalization scope

Lean represents the two vertex copies by H ⊕ H and the dihedral group by Mathlib's DihedralGroup n. The graph is a SimpleGraph built from the three edge rules above, with the inverse-closure and loop-exclusion conditions retained in BiCayData. Semisymmetry uses the full graph automorphism group, unoriented edge transitivity, and constant graph degree. The source's standing finite, connected, simple, and undirected conventions are explicit through finite groups, simple graphs, and connectedness hypotheses. The goal and Theorem 6.1 use n≥3n\ge3n≥3, the range of dihedral groups in §6; the existence examples are also required to be connected.

The definitions do not set valency equal to ∣S∣|S|∣S∣, replace a graph degree by a formula, or make an edge-transitivity claim true by ranging over an empty edge set. The setting requires valid inverse-closed side sets rather than silently symmetrising an invalid one. The arithmetic hypothesis uses the multiplicative order in Z/nZ\mathbb Z/n\mathbb ZZ/nZ and enumerates the indices of Zk\mathbb Z_kZk​ by 0,…,k−10,\ldots,k-10,…,k−1. Normal edge-transitivity in Example 6.2 is stated for the subgroup XXX and semisymmetry for the full group Aut⁡(Γ)\operatorname{Aut}(\Gamma)Aut(Γ); the two are not interchanged. The restriction n≥3n\ge3n≥3 in the lower bound is that of Theorem 6.1 and is a disclosed hypothesis of the goal: D1D_1D1​ and D2D_2D2​ are abelian, and the paper treats bi-Cayley graphs over abelian groups separately. The auxiliary parameter statement is not in the paper; it is supplied so that the existence claim has a complete proof path. The third sentence of Theorem 1.5 is not formalized. Contributions establishing the structural lower bound, the family criterion, or the arithmetic existence statement are all within scope.

Selected references

  • Marston Conder, Jin-Xin Zhou, Yan-Quan Feng and Mi-Mi Zhang, Edge-transitive bi-Cayley graphs, preprint arXiv:1606.04625v1, 2016; published in Journal of Combinatorial Theory, Series B, 2020. Preprint · DOI.
9 thms1 active userReviewed
CombinatoricsGroup Theory·Captain: mikedeng1

Edge-Transitive Bi-Cayley Graphs 1: For Connected BiCay(H, ∅, ∅, S), the Normaliser of R(H) Is 2-Arc-Transitive Iff Three Conditions on Aut(H) Hold, and Never 3-Arc-TransitiveResearch Paper

Motivation

A graph is symmetric when its automorphism group moves any vertex, edge or arc to any other, and sss-arc-transitive when it moves any walk of length sss without immediate backtracking to any other. Tutte's work on cubic graphs made sss-arc-transitivity a central measure of symmetry in algebraic graph theory, and most constructions of highly symmetric graphs start from a group acting regularly on vertices (Cayley graphs) or semiregularly with a few orbits.

Bi-Cayley graphs are the case of two orbits: a group HHH acts semiregularly by right multiplication with exactly two vertex orbits. The Petersen graph, the Gray graph (the smallest cubic semisymmetric graph) and the smallest member of Bouwer's family of half-arc-transitive graphs are all bi-Cayley graphs, and all of them are edge-transitive. Deciding edge-transitivity of a bi-Cayley graph is hard in general, but the subgroup of automorphisms that normalise the semiregular group R(H)R(H)R(H) is explicitly known, which makes the "normal" version of the question tractable.

C. H. Li asked whether 3-arc-transitive bi-normal Cayley graphs exist and asked for a description of the 2-arc-transitive ones (Li, Proc. Amer. Math. Soc. 133, Question 1.2(a) and Problem 1.3(b), as cited in the paper, p. 2). M. Conder, J.-X. Zhou, Y.-Q. Feng and M.-M. Zhang answered both through their Theorem 1.1 (arXiv:1606.04625v1; J. Combin. Theory Ser. B, 2020). This mission formalizes that theorem and the lemmas its proof uses.

Setting

Let HHH be a finite group with identity 111. Let R,L,S⊆HR, L, S \subseteq HR,L,S⊆H with ∣R∣=∣L∣|R| = |L|∣R∣=∣L∣, 1∉R=R−11 \notin R = R^{-1}1∈/R=R−1 and 1∉L=L−11 \notin L = L^{-1}1∈/L=L−1. The bi-Cayley graph Γ=BiCay(H,R,L,S)\Gamma = \mathrm{BiCay}(H,R,L,S)Γ=BiCay(H,R,L,S) has vertex set H0∪H1H_0 \cup H_1H0​∪H1​, two disjoint copies of HHH (the copy of hhh in HiH_iHi​ is hih_ihi​), and edges {h0,(xh)0}\{h_0,(xh)_0\}{h0​,(xh)0​}, {h1,(yh)1}\{h_1,(yh)_1\}{h1​,(yh)1​} and {h0,(zh)1}\{h_0,(zh)_1\}{h0​,(zh)1​} for h∈Hh \in Hh∈H, x∈Rx \in Rx∈R, y∈Ly \in Ly∈L, z∈Sz \in Sz∈S. The main theorem concerns R=L=∅R = L = \emptysetR=L=∅, where Γ\GammaΓ is bipartite with parts H0H_0H0​, H1H_1H1​ and every vertex has ∣S∣|S|∣S∣ neighbours.

For g∈Hg \in Hg∈H, R(g):hi↦(hg)iR(g) : h_i \mapsto (hg)_iR(g):hi​↦(hg)i​ is an automorphism of Γ\GammaΓ; these form the subgroup R(H)≤Aut(Γ)R(H) \le \mathrm{Aut}(\Gamma)R(H)≤Aut(Γ), and

X=NAut(Γ)(R(H))X = N_{\mathrm{Aut}(\Gamma)}(R(H))X=NAut(Γ)​(R(H))

is its normaliser. For α∈Aut(H)\alpha \in \mathrm{Aut}(H)α∈Aut(H) (written as an exponent, h↦hαh \mapsto h^\alphah↦hα) and x,y,g∈Hx, y, g \in Hx,y,g∈H the paper defines the permutations

δα,x,y:h0↦(xhα)1, h1↦(yhα)0,σα,g:h0↦(hα)0, h1↦(ghα)1,\delta_{\alpha,x,y} : h_0 \mapsto (xh^\alpha)_1,\ h_1 \mapsto (yh^\alpha)_0, \qquad \sigma_{\alpha,g} : h_0 \mapsto (h^\alpha)_0,\ h_1 \mapsto (gh^\alpha)_1,δα,x,y​:h0​↦(xhα)1​, h1​↦(yhα)0​,σα,g​:h0​↦(hα)0​, h1​↦(ghα)1​,

and the sets I\mathrm{I}I of those δα,x,y\delta_{\alpha,x,y}δα,x,y​ with Rα=x−1LxR^\alpha = x^{-1}LxRα=x−1Lx, Lα=y−1RyL^\alpha = y^{-1}RyLα=y−1Ry, Sα=y−1S−1xS^\alpha = y^{-1}S^{-1}xSα=y−1S−1x, and F\mathrm{F}F of those σα,g\sigma_{\alpha,g}σα,g​ with Rα=RR^\alpha = RRα=R, Lα=g−1LgL^\alpha = g^{-1}LgLα=g−1Lg, Sα=g−1SS^\alpha = g^{-1}SSα=g−1S.

An sss-arc is a sequence (v0,…,vs)(v_0,\dots,v_s)(v0​,…,vs​) of vertices in which consecutive vertices are adjacent and vi−1≠vi+1v_{i-1} \ne v_{i+1}vi−1​=vi+1​. A subgroup K≤Aut(Γ)K \le \mathrm{Aut}(\Gamma)K≤Aut(Γ) is transitive on sss-arcs if any sss-arc can be mapped to any other by an element of KKK; 1-arcs are arcs. Γ\GammaΓ is normal edge-transitive if XXX is transitive on edges, and normal locally arc-transitive if the stabiliser X10X_{1_0}X10​​ is transitive on the neighbourhood Γ(10)\Gamma(1_0)Γ(10​). KKK acts semisymmetrically if it is edge-transitive but not vertex-transitive. Aut(H,Δ)\mathrm{Aut}(H,\Delta)Aut(H,Δ) is the setwise stabiliser of Δ⊆H\Delta \subseteq HΔ⊆H in Aut(H)\mathrm{Aut}(H)Aut(H).

Formalization targets

Goal: Theorem 1.1

Let Γ=BiCay(H,∅,∅,S)\Gamma = \mathrm{BiCay}(H,\emptyset,\emptyset,S)Γ=BiCay(H,∅,∅,S) be connected, with 1∈S1 \in S1∈S and ∣S∣≥2|S| \ge 2∣S∣≥2. Then

X is transitive on the 2-arcs of Γ  ⟺  (a)∧(b)∧(c),X \text{ is transitive on the 2-arcs of } \Gamma \iff \text{(a)} \wedge \text{(b)} \wedge \text{(c)},X is transitive on the 2-arcs of Γ⟺(a)∧(b)∧(c),

where (a) Sα=S−1S^\alpha = S^{-1}Sα=S−1 for some α∈Aut(H)\alpha \in \mathrm{Aut}(H)α∈Aut(H); (b) Aut(H,S∖{1})\mathrm{Aut}(H, S\setminus\{1\})Aut(H,S∖{1}) is transitive on S∖{1}S \setminus \{1\}S∖{1}; (c) Sβ=s−1SS^\beta = s^{-1}SSβ=s−1S for some s∈S∖{1}s \in S\setminus\{1\}s∈S∖{1} and β∈Aut(H)\beta \in \mathrm{Aut}(H)β∈Aut(H). Furthermore, if ∣S∣≥3|S| \ge 3∣S∣≥3, XXX is not transitive on the 3-arcs of Γ\GammaΓ.

Milestones

  • Proposition 2.2 (cited by the paper from Zhou and Feng): for connected Γ\GammaΓ, every element of F\mathrm{F}F and of I\mathrm{I}I is an automorphism in XXX, and X=R(H) (F∪I)X = R(H)\,(\mathrm{F} \cup \mathrm{I})X=R(H)(F∪I).
  • Proposition 2.2(a): for δα,x,y∈I\delta_{\alpha,x,y} \in \mathrm{I}δα,x,y​∈I, ⟨R(H),δα,x,y⟩\langle R(H),\delta_{\alpha,x,y}\rangle⟨R(H),δα,x,y​⟩ is vertex-transitive.
  • Lemma 3.2: X1011={σα,1∣α∈Aut(H,S∖{1})}X_{1_01_1} = \{\sigma_{\alpha,1} \mid \alpha \in \mathrm{Aut}(H,S\setminus\{1\})\}X10​11​​={σα,1​∣α∈Aut(H,S∖{1})}, and, for ∣S∣≥3|S| \ge 3∣S∣≥3, XXX is not 3-arc-transitive.
  • Proposition 3.3: Γ\GammaΓ is normal locally arc-transitive iff Γ(10)\Gamma(1_0)Γ(10​) is an orbit of F\mathrm{F}F; in that case XXX is arc-transitive iff Sα=S−1S^\alpha = S^{-1}Sα=S−1 for some α\alphaα, and semisymmetric otherwise.

A companion item, Lemma 3.1, shows that a connected normal edge-transitive BiCay(H,R,L,S)\mathrm{BiCay}(H,R,L,S)BiCay(H,R,L,S) has R=L=∅R = L = \emptysetR=L=∅; it explains why the goal is stated for BiCay(H,∅,∅,S)\mathrm{BiCay}(H,\emptyset,\emptyset,S)BiCay(H,∅,∅,S).

Significance

Combined with Lemma 3.1, Theorem 1.1 reduces 2-arc-transitivity of the normaliser of R(H)R(H)R(H) to three conditions on Aut(H)\mathrm{Aut}(H)Aut(H) and the single set SSS, all checkable in the group without building the graph. Its second part rules out normal 3-arc-transitive bi-Cayley graphs of valency at least 3, which settles Li's question on 3-arc-transitive bi-normal Cayley graphs negatively. Proposition 3.3 is the arc-transitive versus semisymmetric dichotomy for normal locally arc-transitive bi-Cayley graphs, the starting point of the paper's later constructions of semisymmetric and half-arc-transitive examples.

All results are proved in the paper (Proposition 2.2 by citation). None of them has a machine-checked proof that this mission knows of. Mathlib has simple graphs, graph isomorphisms and normalisers, but no bi-Cayley graphs, no sss-arcs and no transitivity notions for subgroups of Aut(Γ)\mathrm{Aut}(\Gamma)Aut(Γ). A formal proof would supply the first verified description of the normaliser of a semiregular group with two orbits, reusable for any bi-Cayley or bi-circulant result.

Difficulty

The bulk of the work is Proposition 2.2. That σα,g\sigma_{\alpha,g}σα,g​ and δα,x,y\delta_{\alpha,x,y}δα,x,y​ preserve edges exactly under the conditions of (2.2) is a direct computation. The converse is harder: every automorphism normalising R(H)R(H)R(H) must be one of these permutations composed with some R(h)R(h)R(h). That needs the conjugation action of XXX on R(H)R(H)R(H) to induce an automorphism of HHH, and XXX to preserve the partition {H0,H1}\{H_0, H_1\}{H0​,H1​} into R(H)R(H)R(H)-orbits.

The 3-arc clause looks like a stabiliser count, but it holds only once 111_111​ has at least two neighbours besides 101_010​: for ∣S∣=2|S| = 2∣S∣=2 the graph is a cycle of length 2∣H∣2|H|2∣H∣, the normaliser is the whole dihedral automorphism group, and that group is transitive on 3-arcs. Any proof must use ∣S∣≥3|S| \ge 3∣S∣≥3 at exactly this step.

Formalization scope

The vertex set is H ⊕ H (Sum.inl h =h0= h_0=h0​, Sum.inr h =h1= h_1=h1​). The data (R,L,S)(R, L, S)(R,L,S) is a structure of three Finsets carrying the side conditions R−1=RR^{-1} = RR−1=R, L−1=LL^{-1} = LL−1=L, 1∉R1 \notin R1∈/R, 1∉L1 \notin L1∈/L and ∣R∣=∣L∣|R| = |L|∣R∣=∣L∣. Adjacency is given directly by the three edge rules. Aut(Γ)\mathrm{Aut}(\Gamma)Aut(Γ) is the group Γ ≃g Γ, R(H)R(H)R(H) is the subgroup generated by the R(g)R(g)R(g), and XXX is its Subgroup.normalizer. Automorphisms of HHH are MulAut H, and Δα\Delta^\alphaΔα is the image of Δ\DeltaΔ under α\alphaα. σα,g\sigma_{\alpha,g}σα,g​ and δα,x,y\delta_{\alpha,x,y}δα,x,y​ are permutations of H ⊕ H, and "σα,g∈X\sigma_{\alpha,g} \in Xσα,g​∈X" means that some element of XXX has σα,g\sigma_{\alpha,g}σα,g​ as its vertex map. Products are written pointwise to avoid the composition-order convention of the group of graph automorphisms. An sss-arc is a function on Fin (s+1).

The paper assumes throughout that groups and graphs are finite and graphs connected; each theorem carries Fintype H and connectivity. Three hypotheses are added to the printed Theorem 1.1, each needed for the statement to be true or non-vacuous: 1∈S1 \in S1∈S (the paper's normalisation of a bi-Cayley triple, used by conditions (b) and (c) and by the proof); ∣S∣≥2|S| \ge 2∣S∣≥2 (for ∣S∣=1|S| = 1∣S∣=1 the graph is K2K_2K2​, which has no 2-arcs, so 2-arc-transitivity would hold vacuously while (c) fails); and ∣S∣≥3|S| \ge 3∣S∣≥3 for the 3-arc clause (false for cycles). Lemma 3.2 carries the same 1∈S1 \in S1∈S and ∣S∣≥3|S| \ge 3∣S∣≥3, and Proposition 3.3 carries 1∈S1 \in S1∈S.

Trivializing formalizations are ruled out in four ways. 2-arc-transitivity is never vacuous, since ∣S∣≥2|S| \ge 2∣S∣≥2 guarantees 2-arcs. RRR and LLL are never symmetrised silently, since the side conditions are data rather than repairs. Every transitivity statement is about XXX, not the full Aut(Γ)\mathrm{Aut}(\Gamma)Aut(Γ). Valency is the degree in the graph, never a formula in ∣S∣|S|∣S∣.

Contributions are welcome at every level: proofs of the milestones, a general library of bi-Cayley graphs and their normalisers, and lemmas on sss-arcs and stabilisers in SimpleGraph. The definitions are repeated identically in the companion missions on bi-abelian graphs and bi-dihedrants.

Selected references

  • M. Conder, J.-X. Zhou, Y.-Q. Feng and M.-M. Zhang, Edge-transitive bi-Cayley graphs, arXiv:1606.04625v1, 2016; J. Combin. Theory Ser. B, 2020. https://arxiv.org/abs/1606.04625v1, https://doi.org/10.1016/j.jctb.2020.05.006
  • J.-X. Zhou and Y.-Q. Feng, The automorphisms of bi-Cayley graphs, J. Combin. Theory Ser. B 116 (2016) 504–532 (reference [55] of the paper; source of Proposition 2.2). https://doi.org/10.1016/j.jctb.2015.10.004
  • C. H. Li, Finite s-arc transitive Cayley graphs and flag transitive projective planes, Proc. Amer. Math. Soc. 133 (2005) 31–41 (reference [29] of the paper; Question 1.2(a) and Problem 1.3(b)). https://doi.org/10.1090/S0002-9939-04-07549-5
6 thms1 active userReviewed
CombinatoricsGroup Theory·Captain: mikedeng1

Edge-Transitive Bi-Cayley Graphs 2: Edge-Transitive Bi-Cayley Graphs over Abelian Groups Are Vertex-Transitive, and Half-Arc-Transitive Ones Have Valency at Least 6Research Paper

Why bi-Cayley graphs over abelian groups matter

A Cayley graph turns multiplication in a group into adjacency between vertices. Every Cayley graph is vertex transitive: multiplying every vertex by the same group element moves any chosen vertex to any other. A bi-Cayley graph starts with two copies of a group and permits edges both within and between the copies. Its built-in group action is transitive on each copy, but it need not exchange them. This makes bi-Cayley graphs a setting for distinguishing edge, vertex, and arc transitivity. Conder, Zhou, Feng and Zhang study these distinctions and their consequences for graph constructions in Edge-transitive bi-Cayley graphs.

For abelian groups, the paper proves a strong restriction: once a connected bi-Cayley graph is edge transitive, it is also vertex transitive. Its half-arc-transitive examples face additional local restrictions, including a minimum valency of six. The paper exhibits a graph of valency six in Example 4.3, so the numerical bound is attained in its class. The result separates these graphs from semisymmetric graphs, which are edge transitive but not vertex transitive. The claims and example appear in Proposition 1.3 and §4.

Graphs, actions, and terminology

Let HHH be a finite group and let R,L,SR,L,SR,L,S be subsets of HHH. Require R=R−1R=R^{-1}R=R−1, L=L−1L=L^{-1}L=L−1, 1∉R∪L1\notin R\cup L1∈/R∪L, and ∣R∣=∣L∣|R|=|L|∣R∣=∣L∣. The graph Γ=BiCay⁡(H,R,L,S)\Gamma=\operatorname{BiCay}(H,R,L,S)Γ=BiCay(H,R,L,S) has vertices H0∪H1H_0\cup H_1H0​∪H1​, two labelled copies of HHH. For h∈Hh\in Hh∈H, an element x∈Rx\in Rx∈R gives an edge from h0h_0h0​ to (xh)0(xh)_0(xh)0​; an element y∈Ly\in Ly∈L gives an edge from h1h_1h1​ to (yh)1(yh)_1(yh)1​; and an element z∈Sz\in Sz∈S gives an edge from h0h_0h0​ to (zh)1(zh)_1(zh)1​. The inverse closure of RRR and LLL makes the within-copy edges undirected. Excluding 111 from those two sets excludes loops. This is the construction on p. 1.

Every g∈Hg\in Hg∈H acts by right multiplication R(g):hi↦(hg)iR(g):h_i\mapsto(hg)_iR(g):hi​↦(hg)i​. The subgroup R(H)R(H)R(H) of graph automorphisms has two vertex orbits, H0H_0H0​ and H1H_1H1​. Its normaliser in Aut⁡(Γ)\operatorname{Aut}(\Gamma)Aut(Γ) is useful for describing further symmetries. For a group automorphism α\alphaα and group elements g,x,yg,x,yg,x,y, the paper defines permutations σα,g\sigma_{\alpha,g}σα,g​, which preserve the two copies, and δα,x,y\delta_{\alpha,x,y}δα,x,y​, which exchange them. Their precise edge-preservation conditions are the sets FFF and III of equation (2.2). The mission includes those definitions and the relevant part of the normaliser characterization, Proposition 2.2.

A graph is vertex transitive if its full automorphism group can move any vertex to any other, and edge transitive if it can move any undirected edge to any other. Arc transitivity asks the same for ordered adjacent pairs. A graph is half-arc-transitive when it is vertex and edge transitive but not arc transitive. A semisymmetric graph has constant valency, is edge transitive, and is not vertex transitive. These are the paper's definitions on p. 4. Valency means the actual number of neighbours of a vertex.

Formalization targets

The central target is the corrected reading of Proposition 1.3: for connected Γ\GammaΓ over an abelian HHH,

Γ edge transitive  ⟹  Γ vertex transitive,\Gamma\text{ edge transitive}\implies\Gamma\text{ vertex transitive},Γ edge transitive⟹Γ vertex transitive,

and, if Γ\GammaΓ is half-arc-transitive,

R∪L≠∅,(∀x∈R∪L) ord⁡(x)≠2,∣R∣=∣L∣ is even,∣S∣>2,deg⁡Γ(v)≥6(v∈V(Γ)).R\cup L\ne\varnothing,\quad (\forall x\in R\cup L)\ \operatorname{ord}(x)\ne2, \quad |R|=|L|\text{ is even},\quad |S|>2,\quad \deg_{\Gamma}(v)\ge6\quad(v\in V(\Gamma)).R∪L=∅,(∀x∈R∪L) ord(x)=2,∣R∣=∣L∣ is even,∣S∣>2,degΓ​(v)≥6(v∈V(Γ)).

The attack path records the normaliser characterization in Proposition 2.2, both parts of Proposition 4.1, its explicit bipartition step, and Proposition 4.2. The companion Corollary 4.4 expresses the absence of connected semisymmetric examples over abelian groups.

What the result establishes

The vertex-transitivity conclusion collapses one apparent possibility in the symmetry classification: an edge-transitive bi-Cayley graph over an abelian group cannot be semisymmetric. For half-arc-transitive graphs, the theorem rules out involutions in the within-copy connection sets, forces even set sizes, and excludes valencies below six. Example 4.3 supplies an attained bound, with H=C28H=C_{28}H=C28​ and valency six. These are concrete restrictions on any later classification or construction in this family; they follow from §4.

The mathematical results are proved in the cited paper. The mission asks for machine-checked Lean proofs of their statements and for a reusable finite bi-Cayley graph interface. The setting includes the right regular action, its normaliser, the two families of candidate automorphisms, and the transitivity predicates. The proof obligations are open in this local proposal; a successful development would make the abelian argument and its underlying graph-action lemmas available for later formalizations of bi-Cayley graphs.

The main obstacle

The right multiplication action gives only two known vertex orbits, one for each copy of HHH. Edge transitivity alone does not name an automorphism that swaps them. Inversion is an automorphism of an abelian group, but its induced copy-swapping permutation preserves graph edges directly only in the relevant bipartite case R=L=∅R=L=\varnothingR=L=∅. The paper must account for edge-transitive graphs whose connection sets initially allow edges inside the copies. For the valency bound, the absence of arc transitivity restricts how a vertex stabiliser can permute its neighbours; degree counting without those orbit constraints cannot yield ∣S∣>2|S|>2∣S∣>2. These are the points at which a proof must use the full automorphism action rather than just the visible right regular subgroup.

Formalization scope

Lean uses the vertex type H⊔HH\sqcup HH⊔H, with Sum.inl h for h0h_0h0​ and Sum.inr h for h1h_1h1​. The group and graph are finite; the graph is simple and undirected by construction, and every theorem assumes it is connected. The finite sets R,L,SR,L,SR,L,S retain the conditions R=R−1R=R^{-1}R=R−1, L=L−1L=L^{-1}L=L−1, 1∉R∪L1\notin R\cup L1∈/R∪L, and ∣R∣=∣L∣|R|=|L|∣R∣=∣L∣. The graph's degree, rather than a separately defined cardinality formula, represents valency. “Involution” means an element of order exactly two. All graph-level transitivity predicates quantify over the full automorphism group; the normaliser appears only in the separate Proposition 2.2 milestone. The paper's exponent notation for group automorphisms is rendered as function application, while R(g)R(g)R(g) remains right multiplication.

The printed first sentence of Proposition 1.3 says every connected bi-Cayley graph over an abelian group is vertex transitive. That sentence is false: the connected generalized Petersen graph GP(7,2)GP(7,2)GP(7,2) is a bi-Cayley graph over C7C_7C7​ and is not vertex transitive. Proposition 4.1(b), which the paper uses to establish Proposition 1.3, includes the missing edge-transitivity condition. The Lean goal states that corrected condition explicitly. It does not add an assumption to the half-arc-transitive clause, since half-arc-transitivity already entails edge transitivity. Proposition 4.2 prints two clauses labelled “(a)”; they are treated as distinct clauses while its quotation retains the printed labels.

The definition of the graph keeps RRR and LLL inverse closed as data; no graph constructor silently symmetrises arbitrary sets. Transitivity is taken over actual edges or arcs, and the degree bound is checked at every vertex. Contributions to the graph-action lemmas, the normaliser result, Proposition 4.1, and the vertex-stabiliser analysis behind Proposition 4.2 are all within scope. The Cayley-graph isomorphism in Proposition 2.2(b) is outside this mission's stated milestone because it is not used by the target.

Selected references

  • M. Conder, J.-X. Zhou, Y.-Q. Feng, and M.-M. Zhang, Edge-transitive bi-Cayley graphs, Journal of Combinatorial Theory, Series B, 2020; arXiv:1606.04625v1; DOI:10.1016/j.jctb.2020.05.006.
7 thms1 active userReviewed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Twin-width I: Tractable FO Model Checking 4: Every K_t-Minor-Free Graph Has Twin-width at Most 4c_g(t)·2^(4c_g(t)+2), a Triple-Exponential Bound in tResearch Paper

Why excluded minors matter here

Many graph classes are defined by excluding a fixed graph as a minor. A minor can be formed by deleting vertices, deleting edges, and contracting edges. Planar graphs, for example, exclude K5K_5K5​ as a minor. Édouard Bonnet, Eun Jung Kim, Stéphan Thomassé, and Rémi Watrigant proved that every class excluding a fixed clique minor has bounded twin-width, a graph parameter that controls how much adjacency information is lost while vertices are successively combined. Their result connects a classical graph-minor restriction with the contraction parameter used in their first-order model-checking framework. The explicit bound is very large, but it establishes that the width depends only on the excluded clique size, rather than on the graph's number of vertices. Bonnet et al., Twin-width I, Theorem 6.3.

Graphs, contractions, and matrix orders

Let GGG be a finite simple graph. A partition of its vertices is a collection of nonempty, disjoint parts covering all vertices. Two distinct parts are homogeneous if every vertex pair across them is adjacent, or every such pair is nonadjacent. Otherwise the pair of parts is marked red. A ddd-partition has at most ddd red neighbours at every part. A contraction sequence begins with singleton parts, merges two parts at each step, and finishes with at most one part. The statement tww⁡(G)≤d\operatorname{tww}(G)\le dtww(G)≤d says that such a sequence exists with every intermediate partition a ddd-partition. This partition formulation is equivalent to the paper's original description in terms of trigraph contractions. Bonnet et al., pp. 3:11–3:12, 3:32.

Write KtK_tKt​ for the complete graph on ttt vertices. The graph GGG is KtK_tKt​-minor free when KtK_tKt​ cannot be obtained from GGG by the minor operations above. The Lean statement uses a branch-set formulation: a minor model consists of disjoint, nonempty, connected vertex sets, with a connecting edge for each edge of the minor. Bonnet et al., p. 3:9.

An ordering σ\sigmaσ of the vertices gives a Boolean adjacency matrix Aσ(G)A_\sigma(G)Aσ​(G), with entry 111 at an edge and 000 otherwise. A division cuts its rows and columns into nonempty consecutive intervals. A zone is mixed if it is neither constant down its columns nor constant across its rows. A kkk-mixed minor is a division with kkk row intervals and kkk column intervals in which every zone is mixed; a kkk-grid minor requires a 111 in every zone. A matrix is kkk-mixed free or kkk-grid free if the corresponding minor does not exist. A mixed Boolean zone always contains a 111, so grid freeness implies mixed freeness. Bonnet et al., pp. 3:18–3:19.

The order used in the proof is a depth-first search (DFS). A discovery order v1,…,vnv_1,\dots,v_nv1​,…,vn​ lists the vertices; at each time the active vertex is the lastly discovered vertex that still has an undiscovered neighbour, and the next vertex is a neighbour of the active vertex, which becomes its parent in the DFS tree T\mathcal TT. The Lex-DFS breaks ties: among the connected components of the undiscovered vertices that meet the active vertex's neighbourhood, it enters one whose word (the 0/10/10/1 vector recording which discovered vertices have a neighbour in the component) is lexicographically largest. The proof also uses minimal subtrees of T\mathcal TT: the union of the tree paths between the vertices of a set. Bonnet et al., pp. 3:26–3:27.

Formalization targets

The main target is the explicit excluded-clique-minor bound. With

ck=83(k+1)224k,g(t)=2(24t+1+2)2,f(t)=4cg(t)24cg(t)+2,c_k=\frac83(k+1)^2 2^{4k},\qquad g(t)=2(2^{4t+1}+2)^2,\qquad f(t)=4c_{g(t)}2^{4c_{g(t)}+2},ck​=38​(k+1)224k,g(t)=2(24t+1+2)2,f(t)=4cg(t)​24cg(t)​+2,

the goal says

Kt⪯̸G⟹tww⁡(G)≤f(t).K_t\not\preceq G\quad\Longrightarrow\quad \operatorname{tww}(G)\le f(t).Kt​⪯G⟹tww(G)≤f(t).

The +2+2+2 in g(t)g(t)g(t) follows the value defined and used in the proof of Theorem 6.3 on p. 3:26; the printed theorem statement has +1+1+1. The formal target uses the proof-supported bound and records the printed statement verbatim in its provenance. It does not assert the paper's asymptotic triple-exponential notation as a separate theorem. Bonnet et al., pp. 3:26–3:29.

The milestones follow the proof in order:

  1. twin-width is the maximum over connected components;
  2. a Lex-DFS exists from any vertex of a connected graph;
  3. Lemma 6.4: in any DFS, the vertices ai,ja_{i,j}ai,j​ of the proof lie on one branch of T\mathcal TT;
  4. Lemma 6.5: the enhancements Bj∗B^*_jBj∗​ avoid the paths Ai′A'_iAi′​, except the last ones;
  5. ∣V(H)∣≤α(H)ω(H)|V(H)|\le\alpha(H)\omega(H)∣V(H)∣≤α(H)ω(H) for intersection graphs of subtrees of a tree;
  6. the Helly property for subtrees of a tree;
  7. Lemma 6.6: each component of T[v]−{v}\mathcal T[v]-\{v\}T[v]−{v} meets at most two enhancements of the clique;
  8. for a connected KtK_tKt​-minor-free graph, the adjacency matrix in any Lex-DFS order is g(t)g(t)g(t)-grid free;
  9. grid free implies mixed free for 0/10/10/1 matrices;
  10. the graph form of Theorem 5.8 with the explicit constant applied on p. 3:29.

The Helly theorem is a reference to an existing open platform statement rather than a duplicate draft. Bonnet et al., pp. 3:22–3:29.

What the bound supplies

The theorem supplies a uniform contraction-width bound for every fixed excluded clique minor. In particular, the paper applies it to planar graphs through exclusion of K5K_5K5​. Its numerical value is far from small; the article remarks that the planar specialization has billions of digits. Without the result, excluding a clique minor alone would not yield the particular twin-width certificate required for those applications. Bonnet et al., p. 3:29.

The mathematical result is proved in the cited paper. This mission asks for machine-checked statements and, eventually, proofs of its graph and matrix claims. The reusable output includes a precise partition model of graph twin-width, interval-division models of mixed and grid minors, and a graph statement that converts a mixed-free adjacency order into a contraction-width bound. No machine-checked proof of these new target items is claimed here.

The main obstruction

A dense pattern of nonzero matrix zones does not directly give a graph minor. A matrix zone merely contains an edge somewhere between two vertex blocks. Those blocks need not induce connected subgraphs, and contracting them as if they were connected would use operations unavailable to graph minors. The proof also needs a vertex order with strong enough structure to control how the relevant connected sets overlap. These are the points that separate the grid-free milestone from a direct application of the matrix grid theorem. Bonnet et al., pp. 3:26–3:29.

The conversion from a mixed-free adjacency matrix to graph twin-width has a second issue: matrix twin-width permits separate row and column contractions, while a graph contraction merges one vertex partition and must keep the two matrix axes synchronized. Theorem 5.8 addresses the graph case. The milestone uses the explicit constant that the paper applies at the end of Theorem 6.3, while the printed Theorem 5.8 states an asymptotic bound. Bonnet et al., pp. 3:22–3:23, 3:29.

Formalization scope

Graphs are finite simple graphs, and twin-width and the number of division parts are natural numbers. The real constants ckc_kck​ and f(t)f(t)f(t) are compared with natural-number twin-width through ⌊f(t)⌋\lfloor f(t)\rfloor⌊f(t)⌋; the explicit value is positive. Vertex orders are equivalences with Fin⁡(n)\operatorname{Fin}(n)Fin(n). Matrix divisions use strictly increasing cut points, so every part is a nonempty consecutive interval. The adjacency matrix has a zero diagonal because simple graphs have no loops. The graph-minor predicate and the subtree representation come from existing published definitions with the same conventions; graph twin-width and the matrix objects are local draft definitions shared across this mission. Bonnet et al., pp. 3:9, 3:12, 3:17–3:19.

The grid-free milestone holds for every Lex-DFS order, as in the paper, and the existence of such an order is a separate milestone. Lemmas 6.4–6.6 are stated in the proof's context: a DFS, and vertices ai,ja_{i,j}ai,j​, bi,jb_{i,j}bi,j​ and blocks BjB_jBj​ with exactly the order and adjacency properties the page uses. Lemma 6.5 is posed for every block except the last. As printed it also covers the last block Bg(t)/2B_{g(t)/2}Bg(t)/2​, and for that block it can fail. A five-vertex example satisfies every hypothesis of the lemma and violates the printed conclusion; it is checked in Lean in the sanity file, and the formal item records the change. The theorem's bound is uniform over all finite graphs on the excluded-minor hypothesis. A twin-width definition using an unrestricted real infimum, a sequence that skips directly to one part, red degree counted over homogeneous pairs, or a grid minor over arbitrary noninterval partitions would change the claim. Components and empty graphs use the endpoint conventions stated above. The k≥1k\ge1k≥1 hypothesis in the graph conversion milestone rules out the degenerate zero-division case.

Contributions toward the core grid-free order claim, the symmetric contraction bridge, the component theorem, and the reusable minor and division lemmas are within scope. The running-time claims elsewhere in the paper and the pointer construction and the two Kt,tK_{t,t}Kt,t​-minor constructions inside the proof of Theorem 6.3 are not separate milestones: they are unnumbered and their statements need the pointers. They remain part of the grid-free milestone.

Selected references

  • É. Bonnet, E. J. Kim, S. Thomassé, and R. Watrigant, Twin-width I: Tractable FO Model Checking, Journal of the ACM 69(1), Article 3, 2021. DOI:10.1145/3486655.
  • A. Marcus and G. Tardos, Excluded permutation matrices and the Stanley–Wilf conjecture, Journal of Combinatorial Theory, Series A 107(1), 153–160, 2004. DOI:10.1016/j.jcta.2004.04.002.
  • F. Gavril, The intersection graphs of subtrees in trees are exactly the chordal graphs, Journal of Combinatorial Theory, Series B 16(1), 47–56, 1974. DOI:10.1016/0095-8956(74)90094-X.
17 thms1 active userReviewed
Combinatorics·Captain: mikedeng1

Three-Coloring and List Three-Coloring of Graphs Without Induced Paths on Seven Vertices 3: In a Lemma 11 Instance, Two Vertices of Xᵢ Have Neighbor Triples Joined by an EdgeResearch Paper

Motivation

A graph is 3-colorable if its vertices can be colored with three colors so that adjacent vertices get different colors. Deciding 3-colorability is NP-complete in general, and one of the standard ways to map the boundary between easy and hard cases is to forbid an induced subgraph. For the path PtP_tPt​ on ttt vertices, 3-colorability of P6P_6P6​-free graphs was known to be polynomial (Randerath and Schiermeyer, 2004, reference [21] of the paper), and Randerath, Schiermeyer and Tewes asked in 2002 whether the same holds for t=7t = 7t=7.

Bonomo, Chudnovsky, Maceli, Schaudt, Stein and Zhong (Combinatorica, 2017) answered this question. Their Theorem 1 states that the list version, in which every vertex vvv may only receive a color from a prescribed list L(v)⊆{1,2,3}L(v) \subseteq \{1,2,3\}L(v)⊆{1,2,3}, can be decided for P7P_7P7​-free graphs in time O(∣V(G)∣21(∣V(G)∣+∣E(G)∣))O(|V(G)|^{21}(|V(G)|+|E(G)|))O(∣V(G)∣21(∣V(G)∣+∣E(G)∣)). The algorithm first reduces the instance to a polynomial number of configurations described by a seed, and then calls Lemma 11 (p. 14) on each configuration. Lemma 11 rests on one structural fact about P7P_7P7​-free graphs, Claim 13 (p. 16), which is the goal of this mission.

Timeline:

  • 2002: Randerath, Schiermeyer and Tewes ask whether 3-colorability of P7P_7P7​-free graphs is polynomial.
  • 2004: Randerath and Schiermeyer give a polynomial algorithm for P6P_6P6​-free graphs.
  • 2017: Bonomo et al. give a polynomial algorithm for list 3-coloring of P7P_7P7​-free graphs; Claims 12 and 13 are its structural core.
  • As of the paper, no ttt is known for which 3-coloring PtP_tPt​-free graphs is NP-complete (p. 2); the case t=8t = 8t=8 is open.

Setting

All graphs are finite and simple. For a set SSS of vertices of a graph GGG, N(S)N(S)N(S) is the set of vertices outside SSS with a neighbor in SSS, and S‾=S∪N(S)\overline{S} = S \cup N(S)S=S∪N(S). A graph is PtP_tPt​-free if it has no induced subgraph isomorphic to PtP_tPt​.

A palette LLL assigns each vertex a list L(v)⊆{1,2,3}L(v) \subseteq \{1,2,3\}L(v)⊆{1,2,3}. A seed of (G,L)(G,L)(G,L) is a nonempty set SSS that induces a connected subgraph, is 2-dominating (every vertex is at distance at most 222 from SSS), and satisfies ∣L(v)∣=1|L(v)| = 1∣L(v)∣=1 for v∈Sv \in Sv∈S and ∣L(v)∣=2|L(v)| = 2∣L(v)∣=2 for v∈N(S)v \in N(S)v∈N(S).

The Lemma 11 setting consists of a graph GGG, a palette LLL and a seed SSS such that

  1. GGG is connected and P7P_7P7​-free;
  2. adjacent v∈Sv \in Sv∈S and w∈N(S)w \in N(S)w∈N(S) have disjoint lists;
  3. the set XXX of vertices with ∣L(v)∣=3|L(v)| = 3∣L(v)∣=3 is stable and anticomplete to V(G)∖(S‾∪X)V(G) \setminus (\overline{S} \cup X)V(G)∖(S∪X);
  4. no vertex of XXX has a connected neighborhood.

For i∈{1,2,3}i \in \{1,2,3\}i∈{1,2,3}, DiD_iDi​ is the set of v∈N(S)v \in N(S)v∈N(S) with L(v)={1,2,3}∖{i}L(v) = \{1,2,3\} \setminus \{i\}L(v)={1,2,3}∖{i}, and Ni(x)=N(x)∩DiN_i(x) = N(x) \cap D_iNi​(x)=N(x)∩Di​. For {i,j,k}={1,2,3}\{i,j,k\} = \{1,2,3\}{i,j,k}={1,2,3}, XiX_iXi​ is the set of x∈Xx \in Xx∈X for which Nj(x)N_j(x)Nj​(x) is not complete to Nk(x)N_k(x)Nk​(x), and for x∈Xix \in X_ix∈Xi​ the paper fixes non-adjacent nj(x)∈Nj(x)n_j(x) \in N_j(x)nj​(x)∈Nj​(x), nk(x)∈Nk(x)n_k(x) \in N_k(x)nk​(x)∈Nk​(x).

Formalization targets

Goal: Claim 13 (p. 16)

In the Lemma 11 setting, let {i,j,k}={1,2,3}\{i,j,k\} = \{1,2,3\}{i,j,k}={1,2,3}, x,y∈Xix, y \in X_ix,y∈Xi​, and let nj∈Nj(x)n_j \in N_j(x)nj​∈Nj​(x), nk∈Nk(x)n_k \in N_k(x)nk​∈Nk​(x) be non-adjacent. Then

∃ a∈{x,nj,nk}, ∃ b∈{y,nj(y),nk(y)}:ab∈E(G).\exists\, a \in \{x, n_j, n_k\},\ \exists\, b \in \{y, n_j(y), n_k(y)\}:\quad ab \in E(G).∃a∈{x,nj​,nk​}, ∃b∈{y,nj​(y),nk​(y)}:ab∈E(G).

Milestones

  1. Display (1), p. 6. For a seed SSS and two non-adjacent v,w∈N(S)v, w \in N(S)v,w∈N(S) there is an induced vvv–www path on at least 333 vertices whose inner vertices lie in SSS.
  2. p. 15. For d∈Did \in D_id∈Di​ and s∈S∩N(d)s \in S \cap N(d)s∈S∩N(d), L(s)={i}L(s) = \{i\}L(s)={i}.
  3. p. 15. No vertex of XXX has a neighbor in SSS.
  4. Claim 12, p. 15. For i≠ji \ne ji=j, ui,vi∈Diu_i, v_i \in D_iui​,vi​∈Di​, uj,vj∈Dju_j, v_j \in D_juj​,vj​∈Dj​ with {ui,vi,uj,vj}\{u_i, v_i, u_j, v_j\}{ui​,vi​,uj​,vj​} stable, there is an induced path PPP with ends a,ba, ba,b among them such that {a,b}≠{ui,uj}\{a,b\} \ne \{u_i,u_j\}{a,b}={ui​,uj​}, {a,b}≠{vi,vj}\{a,b\} \ne \{v_i,v_j\}{a,b}={vi​,vj​}, the interior P∗P^*P∗ lies in SSS, and P∗P^*P∗ is anticomplete to the other two vertices.

Significance

Claim 13 says that the neighbor triples {x,nj(x),nk(x)}\{x, n_j(x), n_k(x)\}{x,nj​(x),nk​(x)} of the vertices of XiX_iXi​ pairwise touch. This is what allows the proof of Lemma 11 (Claims 14 to 17) to refine the palette of all of XiX_iXi​ simultaneously with only polynomially many guesses: without it, the vertices of XiX_iXi​ could interact independently and the number of cases would be exponential. Claim 12 and display (1) are the general tools for building long induced paths through a seed, used throughout §3.1.

The paper's result is proved and published. None of it is formalized in a proof assistant as far as the platform's index shows. This mission produces a machine-checked version of the structural core of Lemma 11, and the definitions it needs (seeds, the Lemma 11 setting, the sets DiD_iDi​ and Ni(x)N_i(x)Ni​(x)) are reusable for the rest of the paper's argument.

Difficulty

The proof of Claim 13 is short on paper, but it assembles an induced P7P_7P7​ from three pieces: the path through SSS given by Claim 12 and the two short paths nj−x−nkn_j - x - n_knj​−x−nk​ and nj(y)−y−nk(y)n_j(y) - y - n_k(y)nj​(y)−y−nk​(y). Showing that the union is induced requires every non-adjacency between the pieces, and these come from different hypotheses: the stability of XXX, the fact that XXX has no neighbor in SSS, the list structure of the vertices of SSS next to DiD_iDi​, and part (c) of Claim 12. Claim 12 itself is an extremal argument (choose the neighbor closest to an end of a path), which in a formal setting requires explicit manipulation of paths as sequences of vertices: splitting, reversing and concatenating them while keeping them induced. The obvious approach of taking any shortest path through SSS between two of the four vertices fails, because such a path can have interior vertices adjacent to the other two.

Formalization scope

Graphs are SimpleGraph V on a finite type V with decidable equality; vertex sets are Finset V. The colors 1,2,31,2,31,2,3 are 0,1,2 in Fin 3, and a palette is V → Finset (Fin 3). P7P_7P7​-freeness is the negation of Mathlib's induced containment IsIndContained of pathGraph 7; ordinary subgraph containment would be a different, much stronger condition. Connectedness uses Mathlib's Connected, which includes nonemptiness. An induced path is the published platform definition StrongPerfectGraph.Main.IsInducedPath (a nonempty list of distinct vertices in which two vertices are adjacent exactly when consecutive); its ends are the first and last list entries and its interior P∗P^*P∗ is the list with both ends removed.

The Lemma 11 setting is a structure Lemma11Hyp G L S bundling its hypotheses. Lemma 11's own conclusion, a decision procedure with running time O(∣V(G)∣9(∣V(G)∣+∣E(G)∣))O(|V(G)|^9(|V(G)|+|E(G)|))O(∣V(G)∣9(∣V(G)∣+∣E(G)∣)), and the running time of Theorem 1 are not formalized, since there is no machine model. The paper's arbitrary choice of nj(y),nk(y)n_j(y), n_k(y)nj​(y),nk​(y) is replaced by quantification over every non-adjacent pair mj∈Nj(y)m_j \in N_j(y)mj​∈Nj​(y), mk∈Nk(y)m_k \in N_k(y)mk​∈Nk​(y), which is Claim 13 for every possible choice. The case x=yx = yx=y is allowed.

One hypothesis is added to Claim 12: the four vertices ui,vi,uj,vju_i, v_i, u_j, v_jui​,vi​,uj​,vj​ are pairwise distinct, and the ends a,ba, ba,b of the path are required to be distinct. The paper's proof and its use in Claim 13 take the four vertices distinct; without the assumption the claim can fail, and without a≠ba \ne ba=b a one-vertex path would satisfy it vacuously.

The following trivializing readings are excluded: dropping the overline in "anticomplete to V(G)∖(S‾∪X)V(G) \setminus (\overline{S} \cup X)V(G)∖(S∪X)" (which, together with milestone 3, would force XXX to be isolated); a one-vertex path or repeated vertices in Claim 12; and subgraph instead of induced-subgraph containment for P7P_7P7​-freeness. The setting is satisfiable with X≠∅X \neq \emptysetX=∅: the 5-cycle y−a−s−t−b−yy - a - s - t - b - yy−a−s−t−b−y with S={s,t}S = \{s,t\}S={s,t}, L(s)={2}L(s) = \{2\}L(s)={2}, L(t)={3}L(t) = \{3\}L(t)={3}, L(a)={1,3}L(a) = \{1,3\}L(a)={1,3}, L(b)={1,2}L(b) = \{1,2\}L(b)={1,2}, L(y)={1,2,3}L(y) = \{1,2,3\}L(y)={1,2,3} satisfies every hypothesis, with y∈X1y \in X_1y∈X1​.

Proofs of the milestones and the goal are welcome, as are general lemmas about induced paths as lists (subpaths, concatenation through a connected set) that would serve other graph-theory missions.

Selected references

  • F. Bonomo, M. Chudnovsky, P. Maceli, O. Schaudt, M. Stein and M. Zhong, Three-coloring and list three-coloring of graphs without induced paths on seven vertices, Combinatorica 38 (2018) 779–801 (OnlineFirst 2017). https://doi.org/10.1007/s00493-017-3553-8
  • B. Randerath and I. Schiermeyer, 3-Colorability ∈ P for P₆-free graphs, Discrete Applied Mathematics 136 (2004) 299–313 (reference [21] of the paper).
  • B. Randerath, I. Schiermeyer and M. Tewes, Three-colorability and forbidden subgraphs. II: polynomial algorithms, Discrete Mathematics 251 (2002) 137–153 (reference [22] of the paper).
  • E. Camby and O. Schaudt, A new characterization of PkP_kPk​-free graphs, Algorithmica (reference [2] of the paper; the source of Theorem 4, which the other missions of this series use).
9 thms1 active userReviewed
CombinatoricsMathematical Logic·Captain: mikedeng1

Twin-width I: Tractable FO Model Checking 5: First-Order Interpretations of a Class of Graphs of Bounded Twin-width Have Bounded Twin-widthResearch Paper

Why interpretations

Twin-width is a graph parameter introduced by Bonnet, Kim, Thomassé and Watrigant in Twin-width I: Tractable FO Model Checking (J. ACM 69(1), Article 3, 2021; arXiv:2004.14789). Its main algorithmic consequence is that first-order model checking is fixed-parameter tractable on graphs given with a contraction sequence of bounded width. A parameter of this kind is useful only if the classes it bounds are closed under the constructions that arise in practice. Many such constructions are first-order definable: the complement of a graph, its square (join vertices at distance at most two), the map graph of a planar map. Section 8 of the paper shows that twin-width is stable under every such construction: a first-order interpretation of a class of bounded twin-width has bounded twin-width. Bounded VC-dimension, by contrast, is not preserved (interval graphs have VC-dimension at most two but interpret every graph, p. 3:41).

The result has since become the basis of the "model-theoretic" view of twin-width: bounded twin-width classes are closed under first-order transductions (Theorem 8.1 of the paper), and later work identifies them, among ordered structures, with the monadically dependent classes (Twin-width IV, Bonnet, Giocanti, Ossona de Mendez, Simon, Thomassé and Toruńczyk).

Setting

Let GGG be a finite simple graph on a vertex set VVV.

Two vertex sets X,YX,YX,Y are homogeneous if all pairs in X×YX\times YX×Y are edges or none is. For a partition P\mathcal PP of VVV, the red graph GPG_{\mathcal P}GP​ has the parts as vertices, two distinct parts being adjacent when they are not homogeneous. P\mathcal PP is a ddd-partition if GPG_{\mathcal P}GP​ has maximum degree at most ddd. The graph GGG has twin-width at most ddd, tww⁡(G)≤d\operatorname{tww}(G)\le dtww(G)≤d, if there is a sequence of ddd-partitions from the partition into singletons to a partition with at most one part, each obtained from the previous one by merging two parts (p. 3:32).

A prenex formula of depth ℓ\ellℓ with two free variables is

φ(x,y)=Q1x1 Q2x2⋯Qℓxℓ φ∗,Qi∈{∀,∃},\varphi(x,y)=Q_1x_1\,Q_2x_2\cdots Q_\ell x_\ell\ \varphi^*,\qquad Q_i\in\{\forall,\exists\},φ(x,y)=Q1​x1​Q2​x2​⋯Qℓ​xℓ​ φ∗,Qi​∈{∀,∃},

with φ∗\varphi^*φ∗ a Boolean combination of atoms u=vu=vu=v and E(u,v)E(u,v)E(u,v), u,v∈{x1,…,xℓ,x,y}u,v\in\{x_1,\dots,x_\ell,x,y\}u,v∈{x1​,…,xℓ​,x,y}. The interpretation φ(G)\varphi(G)φ(G) is the graph on VVV in which distinct u,vu,vu,v are adjacent iff G⊨φ(u,v)∧φ(v,u)G\models\varphi(u,v)\wedge\varphi(v,u)G⊨φ(u,v)∧φ(v,u). For a class G\mathcal GG, φ(G)\varphi(\mathcal G)φ(G) is the class of all induced subgraphs of the graphs φ(G)\varphi(G)φ(G), G∈GG\in\mathcal GG∈G.

The proof works with morphism-trees. The complete tree MTk(V)MT_k(V)MTk​(V) has as nodes all tuples (v1,…,vi)(v_1,\dots,v_i)(v1​,…,vi​) of vertices with i≤ki\le ki≤k, the parent of a tuple being its prefix. Two sibling nodes are equivalent in (G,P)(G,\mathcal P)(G,P) if an automorphism of the tree swaps them, preserving equalities and adjacencies among the entries of every tuple and the part of P\mathcal PP of every entry. A reduction deletes one of two equivalent siblings with its subtree, repeatedly; a reduct admits no further reduction. Two vertices u,u′u,u'u,u′ are adjacent in Eℓ+2(G,P)E_{\ell+2}(G,\mathcal P)Eℓ+2​(G,P) when (u)(u)(u) and (u′)(u')(u′) are equivalent in some reduction of MTℓ+2(V)MT_{\ell+2}(V)MTℓ+2​(V), and Iℓ+2(G,P)I_{\ell+2}(G,\mathcal P)Iℓ+2​(G,P) is the partition into the connected components of Eℓ+2(G,P)E_{\ell+2}(G,\mathcal P)Eℓ+2​(G,P).

Formalization targets

Goal: Theorem 8.3 (graphs)

For every prenex formula φ(x,y)\varphi(x,y)φ(x,y) and every ddd there is DDD such that

tww⁡(G)≤d ⟹ tww⁡(φ(G)[S])≤Dfor every graph G and every S⊆V(G).\operatorname{tww}(G)\le d\ \Longrightarrow\ \operatorname{tww}\bigl(\varphi(G)[S]\bigr)\le D\quad\text{for every graph } G \text{ and every } S\subseteq V(G).tww(G)≤d ⟹ tww(φ(G)[S])≤Dfor every graph G and every S⊆V(G).

The bound DDD depends on φ\varphiφ and ddd only. No explicit function is asserted: the paper gives none (the proof goes through a tower-type bound on reducts).

Milestones

  1. Lemma 5.2 in graph form (p. 3:18, applied on p. 3:43): rrr-refining hhh-partitions from the finest to the coarsest partition give tww⁡(G)≤r(h+1)\operatorname{tww}(G)\le r(h+1)tww(G)≤r(h+1).
  2. Lemma 7.3 (p. 3:32): reducts of MTℓ(G)MT_\ell(G)MTℓ​(G) have size bounded by a function of ℓ\ellℓ.
  3. Lemmas 7.11–7.14 (pp. 3:35–3:38): restriction to connected tuples rooted at a part and the pruned shuffle commute with reductions; the pruned shuffle of all MTℓ(G,P,X)MT_\ell(G,\mathcal P,X)MTℓ​(G,P,X) is MTℓ(G,P)MT_\ell(G,\mathcal P)MTℓ​(G,P), and that of reducts is a reduction of it.
  4. Lemma 8.4 (p. 3:41) and its extension to reductions (p. 3:42): equivalent nodes (u,v),(u,v′)(u,v),(u,v')(u,v),(u,v′) of MTℓ+2MT_{\ell+2}MTℓ+2​ satisfy the same prenex formulas of depth ℓ\ellℓ.
  5. Iℓ+2I_{\ell+2}Iℓ+2​ refines P\mathcal PP and is monotone under coarsening (p. 3:42).
  6. Lemma 8.5 (p. 3:42): a part of a ddd-partition meets boundedly many components of Eℓ+2(G,P)E_{\ell+2}(G,\mathcal P)Eℓ+2​(G,P).
  7. Lemma 8.6 (p. 3:43): parts of Iℓ+2(G,P)I_{\ell+2}(G,\mathcal P)Iℓ+2​(G,P) inside parts at distance at least 3ℓ+23^{\ell+2}3ℓ+2 in GPG_{\mathcal P}GP​ are homogeneous in φ(G)\varphi(G)φ(G).

Significance

Theorem 8.3 turns every first-order definable graph construction into a source of bounded-twin-width classes: squares and fixed powers of planar graphs, complements, map graphs (as transductions, via Theorem 8.1), kkk-planar graphs and bounded-degree string graphs. Combined with the model-checking algorithm it gives FPT first-order model checking on all these classes, provided a contraction sequence is available. The theorem is also the first step of the model-theoretic characterisations of bounded twin-width in the later papers of the series.

The result is proved in the paper; no machine-checked proof of it, or of any twin-width bound, exists to our knowledge, and Mathlib has no notion of twin-width or of first-order interpretations of graphs. The mission asks for a formal proof of the graph case. The morphism-tree milestones are also the combinatorial core of the paper's linear-time model-checking algorithm (Theorems 7.5, 7.15), whose running-time statements are not part of this mission.

Difficulty

The obvious approach refines the ddd-partitions Pi\mathcal P_iPi​ of GGG by the "type" of each vertex with respect to φ\varphiφ. The number of such types is not bounded: a vertex's behaviour depends on the whole graph through the quantifiers. The paper replaces types by the components of Eℓ+2(G,Pi)E_{\ell+2}(G,\mathcal P_i)Eℓ+2​(G,Pi​), defined through reductions of morphism-trees that respect Pi\mathcal P_iPi​. Two facts must then be shown, and neither is local in an obvious way. First, the number of components inside a part is bounded (Lemma 8.5); this needs the pruned shuffle, which assembles a reduction of the whole tree MTℓ+2(G,P)MT_{\ell+2}(G,\mathcal P)MTℓ+2​(G,P) from reductions of the local trees MTℓ+2(G,P,X)MT_{\ell+2}(G,\mathcal P,X)MTℓ+2​(G,P,X), and the fact that it commutes with reductions (Lemma 7.12), where an automorphism of one local tree must be extended to the shuffle without breaking adjacency between non-homogeneous parts. Second, far-apart components are homogeneous in φ(G)\varphi(G)φ(G) (Lemma 8.6), which needs the evaluation argument of Lemma 8.4 inside a reduction of the partitioned tree.

Formalization scope

  • Graphs are SimpleGraph V on a Fintype with decidable equality. Twin-width is the predicate TwinWidthLE G d in the paper's partition form (pp. 3:12, 3:32); trigraphs are not formalized, and twin-width is never an infimum, so no junk value of sInf can make a bound vacuous. Merge steps merge exactly two distinct parts, so a sequence cannot jump to one part.
  • Formulas are Mathlib's FirstOrder.Language.graph.Formula (Fin 2), evaluated in G.structure; variable 0 is xxx, variable 1 is yyy. Prenex is Mathlib's IsPrenex; quantifier depth qdepth is defined in the mission (∃=¬∀¬\exists=\neg\forall\neg∃=¬∀¬ counts once).
  • The interpretation requires u≠vu\neq vu=v for an edge (a simple graph has no loops). The goal is the graph case of Theorem 8.3, which the page states for augmented binary structures; Section 8 itself says it works with undirected graphs (p. 3:40).
  • In the goal the bound DDD is chosen after φ\varphiφ and ddd and before the graph: a bound chosen after the graph would make the statement trivial. The hereditary closure (induced subgraphs on every SSS) is part of the statement.
  • Morphism-trees are sets of tuples (Set (List V)), the paper's identification of a node with its current path (p. 3:30), which is exact for reductions of complete trees. Equivalent siblings are distinct. Distances in GPG_{\mathcal P}GP​ are in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞}, so disconnected parts are far apart. The sequence graph uses 0-based indices, hence the exponent 3ℓ−(k+1)3^{\ell-(k+1)}3ℓ−(k+1). The pruned shuffle is defined by its characterisation (each component of the sequence graph of a tuple lies in the tree of its local root). Lemmas of §7.3 carry its standing assumption ℓ>0\ell>0ℓ>0.
  • Lemma 5.2 is posed for graphs with the bound r(h+1)r(h+1)r(h+1) (the matrix lemma's rtrtrt with t=h+1t=h+1t=h+1).
  • Not in the mission: Theorem 8.1 (transductions), Lemma 8.2 and Lemma 5.1 (augmented binary structures), and all running-time statements (Theorems 1.1, 7.1, 7.5, 7.15, Corollaries 7.16–7.17).

Reusable beyond this mission: the partition form of twin-width, quantifier depth of first-order formulas, and the morphism-tree machinery (reductions, pruned shuffles), which formalizes the algorithmic core of first-order model checking on bounded twin-width. Contributions welcome: proofs of the §7 lemmas, of Lemma 8.4 by induction on the quantifier prefix, and the closure of twin-width under induced subgraphs.

Selected references

  • É. Bonnet, E. J. Kim, S. Thomassé, R. Watrigant, Twin-width I: Tractable FO Model Checking, J. ACM 69(1), Article 3, 2021. https://doi.org/10.1145/3486655 (arXiv:2004.14789, https://arxiv.org/abs/2004.14789)
  • É. Bonnet, U. Giocanti, P. Ossona de Mendez, P. Simon, S. Thomassé, S. Toruńczyk, Twin-width IV: Ordered Graphs and Matrices, J. ACM, 2024. https://arxiv.org/abs/2102.03117
15 thms1 active userReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems V: On Percolated Graphs Sampled From a Lipschitz Graphon the Mean-Square Error Is at Most κ(q)/(nβₙ)^{1/q}Research Paper

Motivation

Large systems of interacting diffusions — neurons, oscillators, agents in a network game, particles in a queueing or load-balancing model — are usually analysed through their mean-field limit: as the population grows, each particle interacts with the empirical law of the others, and the system is replaced by a single nonlinear (McKean–Vlasov) equation. The classical theory assumes that every particle interacts with every other one with the same strength. Real networks are heterogeneous and sparse: particle iii interacts with particle jjj only through an edge of a random graph, and the number of neighbours of a particle is much smaller than the population size.

Bayraktar, Chakraborty and Wu (Ann. Appl. Probab. 33(5), 2023) describe the heterogeneity by a graphon and prove laws of large numbers and rates of convergence for such systems. This mission is the fifth of a series formalizing that paper. Its subject is the paper's quantitative result for not-so-dense graphs: when the edges are sampled from a Lipschitz graphon and then thinned by a sparsity factor βn\beta_nβn​, the average mean-square distance between each particle and its graphon limit is of order (nβn)−1/q(n\beta_n)^{-1/q}(nβn​)−1/q for every q>1q>1q>1.

Setting

Let I=[0,1]I=[0,1]I=[0,1] with Lebesgue measure. A graphon is a measurable symmetric function G:I×I→[0,1]G:I\times I\to[0,1]G:I×I→[0,1]. Fix a horizon T>0T>0T>0 and write Cd=C([0,T]:Rd)\mathcal C_d=C([0,T]:\mathbb R^d)Cd​=C([0,T]:Rd) with the norm ∥x∥∗,T=sup⁡0≤s≤T∣xs∣\|x\|_{*,T}=\sup_{0\le s\le T}|x_s|∥x∥∗,T​=sup0≤s≤T​∣xs​∣.

One probability space carries, for every label u∈Iu\in Iu∈I, an initial state Xu(0)X_u(0)Xu​(0) with law μu(0)\mu_u(0)μu​(0) and a standard ddd-dimensional Brownian motion BuB_uBu​, the whole family {Xu(0),Bu:u∈I}\{X_u(0),B_u:u\in I\}{Xu​(0),Bu​:u∈I} being independent. The limit system (4.2) is a continuum of diffusions, one per label:

Xu(t)=Xu(0)+∫0t ⁣∫I ⁣∫Rdb(Xu(s),x) G(u,v) μv,s(dx) dv ds+∫0tσ(Xu(s)) dBu(s),μu,t=L(Xu(t)).X_u(t)=X_u(0)+\int_0^t\!\int_I\!\int_{\mathbb R^d}b(X_u(s),x)\,G(u,v)\,\mu_{v,s}(dx)\,dv\,ds+\int_0^t\sigma(X_u(s))\,dB_u(s),\qquad \mu_{u,t}=\mathcal L(X_u(t)).Xu​(t)=Xu​(0)+∫0t​∫I​∫Rd​b(Xu​(s),x)G(u,v)μv,s​(dx)dvds+∫0t​σ(Xu​(s))dBu​(s),μu,t​=L(Xu​(t)).

The particles are independent but not identically distributed; their laws are coupled through the dvdvdv-integral.

The not-so-dense nnn-particle system (4.1) uses the same initial states and Brownian motions at the labels i/ni/ni/n:

Xin(t)=Xi/n(0)+∫0t1nβn∑j=1nξijn b(Xin(s),Xjn(s)) ds+∫0tσ(Xin(s)) dBi/n(s),i=1,…,n.X^n_i(t)=X_{i/n}(0)+\int_0^t\frac{1}{n\beta_n}\sum_{j=1}^n\xi^n_{ij}\,b(X^n_i(s),X^n_j(s))\,ds+\int_0^t\sigma(X^n_i(s))\,dB_{i/n}(s),\qquad i=1,\dots,n.Xin​(t)=Xi/n​(0)+∫0t​nβn​1​j=1∑n​ξijn​b(Xin​(s),Xjn​(s))ds+∫0t​σ(Xin​(s))dBi/n​(s),i=1,…,n.

The edges are percolated samples of GGG (Condition 4.3): ξijn=ξjin∼Bernoulli(βnG(i/n,j/n))\xi^n_{ij}=\xi^n_{ji}\sim\mathrm{Bernoulli}(\beta_nG(i/n,j/n))ξijn​=ξjin​∼Bernoulli(βn​G(i/n,j/n)), independent for 1≤i≤j≤n1\le i\le j\le n1≤i≤j≤n and independent of the noise. A particle has about nβnn\beta_nnβn​ neighbours, which is why the interaction is scaled by 1/(nβn)1/(n\beta_n)1/(nβn​).

The standing hypotheses are Condition 4.1 — u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) measurable with bounded second moments; bbb bounded and Lipschitz; σ\sigmaσ bounded, Lipschitz and invertible with bounded inverse; βn∈(0,1]\beta_n\in(0,1]βn​∈(0,1] and nβn→∞n\beta_n\to\inftynβn​→∞ — and Condition 2.3: on finitely many intervals I1,…,INI_1,\dots,I_NI1​,…,IN​ covering III, the map u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) is Lipschitz in the Wasserstein distance W2W_2W2​, and GGG is Lipschitz on every block Ii×IjI_i\times I_jIi​×Ij​.

Formalization targets

Goal: Theorem 4.2

Under Conditions 2.3, 4.1 and 4.3, for each q∈(1,∞)q\in(1,\infty)q∈(1,∞) there is κ(q)∈(0,∞)\kappa(q)\in(0,\infty)κ(q)∈(0,∞) with

1n∑i=1nE∥Xin−Xi/n∥∗,T2≤κ(q)(nβn)1/qfor all n∈N.\frac1n\sum_{i=1}^n\mathbb E\big\|X^n_i-X_{i/n}\big\|_{*,T}^2\le\frac{\kappa(q)}{(n\beta_n)^{1/q}}\qquad\text{for all }n\in\mathbb N .n1​i=1∑n​E​Xin​−Xi/n​​∗,T2​≤(nβn​)1/qκ(q)​for all n∈N.

The constant depends on qqq and on the data, never on nnn. The bound is on the average over particles, not on the maximum.

Milestones

  1. Lemma 7.1: all moments of the coupling error are bounded uniformly in nnn and iii.
  2. (7.10): E[ξijn∣Xjn(s)−Xj/n(s)∣2]≤(2E∣Xjn(s)−Xj/n(s)∣2+κ(q)(nβn)−1/q)βnG(i/n,j/n)\mathbb E[\xi^n_{ij}|X^n_j(s)-X_{j/n}(s)|^2]\le\big(2\mathbb E|X^n_j(s)-X_{j/n}(s)|^2+\kappa(q)(n\beta_n)^{-1/q}\big)\beta_nG(i/n,j/n)E[ξijn​∣Xjn​(s)−Xj/n​(s)∣2]≤(2E∣Xjn​(s)−Xj/n​(s)∣2+κ(q)(nβn​)−1/q)βn​G(i/n,j/n).
  3. (7.15): the interaction error Rsn,2R^{n,2}_sRsn,2​ is at most (κ/n)∑jE∣Xjn(s)−Xj/n(s)∣2+κ(q)(nβn)−1/q(\kappa/n)\sum_j\mathbb E|X^n_j(s)-X_{j/n}(s)|^2+\kappa(q)(n\beta_n)^{-1/q}(κ/n)∑j​E∣Xjn​(s)−Xj/n​(s)∣2+κ(q)(nβn​)−1/q.
  4. Theorem 2.1(b) for (4.2): W2,T(μu,μv)≤κ∣u−v∣W_{2,T}(\mu_u,\mu_v)\le\kappa|u-v|W2,T​(μu​,μv​)≤κ∣u−v∣ for u,vu,vu,v in one interval IiI_iIi​.
  5. §7.3 display: the discretization error Rsn,4R^{n,4}_sRsn,4​ is at most κ/n2\kappa/n^2κ/n2.

Significance

Theorem 4.2 gives an explicit rate for the convergence of a particle system on a sparse random graph to its graphon limit. Its companion, Theorem 4.1 (mission IV), proves the law of large numbers under weaker hypotheses but without a rate. The rate says that the effective sample size is nβnn\beta_nnβn​, the typical degree, not nnn: graphs with βn→0\beta_n\to0βn​→0 still approximate the graphon system as long as the degree grows. Results of this type justify replacing a large heterogeneous network by its graphon limit in control and game problems, such as graphon mean-field games, where the error of the replacement must be quantified.

The paper proves the theorem by hand; none of its results has a machine-checked proof. The formalization produces precise statements of the not-so-dense model, of the estimates of §7.1 that hold for any sparsity, and of the discretization bound, each usable on its own. The milestones (7.10) and (7.15) are the paper's key estimates for sparse interactions and are shared with mission IV.

Difficulty

The obvious approach is to compare each XinX^n_iXin​ with Xi/nX_{i/n}Xi/n​ by Gronwall's inequality, after bounding the difference of drifts. That step fails at the interaction term: the edge ξijn\xi^n_{ij}ξijn​ and the error Xjn−Xj/nX^n_j-X_{j/n}Xjn​−Xj/n​ are dependent, because the edge enters the dynamics of particle jjj. A plain Cauchy–Schwarz bound loses a factor 1/βn1/\beta_n1/βn​ and gives nothing when βn→0\beta_n\to0βn​→0. Estimates (7.10) and (7.15) quantify how weak this dependence is, and they are where the hypotheses that σ\sigmaσ is invertible with bounded inverse and that nβn→∞n\beta_n\to\inftynβn​→∞ are needed. The discretization term needs the Lipschitz regularity of the limit laws in the label (Theorem 2.1(b)), itself a fixed-point estimate on a space of measure-valued maps.

Formalization scope

The Lean development uses Mathlib's unitInterval for III, Fin d → ℝ with its sup norm for Rd\mathbb R^dRd (all constants are existential, so the choice of norm does not matter), time ℝ≥0, and the path space Cd\mathcal C_dCd​ as continuous maps on [0,T][0,T][0,T] with the sup norm. Brownian motions, Itô integrals and Itô processes are those of the published Peng1990.SMP.Stochastic; Wasserstein distances are the published WassersteinDRO.Duality.wassersteinDistance, valued in [0,∞][0,\infty][0,∞]. Particle iii is Lean's i : Fin n with label (i+1)/n(i+1)/n(i+1)/n.

Committed conventions:

  • A solution of (4.2) has continuous paths, AE-measurable path maps, and path laws in the class M\mathcal MM (measurable in uuu, bounded second moments); each XuX_uXu​ is a strong solution for the filtration of (Xu(0),Bu)(X_u(0),B_u)(Xu​(0),Bu​). A solution of (4.1) is a strong solution for the filtration of the edges, the initial states and the Brownian motions of the nnn labels.
  • The edges are {0,1}\{0,1\}{0,1}-valued, symmetric, with the diagonal included, mutually independent over i≤ji\le ji≤j, and independent of the σ-algebra of the noise.
  • Condition 4.1(d) is required for n≥1n\ge1n≥1; the invertibility of σ\sigmaσ is that of a d×dd\times dd×d matrix; the intervals of Condition 2.3 are order-connected subsets of III.
  • Expectations of nonnegative quantities are lower Lebesgue integrals, so no expectation is silently 000 for a non-integrable integrand.

The solution predicates are not vacuous: the zero-coefficient system Xu(t)=Xu(0)X_u(t)=X_u(0)Xu​(t)=Xu​(0) satisfies both, Condition 4.1 and Condition 4.3 are satisfiable, and every statement quantifies over solutions rather than assuming the conclusion. The constant κ(q)\kappa(q)κ(q) is chosen after qqq and before nnn, so it cannot depend on nnn, and the bound is the paper's average over iii, not a weaker or a stronger form. At n=0n=0n=0 both sides of every rate are 000.

The paper's proofs are in its §5 (Theorem 2.1), §7.1 (Lemma 7.1, (7.10), (7.15)) and §7.3 (Theorem 4.2). Stochastic analysis at that level — Burkholder–Davis–Gundy and Rosenthal inequalities, Girsanov's theorem for these SDEs, Gronwall arguments for measure-valued fixed points — is mostly absent from Mathlib, and contributions of such infrastructure are welcome, as are proofs of any milestone.

Selected references

  • E. Bayraktar, S. Chakraborty, R. Wu, Graphon mean field systems, Ann. Appl. Probab. 33(5):3587–3619, 2023. https://doi.org/10.1214/22-AAP1901
  • L. Lovász, Large Networks and Graph Limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
  • B. Bollobás, C. Borgs, J. Chayes, O. Riordan, Percolation on dense graph sequences, Ann. Probab. 38(1):150–183, 2010. https://doi.org/10.1214/09-AOP478
  • S. Peng, A general stochastic maximum principle for optimal control problems, SIAM J. Control Optim. 28(4):966–979, 1990. https://doi.org/10.1137/0328054
11 thms1 active userReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems IV: Law of Large Numbers for Particle Systems on Percolated Not-So-Dense Graphs With nβₙ → ∞Research Paper

Motivation

Large networks are often modeled through a graphon, a measurable function that records how likely two labeled vertices are to interact. For interacting diffusions, this lets one ask whether a finite random network has a predictable population limit. Bayraktar, Chakraborty and Wu study that question for graphon weighted systems, including graphs obtained by independently retaining potential edges with a sparsity parameter βn\beta_nβn​. Their not-so-dense regime allows βn\beta_nβn​ to decrease, provided the typical interaction scale nβnn\beta_nnβn​ diverges. Theorem 4.1 is the paper's law of large numbers for this regime; its two conclusions distinguish convergence of the population's empirical law from mean-square closeness of each particle to a coupled graphon limit. Bayraktar, Chakraborty and Wu, 2023.

The authors present three main convergence results: continuity of the continuum graphon system in Section 2, a dense-system law of large numbers in Section 3, and the percolated-system law of large numbers in Section 4. Section 7 establishes the latter. The graphon and cut-metric framework used here follows the graph-limit treatment cited by the paper. Bayraktar, Chakraborty and Wu, 2023; Lovász, 2012.

Setting

Let I=[0,1]I=[0,1]I=[0,1]. A graphon is a symmetric measurable G:I2→[0,1]G:I^2\to[0,1]G:I2→[0,1]. The cut norm ∥W∥□\|W\|_\square∥W∥□​ measures the largest absolute integral of WWW over a measurable rectangle S×US\times US×U. A sequence of step graphons GnG_nGn​ approaches GGG when ∥Gn−G∥□→0\|G_n-G\|_\square\to0∥Gn​−G∥□​→0. The related operator norm ∥W∥∞→1\|W\|_{\infty\to1}∥W∥∞→1​ measures the effect of WWW on bounded measurable functions. Remark 2.1 says cut convergence implies convergence in this operator norm. Bayraktar, Chakraborty and Wu, 2023, Remark 2.1.

Fix a dimension d≥1d\ge1d≥1 and horizon T>0T>0T>0. A label u∈Iu\in Iu∈I has an initial state Xu(0)∈RdX_u(0)\in\mathbb R^dXu​(0)∈Rd and an independent ddd-dimensional Brownian motion BuB_uBu​. The initial states are independent of each other and of the Brownian motions. The limiting graphon particle XuX_uXu​ is a continuous path whose drift averages b(Xu(s),x)G(u,v)b(X_u(s),x)G(u,v)b(Xu​(s),x)G(u,v) against the time-sss laws of all labels vvv, and whose diffusion matrix is σ(Xu(s))\sigma(X_u(s))σ(Xu​(s)). The time-sss law at label vvv is μv,s=L(Xv(s))\mu_{v,s}=\mathcal L(X_v(s))μv,s​=L(Xv​(s)). The family of path laws μu=L(Xu)\mu_u=\mathcal L(X_u)μu​=L(Xu​) is measurable in uuu and has uniformly bounded second moments. Bayraktar, Chakraborty and Wu, 2023, (4.2).

For nnn particles, label iii means i/ni/ni/n, with 1≤i≤n1\le i\le n1≤i≤n. Each undirected edge weight ξijn=ξjin\xi^n_{ij}=\xi^n_{ji}ξijn​=ξjin​ is Bernoulli with mean βnGn(i/n,j/n)\beta_nG_n(i/n,j/n)βn​Gn​(i/n,j/n); diagonal weights are included. These edge variables are independent of the Brownian motions and initial states. The finite particle XinX_i^nXin​ uses the same initial state and Brownian motion as Xi/nX_{i/n}Xi/n​. Its drift is the sum of ξijnb(Xin,Xjn)\xi^n_{ij}b(X_i^n,X_j^n)ξijn​b(Xin​,Xjn​) divided by nβnn\beta_nnβn​; its diffusion is σ(Xin)\sigma(X_i^n)σ(Xin​). The coefficients bbb and σ\sigmaσ are bounded and Lipschitz, and σ(x)\sigma(x)σ(x) is invertible with uniformly bounded inverse. Initial laws have uniformly bounded second moments, 0<βn≤10<\beta_n\le10<βn​≤1, and nβn→∞n\beta_n\to\inftynβn​→∞. Bayraktar, Chakraborty and Wu, 2023, Conditions 4.1–4.2.

Formalization targets

Let Cd=C([0,T];Rd)\mathcal C_d=C([0,T];\mathbb R^d)Cd​=C([0,T];Rd), μn=n−1∑i=1nδXin\mu^n=n^{-1}\sum_{i=1}^{n}\delta_{X_i^n}μn=n−1∑i=1n​δXin​​, and μˉ=∫Iμu du\bar\mu=\int_I\mu_u\,duμˉ​=∫I​μu​du. The main target, Theorem 4.1, first asserts convergence in probability in the weak topology on P(Cd)\mathcal P(\mathcal C_d)P(Cd​):

μn→Pμˉ.\mu^n\xrightarrow{\mathbb P}\bar\mu.μnP​μˉ​.

This uses Condition 2.2(a): on each of finitely many intervals covering III, the initial law varies continuously in Wasserstein distance W2W_2W2​. If the graphon also satisfies Condition 2.2(b), its stated sectional continuity outside null sets, the theorem adds

1n∑i=1nE∥Xin−Xi/n∥∗,T2⟶0,∥x∥∗,T=sup⁡0≤s≤T∣xs∣.\frac1n\sum_{i=1}^{n}\mathbb E\|X_i^n-X_{i/n}\|_{*,T}^2\longrightarrow0, \qquad \|x\|_{*,T}=\sup_{0\le s\le T}|x_s|.n1​i=1∑n​E∥Xin​−Xi/n​∥∗,T2​⟶0,∥x∥∗,T​=0≤s≤Tsup​∣xs​∣.

The milestones record the paper's cut-to-operator observation, uniform coupling moments, the one-edge and two-edge estimates, the bound on the second interaction-error term, the quantitative coupling lemma, and the empirical law of large numbers for the independent limiting particles. The second displayed limit is conditional on 2.2(b); the first is not. Bayraktar, Chakraborty and Wu, 2023, Theorem 4.1.

Significance

The first conclusion identifies the population-level limit even when the sampled network has a decreasing edge-retention probability. It gives a deterministic law on path space, rather than only a limit at one time. The second conclusion is stronger: it compares each finite particle with the limiting particle built from the same initial state and Brownian motion, averaged over all labels. This comparison is available when the extra graphon continuity assumption holds. Bayraktar, Chakraborty and Wu, 2023, Theorem 4.1 and Remark 4.1.

The mathematical results were proved in the published paper. The present mission asks for machine-checked Lean proofs of its specified statements; the draft theorems themselves currently contain sorry. Reusable outputs would include the percolated-edge model, a path-law formulation of graphon diffusions, and estimates for random interactions that depend on the particle states they influence. The paper's proofs appear in Sections 5–7, with the sparse-system argument in Section 7. Bayraktar, Chakraborty and Wu, 2023.

Difficulty

The edge variables are independent of the driving noise, but a particle's state depends on its incident edges. Therefore an expression such as E[ξijn∣Xjn−Xj/n∣2]\mathbb E[\xi^n_{ij}|X_j^n-X_{j/n}|^2]E[ξijn​∣Xjn​−Xj/n​∣2] cannot be factored into the edge mean times the error moment. The same dependence affects pairs of edges in the squared interaction sum. Equations (7.10), (7.14), and (7.15) quantify the discrepancy. A second difficulty is that cut convergence controls integrals against fixed bounded tests, while the dynamics and their path laws supply tests that vary with the system. The paper's regularity conditions and quantitative estimates are needed to connect those forms of control. Bayraktar, Chakraborty and Wu, 2023, Sections 7.1–7.2.

Formalization scope

Lean represents III by Mathlib's unit interval, Rd\mathbb R^dRd by Fin d → ℝ, and continuous paths by maps from [0,T][0,T][0,T] into that vector space. The chosen vector and matrix norms are fixed sup-type norms; finite-dimensional norm equivalence preserves the limit statements and the estimates with existential constants. The horizon and dimension are explicitly positive. Particle index i : Fin n denotes the paper's (i+1)(i+1)(i+1)st particle, and empirical laws use positive sizes n+1n+1n+1 to avoid division by zero. Wasserstein distances use extended nonnegative reals, while nonnegative expectations use lower integrals.

The noise lives on one probability space. The continuum equation uses the natural filtration of its own initial state and Brownian motion. The finite equation uses the natural filtration of its edge array, sampled initial states, and Brownian motions. A solution carries a measurable continuous-path law family in the paper's class M\mathcal MM; this prevents a nonmeasurable pushforward or a default-zero integral from silently satisfying an equation. Edge symmetry, Bernoulli marginals, independence, and diagonal self-loops are explicit. The zero-coefficient system has a separate sorry-free Lean witness for both solution predicates, ruling out a vacuous solution concept.

The formalization uses the published Peng1990.SMP.Stochastic Itô interface and WassersteinDRO.Duality.wassersteinDistance definition. Contributions can prove the listed estimates, strengthen their reusable stochastic and measure-theoretic infrastructure, or connect this mission's local setting to the other graphon missions after their definitions are published. The auxiliary changed-measure particles of Section 7.1 are outside the initial milestone list. Bayraktar, Chakraborty and Wu, 2023, Sections 4 and 7.

Selected references

  • E. Bayraktar, S. Chakraborty and R. Wu, Graphon mean field systems, Annals of Applied Probability 33(5):3587–3619, 2023. DOI.
  • L. Lovász, Large Networks and Graph Limits, American Mathematical Society Colloquium Publications 60, 2012. DOI.
13 thms1 active userReviewed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Twin-width I: Tractable FO Model Checking 3: Grid Minor Theorem for Twin-width — t-Twin-Ordered Matrices Are (2t+2)-Mixed Free, and t-Mixed Free Matrices Have Twin-width at Most 4c_t·α^(4c_t+2)Research Paper

Motivation

Twin-width is a graph and matrix invariant introduced by Bonnet, Kim, Thomassé and Watrigant (J. ACM 69(1), Article 3, 2021). Classes of bounded twin-width include bounded tree-width and clique-width graphs, proper minor-closed classes, posets of bounded width, and permutations avoiding a fixed pattern; on all of them, first-order model checking is fixed-parameter tractable once a contraction sequence is given. The difficulty in applying this framework is to bound twin-width for a concrete class, because a contraction sequence is an intricate object.

Section 5 of the paper supplies the tool that makes such bounds routine: the Grid Minor Theorem for twin-width (Theorem 5.4). It reduces bounding twin-width to finding an order of the rows and columns of a matrix (for a graph, an order of its vertices) in which no large mixed minor appears. The paper then uses it for pattern-avoiding permutations, posets of bounded width and KtK_tKt​-minor free graphs (Section 6, Table 1). The argument is modelled on the proof of the Stanley–Wilf conjecture by Marcus and Tardos (JCTA 107, 2004) and on its use by Guillemot and Marx for permutation pattern matching (SODA 2014).

Timeline. Füredi and Hajnal (1992) conjectured a linear bound on the number of 111s in a 0,10,10,1-matrix avoiding a permutation pattern. Marcus and Tardos (2004) proved it, with constant ct=2t4(t2t)c_t=2t^4\binom{t^2}{t}ct​=2t4(tt2​), via ttt-grid minors. Fox (arXiv:1310.8378) improved the constant to 3t 28t3t\,2^{8t}3t28t, and Cibulka and Kynčl (arXiv:1607.07491) to 83(t+1)224t\tfrac83(t+1)^22^{4t}38​(t+1)224t. Guillemot and Marx (2014) built from the grid structure a decomposition of permutations; Bonnet et al. (2020 preprint, 2021 journal) generalized it to arbitrary matrices over a finite alphabet, which is Theorem 5.4.

Setting

Let M=(mi,j)M=(m_{i,j})M=(mi,j​) be an n×mn\times mn×m matrix with entries in a finite alphabet AAA of size α\alphaα. For sets RRR of rows and CCC of columns, the zone R∩CR\cap CR∩C is the submatrix on those rows and columns; it is constant if all its entries are equal.

A partition of MMM is a pair (R,C)(\mathcal R,\mathcal C)(R,C) of a partition of the rows and a partition of the columns. Its error value is the largest number, over all parts RiR_iRi​ (respectively CjC_jCj​), of parts of the other side forming a non-constant zone with it. A contraction sequence is a sequence of partitions starting from the finest one (all singletons), ending at the coarsest one (one row part, one column part), each obtained from the previous by merging two row parts or two column parts. MMM has twin-width at most ttt, tww⁡(M)≤t\operatorname{tww}(M)\le ttww(M)≤t, if some contraction sequence has error value at most ttt throughout.

A division is a partition whose parts consist of consecutive indices. MMM is ttt-twin-ordered if it has a contraction sequence of error value at most ttt consisting of divisions.

A zone is vertical if its rows are all equal, horizontal if its columns are all equal, and mixed otherwise. A ttt-mixed minor is a division into ttt row intervals and ttt column intervals all of whose t2t^2t2 zones are mixed; MMM is ttt-mixed free if it has none. For 0,10,10,1-matrices, a ttt-grid minor is such a division in which every zone contains a 111. Finally ct=83(t+1)224tc_t=\tfrac83(t+1)^22^{4t}ct​=38​(t+1)224t.

Formalization targets

Goal: Theorem 5.4

For n,m≥1n,m\ge1n,m≥1:

$$M\ \text{ttt-twin-ordered}\ \Longrightarrow\ M\ \text{(2t+2)(2t+2)(2t+2)-mixed free},

t\ge1,\ M\ \text{ttt-mixed free}\ \Longrightarrow\ \operatorname{tww}(M)\le 4c_t,\alpha^{4c_t+2}.$$

The constant is the one the paper's proof produces; the paper's summary 22O(t)2^{2^{O(t)}}22O(t) is not formalized.

Milestones, in the order of the proof

  1. Lemma 5.5: a matrix is mixed iff it contains a corner, a mixed 2×22\times22×2 submatrix on consecutive rows and columns.
  2. Lemma 5.6: fusing two consecutive row parts does not increase the mixed value of a column interval (mixed zones plus mixed cuts).
  3. Lemma 5.2: a sequence of partitions, each rrr-refining the next, with error value at most ttt, certifies tww⁡(M)≤rt\operatorname{tww}(M)\le rttww(M)≤rt.
  4. Theorem 5.3 (Marcus–Tardos, constant of Cibulka–Kynčl): ctmax⁡(n,m)c_t\max(n,m)ct​max(n,m) entries 111 force a ttt-grid minor.
  5. Lemma 5.7: a ttt-mixed free matrix has a division sequence of mixed value at most 2ct2c_t2ct​.
  6. The first item of Theorem 5.4.
  7. The row-type count on p. 3:22: within a row part, rows restricted to the non-mixed zones take at most αt′+1\alpha^{t'+1}αt′+1 values.

Significance

Theorem 5.4 characterizes bounded twin-width of matrices, up to the order of rows and columns, by the exclusion of a mixed minor. Its second item is used in Section 6 of the same paper to show that KtK_tKt​-minor free graphs, posets of bounded width and pattern-avoiding permutations have bounded twin-width, and in later papers of the twin-width series as the standard route to twin-width upper bounds. Without it, each such bound needs its own explicit contraction sequence.

The result is proved in the paper; no machine-checked proof of it exists in Mathlib or on the platform. The mission produces a machine-checked statement and, once solved, a proof of the grid theorem with its explicit constant, a formal treatment of contraction sequences of matrices, and a statement of the Marcus–Tardos theorem with the Cibulka–Kynčl constant, which is itself absent from Mathlib.

Difficulty

The first item is comparatively short. The second item is where the work lies. The natural first idea — contract greedily as long as the error value stays below a bound depending on ttt — does not work: the error value is not monotone under fusions, and nothing in the absence of a mixed minor directly prevents a greedy sequence from getting stuck with large error. Any proof has to identify a quantity that is monotone under fusion and is bounded by a density argument, and then turn a sequence controlled by that quantity into a contraction sequence of bounded error value; the paper's choices (mixed value, Lemmas 5.6 and 5.7) are the milestones. The density argument rests on the Marcus–Tardos theorem with an explicit single-exponential constant, a substantial combinatorial result in its own right and not available in Mathlib.

Formalization scope

Matrices are Matrix (Fin n) (Fin m) A with [Fintype A], and α\alphaα is Fintype.card A. Partitions are Mathlib Finpartitions of univ; the finest is ⊥, and "coarsest" means at most one part. Twin-width is the partition form of p. 3:18, as a predicate MatTwinWidthLE M t with exactly one merge, of rows or of columns, per step; it is never an infimum over a possibly empty set. Divisions in minors are given by strictly increasing cut points 0=r0<⋯<rt=n0=r_0<\dots<r_t=n0=r0​<⋯<rt​=n, so all parts are non-empty. Vertical and horizontal mean "all rows (columns) of the zone agree". Grid minors use Bool entries. Since ctc_tct​ is not an integer, α4ct+2\alpha^{4c_t+2}α4ct​+2 is a real power, and the bound is ⌊4ctα4ct+2⌋\lfloor 4c_t\alpha^{4c_t+2}\rfloor⌊4ct​α4ct​+2⌋; the mixed-value bound of Lemma 5.7 is ⌊2ct⌋\lfloor 2c_t\rfloor⌊2ct​⌋.

Hypotheses added relative to the page: n,m≥1n,m\ge1n,m≥1 in Theorems 5.3, 5.4 and Lemma 5.7; t≥1t\ge1t≥1 in the second item of Theorem 5.4, in Theorem 5.3 and in Lemma 5.7, since for t=0t=0t=0 every non-empty matrix is vacuously 000-mixed free and the bound would be false; m≥1m\ge1m≥1 in the row-type count. Theorem 5.3 is stated with the explicit constant ctc_tct​ that the paper quotes and uses, not as a bare existence statement.

The statements are not trivializable by the usual routes: twin-width is not a junk infimum, a contraction sequence cannot jump to one part in a single step, the error value counts non-constant zones per part on both sides, minors range over divisions with non-empty interval parts rather than arbitrary partitions, and all bounds depend only on ttt and α\alphaα, never on MMM.

A complete development needs: finite partitions and their merges, divisions as interval partitions, the Marcus–Tardos theorem, and bookkeeping for sequences of partitions. The partition machinery and Theorem 5.3 are reusable beyond this mission (graph twin-width in the other missions of this series, pattern avoidance). Proofs of individual milestones, and of the Marcus–Tardos theorem from its original argument, are welcome.

Selected references

  • É. Bonnet, E. J. Kim, S. Thomassé, R. Watrigant, Twin-width I: Tractable FO Model Checking, J. ACM 69(1), Article 3, 2021. https://doi.org/10.1145/3486655
  • A. Marcus, G. Tardos, Excluded permutation matrices and the Stanley–Wilf conjecture, J. Combin. Theory Ser. A 107(1), 153–160, 2004. https://doi.org/10.1016/j.jcta.2004.04.002
  • J. Fox, Stanley–Wilf limits are typically exponential, arXiv:1310.8378, 2013. https://arxiv.org/abs/1310.8378
  • J. Cibulka, J. Kynčl, Füredi–Hajnal limits are typically subexponential, arXiv:1607.07491, 2016. https://arxiv.org/abs/1607.07491
  • S. Guillemot, D. Marx, Finding small patterns in permutations in linear time, SODA 2014, 82–101. https://doi.org/10.1137/1.9781611973402.7
10 thms1 active userReviewed
Markov ChainOperations ResearchProbability·Captain: mikedeng1

Routing Betweenness Centrality: Under Source-Oblivious Loop-Free Routing, the Betweenness of a Node Sequence Is a Sum over Targets of Chained Pairwise DependenciesResearch Paper

Motivation

Betweenness centrality measures how much of the communication in a network passes through a node. Freeman's original measure (Freeman 1977) counts the fraction of shortest paths between each pair of nodes that pass through it, and Brandes' algorithm (Brandes 2001) made it computable on large graphs. Variants replace shortest paths by maximum flows (Freeman, Borgatti and White 1991), random walks (Newman 2005) or load-splitting rules (Goh, Kahng and Kim 2001), and group variants measure the traffic that passes through a set of nodes (Everett and Borgatti 1999).

Dolev, Elovici and Puzis (technical report 2009, J. ACM 2010) define routing betweenness centrality (RBC): given the actual routing scheme of a communication network and its traffic matrix, the RBC of a node is the expected number of packets that pass through it. Each of the measures above is a special choice of routing scheme. The application is the placement of traffic monitors: the RBC of an ordered sequence of monitors is the expected number of packets sampled by all of them in order, which measures, for example, redundant inspection along a path.

Setting

Let VVV be a finite set of nodes. A routing scheme assigns to every source sss, current node uuu, next node vvv and target ttt a number R(s,u,v,t)R(s,u,v,t)R(s,u,v,t), the probability that uuu forwards to vvv a packet with source address sss and target address ttt. Routing decisions are independent. The scheme is assumed to be

  • nonnegative;
  • total: every node u≠tu\ne tu=t forwards a packet targeted at ttt with total probability one, ∑vR(s,u,v,t)=1\sum_v R(s,u,v,t)=1∑v​R(s,u,v,t)=1, and the target forwards nothing;
  • loop-free: for every pair (s,t)(s,t)(s,t) the directed graph with an arc a→ba\to ba→b whenever R(s,a,b,t)>0R(s,a,b,t)>0R(s,a,b,t)>0 is acyclic;
  • source-oblivious where stated: R(s,u,v,t)R(s,u,v,t)R(s,u,v,t) does not depend on sss.

A route from sss to ttt is a sequence of nodes from sss to ttt whose every hop has positive forwarding probability; its probability is the product of the hop probabilities. For nodes s,t,vs,t,vs,t,v and a sequence S=(s1,…,sk)S=(s_1,\dots,s_k)S=(s1​,…,sk​):

  • the pairwise dependency δs,t(v)\delta_{s,t}(v)δs,t​(v) is the total probability of the routes from sss to ttt that visit vvv;
  • the sequence dependency δ~s,t(S)\tilde\delta_{s,t}(S)δ~s,t​(S) is the total probability of the routes from sss to ttt that visit s1s_1s1​, then s2s_2s2​, …, then sks_ksk​ (consecutive repetitions in SSS collapsed).

A traffic matrix T(s,t)T(s,t)T(s,t) gives the number of packets sent from sss to ttt, and sampling rates ρv∈[0,1]\rho_v\in[0,1]ρv​∈[0,1] give the fraction of passing packets a monitor at vvv samples. The target dependency is δ∙,t(v)=∑sδs,t(v)T(s,t)\delta_{\bullet,t}(v)=\sum_{s}\delta_{s,t}(v)T(s,t)δ∙,t​(v)=∑s​δs,t​(v)T(s,t), and the RBC of the sequence SρS_\rhoSρ​ of distinct nodes is

δ~∙,∙(Sρ)=∏r∈Sρr⋅∑s,t∈Vδ~s,t(S) T(s,t).\tilde\delta_{\bullet,\bullet}(S_\rho)=\prod_{r\in S}\rho_r\cdot\sum_{s,t\in V}\tilde\delta_{s,t}(S)\,T(s,t).δ~∙,∙​(Sρ​)=r∈S∏​ρr​⋅s,t∈V∑​δ~s,t​(S)T(s,t).

Formalization targets

Goal: sequence RBC as a sum over targets of chained dependencies

For a source-oblivious scheme and a sequence S=(s1,…,sk)S=(s_1,\dots,s_k)S=(s1​,…,sk​) of distinct nodes, k≥1k\ge1k≥1,

∏r∈Sρr⋅∑s,t∈Vδ~s,t(S) T(s,t)=∑t∈Vδ∙,t(s1) ρs1∏i=1k−1δsi,t(si+1) ρsi+1.\prod_{r\in S}\rho_r\cdot\sum_{s,t\in V}\tilde\delta_{s,t}(S)\,T(s,t)=\sum_{t\in V}\delta_{\bullet,t}(s_1)\,\rho_{s_1}\prod_{i=1}^{k-1}\delta_{s_i,t}(s_{i+1})\,\rho_{s_{i+1}}.r∈S∏​ρr​⋅s,t∈V∑​δ~s,t​(S)T(s,t)=t∈V∑​δ∙,t​(s1​)ρs1​​i=1∏k−1​δsi​,t​(si+1​)ρsi+1​​.

The right side is the value returned by the paper's Algorithm 6. The identity contains no running time and no constant: it states the exact equality the algorithm relies on.

Milestones

  1. Proposition 1: for a loop-free scheme (not necessarily source-oblivious), two different orderings of the same set of nodes cannot both have positive sequence dependency.
  2. Lemma 1 (Dependency chaining): for s≠ts\ne ts=t, δ~s,t((s1,…,sk))=δs,t(s1)⋅δ~s1,t((s2,…,sk))\tilde\delta_{s,t}((s_1,\dots,s_k))=\delta_{s,t}(s_1)\cdot\tilde\delta_{s_1,t}((s_2,\dots,s_k))δ~s,t​((s1​,…,sk​))=δs,t​(s1​)⋅δ~s1​,t​((s2​,…,sk​)).
  3. Eq. (16): for s≠ts\ne ts=t, δ~s,t((s1,…,sk))=δs,t(s1)∏i=2kδsi−1,t(si)\tilde\delta_{s,t}((s_1,\dots,s_k))=\delta_{s,t}(s_1)\prod_{i=2}^k\delta_{s_{i-1},t}(s_i)δ~s,t​((s1​,…,sk​))=δs,t​(s1​)∏i=2k​δsi−1​,t​(si​).
  4. Eq. (17): δ~∙,t((s1,…,sk))=δ∙,t(s1)∏i=2kδsi−1,t(si)\tilde\delta_{\bullet,t}((s_1,\dots,s_k))=\delta_{\bullet,t}(s_1)\prod_{i=2}^k\delta_{s_{i-1},t}(s_i)δ~∙,t​((s1​,…,sk​))=δ∙,t​(s1​)∏i=2k​δsi−1​,t​(si​).
  5. Eq. (1) (off the goal's path): δs,t(s)=1\delta_{s,t}(s)=1δs,t​(s)=1 and δs,t(v)=∑uδs,t(u)R(s,u,v,t)\delta_{s,t}(v)=\sum_{u}\delta_{s,t}(u)R(s,u,v,t)δs,t​(v)=∑u​δs,t​(u)R(s,u,v,t) for v≠sv\ne sv=s, the sum over predecessors uuu of vvv.
  6. Eq. (9) (off the goal's path): δ∙,t(v)=T(v,t)+∑uδ∙,t(u)R(⊘,u,v,t)\delta_{\bullet,t}(v)=T(v,t)+\sum_u\delta_{\bullet,t}(u)R(\oslash,u,v,t)δ∙,t​(v)=T(v,t)+∑u​δ∙,t​(u)R(⊘,u,v,t) for source-oblivious schemes.

Significance

The goal reduces the RBC of a sequence of kkk monitors to O(nk)O(nk)O(nk) arithmetic on precomputed single-node quantities, instead of a pass over the routing graph of every source–target pair. This is what makes repeated evaluation of candidate monitor sequences, as in greedy or search-based placement, affordable on large networks. Eqs. (1) and (9) are the recursions behind the paper's algorithms for the RBC of single nodes; Eq. (9) replaces the loop over all source–target pairs by a loop over targets.

The paper's proof of Lemma 1 is an argument about conditional probabilities of informally described events. A formalization fixes what "the probability that a packet passes through SSS" means, as a sum over a finite set of routes, and checks the factorization against that definition, including the corner cases the prose passes over: sequences with repeated nodes, a first node equal to the source or to the target, and the source s=ts=ts=t in the sum over all sources. No machine-checked version of these identities is known.

Difficulty

Dependency chaining is a Markov property: after the packet reaches s1s_1s1​, its future does not depend on its past. The obvious argument splits every route through s1s_1s1​ into a prefix ending at s1s_1s1​ and a suffix starting there. Two points need care. First, the suffix must be a route of a packet with source s1s_1s1​, which is where source-obliviousness enters; for source-dependent routing the lemma is false. Second, the prefix probabilities must sum to δs,t(s1)\delta_{s,t}(s_1)δs,t​(s1​) while the suffix probabilities sum to one, which requires that no packet is lost: this is where totality of the forwarding probabilities and loop-freeness (every walk reaches the target in at most ∣V∣|V|∣V∣ steps) are needed. The passage from Lemma 1 to the goal then requires the sum over all sources, including s=ts=ts=t, where Lemma 1 does not apply.

Formalization scope

Nodes form a type V with [Fintype V] [DecidableEq V]; probabilities, traffic and sampling rates are real numbers; sequences are Lean lists, indexed from 000. Routes are repetition-free lists, and every dependency is a finite sum over them; no infinite sums, suprema or measure theory are used. The definitions file bundles the standing assumptions in IsRoutingScheme and source-obliviousness in IsSourceOblivious.

The paper's probabilities are informal. The following readings are explicit:

  • every non-target node forwards with total probability one (implicit in the paper; without it Lemma 1 fails);
  • the target forwards nothing;
  • the paper's convention R(⊘,v,v,⊘)=1R(\oslash,v,v,\oslash)=1R(⊘,v,v,⊘)=1 is not adopted (it contradicts loop-freeness); collapsing consecutive repetitions in SSS reproduces its only use;
  • the "don't care" source ⊘\oslash⊘ of Eq. (9) is an arbitrary fixed node, harmless under source-obliviousness;
  • the printed index slips in Eqs. (16) and (17) are corrected;
  • no sign condition is imposed on TTT; 0≤ρv≤10\le\rho_v\le10≤ρv​≤1 is kept in the goal as on the page.

Defining δ~\tilde\deltaδ~ by the recursion (3) or the product (16), the target sequence dependency by (17), or the sequence RBC as the sum over targets of the chained product would make the goal true by definition; every dependency here is a sum of route probabilities, and the RBC of a sequence is the double sum of Eq. (4).

The running-time claims of the paper are out of scope. Lemma 2 of the paper (p. 17) is false as printed for a nonempty monitor set and is not posed. Contributions welcome: proofs of the milestones, a general lemma that the total route probability from any node is one under a total loop-free scheme (reusable for absorbing Markov chains on finite DAGs), and the prefix–suffix decomposition of routes.

Selected references

  • S. Dolev, Y. Elovici, R. Puzis, Routing Betweenness Centrality, Technical Report #2009-09, Ben-Gurion University of the Negev, 2009; J. ACM 57(4), 2010. https://doi.org/10.1145/1734213.1734219
  • L. C. Freeman, A set of measures of centrality based on betweenness, Sociometry 40(1), 1977. https://doi.org/10.2307/3033543
  • U. Brandes, A faster algorithm for betweenness centrality, J. Math. Sociology 25(2), 2001. https://doi.org/10.1080/0022250X.2001.9990249
  • L. C. Freeman, S. P. Borgatti, D. R. White, Centrality in valued graphs: a measure of betweenness based on network flow, Social Networks 13, 1991. https://doi.org/10.1016/0378-8733(91)90017-N
  • M. E. J. Newman, A measure of betweenness centrality based on random walks, Social Networks 27, 2005. https://doi.org/10.1016/j.socnet.2004.11.009
  • K.-I. Goh, B. Kahng, D. Kim, Universal behavior of load distribution in scale-free networks, Phys. Rev. Lett. 87, 2001. https://doi.org/10.1103/PhysRevLett.87.278701
  • M. G. Everett, S. P. Borgatti, The centrality of groups and classes, J. Math. Sociology 23(3), 1999. https://doi.org/10.1080/0022250X.1999.9990219
8 thms1 active userReviewed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

A c^k n 5-Approximation Algorithm for Treewidth 1: If tw(G) ≤ k, Pushed Terminal Separations Around a 3/4-Balanced Separation Glue to a Separation of Order k+1 with Sides ≤ 7/8|V|+(|B|+k+1)/2Research Paper

Motivation

Treewidth measures how closely a graph resembles a tree. Many problems that are NP-hard on general graphs can be solved in time linear in the number of vertices once a tree decomposition of bounded width is available, so computing such decompositions quickly is a basic task in algorithmic graph theory. Bodlaender, Drange, Dregi, Fomin, Lokshtanov and Pilipczuk (SIAM J. Comput. 45(2), 2016) gave an algorithm that, for a graph GGG on nnn vertices and an integer kkk, either reports that the treewidth of GGG exceeds kkk or returns a tree decomposition of width at most 5k+45k+45k+4, in time 2O(k)n2^{O(k)} n2O(k)n.

The paper's headline is that algorithm, a statement about pseudo-code and running times. This mission does not formalize it. It formalizes one of the paper's combinatorial lemmas, Lemma 6.5 (p. 369), which the algorithm's data structure uses to answer its balanced-separator query: it turns the search for a small balanced separator into two maximization problems that a dynamic program can solve.

Timeline. Robertson and Seymour (Graph Minors II, J. Algorithms 7, 1986) proved that graphs of treewidth kkk have balanced separators of size k+1k+1k+1 (Lemma 2.1 of the paper). Robertson and Seymour (Graph Minors XIII, 1995) used it in an O(33kn2)O(3^{3k} n^2)O(33kn2) approximation algorithm, Reed (1992) improved the dependence on nnn to nlog⁡nn \log nnlogn, and Bodlaender et al. (2016) reached linear time with a single-exponential dependence on kkk.

Setting

All graphs are finite and simple. For X⊆V(G)X \subseteq V(G)X⊆V(G), G∖XG \setminus XG∖X is the graph obtained by deleting XXX; its components are the vertex sets of its connected components.

A separation of GGG is a partition (L,X,R)(L, X, R)(L,X,R) of V(G)V(G)V(G) (parts may be empty) such that no edge of GGG joins LLL to RRR. Its order is ∣X∣|X|∣X∣. It is α\alphaα-balanced if ∣L∣,∣R∣≤α∣V(G)∣|L|, |R| \le \alpha|V(G)|∣L∣,∣R∣≤α∣V(G)∣.

For W⊆V(G)W \subseteq V(G)W⊆V(G) and sets of terminals TL,TRT_L, T_RTL​,TR​, a terminal separation of the induced subgraph G[W]G[W]G[W] of order ℓ\ellℓ (Definition 6.4) is a partition (L,X,R)(L, X, R)(L,X,R) of WWW with

TL⊆L,TR⊆R,no edge between L and R,∣X∣≤ℓ.T_L \subseteq L,\qquad T_R \subseteq R,\qquad \text{no edge between } L \text{ and } R,\qquad |X| \le \ell.TL​⊆L,TR​⊆R,no edge between L and R,∣X∣≤ℓ.

It is left-pushed if ∣L∣|L|∣L∣ is maximum, and right-pushed if ∣R∣|R|∣R∣ is maximum, among all terminal separations of G[W]G[W]G[W] of order ℓ\ellℓ with the same terminals.

The treewidth tw(G)\mathrm{tw}(G)tw(G) is the least width of a tree decomposition of GGG (Definition 1.1), where the width is the largest bag size minus one.

Formalization targets

Goal: Lemma 6.5

Let tw(G)≤k\mathrm{tw}(G) \le ktw(G)≤k and let (A1,B,A2)(A_1, B, A_2)(A1​,B,A2​) be a separation of GGG with ∣A1∣,∣A2∣≤34∣V(G)∣|A_1|, |A_2| \le \tfrac34|V(G)|∣A1​∣,∣A2​∣≤43​∣V(G)∣. Then there are a partition (TL,XB,TR)(T_L, X_B, T_R)(TL​,XB​,TR​) of BBB and integers k1,k2k_1, k_2k1​,k2​ with k1+k2+∣XB∣≤k+1k_1 + k_2 + |X_B| \le k+1k1​+k2​+∣XB​∣≤k+1 such that, for Gi=G[Ai∪(B∖XB)]G_i = G[A_i \cup (B \setminus X_B)]Gi​=G[Ai​∪(B∖XB​)] with terminals TL,TRT_L, T_RTL​,TR​,

  1. G1G_1G1​ and G2G_2G2​ have terminal separations of orders k1k_1k1​ and k2k_2k2​;
  2. for every left-pushed (L1,X1,R1)(L_1, X_1, R_1)(L1​,X1​,R1​) of order k1k_1k1​ in G1G_1G1​ and every right-pushed (L2,X2,R2)(L_2, X_2, R_2)(L2​,X2​,R2​) of order k2k_2k2​ in G2G_2G2​, the triple (L1∪TL∪L2, X1∪XB∪X2, R1∪TR∪R2)(L_1 \cup T_L \cup L_2,\ X_1 \cup X_B \cup X_2,\ R_1 \cup T_R \cup R_2)(L1​∪TL​∪L2​, X1​∪XB​∪X2​, R1​∪TR​∪R2​) is a terminal separation of GGG of order at most k+1k+1k+1 with
∣L1∪TL∪L2∣, ∣R1∪TR∪R2∣≤78∣V(G)∣+∣B∣+(k+1)2.|L_1 \cup T_L \cup L_2|,\ |R_1 \cup T_R \cup R_2| \le \tfrac78|V(G)| + \frac{|B| + (k+1)}{2}.∣L1​∪TL​∪L2​∣, ∣R1​∪TR​∪R2​∣≤87​∣V(G)∣+2∣B∣+(k+1)​.

Milestones

In the order of the proof, the milestones are: Lemma 2.1 (balanced SSS-separators of size k+1k+1k+1, cited from Graph Minors II); the packing step of the proof sketch of Lemma 2.2 (components holding at most 12∣S∣\tfrac12|S|21​∣S∣ of SSS each split into two sides holding at most 23∣S∣\tfrac23|S|32​∣S∣ each); the folklore 23\tfrac2332​-balanced separation of order k+1k+1k+1; the restriction of that separation to G1,G2G_1, G_2G1​,G2​, which gives part (i); four lower bounds on ∣L∩Ai∣|L \cap A_i|∣L∩Ai​∣, ∣R∩Ai∣|R \cap A_i|∣R∩Ai​∣; the dichotomy that two opposite "quadrants" have at least 18∣V(G)∣−∣B∣+(k+1)2\tfrac18|V(G)| - \tfrac{|B|+(k+1)}{2}81​∣V(G)∣−2∣B∣+(k+1)​ vertices each; the comparison of pushed separations with the restricted one; and the gluing step of part (ii).

Significance

Lemma 6.5 is what makes the paper's balanced-separator query computable: given a bag BBB that splits a component in a balanced way, a left-pushed and a right-pushed terminal separation are each the solution of a maximization problem, and the lemma guarantees that gluing them gives a separator of order k+1k+1k+1 whose sides are bounded away from ∣V(G)∣|V(G)|∣V(G)∣. The 78\tfrac7887​ constant feeds into the approximation ratio of the whole algorithm.

The proof is complete in the paper and elementary given Lemma 2.1. No part of it is formalized in Lean's Mathlib or, as far as a search of the Prove2Me catalog shows, on the platform. Lemma 2.1 itself, Robertson and Seymour's balanced-separator theorem for bounded treewidth, is not in Mathlib; a machine-checked proof of it would be reusable wherever treewidth separators are used (divide-and-conquer on tree decompositions, Lipton–Tarjan style arguments, grid-minor theorems).

Difficulty

Everything after Lemma 2.1 is counting with partitions, but the bookkeeping is the obstacle: the two separations (A1,B,A2)(A_1, B, A_2)(A1​,B,A2​) and (L,X,R)(L, X, R)(L,X,R) cut V(G)V(G)V(G) into nine pieces, the induced subgraphs G1,G2G_1, G_2G1​,G2​ overlap on B∖XBB \setminus X_BB∖XB​, and pushedness is a maximum over a family of partitions of a subset, not of V(G)V(G)V(G).

Lemma 2.1 is the hard milestone. It is a statement about tree decompositions, and any proof has to relate subtrees of the decomposition tree to components of G∖XG \setminus XG∖X for a bag XXX, infrastructure that does not yet exist in Mathlib.

Formalization scope

  • Vertices. A finite type V : Type with decidable equality; vertex sets are Finset V, V(G)V(G)V(G) is Finset.univ, ∣V(G)∣|V(G)|∣V(G)∣ is Fintype.card V. Universe 0 matches the published definition of treewidth.
  • Treewidth. "tw(G)≤k\mathrm{tw}(G) \le ktw(G)≤k" is the published RobertsonSeymour1986.GM5.TreewidthLE G k: a decomposition tree on a finite type whose bags have at most k+1k+1k+1 vertices. Definition 1.1 roots the tree; rooting does not change the width.
  • Components of G∖XG \setminus XG∖X. The component of u∉Xu \notin Xu∈/X is avoidComp G X u: the vertices reachable from uuu by a walk of GGG none of whose vertices, endpoints included, lies in XXX.
  • Partitions are pairwise disjoint triples with the given union; empty parts are allowed.
  • Induced subgraphs G[W]G[W]G[W] are encoded by their vertex set WWW: a terminal separation of G[W]G[W]G[W] is a partition of WWW, with GGG's adjacency restricted to WWW. Pushedness maximizes over partitions of WWW only.
  • "Of order ℓ\ellℓ" means ∣X∣≤ℓ|X| \le \ell∣X∣≤ℓ, as in Definition 6.4 (iii).
  • Arithmetic. Bounds with fractions or subtraction (34\tfrac3443​, 23\tfrac2332​, 78\tfrac7887​, 18∣V(G)∣−∣B∣+(k+1)2\tfrac18|V(G)| - \tfrac{|B|+(k+1)}{2}81​∣V(G)∣−2∣B∣+(k+1)​) are stated over R\mathbb RR; bounds without subtraction may be cleared of denominators in N\mathbb NN (2∣C∩S∣≤∣S∣2|C \cap S| \le |S|2∣C∩S∣≤∣S∣, 3∣L′∩S∣≤2∣S∣3|L' \cap S| \le 2|S|3∣L′∩S∣≤2∣S∣).
  • Corrections of the page. The printed bound of Lemma 6.5 has ∣X∣|X|∣X∣ where the proof derives ∣B∣|B|∣B∣; the goal states ∣B∣|B|∣B∣. The proof's "14\tfrac1441​- and 13\tfrac1331​-balanced" means 34\tfrac3443​ and 23\tfrac2332​; the milestone uses the latter.

Part (i) is part of the goal: without it, left- and right-pushed separations of the chosen orders might not exist and part (ii) would hold vacuously. Choosing XB=BX_B = BXB​=B does not trivialize the goal, because the size bound and the budget k1+k2+∣XB∣≤k+1k_1 + k_2 + |X_B| \le k+1k1​+k2​+∣XB​∣≤k+1 must hold simultaneously.

Welcome contributions: a proof of Lemma 2.1 from tree decompositions (reusable on its own), the packing lemma, and the counting milestones of the proof of Lemma 6.5.

Selected references

  • H. L. Bodlaender, P. G. Drange, M. S. Dregi, F. V. Fomin, D. Lokshtanov, M. Pilipczuk, A cknc^k nckn 5-approximation algorithm for treewidth, SIAM J. Comput. 45(2):317–378, 2016. https://doi.org/10.1137/130947374
  • N. Robertson, P. D. Seymour, Graph minors. II. Algorithmic aspects of tree-width, J. Algorithms 7(3):309–322, 1986. https://doi.org/10.1016/0196-6774(86)90023-4
  • N. Robertson, P. D. Seymour, Graph minors. XIII. The disjoint paths problem, J. Combin. Theory Ser. B 63(1):65–110, 1995. https://doi.org/10.1006/jctb.1995.1006
  • B. A. Reed, Finding approximate separators and computing tree width quickly, STOC 1992, 221–228. https://doi.org/10.1145/129712.129734
12 thms1 active userReviewed
Combinatorics·Captain: mikedeng1

Three-Coloring and List Three-Coloring of Graphs Without Induced Paths on Seven Vertices 1: Every Connected P₇-Free Graph Has a Connected 2-Dominating Set of at Most Three Vertices or a K₄Research Paper

Motivation

Graph coloring asks whether vertices can be assigned colors so that adjacent vertices receive different colors. Restricting the graphs by forbidden induced paths can expose structure that is absent from general graphs. Bonomo, Chudnovsky, Maceli, Schaudt, Stein and Zhong study list three-coloring when the graph has no induced path on seven vertices, building on earlier work on coloring graphs with shorter forbidden induced paths. Their paper proves an algorithmic result, but its treatment also uses a purely graph-theoretic statement: a connected P7P_7P7​-free graph has a small connected set that reaches every vertex within two steps, unless the graph contains a four-vertex clique. This mission isolates that statement as Corollary 5 of the paper.

The structural input to the corollary is a theorem of Camby and Schaudt about connected dominating sets in PtP_tPt​-free graphs, quoted as Theorem 4 in the coloring paper. That result applies to every allowed path length ttt, rather than only to the seven-vertex case. It is a milestone here in the precise form printed in the source article. The Camby–Schaudt paper develops connected domination as a characterization of graphs without long induced paths.

Setting

A simple graph GGG has a vertex set V(G)V(G)V(G) and an undirected edge relation with no loops. The graphs in this mission are finite. The path PtP_tPt​ has ttt vertices and consecutive vertices adjacent. A graph is PtP_tPt​-free when it has no induced subgraph isomorphic to PtP_tPt​: the graph may contain a path with ttt vertices as a non-induced subgraph if extra edges join vertices of that path. This distinction matters throughout the source paper.

For a vertex set SSS, its open neighborhood N(S)N(S)N(S) consists of vertices outside SSS with a neighbor in SSS. Write Sˉ=S∪N(S)\bar S=S\cup N(S)Sˉ=S∪N(S). The set SSS is dominating when Sˉ=V(G)\bar S=V(G)Sˉ=V(G), so each vertex lies in SSS or has an edge to SSS. It is 2-dominating when Sˉ‾=V(G)\overline{\bar S}=V(G)Sˉ=V(G), meaning each vertex is at graph distance at most two from some vertex of SSS. A dominating or 2-dominating set is connected when the subgraph G[S]G[S]G[S] induced by SSS is connected. In this mission, connectedness includes nonemptiness, as in the paper's use of such sets.

A clique of size four is a set of four pairwise adjacent vertices. The notation K4⊆GK_4\subseteq GK4​⊆G below asserts the existence of such a set. Since all edges between those four vertices are present, it is also an induced complete subgraph. The notation G[S]≅Pt−2G[S]\cong P_{t-2}G[S]≅Pt−2​ means there is a graph isomorphism from the induced graph on SSS to that path.

Formalization targets

Corollary 5: a small connected 2-dominating set or a clique

For a finite, connected, P7P_7P7​-free graph GGG, the goal is

(∃S⊆V(G):∣S∣≤3,  G[S] connected,  Sˉ‾=V(G))∨K4⊆G.\left(\exists S\subseteq V(G): |S|\le 3,\;G[S]\text{ connected},\; \overline{\bar S}=V(G)\right)\quad\lor\quad K_4\subseteq G.(∃S⊆V(G):∣S∣≤3,G[S] connected,Sˉ=V(G))∨K4​⊆G.

The four-vertex clique is a genuine alternative in the claim; the mission does not assume its absence. The result is the first, structural sentence of Corollary 5. Its following sentence gives a running time for finding the set or clique and is outside this mission.

Theorem 4: connected domination across forbidden path lengths

For every t≥3t\ge 3t≥3, the cited theorem asserts that each finite connected PtP_tPt​-free graph has a connected dominating set SSS for which

G[S] is Pt−2-free∨G[S]≅Pt−2.G[S]\text{ is }P_{t-2}\text{-free}\quad\lor\quad G[S]\cong P_{t-2}.G[S] is Pt−2​-free∨G[S]≅Pt−2​.

The other milestones are statements identified in the proof of Corollary 5: a connected P3P_3P3​-free graph is complete, domination inside a dominating set yields 2-domination in the original graph, and the three non-leaf vertices work when the induced dominating graph is a P5P_5P5​. These are listed in the paper's order of use, with the quoted theorem also available as a separate target.

Significance

Corollary 5 gives a bounded set of vertices from which the rest of a connected P7P_7P7​-free graph can be reached in two steps, except when a K4K_4K4​ is present. Its bound is the concrete number three, independent of graph size. In the coloring paper this supplies a small structural seed for later restrictions on vertex palettes. A formal statement of the corollary makes that finite structural claim available independently of any coloring algorithm.

Theorem 4 has broader reuse: it links forbidden induced paths to connected dominating sets for general ttt. The remaining milestones isolate the graph facts needed to connect this general result to the seven-vertex corollary. The article and the cited theorem provide mathematical proofs; these Lean declarations are open theorem statements, with machine-checked proofs still to be supplied. This mission therefore asks for formal proofs of known graph-theoretic results, not for a new resolution of an open mathematical problem.

Difficulty

The difficult step is finding a connected dominating set whose induced graph has the much stronger path restriction in Theorem 4. Ordinary domination by itself does not impose that restriction, and replacing induced-path exclusion with exclusion of all paths would change the class of graphs. The later cases require careful handling of connectedness and the passage from domination within G[S]G[S]G[S] to distance-two domination in GGG. These issues are mathematical as well as formal: an empty set or a merely preconnected induced graph would erase the intended content.

Formalization scope

Lean represents GGG as SimpleGraph V with a finite vertex type. Vertex sets are Finset V; PtP_tPt​ is Mathlib's pathGraph t; PtP_tPt​-freeness is the absence of an induced copy using IsIndContained. The definitions of N(S)N(S)N(S), Sˉ\bar SSˉ, domination and 2-domination are shared within P7ThreeColor.Seed. A connected induced graph uses Connected, which requires at least one vertex. Four-clique existence is ¬ G.CliqueFree 4. The only added standing convention is finiteness, implicit in the source's graph algorithm and finite cardinalities. No extra colorability or clique-free hypothesis is attached to the goal.

The running-time claim in the second sentence of Corollary 5 and the paper's headline Theorem 1 are not formalized because this mission has no machine model. A statement that merely postulates a small set as a hypothesis, or substitutes ordinary domination for 2-domination, would miss the result. Contributions are welcome for proofs of the goal and milestones, reusable facts about induced paths and connected domination, and auxiliary Mathlib lemmas needed to cross between graph isomorphisms, induced subgraphs and finite sets.

Selected references

  • F. Bonomo, M. Chudnovsky, P. Maceli, O. Schaudt, M. Stein and M. Zhong, Three-coloring and list three-coloring of graphs without induced paths on seven vertices, Combinatorica, OnlineFirst 2017, DOI: 10.1007/s00493-017-3553-8.
  • E. Camby and O. Schaudt, A new characterization of PkP_kPk​-free graphs, Algorithmica 75 (2016), 205–217, DOI: 10.1007/s00453-015-9989-6; arXiv:1402.7213.
6 thms1 active userReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems II: The Empirical Measure of the Dense Graphon-Weighted Particle System Converges in Probability to the Averaged Graphon LawResearch Paper

Interacting diffusions on large weighted graphs

Classical mean field theory describes nnn exchangeable diffusions, each interacting with the empirical average of all the others. As n→∞n\to\inftyn→∞ the particles become independent copies of a McKean–Vlasov process (propagation of chaos). Many systems of interest are not exchangeable. Particles sit at the vertices of a network and interact only with their neighbours, as in models of opinion dynamics, epidemics on contact networks, synchronisation of oscillators and systemic risk in financial networks. When the network is dense and its adjacency structure converges, the natural limit object is a graphon, a symmetric measurable function G:[0,1]2→[0,1]G:[0,1]^2\to[0,1]G:[0,1]2→[0,1] that encodes the limiting edge density between "positions" u,v∈[0,1]u,v\in[0,1]u,v∈[0,1] (Lovász, Large Networks and Graph Limits, 2012).

Bayraktar, Chakraborty and Wu (Ann. Appl. Probab. 33(5), 2023) develop a law of large numbers for such systems. This mission formalizes their second main result, Theorem 3.1. Under cut-metric convergence of the interaction kernels, the empirical measure of the nnn-particle system converges in probability to the averaged law of a continuum graphon particle system. Related work includes Bhamidi, Budhiraja and Wu (Stochastic Process. Appl., 2019) on weakly interacting particles on inhomogeneous random graphs, Coppini, Dietert and Giacomin (Stoch. Dyn., 2020) on Erdős–Rényi graphs, and Bet, Coppini and Nardi (arXiv:2006.07670) on oscillators on dense random graphs.

Setting

Write I=[0,1]I=[0,1]I=[0,1] with Lebesgue measure and fix a horizon T>0T>0T>0. Let Cd=C([0,T]:Rd)\mathcal C_d=C([0,T]:\mathbb R^d)Cd​=C([0,T]:Rd) with the uniform norm ∥x∥∗,T=sup⁡s≤T∣xs∣\|x\|_{*,T}=\sup_{s\le T}|x_s|∥x∥∗,T​=sups≤T​∣xs​∣. On one probability space take independent initial states Xu(0)X_u(0)Xu​(0) with laws μu(0)\mu_u(0)μu​(0) and independent ddd-dimensional Brownian motions BuB_uBu​, for all u∈Iu\in Iu∈I. The graphon particle system (2.1) for a graphon GGG and Lipschitz coefficients b,σb,\sigmab,σ is the continuum of nonlinear diffusions

Xu(t)=Xu(0)+∫0t ⁣ ⁣∫I ⁣∫Rdb(Xu(s),x)G(u,v) μv,s(dx) dv ds+∫0t ⁣ ⁣∫I ⁣∫Rdσ(Xu(s),x)G(u,v) μv,s(dx) dv dBu(s),X_u(t)=X_u(0)+\int_0^t\!\!\int_I\!\int_{\mathbb R^d}b(X_u(s),x)G(u,v)\,\mu_{v,s}(dx)\,dv\,ds+\int_0^t\!\!\int_I\!\int_{\mathbb R^d}\sigma(X_u(s),x)G(u,v)\,\mu_{v,s}(dx)\,dv\,dB_u(s),Xu​(t)=Xu​(0)+∫0t​∫I​∫Rd​b(Xu​(s),x)G(u,v)μv,s​(dx)dvds+∫0t​∫I​∫Rd​σ(Xu​(s),x)G(u,v)μv,s​(dx)dvdBu​(s),

where μv,s\mu_{v,s}μv,s​ is the law of Xv(s)X_v(s)Xv​(s). Write μu∈P(Cd)\mu_u\in\mathcal P(\mathcal C_d)μu​∈P(Cd​) for the path law of XuX_uXu​ and μˉ=∫Iμu du\bar\mu=\int_I\mu_u\,duμˉ​=∫I​μu​du for the averaged law.

The nnn-particle system (3.1) places particle iii at label i/ni/ni/n, starts it at Xi/n(0)X_{i/n}(0)Xi/n​(0), drives it by Bi/nB_{i/n}Bi/n​, and weights its interaction with particle jjj by ξijn\xi^n_{ij}ξijn​:

Xin(t)=Xi/n(0)+∫0t1n∑j=1nξijn b(Xin(s),Xjn(s)) ds+∫0t1n∑j=1nξijn σ(Xin(s),Xjn(s)) dBi/n(s).X^n_i(t)=X_{i/n}(0)+\int_0^t\frac1n\sum_{j=1}^n\xi^n_{ij}\,b(X^n_i(s),X^n_j(s))\,ds+\int_0^t\frac1n\sum_{j=1}^n\xi^n_{ij}\,\sigma(X^n_i(s),X^n_j(s))\,dB_{i/n}(s).Xin​(t)=Xi/n​(0)+∫0t​n1​j=1∑n​ξijn​b(Xin​(s),Xjn​(s))ds+∫0t​n1​j=1∑n​ξijn​σ(Xin​(s),Xjn​(s))dBi/n​(s).

The weights come from a step graphon GnG_nGn​, constant on the squares of the nnn-grid. Either ξijn=Gn(i/n,j/n)\xi^n_{ij}=G_n(i/n,j/n)ξijn​=Gn​(i/n,j/n), or ξijn=ξjin\xi^n_{ij}=\xi^n_{ji}ξijn​=ξjin​ are independent Bernoulli(Gn(i/n,j/n))(G_n(i/n,j/n))(Gn​(i/n,j/n)) variables for i≤ji\le ji≤j, independent of the noise (Condition 3.1). The kernels converge in the cut norm ∥W∥□=sup⁡S,T∣∫S×TW∣\|W\|_\square=\sup_{S,T}\big|\int_{S\times T}W\big|∥W∥□​=supS,T​​∫S×T​W​. Two further norms appear: the operator norm ∥W∥=sup⁡∥g∥∞≤1∫I∣∫IW(u,v)g(v) dv∣ du\|W\|=\sup_{\|g\|_\infty\le1}\int_I|\int_I W(u,v)g(v)\,dv|\,du∥W∥=sup∥g∥∞​≤1​∫I​∣∫I​W(u,v)g(v)dv∣du, and the Wasserstein-2 distance W2,TW_{2,T}W2,T​ on P(Cd)\mathcal P(\mathcal C_d)P(Cd​).

Formalization targets

Goal: Theorem 3.1

Assume the initial laws are measurable in uuu with a uniform (2+ε)(2+\varepsilon)(2+ε)-moment and b,σb,\sigmab,σ are Lipschitz (Condition 2.1). Assume u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) is W2W_2W2​-continuous on each interval of a finite cover of III (Condition 2.2(a)), and that Condition 3.1 holds with Gn→GG_n\to GGn​→G in the cut metric. Then

μn:=1n∑i=1nδXin⟶μˉin P(Cd) in probability.\mu^n:=\frac1n\sum_{i=1}^n\delta_{X^n_i}\longrightarrow\bar\mu\quad\text{in }\mathcal P(\mathcal C_d)\text{ in probability}.μn:=n1​i=1∑n​δXin​​⟶μˉ​in P(Cd​) in probability.

If, in addition, GGG is continuous at (u,v)(u,v)(u,v) for a.e. vvv at each interior point uuu of the cover (Condition 2.2(b)), then

1n∑i=1nE∥Xin−Xi/n∥∗,T2⟶0.\frac1n\sum_{i=1}^n\mathbb E\|X^n_i-X_{i/n}\|^2_{*,T}\longrightarrow0 .n1​i=1∑n​E∥Xin​−Xi/n​∥∗,T2​⟶0.

Milestones

In order: Remark 2.1, where cut-norm convergence implies operator-norm convergence; Theorem 2.1(a), continuity of u↦μuu\mapsto\mu_uu↦μu​ in W2,TW_{2,T}W2,T​; the weak-LLN bound (6.4); the discretization bound (6.11); Lemma 6.1,

lim sup⁡n1n∑iE∥Xin−Xi/n∥∗,T2≤κ(M)lim sup⁡n∥Gn−G∥+κM−ε;\limsup_{n}\frac1n\sum_i\mathbb E\|X^n_i-X_{i/n}\|^2_{*,T}\le\kappa(M)\limsup_n\|G_n-G\|+\kappa M^{-\varepsilon};nlimsup​n1​i∑​E∥Xin​−Xi/n​∥∗,T2​≤κ(M)nlimsup​∥Gn​−G∥+κM−ε;

Lemma 6.2, a law of large numbers for the continuum particles X1/n,…,Xn/nX_{1/n},\dots,X_{n/n}X1/n​,…,Xn/n​; the approximation of GGG by continuous graphons; the continuity (6.12) of μˉ\bar\muμˉ​ in the graphon, in a corrected qualitative form; the empirical coupling bound for W2,TW_{2,T}W2,T​; and (6.14).

Significance

Theorem 3.1 shows that the averaged graphon law μˉ\bar\muμˉ​ is the correct macroscopic description of a dense heterogeneous particle system. This holds for deterministic weights and for Erdős–Rényi-type random weights, and the graphon may be discontinuous. It is the foundation for graphon mean field games and for control of large networked systems, where μˉ\bar\muμˉ​ replaces the intractable nnn-particle law. Under the regularity Condition 2.2(b) the convergence holds in L2L^2L2, particle by particle, along the natural coupling.

The result is proved in the paper (§6.1–§6.2). Nothing in it is formalized. A formal proof would produce reusable infrastructure: graphons, the cut and operator norms, Lusin approximation of graphons, Wasserstein distances between empirical measures, and laws of large numbers for random probability measures on path space. It also needs pathwise stability estimates for SDE systems built on Itô integrals, which goes well beyond what Mathlib has today.

Difficulty

The particles are neither exchangeable nor identically distributed. The classical propagation-of-chaos argument compares each particle with an i.i.d. copy of one McKean–Vlasov process, and it does not apply here. Each XinX^n_iXin​ must be compared with its own limit Xi/nX_{i/n}Xi/n​, and the error has two sources: the random or discretized weights ξijn\xi^n_{ij}ξijn​, and the replacement of GGG by GnG_nGn​. Cut-metric convergence controls Gn−GG_n-GGn​−G only through integrals against bounded test functions. The interaction term, however, integrates the unbounded function b(Xi/n(s),⋅)b(X_{i/n}(s),\cdot)b(Xi/n​(s),⋅) against varying laws μj/n,s\mu_{j/n,s}μj/n,s​. When GGG is merely measurable, the map u↦μuu\mapsto\mu_uu↦μu​ need not be continuous, and the error of replacing ∫I⋅ dv\int_I\cdot\,dv∫I​⋅dv by the average over the labels j/nj/nj/n need not vanish. The direct estimate therefore fails for (3.3) without Condition 2.2(b), and this is the case the goal must cover.

Formalization scope

The setting file GraphonMF.DenseLLN.Setting fixes these conventions:

  • Rd\mathbb R^dRd is Fin d → ℝ with the sup norm. All statements are invariant under this choice: constants are existential and the other conclusions are limits.
  • Processes are path-valued maps Ω→Cd\Omega\to\mathcal C_dΩ→Cd​, and Cd\mathcal C_dCd​ carries its Borel σ-algebra.
  • Each SDE is encoded through the published Itô-process predicate Peng1990.SMP.IsItoProcess. The filtrations are the uncompleted natural filtrations σ(Xu(0),Bu(r):r≤t)\sigma(X_u(0),B_u(r):r\le t)σ(Xu​(0),Bu​(r):r≤t) and σ(ξn,Xj/n(0),Bj/n(r):j≤n,r≤t)\sigma(\xi^n,X_{j/n}(0),B_{j/n}(r):j\le n,r\le t)σ(ξn,Xj/n​(0),Bj/n​(r):j≤n,r≤t).
  • A solution of (2.1) carries the membership of its law family in the class M\mathcal MM (measurable in uuu, uniformly bounded second moments), which makes every Bochner integral of the drift meaningful.
  • W2W_2W2​ and W2,TW_{2,T}W2,T​ are the published WassersteinDRO.Duality.wassersteinDistance with values in [0,∞][0,\infty][0,∞].
  • Expectations of nonnegative quantities are lower Lebesgue integrals.
  • Convergence in probability is topological: P(μn∉U)→0\mathbb P(\mu^n\notin U)\to0P(μn∈/U)→0 for every weak-topology neighbourhood UUU of μˉ\bar\muμˉ​.
  • The nnn-particle system shares the noise of the graphon system: particle iii (Lean index iii, label (i+1)/n(i+1)/n(i+1)/n) uses X(i+1)/n(0)X_{(i+1)/n}(0)X(i+1)/n​(0) and B(i+1)/nB_{(i+1)/n}B(i+1)/n​.
  • Each ξn\xi^nξn is assumed measurable.

Two encodings would trivialize the targets, and both are ruled out. A pointwise or almost-sure reading of (3.3), or convergence of ∫f dμn\int f\,d\mu^n∫fdμn for a single test function, is not the claim. Choosing the constants of Lemma 6.1 and (6.11) after the sequence (Gn)(G_n)(Gn​) would make them nearly empty, so they are chosen before it. A sorry-free check confirms that the solution predicates are satisfiable: with b=σ=0b=\sigma=0b=σ=0 the constant paths solve both systems.

Solvers will need several pieces of library work: the inequality ∥W∥∞→1≤4∥W∥□\|W\|_{\infty\to1}\le4\|W\|_\square∥W∥∞→1​≤4∥W∥□​, BDG-type moment bounds for the Itô integrals of Peng's library, Gronwall estimates, Lusin's theorem on [0,1]2[0,1]^2[0,1]2, and the metrization of weak convergence on P(Cd)\mathcal P(\mathcal C_d)P(Cd​) (Lévy–Prokhorov). Contributions to any of these are welcome independently of the goal. The paper's proofs are in its §5 (Theorem 2.1) and §6.1–§6.2 (Lemmas 6.1, 6.2 and Theorem 3.1).

Selected references

  • E. Bayraktar, S. Chakraborty, R. Wu, Graphon mean field systems, Ann. Appl. Probab. 33(5):3587–3619, 2023. https://doi.org/10.1214/22-AAP1901
  • L. Lovász, Large Networks and Graph Limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
  • S. Bhamidi, A. Budhiraja, R. Wu, Weakly interacting particle systems on inhomogeneous random graphs, Stochastic Process. Appl. 129:2174–2206, 2019. https://doi.org/10.1016/j.spa.2018.06.014
  • F. Coppini, H. Dietert, G. Giacomin, A law of large numbers and large deviations for interacting diffusions on Erdős–Rényi graphs, Stoch. Dyn. 20:2050010, 2020. https://doi.org/10.1142/S0219493720500100
  • G. Bet, F. Coppini, F. R. Nardi, Weakly interacting oscillators on dense random graphs, arXiv preprint, 2020. https://arxiv.org/abs/2006.07670
  • A.-S. Sznitman, Topics in propagation of chaos, École d'Été de Probabilités de Saint-Flour XIX, LNM 1464, 1991. https://doi.org/10.1007/BFb0085169
15 thms1 active userReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems III: With Edge Weights Sampled From a Lipschitz Graphon, Each Particle Is Within κ/n of Its Graphon Limit in Mean SquareResearch Paper

Interacting diffusions on large heterogeneous networks

Classical mean-field theory describes nnn diffusions that interact symmetrically, each feeling the empirical average of all the others, and shows that as n→∞n\to\inftyn→∞ every particle behaves like an independent copy of a single McKean–Vlasov process. Many systems of interest are not symmetric: neurons, agents in a financial network, or oscillators interact through a graph, and particles in different parts of the graph feel different neighbourhoods. A graphon G:[0,1]2→[0,1]G:[0,1]^2\to[0,1]G:[0,1]2→[0,1] is the limit object of dense graph sequences (Lovász, Large Networks and Graph Limits, AMS 2012), and it is the natural way to describe the interaction structure of such systems in the limit.

Bayraktar, Chakraborty and Wu (Ann. Appl. Probab. 33(5), 2023) study a continuum of diffusions indexed by labels u∈[0,1]u\in[0,1]u∈[0,1] whose interaction is weighted by a graphon, and prove laws of large numbers and rates of convergence for nnn-particle systems on graphs sampled from the graphon. This mission targets their uniform rate of convergence for weights sampled from a Lipschitz graphon (Theorem 3.2).

Setting

Write I=[0,1]I=[0,1]I=[0,1] with Lebesgue measure. Fix a horizon T∈(0,∞)T\in(0,\infty)T∈(0,∞) and let Cd=C([0,T]:Rd)\mathcal C_d=C([0,T]:\mathbb R^d)Cd​=C([0,T]:Rd) with ∥x∥∗,t=sup⁡0≤s≤t∣xs∣\|x\|_{*,t}=\sup_{0\le s\le t}|x_s|∥x∥∗,t​=sup0≤s≤t​∣xs​∣.

On one probability space, {Bu:u∈I}\{B_u:u\in I\}{Bu​:u∈I} are i.i.d. ddd-dimensional Brownian motions and {Xu(0):u∈I}\{X_u(0):u\in I\}{Xu​(0):u∈I} are independent initial states with laws μu(0)\mu_u(0)μu​(0), independent of the Brownian motions. Given a graphon GGG and coefficients b,σb,\sigmab,σ, the graphon particle system (2.1) is

Xu(t)=Xu(0)+∫0t ⁣ ⁣∫I ⁣∫Rdb(Xu(s),x)G(u,v) μv,s(dx) dv ds+∫0t ⁣ ⁣∫I ⁣∫Rdσ(Xu(s),x)G(u,v) μv,s(dx) dv dBu(s),X_u(t)=X_u(0)+\int_0^t\!\!\int_I\!\int_{\mathbb R^d}b(X_u(s),x)G(u,v)\,\mu_{v,s}(dx)\,dv\,ds+\int_0^t\!\!\int_I\!\int_{\mathbb R^d}\sigma(X_u(s),x)G(u,v)\,\mu_{v,s}(dx)\,dv\,dB_u(s),Xu​(t)=Xu​(0)+∫0t​∫I​∫Rd​b(Xu​(s),x)G(u,v)μv,s​(dx)dvds+∫0t​∫I​∫Rd​σ(Xu​(s),x)G(u,v)μv,s​(dx)dvdBu​(s),

with μu,t=L(Xu(t))\mu_{u,t}=\mathcal L(X_u(t))μu,t​=L(Xu​(t)). The particles are independent but not identically distributed, and their laws are coupled through the dvdvdv-integral.

The nnn-particle system (3.1) uses the same randomness at the labels i/ni/ni/n:

Xin(t)=Xi/n(0)+∫0t1n∑j=1nξijn b(Xin(s),Xjn(s)) ds+∫0t1n∑j=1nξijn σ(Xin(s),Xjn(s)) dBi/n(s),X^n_i(t)=X_{i/n}(0)+\int_0^t\frac1n\sum_{j=1}^n\xi^n_{ij}\,b(X^n_i(s),X^n_j(s))\,ds+\int_0^t\frac1n\sum_{j=1}^n\xi^n_{ij}\,\sigma(X^n_i(s),X^n_j(s))\,dB_{i/n}(s),Xin​(t)=Xi/n​(0)+∫0t​n1​j=1∑n​ξijn​b(Xin​(s),Xjn​(s))ds+∫0t​n1​j=1∑n​ξijn​σ(Xin​(s),Xjn​(s))dBi/n​(s),

so particle iii and the continuum particle at u=i/nu=i/nu=i/n share their initial state and Brownian motion.

The hypotheses are:

  • Condition 2.1: u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) is measurable, sup⁡uE∣Xu(0)∣2+ε<∞\sup_u\mathbb E|X_u(0)|^{2+\varepsilon}<\inftysupu​E∣Xu​(0)∣2+ε<∞ for some ε>0\varepsilon>0ε>0, and bbb, σ\sigmaσ are Lipschitz in both arguments.
  • Condition 2.3: there are finitely many intervals I1,…,INI_1,\dots,I_NI1​,…,IN​ covering III and a κ>0\kappa>0κ>0 with W2(μu1(0),μu2(0))≤κ∣u1−u2∣W_2(\mu_{u_1}(0),\mu_{u_2}(0))\le\kappa|u_1-u_2|W2​(μu1​​(0),μu2​​(0))≤κ∣u1​−u2​∣ on each IiI_iIi​ and ∣G(u1,v1)−G(u2,v2)∣≤κ(∣u1−u2∣+∣v1−v2∣)|G(u_1,v_1)-G(u_2,v_2)|\le\kappa(|u_1-u_2|+|v_1-v_2|)∣G(u1​,v1​)−G(u2​,v2​)∣≤κ(∣u1​−u2​∣+∣v1​−v2​∣) on each block Ii×IjI_i\times I_jIi​×Ij​. The graphon may jump across block boundaries.
  • Condition 3.2: the weights are sampled from GGG itself, either deterministically, ξijn=G(i/n,j/n)\xi^n_{ij}=G(i/n,j/n)ξijn​=G(i/n,j/n), or as a random graph, ξijn=ξjin∼Bernoulli(G(i/n,j/n))\xi^n_{ij}=\xi^n_{ji}\sim\mathrm{Bernoulli}(G(i/n,j/n))ξijn​=ξjin​∼Bernoulli(G(i/n,j/n)) independently for i≤ji\le ji≤j and independently of the noise.

Formalization targets

Goal: Theorem 3.2 (p. 3596)

Under Conditions 2.1, 2.3 and 3.2 there is κ∈(0,∞)\kappa\in(0,\infty)κ∈(0,∞) such that

max⁡i=1,…,nE∥Xin−Xi/n∥∗,T2≤κn∀n∈N.\max_{i=1,\dots,n}\mathbb E\|X^n_i-X_{i/n}\|^2_{*,T}\le\frac{\kappa}{n}\qquad\forall n\in\mathbb N.i=1,…,nmax​E∥Xin​−Xi/n​∥∗,T2​≤nκ​∀n∈N.

Milestones (§6.3, in attack order)

  • (6.15): E∥Xin−Xi/n∥∗,t2\mathbb E\|X^n_i-X_{i/n}\|^2_{*,t}E∥Xin​−Xi/n​∥∗,t2​ is bounded by κ\kappaκ times the time-integrated mean-square mismatch between the finite-nnn drift and diffusion and their graphon counterparts.
  • The drift mismatch is split by (6.16) into three terms T~sn,1,T~sn,2,T~sn,3\tilde T^{n,1}_s,\tilde T^{n,2}_s,\tilde T^{n,3}_sT~sn,1​,T~sn,2​,T~sn,3​, bounded by
T~sn,1≤2κmax⁡iE∣Xin(s)−Xi/n(s)∣2  (6.17),T~sn,2≤κn  (6.18),T~sn,3≤κn2  (6.19).\tilde T^{n,1}_s\le2\kappa\max_i\mathbb E|X^n_i(s)-X_{i/n}(s)|^2\ \ (6.17),\qquad\tilde T^{n,2}_s\le\frac{\kappa}{n}\ \ (6.18),\qquad\tilde T^{n,3}_s\le\frac{\kappa}{n^2}\ \ (6.19).T~sn,1​≤2κimax​E∣Xin​(s)−Xi/n​(s)∣2  (6.17),T~sn,2​≤nκ​  (6.18),T~sn,3​≤n2κ​  (6.19).
  • Theorem 2.1(b) (p. 3593): under Condition 2.3, W2,T(μu,μv)≤κ∣u−v∣W_{2,T}(\mu_u,\mu_v)\le\kappa|u-v|W2,T​(μu​,μv​)≤κ∣u−v∣ whenever u,vu,vu,v lie in the same interval IiI_iIi​. It is used in (6.19).

Significance

The result. Theorem 3.2 gives a quantitative propagation of chaos for heterogeneous interactions: with weights sampled from a blockwise Lipschitz graphon, every particle of the finite system, not just an average particle, is within O(n−1/2)O(n^{-1/2})O(n−1/2) in root-mean-square path distance of the continuum particle with the same label. The rate matches the classical mean-field rate (e.g. Sznitman 1991), so heterogeneity of the interaction costs nothing in the exponent. The same estimate covers the deterministic weighted graph and the random graph sampled from GGG (in the annealed sense), and it is the dense counterpart of Theorem 4.2 for percolated graphs.

Formalizing it. The result is proved in the paper (§6.3); it has no machine-checked proof. A formal development would produce a reusable account of graphon particle systems and their finite-nnn approximations: graphon-weighted McKean–Vlasov drifts, couplings through shared noise, and the decomposition of the interaction error into a Lipschitz feedback term, a law-of-large-numbers term and a Riemann-sum discretization term.

Difficulty

The obvious approach is to compare XinX^n_iXin​ with Xi/nX_{i/n}Xi/n​ directly and close a Gronwall inequality. The difficulty is that the comparison involves the empirical, graph-weighted average of the other particles, which is neither the graphon average nor an average of independent terms. Three different errors must be controlled simultaneously and uniformly in iii: the feedback of all particles' errors into particle iii, the fluctuation of a sum of independent but non-identically distributed terms whose weights may themselves be random, and the error of replacing an integral over labels by a sum over the grid {j/n}\{j/n\}{j/n}. The last needs regularity of the law μu\mu_uμu​ in the label uuu, which is itself a theorem about the graphon system (Theorem 2.1(b)) and fails at block boundaries of GGG, where only the small number of boundary cells saves the 1/n21/n^21/n2 rate. Averaging over iii, as in the law of large numbers of Theorem 3.1, is not enough: the goal bounds the maximum.

Formalization scope

  • Rd\mathbb R^dRd is Fin d → ℝ with the sup norm; matrices are measured entrywise. All statements are norm-invariant because every constant is existential. Particle i∈{1,…,n}i\in\{1,\dots,n\}i∈{1,…,n} is Lean's i : Fin n with label (i+1)/n(i+1)/n(i+1)/n; time is ℝ≥0; III is Mathlib's unitInterval.
  • Brownian motions, Itô integrals and Itô processes are the published Peng1990.SMP.Stochastic definitions; an SDE with random initial value X(0)X(0)X(0) is stated as an Itô process for X−X(0)X-X(0)X−X(0). Each continuum particle uses the natural filtration of (Xu(0),Bu)(X_u(0),B_u)(Xu​(0),Bu​); the nnn-particle system uses the natural filtration of the weights, initial states and Brownian motions at the labels.
  • W2W_2W2​ and W2,TW_{2,T}W2,T​ are the published WassersteinDRO.Duality.wassersteinDistance 2, valued in [0,∞][0,\infty][0,∞]; expectations of nonnegative quantities are lower Lebesgue integrals.
  • A solution of (2.1) is required to have continuous paths, measurable path maps, and a law family in the class M\mathcal MM (measurable in uuu, uniformly bounded second moments), as Proposition 2.1 asserts; this makes every inner integral meaningful.
  • The constant κ\kappaκ is chosen after the data and before nnn and iii. A κ\kappaκ chosen after nnn, a bound on the average over iii, or a solution predicate that no process satisfies would trivialize the goal; the first two are excluded by the statement, and a sanity check exhibits the zero-coefficient system as a solution of both systems.
  • Condition 2.2 is not assumed, and the weights come from GGG, not from a step graphon GnG_nGn​; no cut-metric hypothesis appears.
  • Needed infrastructure: Doob/Burkholder–Davis–Gundy for Itô integrals, Gronwall's inequality in integrated form, moment bounds for the graphon system, Kantorovich–Rubinstein-type duality for Lipschitz test functions, and independence of functionals of independent noises. Proofs of the milestones, of Theorem 2.1(b), and of the goal are all welcome; the paper's proofs are in its §5–§6.

Selected references

  • E. Bayraktar, S. Chakraborty, R. Wu, Graphon mean field systems, Ann. Appl. Probab. 33(5):3587–3619, 2023. https://doi.org/10.1214/22-AAP1901
  • L. Lovász, Large Networks and Graph Limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
  • A.-S. Sznitman, Topics in propagation of chaos, École d'Été de Probabilités de Saint-Flour XIX, LNM 1464, Springer, 1991. https://doi.org/10.1007/BFb0085169
  • S. Peng, A general stochastic maximum principle for optimal control problems, SIAM J. Control Optim. 28(4):966–979, 1990. https://doi.org/10.1137/0328054
10 thms1 active userReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems I: The Law of the Graphon Particle System Depends Continuously on the Graphon in the Cut MetricResearch Paper

Motivation

Large interacting diffusion systems often have different interaction strengths between different pairs of agents. A graphon gives a continuum description of such weighted networks: each agent carries a label in [0,1][0,1][0,1], and a kernel assigns the strength of interaction between two labels. The question in this mission is whether the law of the resulting continuum stochastic system changes continuously when its graphon changes. This matters when a large finite network is represented by an approximate graphon, since an approximation of the network should lead to an approximation of the dynamics. Bayraktar, Chakraborty and Wu establish this stability as Theorem 2.1(c) of Graphon mean field systems.

The paper also studies dense and less dense finite-particle systems. Their convergence results use a limiting graphon particle system even when the limiting kernel is only measurable, with no continuity in the labels. Theorem 2.1(c) is the continuity result that makes approximation by simpler graphons useful in those later sections Bayraktar–Chakraborty–Wu, §§3–4.

Setting

Let I=[0,1]I=[0,1]I=[0,1], fix a time horizon T>0T>0T>0, and let Cd=C([0,T];Rd)\mathcal C_d=C([0,T];\mathbb R^d)Cd​=C([0,T];Rd) be the space of continuous paths with the uniform norm ∥x∥∗,t=sup⁡0≤s≤t∣xs∣\|x\|_{*,t}=\sup_{0\le s\le t}|x_s|∥x∥∗,t​=sup0≤s≤t​∣xs​∣. A graphon G:I2→[0,1]G:I^2\to[0,1]G:I2→[0,1] is measurable and symmetric. Its cut norm is the largest absolute integral of GGG over a measurable rectangle S×U⊆I2S\times U\subseteq I^2S×U⊆I2. For two graphons, their cut distance is the cut norm of their difference.

For each label u∈Iu\in Iu∈I, an initial state Xu(0)X_u(0)Xu​(0) and a ddd-dimensional Brownian motion BuB_uBu​ live on one probability space. Initial states are mutually independent, Brownian motions are mutually independent, and the two families are independent of each other. The initial state has law μu(0)\mu_u(0)μu​(0). The process XuX_uXu​ evolves with a drift obtained by averaging b(Xu(s),x)G(u,v)b(X_u(s),x)G(u,v)b(Xu​(s),x)G(u,v) first against the time-sss law of XvX_vXv​ and then over v∈Iv\in Iv∈I. Its diffusion coefficient averages the matrix-valued σ(Xu(s),x)G(u,v)\sigma(X_u(s),x)G(u,v)σ(Xu​(s),x)G(u,v) in the same way and drives BuB_uBu​. This is the graphon particle system (2.1) Bayraktar–Chakraborty–Wu, p. 3591.

Condition 2.1 requires measurable initial laws, a uniform (2+ε)(2+\varepsilon)(2+ε) moment for some ε>0\varepsilon>0ε>0, and a common global Lipschitz bound for bbb and σ\sigmaσ. The path-law family μ=(L(Xu))u∈I\mu=(\mathcal L(X_u))_{u\in I}μ=(L(Xu​))u∈I​ belongs to M\mathcal MM when it is measurable as a family of probability measures on Cd\mathcal C_dCd​ and has uniformly bounded second moments. The distance W2,tW_{2,t}W2,t​ between path laws uses the supremum norm through time ttt; W2,tMW^\mathcal M_{2,t}W2,tM​ is the supremum of these distances over labels Bayraktar–Chakraborty–Wu, pp. 3591–3592.

Formalization targets

Stability of path laws

For graphons GnG_nGn​ converging to GGG in cut metric, let μuGn\mu^{G_n}_uμuGn​​ and μuG\mu^G_uμuG​ be the laws of the associated solutions, with the same initial laws and coefficients. The target is

∥Gn−G∥□⟶0⟹∫I[W2,T(μuGn,μuG)]2 du⟶0.\|G_n-G\|_\square\longrightarrow0 \quad\Longrightarrow\quad \int_I[W_{2,T}(\mu^{G_n}_u,\mu^G_u)]^2\,du\longrightarrow0.∥Gn​−G∥□​⟶0⟹∫I​[W2,T​(μuGn​​,μuG​)]2du⟶0.

The conclusion is an integrated law convergence statement. It does not assert that the laws converge uniformly in uuu or that paths converge almost surely. This is precisely Theorem 2.1(c) Bayraktar–Chakraborty–Wu, p. 3593.

Supporting results

The mission also states the paper's cut-to-operator convergence (Remark 2.1), existence of frozen-flow solutions (5.2), the squared estimate established in the proof of (5.3), well-posedness and moments for the nonlinear system (Proposition 2.1), and the quantitative stability display on p. 3604. These are the source's claims used to reach the target, with their provenance attached to each milestone Bayraktar–Chakraborty–Wu, §§2 and 5.

Significance

Stability shows that the continuum model is robust to graphon approximation in a metric natural for dense networks. Together with Proposition 2.1, it gives a well-defined solution-law family whose dependence on the network kernel can be controlled. The paper uses this framework in its later laws of large numbers for finite interacting systems Bayraktar–Chakraborty–Wu, §§3–4.

The theorem is proved in the paper; this mission asks for a machine-checked proof of the faithful statements. A complete development would also supply reusable Lean interfaces for graphons, cut and operator norms, measurable families of path laws, and stochastic systems whose coefficients depend on those families. The Itô-process and general Wasserstein definitions already available as published modules are reused. The remaining graphon-specific definitions and the estimates in §5 are open proof targets.

Difficulty

A direct stability argument for stochastic differential equations would compare the two coefficients point by point. Cut convergence does not provide pointwise convergence of Gn(u,v)G_n(u,v)Gn​(u,v) to G(u,v)G(u,v)G(u,v), nor a uniform bound on their difference that tends to zero. The drift and diffusion also contain laws generated by the processes being compared, so a change in the graphon feeds back into the law family. The central challenge is to control these coupled differences using only a weak kernel metric while retaining the measurable structure needed to integrate over labels Bayraktar–Chakraborty–Wu, §5.2.

The result also depends on a well-posedness theory for a continuum of stochastic equations. The laws need to vary measurably with uuu even though the paper does not assume the sample paths themselves vary measurably with uuu Bayraktar–Chakraborty–Wu, Remark 2.3(b).

Formalization scope

Lean uses unitInterval for III, nonnegative real time, and continuous maps from [0,T][0,T][0,T] to Fin d → ℝ for Cd\mathcal C_dCd​. The norm on Rd\mathbb R^dRd is the sup norm; the paper's Euclidean norm is equivalent in finite dimension, so existential Lipschitz and stability constants may change while the qualitative claim remains the same. The diffusion coefficient is a d×dd\times dd×d matrix, measured through its entries. Wasserstein quantities and nonnegative expectations are valued in [0,∞][0,\infty][0,∞], avoiding defaults for infinite moments. The path-law family uses the measurable space of measures, with the Borel space on paths. A zero-time or zero-dimensional state space is admitted where the formulas still have their stated meaning; T>0T>0T>0 and ε>0\varepsilon>0ε>0 are explicit.

For each graphon, a solution is a continuous-path process driven by the same given initial variables and Brownian motions, adapted to the natural filtration of those variables for each label. Its path map is almost everywhere measurable, its law family belongs to M\mathcal MM, and the nested interaction integrals are required to be integrable. These clauses prevent a zero measure from arising through a nonmeasurable pushforward and prevent a nonintegrable Bochner integral from silently evaluating to zero. The graphons are explicitly measurable, symmetric, and [0,1][0,1][0,1]-valued, so the real cut and operator norms are used only on bounded kernels. The theorem also asserts measurability of each Wasserstein integrand, reflecting the paper's treatment of the displayed integral.

The paper's proofs are in §§5–7; this mission concerns §5. Contributions to any stated milestone and to the reusable measure and Itô lemmas needed for them are welcome. A formalization that assumes the stability limit, that treats arbitrary unbounded kernels as graphons, or that uses an empty solution class would not establish this target.

Selected references

  • E. Bayraktar, S. Chakraborty and R. Wu, Graphon mean field systems, Annals of Applied Probability 33(5), 3587–3619, 2023. DOI 10.1214/22-AAP1901.
9 thms1 active userReviewed
Combinatorics·Captain: mikedeng1

Turán Graphs with Bounded Matching Number 1: An n-Vertex Graph with Clique Number at Most k and Matching Number at Most s Has at Most max{t(2s+1,k), g(n,k,s)} Edges, and This Is AttainedResearch Paper

Two classical extremal problems and their common generalization

Two of the oldest results in extremal graph theory bound the number of edges of a graph under a single structural restriction. Turán's theorem (Turán 1941) determines the maximum number of edges t(n,k)t(n,k)t(n,k) of an nnn-vertex graph with no complete subgraph on k+1k+1k+1 vertices. The Erdős–Gallai theorem (Erdős and Gallai 1959) determines the maximum number of edges of an nnn-vertex graph with no s+1s+1s+1 pairwise disjoint edges; for n≥2s+1n\ge 2s+1n≥2s+1 it is max⁡{(2s+12),(s2)+s(n−s)}\max\{\binom{2s+1}{2},\binom{s}{2}+s(n-s)\}max{(22s+1​),(2s​)+s(n−s)}, attained by a clique on 2s+12s+12s+1 vertices plus isolated vertices, or by sss vertices joined to everything.

N. Alon and P. Frankl (arXiv:2210.15076, published in J. Combin. Theory Ser. B, 2024, DOI 10.1016/j.jctb.2023.12.002) impose both restrictions at once and determine the maximum for all values of the parameters. Their answer is again the larger of two explicit constructions, and each of the two classical theorems is recovered as a degenerate range of theirs. The question belongs to a line of work on generalized Turán problems with an additional bound on the matching number, which has continued since this paper with other forbidden subgraphs in place of the clique.

Setting

All graphs are finite and simple, on the vertex set {0,1,…,n−1}\{0,1,\dots,n-1\}{0,1,…,n−1}.

  • The clique number of GGG is the largest number of vertices of a complete subgraph. "Clique number at most kkk" means that GGG contains no complete subgraph on k+1k+1k+1 vertices.
  • A matching is a set of pairwise disjoint edges. The matching number ν(G)\nu(G)ν(G) is the largest size of a matching of GGG.
  • A graph is complete kkk-partite if its vertex set is split into kkk classes (some possibly empty) and two vertices are adjacent exactly when they lie in different classes.
  • T(n,k)T(n,k)T(n,k) is the complete kkk-partite graph on nnn vertices whose classes have sizes as equal as possible, and t(n,k)t(n,k)t(n,k) is its number of edges, the Turán number.
  • G(n,k,s)G(n,k,s)G(n,k,s) is the complete kkk-partite graph on nnn vertices with k−1k-1k−1 classes of sizes as equal as possible and total size sss, and one further class of size n−sn-sn−s. Its number of edges is g(n,k,s)g(n,k,s)g(n,k,s).

A graph on nnn vertices is admissible for (k,s)(k,s)(k,s) when its clique number is at most kkk and its matching number is at most sss. Both constructions are admissible when n≥2s+1n\ge 2s+1n≥2s+1: G(n,k,s)G(n,k,s)G(n,k,s) is kkk-colourable, and every edge of it meets the sss vertices of the small classes; the Turán graph T(2s+1,k)T(2s+1,k)T(2s+1,k) on 2s+12s+12s+1 of the vertices, with the other n−2s−1n-2s-1n−2s−1 vertices isolated, has no s+1s+1s+1 disjoint edges.

Formalization targets

Goal: Theorem 1.1

For k≥2k\ge2k≥2 and n≥2s+1n\ge 2s+1n≥2s+1,

max⁡{∣E(G)∣:G on n vertices, ω(G)≤k, ν(G)≤s}  =  max⁡{ t(2s+1,k), g(n,k,s) },\max\bigl\{|E(G)| : G \text{ on } n \text{ vertices},\ \omega(G)\le k,\ \nu(G)\le s\bigr\} \;=\; \max\{\,t(2s+1,k),\ g(n,k,s)\,\},max{∣E(G)∣:G on n vertices, ω(G)≤k, ν(G)≤s}=max{t(2s+1,k), g(n,k,s)},

stated as an upper bound for every admissible graph together with an admissible graph attaining it.

Milestones

The milestones follow the proof of §2 of the paper, in attack order:

  1. the two constructions are admissible (the lower bound);
  2. the barrier: every graph has a vertex set BBB such that every component of G−BG-BG−B is odd and ∣B∣+∑i(∣Ai∣−1)/2=ν(G)|B|+\sum_i(|A_i|-1)/2=\nu(G)∣B∣+∑i​(∣Ai​∣−1)/2=ν(G);
  3. Lemma 2.1: in the §2 choice of an extremal graph and barrier maximizing ∑i∣Ai∣2\sum_i|A_i|^2∑i​∣Ai​∣2, non-adjacent vertices of BBB have the same neighbourhood;
  4. the graph induced on BBB is then complete kkk-partite;
  5. Claim 2.3: for that choice of graph and barrier, without loss of generality every component of G−BG-BG−B has a vertex with no neighbour in the smallest class of BBB;
  6. Lemma 2.2: for an extremal graph maximizing ∑i∣Ai∣2\sum_i |A_i|^2∑i​∣Ai​∣2, at most one component of G−BG-BG−B has more than one vertex;
  7. the case analysis on b=∣B∣b=|B|b=∣B∣ that bounds the remaining structure by max⁡{t(2s+1,k),g(n,k,s)}\max\{t(2s+1,k),g(n,k,s)\}max{t(2s+1,k),g(n,k,s)};
  8. the identity t(s+⌊s/(k−1)⌋,k)+(n−s−⌊s/(k−1)⌋) s=g(n,k,s)t(s+\lfloor s/(k-1)\rfloor,k)+(n-s-\lfloor s/(k-1)\rfloor)\,s=g(n,k,s)t(s+⌊s/(k−1)⌋,k)+(n−s−⌊s/(k−1)⌋)s=g(n,k,s);
  9. the Turán increment t(N+1,k)−t(N,k)=N−⌊N/k⌋t(N+1,k)-t(N,k)=N-\lfloor N/k\rfloort(N+1,k)−t(N,k)=N−⌊N/k⌋ is non-decreasing in NNN, with steps at most 111;
  10. the function f(b)=t(2s−b+1,k)+b(n−2s+b−1)f(b)=t(2s-b+1,k)+b(n-2s+b-1)f(b)=t(2s−b+1,k)+b(n−2s+b−1) has strictly increasing differences on the range of Case 4.

A companion item records the paper's parenthetical remark: for n≤2s+1n\le 2s+1n≤2s+1 the maximum is t(n,k)t(n,k)t(n,k).

Significance

The theorem settles the edge-extremal problem for the two most basic graph parameters, clique number and matching number, jointly and for every admissible parameter value. It contains both classical theorems: when n≤2s+1n\le 2s+1n≤2s+1 the matching condition is void and the answer is Turán's t(n,k)t(n,k)t(n,k); when k≥2s+1k\ge 2s+1k≥2s+1 the clique condition is void and the answer is the Erdős–Gallai bound. The two extremal constructions show that both regimes genuinely occur, and the threshold between them depends on nnn, kkk and sss together. The same paper's Proposition 3.1 extends the answer g(n,k,s)g(n,k,s)g(n,k,s) to every colour-critical forbidden graph for large sss and nnn, and the present theorem is the base case of that programme.

The theorem is proved in the paper; it has not been formalized. Mathlib has Turán's theorem (SimpleGraph.CliqueFree.card_edgeFinset_le), Tutte's theorem on perfect matchings, and extremal-graph infrastructure (SimpleGraph.IsExtremal), but no Tutte–Berge formula, no Gallai–Edmonds decomposition, and no Erdős–Gallai theorem for matchings. A formal proof would add the Tutte–Berge formula in the barrier form used here, a reusable Zykov symmetrization argument, and the first machine-checked proof of the Erdős–Gallai matching bound as a special case.

Difficulty

The arithmetic of the proof (milestones 8–10) is short. The difficulty is concentrated elsewhere.

First, the proof starts from a barrier BBB whose removal leaves only odd components with ∣B∣+∑i(∣Ai∣−1)/2|B|+\sum_i(|A_i|-1)/2∣B∣+∑i​(∣Ai​∣−1)/2 equal to the matching number. This is the Tutte–Berge formula sharpened to a maximal barrier, all of whose components are factor-critical. Mathlib's Tutte theorem covers perfect matchings only, so the deficiency version and the structure of maximal barriers have to be built.

Second, Lemma 2.1, Claim 2.3 and Lemma 2.2 are surgery on an extremal graph. Each modifies the graph (copying a neighbourhood, swapping two classes of BBB in the neighbourhood of a component, merging two components) and must re-verify three things: the number of edges does not decrease, no clique on k+1k+1k+1 vertices appears, and the matching number stays at most sss. The matching bound is re-verified through the barrier, so the components of the modified graph have to be tracked, and the extremality hypothesis has to be used in the right form. In particular, the paper's printed proof of Lemma 2.1 shows only that non-adjacency is transitive on BBB and that non-adjacent vertices of BBB have equal degree; the full "same neighbourhood" conclusion needs a further argument from extremality.

Third, the barrier equation of p. 2 is written with sss, which presumes that an extremal graph has matching number exactly sss; the paper does not prove this, and a solver who wants to chain the milestones needs it.

Formalization scope

Graphs are SimpleGraph (Fin n); edge counts are #G.edgeFinset. Clique number at most kkk is G.CliqueFree (k + 1). The matching number is the supremum, in ℕ, of the edge counts of matching subgraphs (Subgraph.IsMatching); the set is nonempty and bounded, so the supremum is a maximum, and the companion sanity file proves ν(G)≤s\nu(G)\le sν(G)≤s iff every matching has at most sss edges. T(n,k)T(n,k)T(n,k) is Mathlib's turanGraph n k. G(n,k,s)G(n,k,s)G(n,k,s) is an explicit graph on Fin n (vertex v<sv<sv<s in class v mod (k−1)v \bmod (k-1)vmod(k−1), all others in class k−1k-1k−1), and g(n,k,s)g(n,k,s)g(n,k,s) is defined as its edge count, not by a closed formula, so the identity of milestone 8 is a statement about the graph. Extremality is Mathlib's IsExtremal for the admissibility property. The maximality of ∑iai2\sum_i a_i^2∑i​ai2​ in Lemmas 2.1 and 2.2 and Claim 2.3 ranges over all extremal graphs together with all their odd barriers of value sss, because the paper chooses BBB together with GGG.

One hypothesis is added to the paper's statement: k≥2k\ge2k≥2. The paper says "every kkk", but G(n,k,s)G(n,k,s)G(n,k,s) has k−1k-1k−1 classes of total size sss, which is impossible for k=1k=1k=1 and s>0s>0s>0, and no nonempty graph has clique number at most 000. Lemmas 2.1, 2.2 and Claim 2.3 keep the paper's barrier value sss as a hypothesis, while the barrier milestone itself is stated with ν(G)\nu(G)ν(G).

The goal is an upper bound together with an attaining graph; an upper bound alone would be a weaker theorem. The surgery lemmas carry extremality as a hypothesis, never a hypothesis on the edge count that would presuppose the answer.

Contributions welcome beyond this mission: the Tutte–Berge formula and Gallai–Edmonds decomposition for SimpleGraph, the matching number as a reusable definition, and the Erdős–Gallai matching theorem.

Selected references

  • N. Alon and P. Frankl, Turán graphs with bounded matching number, arXiv:2210.15076v1, 2022; J. Combin. Theory Ser. B, 2024. https://arxiv.org/abs/2210.15076, https://doi.org/10.1016/j.jctb.2023.12.002
  • P. Turán, On an extremal problem in graph theory (in Hungarian), Mat. Fiz. Lapok 48 (1941), 436–452.
  • P. Erdős and T. Gallai, On maximal paths and circuits of graphs, Acta Math. Acad. Sci. Hungar. 10 (1959), 337–356. https://doi.org/10.1007/BF02024498
  • L. Lovász and M. D. Plummer, Matching Theory, North-Holland Mathematics Studies 121, 1986. https://doi.org/10.1090/chel/367
  • A. A. Zykov, On some properties of linear complexes, Mat. Sbornik N.S. 24(66) (1949), 163–188.
12 thms1 active userReviewed
🏆Completed
CombinatoricsLinear OptimizationOptimal Transport·Captain: mikedeng1

Computational Optimal Transport IV: The Support Graph of an Extreme Point of the Transportation Polytope U(a, b) Has No Cycle, Hence at Most n + m − 1 Nonzero EntriesTextbook

Motivation

Discrete optimal transport between two histograms is a linear program. Its feasible set, the transportation polytope, has been studied since Hitchcock (1941) and Kantorovich, and the structure of its vertices underlies the classical algorithms for the problem: the north-west corner rule produces a vertex, and the network simplex moves from vertex to vertex. A linear program with a nonempty bounded feasible set attains its minimum at a vertex, so knowing what vertices look like tells us what optimal transport plans can be assumed to look like.

This mission is the fourth of a series formalizing G. Peyré and M. Cuturi, Computational Optimal Transport (Foundations and Trends in Machine Learning, 2019). It covers §3.4.1 of the book, Tree Structure of the Support of All Vertices of U(a, b) (pp. 405–407), whose single numbered result, Proposition 3.4, states that the support of a vertex is a forest. The book credits the result to Brualdi, Combinatorial Matrix Classes (2006, Theorem 8.1.2).

Setting

Fix integers n,mn, mn,m and histograms a∈Σna \in \Sigma_na∈Σn​, b∈Σmb \in \Sigma_mb∈Σm​. The transportation polytope is

U(a,b)={P∈R+n×m:P1m=a, P⊤1n=b},U(a,b) = \{ P \in \mathbb R_+^{n\times m} : P\mathbb 1_m = a,\ P^\top\mathbb 1_n = b \},U(a,b)={P∈R+n×m​:P1m​=a, P⊤1n​=b},

the set of nonnegative n×mn\times mn×m matrices whose row sums are a1,…,ana_1,\dots,a_na1​,…,an​ and whose column sums are b1,…,bmb_1,\dots,b_mb1​,…,bm​. Histograms are nonnegative and sum to one; U(a,b)U(a,b)U(a,b) is the set of couplings between them. The cost of a plan PPP for a cost matrix CCC is ⟨C,P⟩=∑i,jCijPij\langle C, P\rangle = \sum_{i,j} C_{ij}P_{ij}⟨C,P⟩=∑i,j​Cij​Pij​.

A point xxx of a set SSS is an extremal point (vertex) of SSS if, whenever y,z∈Sy, z \in Sy,z∈S and x=(y+z)/2x = (y+z)/2x=(y+z)/2, necessarily x=y=zx = y = zx=y=z.

Take nnn source nodes V={1,…,n}V = \{1,\dots,n\}V={1,…,n} and mmm target nodes V′={1′,…,m′}V' = \{1',\dots,m'\}V′={1′,…,m′}. The complete bipartite graph between them has the nmnmnm edges (i,j′)(i, j')(i,j′). For a matrix PPP, the support S(P)S(P)S(P) is the set of edges (i,j′)(i,j')(i,j′) with Pij>0P_{ij} > 0Pij​>0, and the support graph is G(P)=(V∪V′,S(P))G(P) = (V\cup V', S(P))G(P)=(V∪V′,S(P)). A transport plan is a flow on the complete bipartite graph, with aia_iai​ leaving node iii and bjb_jbj​ entering node j′j'j′; G(P)G(P)G(P) records the edges the flow actually uses.

Formalization targets

Goal: Proposition 3.4 (p. 406)

If PPP is an extremal point of U(a,b)U(a,b)U(a,b), then

G(P) has no cyclesand#{(i,j):Pij≠0}≤n+m−1.G(P) \text{ has no cycles}\qquad\text{and}\qquad \#\{(i,j) : P_{ij}\ne 0\} \le n + m - 1 .G(P) has no cyclesand#{(i,j):Pij​=0}≤n+m−1.

Milestones (proof of Proposition 3.4, p. 407)

  1. If G(P)G(P)G(P) has a cycle, there is a nonzero matrix EEE with E1m=0E\mathbb 1_m = 0E1m​=0, E⊤1n=0E^\top\mathbb 1_n = 0E⊤1n​=0, and Eij≠0⇒Pij>0E_{ij} \ne 0 \Rightarrow P_{ij} > 0Eij​=0⇒Pij​>0.
  2. For such an EEE and P∈U(a,b)P \in U(a,b)P∈U(a,b), both P+tEP + tEP+tE and P−tEP - tEP−tE lie in U(a,b)U(a,b)U(a,b) for some t>0t > 0t>0.
  3. A graph with kkk nodes and no cycles has at most k−1k-1k−1 edges.
  4. For P≥0P \ge 0P≥0, the edges of G(P)G(P)G(P) correspond one-to-one to the nonzero entries of PPP.

Companion (§3.4, p. 405)

For a∈Σna \in \Sigma_na∈Σn​, b∈Σmb \in \Sigma_mb∈Σm​ and any cost CCC, some minimizer of ⟨C,P⟩\langle C, P\rangle⟨C,P⟩ over U(a,b)U(a,b)U(a,b) is an extremal point of U(a,b)U(a,b)U(a,b).

Significance

Proposition 3.4 is the sparsity statement of discrete optimal transport: a vertex of U(a,b)U(a,b)U(a,b) moves mass along at most n+m−1n+m-1n+m−1 source–target pairs, out of nmnmnm possible. With the companion, some optimal plan is that sparse. This is what makes the network simplex method a combinatorial algorithm on spanning trees of the bipartite graph (§3.5 of the book), and it is the reason the north-west corner rule (§3.4.2) can produce a feasible vertex with a tree support. The same fact, specialized to n=mn = mn=m and uniform histograms, is one step on the way from Kantorovich's relaxation back to Monge's assignment problem.

The result is classical and proved; it has not, to our knowledge, been formalized. Mathlib has the transportation polytope only implicitly (no dedicated definition), and it has the edge count of trees (SimpleGraph.IsTree.card_edgeFinset) but not the corresponding bound for forests. The mission produces a machine-checked link between the convex geometry of U(a,b)U(a,b)U(a,b) and the combinatorics of its support graph, reusable by later missions of the series on the network simplex and on the assignment problem.

Difficulty

The convex-geometric half is elementary; the work is in the passage between a cycle of a graph and a matrix. The cycle is an object of Mathlib's graph library (a closed walk with distinct vertices in a graph on V⊔V′V \sqcup V'V⊔V′), while the perturbation is a matrix indexed by V×V′V \times V'V×V′. An informal picture of a cycle as a list of edges hides what the matrix statement needs: each source and each target must be visited at most once, and Mathlib's cycles are walks whose distinctness conditions have to be carried to the matrix level.

The edge bound for forests is not packaged in Mathlib either: it counts the edges of a tree, a connected acyclic graph, while the support graph of a vertex is in general disconnected, and the empty graph must be handled.

Formalization scope

All objects live in the namespace CompOT.Vertices. Indices are Fin n and Fin m (0-based; the book's [ ⁣[n] ⁣]={1,…,n}[\![n]\!] = \{1,\dots,n\}[[n]]={1,…,n}). Matrices are Matrix (Fin n) (Fin m) ℝ, and U(a,b)U(a,b)U(a,b) is the set of matrices with nonnegative entries, row sums aaa and column sums bbb. Extremality is the book's midpoint definition, relative to the set: the points y,zy, zy,z in x=(y+z)/2x = (y+z)/2x=(y+z)/2 range over U(a,b)U(a,b)U(a,b) only. This rules out the trivializing reading in which y,zy, zy,z range over all matrices (no point would then be extremal and the goal would be vacuous); conversely, the sanity file shows a non-extremal coupling, so extremality is not automatic either.

The support graph is a SimpleGraph (Fin n ⊕ Fin m) with Sum.inl i adjacent to Sum.inr j exactly when Pij>0P_{ij} > 0Pij​>0 and no other adjacencies. The book calls the edges (i,j′)(i,j')(i,j′) directed, but the cycles of its proof alternate between sources and targets, so "no cycles" is Mathlib's SimpleGraph.IsAcyclic of this undirected graph. The number of nonzero entries is the cardinality of {(i,j):Pij≠0}\{(i,j) : P_{ij} \ne 0\}{(i,j):Pij​=0}.

The goal, feasibility milestone and companion carry a∈Σna \in \Sigma_na∈Σn​, b∈Σmb \in \Sigma_mb∈Σm​, the book's standing notation for histograms (p. 360). These hypotheses make the index sets nonempty, so natural-number subtraction in the bound n+m−1n+m-1n+m−1 agrees with the book's integer expression.

Useful infrastructure, reusable beyond this mission: the forest edge bound for finite simple graphs, and the correspondence between cycles of a bipartite graph and signed matrices with zero line sums. Contributions of either as standalone lemmas are welcome.

Selected references

  • G. Peyré, M. Cuturi, Computational Optimal Transport, Foundations and Trends in Machine Learning 11(5–6):355–607, 2019. https://doi.org/10.1561/2200000073
  • R. A. Brualdi, Combinatorial Matrix Classes, Cambridge University Press, 2006, §8.1. https://doi.org/10.1017/CBO9780511721182
  • D. Bertsimas, J. N. Tsitsiklis, Introduction to Linear Optimization, Athena Scientific, 1997, Theorem 2.7.
  • F. L. Hitchcock, The distribution of a product from several sources to numerous localities, Journal of Mathematics and Physics 20:224–230, 1941. https://doi.org/10.1002/sapm1941201224
7 thms1 active userReviewed
CombinatoricsComplexity TheoryTheoretical Computer Science·Captain: mikedeng1

Algebraic Approach to Promise Constraint Satisfaction 2: Pol(K_k, K_{2k−1}) Has No Olšák Function, So It Maps by a Minion Homomorphism to Pol(H₂, H_K)Research Paper

Approximate graph colouring

Given a graph that is promised to be kkk-colourable, how many colours are needed to find a proper colouring in polynomial time? The decision version asks to distinguish kkk-colourable graphs from graphs that are not even ccc-colourable. It has been studied since Garey and Johnson (1976) and is conjectured to be NP-hard for all fixed 3≤k≤c3 \le k \le c3≤k≤c. This is the approximate graph colouring problem, and it is one of the motivating examples of promise constraint satisfaction problems (PCSPs).

Timeline of the NP-hardness results without extra assumptions, as surveyed in Example 2.9 and §6.2 of the paper:

  • 2000–2004, Khanna–Linial–Safra and Guruswami–Khanna: 3 versus 4 colours.
  • 2013, Huang (cited as [Hua13] in the paper): kkk versus 2Ω(k1/3)2^{\Omega(k^{1/3})}2Ω(k1/3) colours for large kkk.
  • 2016, Brakensiek–Guruswami: kkk versus 2k−22k - 22k−2 colours for every k≥3k \ge 3k≥3, via polymorphisms.
  • 2018/2019, Barto–Bulín–Krokhin–Opršal (this paper): kkk versus 2k−12k - 12k−1 colours for every k≥3k \ge 3k≥3, in particular 3 versus 5. The argument is an instance of a general algebraic theory of PCSPs.

The method also applies to the question of 5-cycle versus 3-colouring, which Brakensiek and Guruswami had raised as an open case.

Setting

A relational structure A\mathbf AA consists of a finite domain AAA and relations RA⊆Aar(R)R^{\mathbf A} \subseteq A^{\mathrm{ar}(R)}RA⊆Aar(R), one for each symbol RRR of a signature. A PCSP template is a pair (A,B)(\mathbf A, \mathbf B)(A,B) of structures with the same signature and a homomorphism A→B\mathbf A \to \mathbf BA→B. Write Ek={0,…,k−1}E_k = \{0, \dots, k-1\}Ek​={0,…,k−1}.

  • Kk=(Ek;≠k)\mathbf K_k = (E_k; \neq_k)Kk​=(Ek​;=k​) is the complete graph: one binary relation {(a,b)∣a≠b}\{(a, b) \mid a \ne b\}{(a,b)∣a=b}.
  • Hk=(Ek;NAEk)\mathbf H_k = (E_k; \mathrm{NAE}_k)Hk​=(Ek​;NAEk​) carries the ternary not-all-equal relation Ek3∖{(a,a,a)}E_k^3 \setminus \{(a, a, a)\}Ek3​∖{(a,a,a)}; PCSP(Hk,Hc)\mathrm{PCSP}(\mathbf H_k, \mathbf H_c)PCSP(Hk​,Hc​) is approximate colouring of 3-uniform hypergraphs.

An nnn-ary polymorphism from A\mathbf AA to B\mathbf BB is a map f:An→Bf : A^n \to Bf:An→B that, applied coordinatewise to any nnn tuples of RAR^{\mathbf A}RA, returns a tuple of RBR^{\mathbf B}RB. For Kk\mathbf K_kKk​ and Kc\mathbf K_cKc​ these are exactly the ccc-colourings of the power Kkn\mathbf K_k^nKkn​. The set Pol(A,B)\mathrm{Pol}(\mathbf A, \mathbf B)Pol(A,B) of polymorphisms of all arities n≥1n \ge 1n≥1 is closed under minors: f(x1,…,xn)=g(xπ(1),…,xπ(m))f(x_1, \dots, x_n) = g(x_{\pi(1)}, \dots, x_{\pi(m)})f(x1​,…,xn​)=g(xπ(1)​,…,xπ(m)​) for π:[m]→[n]\pi : [m] \to [n]π:[m]→[n]. A minion is any nonempty set of functions of arities n≥1n \ge 1n≥1 that is closed under minors. A minion homomorphism ξ:M→N\xi : \mathscr M \to \mathscr Nξ:M→N preserves arities and commutes with minors. Write HK=Pol(H2,HK)\mathscr H_K = \mathrm{Pol}(\mathbf H_2, \mathbf H_K)HK​=Pol(H2​,HK​).

An Olšák function is a 6-ary function ooo with

o(x,x,y,y,y,x)=o(x,y,x,y,x,y)=o(y,x,x,x,y,y)for all x,y.o(x, x, y, y, y, x) = o(x, y, x, y, x, y) = o(y, x, x, x, y, y) \quad \text{for all } x, y .o(x,x,y,y,y,x)=o(x,y,x,y,x,y)=o(y,x,x,x,y,y)for all x,y.

The free structure FM(A)\mathbf F_{\mathscr M}(\mathbf A)FM​(A) has the ∣A∣|A|∣A∣-ary members of M\mathscr MM as its domain and relations defined by minor identities that follow the tuples of A\mathbf AA.

Formalization targets

Goal: the algebraic core of Theorem 6.5

For every k≥3k \ge 3k≥3,

∃ K≥2:there is a minion homomorphism Pol(Kk,K2k−1)→Pol(H2,HK).\exists\, K \ge 2:\quad \text{there is a minion homomorphism } \mathrm{Pol}(\mathbf K_k, \mathbf K_{2k-1}) \to \mathrm{Pol}(\mathbf H_2, \mathbf H_K).∃K≥2:there is a minion homomorphism Pol(Kk​,K2k−1​)→Pol(H2​,HK​).

KKK is left unspecified, as in the paper. This is the statement from which the paper deduces that kkk versus 2k−12k - 12k−1 colouring is NP-hard.

Milestones

  1. Example 3.5: the Olšák minor condition is not satisfied in Hk\mathscr H_kHk​ for any k≥2k \ge 2k≥2.
  2. Lemma 4.4: homomorphisms FM(A)→B\mathbf F_{\mathscr M}(\mathbf A) \to \mathbf BFM​(A)→B correspond one-to-one to minion homomorphisms M→Pol(A,B)\mathscr M \to \mathrm{Pol}(\mathbf A, \mathbf B)M→Pol(A,B).
  3. Theorem 6.2: for a minion M\mathscr MM on finite sets,
(∃K≥2, ∃ minion homomorphism M→HK)  ⟺  M contains no Olsˇaˊk function.\bigl(\exists K \ge 2,\ \exists\ \text{minion homomorphism } \mathscr M \to \mathscr H_K\bigr) \iff \mathscr M \text{ contains no Olšák function}.(∃K≥2, ∃ minion homomorphism M→HK​)⟺M contains no Olsˇaˊk function.
  1. Lemma 6.4: Pol(Kk,K2k−1)\mathrm{Pol}(\mathbf K_k, \mathbf K_{2k-1})Pol(Kk​,K2k−1​) contains no Olšák function, k≥3k \ge 3k≥3.

One further item, not a milestone, shows the bound is sharp. Proposition 10.1 says that Pol(Kk,K2k)\mathrm{Pol}(\mathbf K_k, \mathbf K_{2k})Pol(Kk​,K2k​) contains an Olšák function.

Significance

Theorem 6.2 is a criterion stated entirely inside the minion. Together with the reduction theory of the paper (a minion homomorphism Pol(A,B)→Pol(A′,B′)\mathrm{Pol}(\mathbf A, \mathbf B) \to \mathrm{Pol}(\mathbf A', \mathbf B')Pol(A,B)→Pol(A′,B′) gives a log-space reduction from PCSP(A′,B′)\mathrm{PCSP}(\mathbf A', \mathbf B')PCSP(A′,B′) to PCSP(A,B)\mathrm{PCSP}(\mathbf A, \mathbf B)PCSP(A,B)) and the NP-hardness of approximate hypergraph colouring (Dinur, Regev and Smyth 2005), it shows that every finite template without an Olšák polymorphism has an NP-hard PCSP (Corollary 6.3). Lemma 6.4 then settles the 3 versus 5 colouring case and improves the 2k−22k - 22k−2 bound to 2k−12k - 12k−1 for every k≥3k \ge 3k≥3. Proposition 10.1 shows that this criterion cannot go beyond 2k−12k - 12k−1.

The results are proved in the paper. None of them has a machine-checked proof that this mission knows of. This mission produces a Lean statement of the minion side of the argument. That includes the free-structure correspondence (Lemma 4.4), which every later application of the theory uses, and Theorem 6.2 for arbitrary finite minions. A finished formalization would include the clique argument behind Lemma 6.4 and would check the key hardness inequality c≥2k−1c \ge 2k - 1c≥2k−1 of the paper by machine.

Difficulty

Lemma 6.4 is a statement about every (2k−1)(2k-1)(2k−1)-colouring of Kk6\mathbf K_k^6Kk6​ that is constant on certain triples of vertices. The power graph has k6k^6k6 vertices, so exhaustive search is out of reach for general kkk, and case checks for k=3k = 3k=3 do not give the general claim. Showing that 2k−12k - 12k−1 colours do not suffice requires an explicit obstruction, uniform in kkk, in the graph obtained by identifying those triples; the identification creates edges that are not present in Kk6\mathbf K_k^6Kk6​, and keeping track of them is the combinatorial core. The bound is tight (Proposition 10.1), so no argument that would also apply to 2k2k2k colours can work.

In Theorem 6.2 the direction from "no Olšák function" to a minion homomorphism must build a minion homomorphism into HK\mathscr H_KHK​ out of nothing but the absence of a single identity. That direction depends on Lemma 4.4, which relates maps defined on functions of a single arity to minion homomorphisms defined on all arities; in Lean the bookkeeping of enumerations, minors and dependent arities is substantial.

Formalization scope

The development reuses the published relational-structure layer PCSPBLPAff_Symmetric_Setting (RelStruct, IsHom, IsPolymorphism, IsPromiseTemplate). Domains EkE_kEk​ are Fin k. Kk\mathbf K_kKk​ and Hk\mathbf H_kHk​ have one-symbol signatures of arity 2 and 3. Polymorphisms apply fff to the columns of an n×ar(R)n \times \mathrm{ar}(R)n×ar(R) matrix whose rows lie in RAR^{\mathbf A}RA, which is the paper's matrix transposed. A minion is a family mem n of nnn-ary functions (Fin n → A) → B with no nullary members, nonempty and closed under the minors fun x => g (x ∘ π). A minion homomorphism is a family of maps that preserves membership and minors on members, with unconstrained values outside the minion, so Lemma 4.4's "one-to-one" is stated up to agreement on the minion. Theorem 6.2 assumes the minion's sets are finite, which is the paper's standing convention (Definition 2.20). The free structure takes an enumeration e:Fin n≃Ae : \mathrm{Fin}\ n \simeq Ae:Fin n≃A of the generating domain and quantifies the enumeration of each relation existentially. In the goal, k≥3k \ge 3k≥3 makes 2k−12k - 12k−1 free of natural-number truncation.

The codomain of the goal is exactly Pol(H2,HK)\mathrm{Pol}(\mathbf H_2, \mathbf H_K)Pol(H2​,HK​) with K≥2K \ge 2K≥2. A minion homomorphism into an arbitrary or existentially chosen minion, or into HK\mathscr H_KHK​ with K≤1K \le 1K≤1, would be trivial and is ruled out by the statement.

Not formalized. The NP-hardness conclusions of Theorem 6.5 and Corollary 6.3 are left out, and so are the log-space reduction of Theorem 3.1 and the cited hardness of approximate hypergraph colouring (Theorem 5.23, [DRS05]). There is no complexity-theoretic substrate in Mathlib for them. Lemma 6.6 (Pol(C5,K3)\mathrm{Pol}(\mathbf C_5, \mathbf K_3)Pol(C5​,K3​)) is a natural follow-up item. Contributions are welcome on every milestone. Lemma 4.4 and the minion layer can be reused beyond this mission, for any minion-homomorphism argument in the PCSP literature.

Selected references

  • L. Barto, J. Bulín, A. Krokhin, J. Opršal, Algebraic approach to promise constraint satisfaction, arXiv:1811.00970v3, 2019; J. ACM 68(4), 2021. https://arxiv.org/abs/1811.00970v3
  • J. Brakensiek, V. Guruswami, New hardness results for graph and hypergraph colorings, CCC 2016. https://doi.org/10.4230/LIPIcs.CCC.2016.14
  • I. Dinur, O. Regev, C. Smyth, The hardness of 3-uniform hypergraph coloring, Combinatorica 25(5), 2005. https://doi.org/10.1007/s00493-005-0032-4
  • M. Olšák, The weakest nontrivial idempotent equations, Bull. London Math. Soc. 49(6), 2017. https://doi.org/10.1112/blms.12097
  • S. Khanna, N. Linial, S. Safra, On the hardness of approximating the chromatic number, Combinatorica 20(3), 2000. https://doi.org/10.1007/s004930070013
9 thms1 active userReviewed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Generation and Properties of Snarks 3: A Snark on 26 Vertices All of Whose 2-Factors Consist of Odd Cycles Refutes the Abreu–Labbate–Sheehan ConjectureResearch Paper

Motivation

A cubic graph has chromatic index 3 or 4 (Vizing). Cubic graphs of chromatic index 4 — the uncolourable ones — are where most of the open problems on cubic graphs live: a minimal counterexample to the cycle double cover conjecture or to Tutte's 5-flow conjecture would be among them. Uncolourable graphs with small edge cuts or short cycles reduce to smaller ones, so attention concentrates on snarks, the uncolourable cubic graphs that are cyclically 4-edge connected and have girth at least 5.

Brinkmann, Goedgebeur, Hägglund and Markström (arXiv:1206.6690v3, J. Combin. Theory Ser. B 103, 2013) generated all snarks on up to 36 vertices and tested a series of published conjectures against them. One of these concerns the 2-factors of a snark, its spanning subgraphs in which every vertex has degree 2. A 3-edge-colouring of a cubic graph gives 2-factors made of even cycles, so a snark is a graph in which some 2-factors are forced to contain odd cycles. At the far end are the snarks in which every 2-factor consists of odd cycles only. Abreu, Labbate and Sheehan (European J. Combin. 33, 2012), studying pseudo and strongly pseudo 2-factor isomorphic graphs, conjectured that this happens only for the Petersen graph, the Blanuša-2 snark and the Flower snarks (paper, p. 11, Conjecture 4.11).

Timeline:

  • 2012: Abreu, Labbate and Sheehan state the conjecture.
  • 2012–2013: Brinkmann, Goedgebeur, Hägglund and Markström refute it with a snark on 26 vertices, printed in Appendix 8.1 of their paper (p. 29); by computer search they report it is the smallest counterexample and that one more exists on 34 vertices (Observation 4.12).

Setting

All graphs are finite and simple. A graph GGG is cubic if every vertex has degree 3. It is colourable if its edges can be coloured with three colours so that edges sharing an endpoint get different colours, and uncolourable otherwise. The girth g(G)g(G)g(G) is the number of vertices of a shortest cycle. GGG is cyclically kkk-edge connected if deleting fewer than kkk edges never leaves two components that both contain a cycle. A snark is an uncolourable, cyclically 4-edge connected cubic graph with g(G)≥5g(G) \ge 5g(G)≥5 (paper, p. 4).

A 2-factor of GGG is a spanning 2-regular subgraph FFF: it uses every vertex, and every vertex has degree exactly 2 in FFF. Each component of FFF is a cycle. Say that all 2-factors of GGG are odd if every component of every 2-factor of GGG has an odd number of vertices.

The witness is the graph G8.1G_{8.1}G8.1​ of Appendix 8.1, given in the paper as a list of higher-numbered neighbours of the vertices 1,…,261, \dots, 261,…,26:

{2,3,4,5,6,7,8,9,10,7,9,8,10,11,12,13,14,15,16,15,17,18,19,18,20,18,21,22,21,23,23,24,22,24,25,26,26,25,26}.\{2, 3, 4, 5, 6, 7, 8, 9, 10, 7, 9, 8, 10, 11, 12, 13, 14, 15, 16, 15, 17, 18, 19, 18, 20, 18, 21, 22, 21, 23, 23, 24, 22, 24, 25, 26, 26, 25, 26\}.{2,3,4,5,6,7,8,9,10,7,9,8,10,11,12,13,14,15,16,15,17,18,19,18,20,18,21,22,21,23,23,24,22,24,25,26,26,25,26}.

It has 26 vertices and 39 edges. In Lean it is G81 : SimpleGraph (Fin 26) with paper vertex iii as i−1i-1i−1.

Formalization targets

Goal: Observation 4.12, existence half

G8.1 is a snarkandevery 2-factor of G8.1 consists of odd cycles.G_{8.1} \text{ is a snark} \quad\text{and}\quad \text{every 2-factor of } G_{8.1} \text{ consists of odd cycles.}G8.1​ is a snarkandevery 2-factor of G8.1​ consists of odd cycles.

Because G8.1G_{8.1}G8.1​ has 26 vertices, it is not the Petersen graph (10 vertices), not the Blanuša-2 snark (18) and not a Flower snark JkJ_kJk​ (4k4k4k vertices). The goal therefore refutes Conjecture 4.11 without formalizing the conjecture's list of exceptions.

Milestones

The goal splits into five properties of the one graph G8.1G_{8.1}G8.1​:

  1. G8.1G_{8.1}G8.1​ is cubic;
  2. g(G8.1)≥5g(G_{8.1}) \ge 5g(G8.1​)≥5;
  3. G8.1G_{8.1}G8.1​ is cyclically 4-edge connected;
  4. every 2-factor of G8.1G_{8.1}G8.1​ consists of odd cycles;
  5. G8.1G_{8.1}G8.1​ is uncolourable.

A supporting statement, not from the paper, asserts that G8.1G_{8.1}G8.1​ has a 2-factor, so that milestone 4 is not vacuous.

Significance

The result. Conjecture 4.11 would have characterised a natural extremal class of snarks by a short explicit list. The counterexample shows that the class is not that list: there are snarks outside the Petersen, Blanuša-2 and Flower families in which no 2-factor contains an even cycle. Such graphs are relevant to questions about 2-factor structure, oddness and cycle double covers, where the parity of 2-factor cycles is exactly what the arguments track.

Formalizing it. The result is proved in the paper by exhibiting the graph and checking it by computer; it has not been machine-checked in a proof assistant. A formal proof certifies every property of the printed graph — including that the printed adjacency list decodes to a graph with these properties — and leaves a reusable library of definitions (snarks, cyclic edge connectivity, 2-factors) that other snark statements can import.

Difficulty

Each property is a finite check, but some are large. Cubicity and girth are local. Cyclic 4-edge connectivity quantifies over all sets of at most three of the 39 edges (about ten thousand) and needs, for each, a component analysis of the remaining graph. The odd-2-factor property quantifies over all 2-factors; in a cubic graph these are the complements of the perfect matchings, and G8.1G_{8.1}G8.1​ has 56 perfect matchings. The obvious route — decide over all subgraphs — is out of reach: there are 2392^{39}239 edge subsets, and the statement is phrased through Mathlib's ConnectedComponent and IsRegularOfDegree, which do not evaluate by kernel computation at this size. Uncolourability has the same problem: the naive search space has 3393^{39}339 colourings.

Formalization scope

Conventions committed to:

  • Vertices are Fin 26, paper vertex iii = Lean vertex i−1i-1i−1; adjacency is the symmetric closure of an explicit edge list G81Edges, so it is decidable.
  • Colourable means G.lineGraph.Colorable 3 (a proper 3-edge-colouring), never G.Colorable 3.
  • Girth is Mathlib's egirth : ℕ∞, so "girth at least 5" is 5 ≤ G.egirth.
  • A 2-factor is a graph F on the same vertex type with F ≤ G and F.IsRegularOfDegree 2 — spanning by construction. Oddness is per component: Odd (Nat.card c.supp) for every component c. It is not the parity of the whole vertex set (26 is even).
  • Cyclic kkk-edge connectivity: for every edge set S with S.card < k, any two components of G.deleteEdges S that both contain a cycle are equal.

The odd-2-factor property is vacuous for a graph without 2-factors. The goal is not trivialised this way: G8.1G_{8.1}G8.1​ has 2-factors, and a supporting theorem states it.

Out of scope, because they rest on exhaustive computer enumeration: that G8.1G_{8.1}G8.1​ is the smallest counterexample, that there is exactly one further counterexample on 34 vertices, Observation 4.10 (all snarks on at most 36 vertices have oddness 2), and Tables 4 and 5. The Petersen graph, the Blanuša-2 snark and the Flower snarks are not defined here; the refutation uses only their orders.

Contributions welcome: proofs of any milestone; general lemmas linking 2-factors of cubic graphs to perfect matchings, and 3-edge-colourings to even 2-factors; decision procedures for cyclic edge connectivity of explicit graphs. The definitions are shared in shape with the other missions of this series and are meant to be reused.

Selected references

  • G. Brinkmann, J. Goedgebeur, J. Hägglund, K. Markström, Generation and properties of snarks, arXiv:1206.6690v3, 2013; J. Combin. Theory Ser. B 103 (2013). https://arxiv.org/abs/1206.6690v3 , https://doi.org/10.1016/j.jctb.2013.05.001
  • M. Abreu, D. Labbate, J. Sheehan, Pseudo and strongly pseudo 2-factor isomorphic regular graphs and digraphs, European J. Combin. 33 (2012), no. 8, 1847–1856. https://mathscinet.ams.org/mathscinet-getitem?mr=2950486
  • G. Brinkmann, K. Coolsaet, J. Goedgebeur, H. Mélot, House of Graphs: a database of interesting graphs, Discrete Appl. Math. 161 (2013), 311–314. https://mathscinet.ams.org/mathscinet-getitem?mr=2973372
12 thms1 active userReviewed
CombinatoricsTheoretical Computer Science·Captain: mikedeng1

Generation and Properties of Snarks 2: A Cyclically 5-Edge-Connected Permutation Snark on 34 Vertices Refutes Zhang's Conjecture That Only the Petersen Graph Is OneResearch Paper

Motivation

A snark is a cubic graph that cannot be properly edge-coloured with three colours and that does not reduce to a smaller such graph. Snarks are the smallest possible counterexamples to several central conjectures about cubic graphs, among them the cycle double cover conjecture, Tutte's 5-flow conjecture and the Berge–Fulkerson conjecture, so structural statements about snarks narrow the search for counterexamples to those conjectures.

Permutation graphs are cubic graphs built from two disjoint cycles of equal length joined by a perfect matching, so that the two cycles are induced. The Petersen graph is the classical example: its outer 5-cycle and inner pentagram are joined by five spokes. In his monograph Integer Flows and Cycle Covers of Graphs (1997), C.-Q. Zhang conjectured that the Petersen graph is the only permutation snark that is also cyclically 5-edge connected.

Brinkmann, Goedgebeur, Hägglund and Markström (arXiv:1206.6690v3, J. Combin. Theory Ser. B 103, 2013) built a generator for snarks and enumerated all snarks up to 36 vertices. Among them they found twelve cyclically 5-edge-connected permutation snarks on 34 vertices, which refute Zhang's conjecture (Observation 4.2), and they printed all twelve in Appendix 8.6 of the paper. This mission formalizes the refutation using the first of those graphs.

Setting

All graphs are finite and simple. A cycle is a 2-regular connected graph. A graph is cubic if every vertex has degree 3. A 2-factor of GGG is a spanning 2-regular subgraph. The girth g(G)g(G)g(G) is the number of vertices in a shortest cycle of GGG.

A cubic graph is colourable if it has a proper 3-edge-colouring: a map c:E(G)→{1,2,3}c : E(G) \to \{1,2,3\}c:E(G)→{1,2,3} giving distinct colours to distinct edges with a common endpoint. Otherwise it is uncolourable.

A graph GGG is cyclically kkk-edge connected if deleting fewer than kkk edges from GGG never leaves two distinct components that both contain a cycle.

A snark is an uncolourable, cyclically 4-edge connected cubic graph with g(G)≥5g(G) \ge 5g(G)≥5.

A cubic graph is a permutation graph if it has a 2-factor consisting of two induced cycles. Equivalently, its vertex set splits as A∪(V∖A)A \cup (V \setminus A)A∪(V∖A) with both induced subgraphs G[A]G[A]G[A] and G[V∖A]G[V \setminus A]G[V∖A] cycles. A permutation snark is a snark that is a permutation graph.

The Petersen graph PPP is the graph on {0,…,9}\{0,\dots,9\}{0,…,9} with the outer cycle i∼i+1(mod5)i \sim i+1 \pmod 5i∼i+1(mod5), the spokes i∼i+5i \sim i+5i∼i+5, and the inner pentagram 5+i∼5+((i+2) mod 5)5+i \sim 5+((i+2) \bmod 5)5+i∼5+((i+2)mod5).

Zhang's Conjecture 4.1. Let GGG be a cubic cyclically 5-edge-connected permutation graph. If GGG is a snark, then GGG must be the Petersen graph.

The witness G1G_1G1​ is the first graph in Appendix 8.6 (p. 36). It is given as a list of higher-numbered neighbours of the vertices 1,…,341,\dots,341,…,34; decoded, it has 51 edges. In the Lean development the paper's vertex iii is i - 1 : Fin 34.

Formalization targets

Goal: Observation 4.2, first sentence

¬ ZhangConjecture,\neg\,\mathrm{ZhangConjecture},¬ZhangConjecture,

that is, there is a finite cubic, cyclically 5-edge connected permutation snark that is not isomorphic to the Petersen graph. This is SnarkGen.Zhang.observation_4_2.

Milestones: the properties of G1G_1G1​

Each milestone is a property of G1G_1G1​ asserted by the heading of Appendix 8.6, "The 12 cyclically 5-edge connected permutation snarks on 34 vertices":

  1. G1G_1G1​ is cubic;
  2. g(G1)≥5g(G_1) \ge 5g(G1​)≥5;
  3. G1G_1G1​ is cyclically 5-edge connected;
  4. G1G_1G1​ is a permutation graph;
  5. G1G_1G1​ is uncolourable;
  6. G1G_1G1​ is a snark;
  7. G1G_1G1​ satisfies all hypotheses of Conjecture 4.1 and G1≇PG_1 \not\cong PG1​≅P.

The goal follows from milestone 7 directly. Two supporting theorems are included: cyclic kkk-edge connectivity is monotone in kkk, and the paper's remark (p. 9) that every permutation snark has order ≡2(mod4)\equiv 2 \pmod 4≡2(mod4).

Significance

Zhang's conjecture predicted that the Petersen graph is isolated among highly connected permutation snarks. Its refutation shows that this is false, and the twelve 34-vertex graphs went on to serve as counterexamples to further conjectures in the same paper (Conjectures 5.17, 5.20 and 5.21, about cycle double covers containing a prescribed 2-factor). They are also listed in the House of Graphs database.

The refutation in the paper is a computer result: the graphs were found by an exhaustive search, and their properties were checked by programs. A formal proof would replace that check with a machine-verified one that is independent of the generator. It would also produce reusable certified tools for finite cubic graphs: a kernel-checked test for uncolourability, cyclic edge connectivity, and girth. No machine-checked proof of any of these properties of G1G_1G1​ is known to exist.

Difficulty

The statements are short, but none of the properties of G1G_1G1​ follows from a general theorem; each must be checked on the specific graph. The obvious approach, to unfold the definitions and let a decision procedure enumerate, does not scale:

  • Uncolourability quantifies over all 3513^{51}351 edge colourings.
  • Cyclic 5-edge connectivity quantifies over all edge sets of size at most 4 (about 270,000 of them among 51 edges), and for each one over the components of the remaining graph and the cycles inside them.
  • The permutation structure requires exhibiting the two induced 17-cycles and proving that each induced subgraph is connected and 2-regular.
  • Girth requires ruling out all cycles of length 3 and 4.

The definitions are stated in terms of Mathlib's general notions (line graphs, walks, connected components), not in terms of a computable representation. Connecting them to an efficient computation is the main work.

Formalization scope

Graphs are SimpleGraph V with [Fintype V]. The conventions are as follows.

  • "Colourable" is G.lineGraph.Colorable 3, a proper 3-colouring of the line graph. It is not a vertex colouring of GGG.
  • Cyclic kkk-edge connectivity quantifies over finite edge sets S ⊆ G.edgeSet with S.card < k. Components and cycles are taken in G.deleteEdges S.
  • Girth is 5 ≤ G.egirth, the extended girth, which is ∞\infty∞ for acyclic graphs.
  • Degrees in the general predicates use classical decidability. Any other decidability instance gives the same degrees.
  • The conjecture quantifies over all finite vertex types in Type and reads "must be the Petersen graph" as Nonempty (H ≃g petersen).

The conjecture keeps every printed hypothesis, including the redundant "cubic". Dropping "permutation graph" or "cyclically 5-edge connected" would make the refutation easier, since other cyclically 5-edge-connected snarks are known. The cyclic connectivity must be computed in the graph after deletion and must require both components to contain a cycle; otherwise it degenerates into ordinary edge connectivity, which is at most 3 for cubic graphs and would make the conjecture vacuous.

Proofs should be checked by the Lean kernel. The decision procedures for colourability and cyclic connectivity are the reusable part of this mission, and contributions that build them for arbitrary finite cubic graphs are welcome.

Out of scope.

  • The second and third sentences of Observation 4.2 ("The smallest counterexamples have 34 vertices, and there are exactly 12 counterexamples of that order"). They rest on an exhaustive enumeration of all snarks up to 36 vertices.
  • The other eleven graphs of Appendix 8.6.
  • Conjectures 5.17, 5.20 and 5.21. Their refutation needs the "no cycle double cover contains this 2-factor" certificate of Observation 5.18.
  • The counts in Table 3.

Selected references

  • G. Brinkmann, J. Goedgebeur, J. Hägglund, K. Markström, Generation and properties of snarks, arXiv:1206.6690v3, 2013; J. Combin. Theory Ser. B 103 (2013). https://arxiv.org/abs/1206.6690v3
  • C.-Q. Zhang, Integer Flows and Cycle Covers of Graphs, Monographs and Textbooks in Pure and Applied Mathematics 205, Marcel Dekker, 1997.
  • G. Brinkmann, K. Coolsaet, J. Goedgebeur, H. Mélot, House of Graphs: a database of interesting graphs, Discrete Appl. Math. 161 (2013), 311–314. https://doi.org/10.1016/j.dam.2012.07.018
15 thms1 active userReviewed
PreviousPage 11 of 14Next

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me