# Determinant-Preserving Sparsification of SDDM Matrices

## Abstract

*Combin. Probab. Comput.,*3 (1994), pp. 97--126]. This leads to a routine that, in quadratic time, sparsifies a graph down to about $n^{1.5}$ edges in a way that preserves both the determinant and the distribution of spanning trees (provided the sparsified graph is viewed as a random object). Extending this algorithm to work with Schur complements and approximate Choleksy factorizations leads to algorithms for counting and sampling spanning trees which are nearly optimal for dense graphs. We give an algorithm that computes a $(1 \pm \delta)$ approximation to the determinant of any SDDM matrix with constant probability in about $n^2 \delta^{-2}$ time. This is the first routine for graphs that outperforms general-purpose routines for computing determinants of arbitrary matrices. We also give an algorithm that generates, in about $n^2 \delta^{-2}$ time, a spanning tree of a weighted undirected graph from a distribution with a total variation distance of $\delta$ from the $\boldsymbol{\mathit{w}}$-uniform distribution.

### Keywords

### MSC codes

## Get full access to this article

View all available purchase options and get full access to this article.

## References

*The random walk construction of uniform spanning trees and uniform labelled trees*, SIAM J. Discrete Math., 3, pp. 450--465, https://doi.org/10.1137/0403039.

*Maintaining information in fully dynamic trees with top trees*, ACM Trans. Algorithms (TALG), 1, pp. 243--264.

