Subcritical Blob Graph Sparsifier
Implementation & benchmark of arXiv:2607.21428 — Subcritical percolation and network archaeology on random recursive tree substrate networks
Source paper: Subcritical percolation and network archaeology on random recursive tree substrate networks arXiv:2607.21428 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Turn stochastic edge thinning into a mass-aware graph sparsifier: first form retained connected blobs, then connect blob summaries using the rank-one shortcut probability from the paper. Choose the retention rate below the explicit subcritical threshold so the contracted graph remains a collection of small decorations around leading blobs, reducing message-passing cost without creating an artificial giant component.
Formulas
Mathematical statement
After edge retention, \(\mathcal{B}_{n,i}\) is a retained backbone blob and \(W_{n,i}=|\mathcal{B}_{n,i}|\) is its vertex mass. Conditional on the blob masses, distinct blobs \(i,j\) are connected by a shortcut with probability \(1-(1-p\lambda/n)^{W_{n,i}W_{n,j}}\), approximately \(p\lambda W_{n,i}W_{n,j}/n\). Starting from a blob of mass \(W_{n,i}=O_{\mathbb P}(n^p)\), the mass in shortcut-distance shell \(r\) satisfies \(A_{n,r}^{(i)}=\rho(p,\lambda)^rW_{n,i}+o_{\mathbb P}(n^p)\), and the total decorated component mass is \(|\mathcal{C}_{n,i}^{\mathrm{full}}|=W_{n,i}/(1-\rho(p,\lambda))+o_{\mathbb P}(n^p)\), where \(\rho(p,\lambda)=p\lambda/(1-2p)<1\). In a GNN, this gives an explicit safe operating region for stochastic graph coarsening: aggregate vertices within blobs, represent each blob by a pooled state, and retain sampled inter-blob edges while monitoring the empirical branching ratio.
Implementation notes
(1) Integration point: replace the full adjacency used by a message-passing GNN with a per-batch coarsened graph. This targets large molecular, citation, interaction, or web graphs where shortcut-like edges create excessive neighborhoods. Estimate \(\lambda\) as \(2|E|/|V|\), or estimate it from the mean degree of a designated shortcut edge layer. Set \(p=\alpha/(\lambda+2)\) with safety factor \(\alpha\in[0.3,0.9]\). (2) Pseudocode: sample each edge independently with probability \(p\); compute retained connected components \(B_i\); replace each \(B_i\) by pooled feature \(h_i=|B_i|^{-1}\sum_{v\in B_i}H_v\) and store \(W_i=|B_i|\); for candidate blob pairs, sample an inter-blob edge with \(q_{ij}=1-(1-p\lambda/N)^{W_iW_j}\), where \(N=|V|\); run the GNN on the blob graph; broadcast blob outputs to member vertices or use a lightweight decoder. Avoid \(O(B^2)\) pair generation by considering only original edges crossing blobs or k-nearest neighbors of pooled features. (3) The exact probability and threshold come from the paper. Estimate the blob masses and \(\lambda\) directly per graph. Monitor empirical shell masses \(A_r\) with BFS and compute \(\hat\rho=\operatorname{median}_r(A_r/A_{r-1})\); reduce \(p\) if \(\hat\rho\geq 0.9\). (4) First experiment: compare a 4-layer GIN against Cluster-GCN and random edge dropout on synthetic power-law and recursive-tree graphs at matched validation accuracy. Sweep \(\alpha\), and record GPU memory, processed edges, wall-clock time, validation accuracy, and largest-component fraction. Success means at least 2x fewer message-passing edges or 30% lower memory with no more than one percentage point accuracy loss, while avoiding a dominant retained component.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.