*An $O({log} \ n/{log} \ {log} \ n)$-approximation algorithm for the asymmetric traveling salesman problem*, in Proceedings of the Twenty-First Annual ACM-SIAM Symposium on Discrete Algorithms (SODA '10), SIAM, Philadelphia, pp. 379--389, https://doi.org/10.1137/1.9781611973075.32.

*The complexity of partial derivatives*, Theoret. Comput. Sci., 22, pp. 317--330.

*Approximating s-t minimum cuts in O͂$(n2)$ time*, in Proceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing (STOC '96), ACM, New York, pp. 47--55.

*A Randomized Algorithm for Approximating the Log Determinant of a Symmetric Positive Definite Matrix*, preprint, https://arxiv.org/abs/1503.00374v1.

*Generating random spanning trees*, in Proceedings of the 30th Annual Symposium on Foundations of Computer Science (FOCS 1989), IEEE, Piscataway, NJ, pp. 442--447.

*Local characteristics, entropy and limit theorems for spanning trees and domino tilings via transfer-impedances*, Ann. Probab., 21, pp. 1329--1371.

*Efficient sampling for Gaussian graphical models via spectral sparsification*, in Proceedings of the 28th Conference on Learning Theory, Proc. Mach. Learn. Res. 40, pp. 364--390, http://jmlr.org/proceedings/papers/v40/Cheng15.pdf.

*Nearly tight oblivious subspace embeddings by trace inequalities*, in Proceedings of the Twenty-Seventh Annual ACM-SIAM Symposium on Discrete Algorithms, SIAM, Philadelphia, pp. 278--287, https://doi.org/10.1137/1.9781611974331.ch21.

*Almost-linear-time algorithms for Markov chains and new spectral primitives for directed graphs*, in Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing (STOC 2017), ACM, New York, pp. 410--419.

*$\ell_p$ row sampling by Lewis weights*, in Proceedings of the Forty-Seventh Annual ACM Symposium on Theory of Computing (STOC '15), ACM, New York, pp. 183--192.

*Unranking and ranking spanning trees of a graph*, J. Algorithms, 10, pp. 271--286.

*Two algorithms for unranking arborescences*, J. Algorithms, 20, pp. 268--281.

*Sampling random spanning trees faster than matrix multiplication*, in Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing (STOC 2017), ACM, New York, pp. 730--742.

*Sparsification---a technique for speeding up dynamic graph algorithms*, J. ACM, 44, pp. 669--696.

*The average impedance of an electrical network*, in Reissner Anniversary Volume, Contributions to Applied Mechanics, J. W. Edwards, Ann Arbor, MI, pp. 333--340.

*A general framework for graph sparsification*, in Proceedings of the Forty-Third Annual ACM Symposium on Theory of Computing, ACM, New York, pp. 71--80.

*Expanders via random spanning trees*, in Proceedings of the Twentieth Annual ACM-SIAM Symposium on Discrete Algorithms (SODA '09), SIAM, Philadelphia, pp. 576--585, https://doi.org/10.1137/1.9781611973068.64.

*Random spanning tree*, J. Algorithms, 4, pp. 214--220.

*Large-scale log-determinant computation through stochastic Chebyshev expansions*, in Proceedings of the 32nd International Conference on Machine Learning, Lille, France, pp. 908--917.

*Generating random spanning trees via fast matrix multiplication*, in LATIN 2016: Theoretical Informatics, Lecture Notes in Comput. Sci. 9644, Springer, Berlin, pp. 522--535.

*Computing the Log-Determinant of Symmetric, Diagonally Dominant Matrices in Near-Linear Time*, preprint, http://arxiv.org/abs/1408.1693.

*Determinant Approximations*, preprint, https://arxiv.org/abs/1105.0437.

*The numbers of spanning trees, Hamilton cycles and perfect matchings in a random graph*, Combin. Probab. Comput., 3, pp. 97--126.

*Density Independent Algorithms for Sparsifying k-step Random Walks*, preprint, https://arxiv.org/abs/1702.06110.

*Faster generation of random spanning trees*, in Proceedings of the 50th Annual IEEE Symposium on Foundations of Computer Science (FOCS 2009), IEEE, Piscataway, NJ, pp. 13--21.

*Ueber die Aufl$\ddot{o}$sung der Gliechungen, auf welche man bei der Untersuchung der linearen Vertheilung galvanischer Str$\ddot{o}$me gef$\ddot{u}$hrt wird*, Ann. Phys. Chem., 148, pp. 497--508.

*Generating random combinatorial objects*, J. Algorithms, 11, pp. 185--207.

*Sparsified Cholesky and multigrid solvers for connection Laplacians*, in Proceedings of the 48th Annual ACM SIGACT Symposium on Theory of Computing (STOC 2016), ACM, New York, pp. 842--850.

*Probability on Trees and Networks*, Cambridge Ser. Statist. Probab. Math. 42, Cambridge University Press, New York.

*Fast generation of random spanning trees and the effective resistance metric*, in Proceedings of the Twenty-Sixth Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2015), pp. 2019--2036, https://doi.org/10.1137/1.9781611973730.134.

*An efficient parallel solver for SDD linear systems*, in Proceedings of the 46th Annual ACM Symposium on Theory of Computing (STOC '14), ACM, New York, pp. 333--342.

*An almost-linear time algorithm for uniform random spanning tree generation*, in Proceedings of the 50th Annual ACM SIGACT Symposium on Theory of Computing (STOC 2018), ACM, New York, pp. 214--227.

*Self-adjusting binary search trees*, J. ACM, 32, pp. 652--686.

*Graph sparsification by effective resistances*, SIAM J. Comput., 40, pp. 1913--1926, https://doi.org/10.1137/080734029.

*Spectral sparsification of graphs*, SIAM J. Comput., 40, pp. 981--1025, https://doi.org/10.1137/08074489X.

*Nearly linear time algorithms for preconditioning and solving symmetric, diagonally dominant linear systems*, SIAM J. Matrix Anal. Appl., 35, pp. 835--885, https://doi.org/10.1137/090771430.

*User-friendly tail bounds for sums of random matrices*, Found. Comput. Math., 12, pp. 389--434.

*$Lx = b$ Laplacian solvers and their algorithmic applications*, Found. Trends Theor. Comput. Sci., 8, pp. 1--141.

*Book review of S. Skiena, “Implementing Discrete Mathematics: Combinatorics and Graph Theory with Mathematica,” Addison-Wesley,*1990, Math. Gaz., 76, pp. 286--288.

*Multiplying matrices faster than Coppersmith-Winograd*, in Proceedings of the Forty-Fourth Annual ACM Symposium on Theory of Computing (STOC '12), ACM, New York, pp. 887--898, https://doi.org/10.1145/2213977.2214056.

## Information & Authors

### Information

#### Published In

#### Copyright

#### History

**Submitted**: 19 January 2018

**Accepted**: 17 October 2019

**Published online**: 10 March 2020

#### Keywords

#### MSC codes

### Authors

#### Funding Information

#### Funding Information

#### Funding Information

## Metrics & Citations

### Metrics

### Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

#### Cited By

- On Randomized Trace Estimates for Indefinite Matrices with an Application to DeterminantsFoundations of Computational Mathematics, Vol. 22, No. 3 | 9 July 2021
- Minor Sparsifiers and the Distributed Laplacian Paradigm2021 IEEE 62nd Annual Symposium on Foundations of Computer Science (FOCS) | 1 Feb 2022