Abstract:The diagonal entries of pseudoinverse of the Laplacian matrix of a graph appear in many important practical applications, since they contain much information of the graph and many relevant quantities can be expressed in terms of them, such as Kirchhoff index and current flow centrality. However, a na\"ıve approach for computing the diagonal of a matrix inverse has cubic computational complexity in terms of the matrix dimension, which is not acceptable for large graphs with millions of nodes. Thus, rigorous solutions to the diagonal of the Laplacian matrices for general graphs, even for particluar graphs are much less. In this paper, we propose a theoretically guaranteed estimation algorithm, which approximates all diagonal entries of the pseudoinverse of a graph Laplacian in nearly linear time with respect to the number of edges in the graph. We execute extensive experiments on real-life networks, which indicate that our algorithm is both efficient and accurate. Also, we determine exact expressions for the diagonal elements of pseudoinverse of the Laplacian matrices for Koch networks and uniform recursive trees, and compare them with those obtained by our approximation algorithm. Finally, we use our algorithm to evaluate the Kirchhoff index of three deterministic model networks, for which the Kirchhoff index can be rigorously determined. These results further show the effectiveness and efficiency of our algorithm.

What problem does this paper attempt to address?

The problem this paper attempts to solve is: how to efficiently and accurately compute the diagonal elements of the pseudo-inverse of the graph Laplacian matrix. ### Background and Motivation In many practical applications, the diagonal elements of the pseudo-inverse \( \mathbf{L}^\dagger \) of the graph Laplacian matrix \( \mathbf{L} \) contain a lot of important information about the graph. For example, these diagonal elements can be used to compute the Kirchhoff index of the graph, which is a measure of the overall connectivity of the network. However, directly computing the diagonal elements of \( \mathbf{L}^\dagger \) has a high computational complexity, especially for large graphs, with a complexity of \( O(N^3) \), which is unacceptable in practical applications. ### Solution To overcome this problem, the paper proposes a theoretically guaranteed approximation algorithm that can estimate all the diagonal elements of the pseudo-inverse of the graph Laplacian matrix in nearly linear time complexity \( O(M) \), where \( M \) is the number of edges in the graph. This algorithm not only improves computational efficiency but also provides error bound guarantees. ### Main Contributions 1. **Approximation Algorithm**: An approximation algorithm is proposed that can compute the diagonal elements of the pseudo-inverse of the graph Laplacian matrix in nearly linear time complexity, and it is proven that the algorithm has high-probability error bound guarantees. 2. **Experimental Validation**: Extensive experiments were conducted on multiple real-world networks, showing that the approximation algorithm outperforms the direct computation method in terms of both computational efficiency and accuracy. 3. **Exact Solutions**: Exact expressions for the diagonal elements of the pseudo-inverse of the Laplacian matrix for two special types of graphs (Koch networks and uniform recursive trees) were derived and compared with the results of the approximation algorithm. 4. **Application Examples**: The algorithm was used to estimate the Kirchhoff index of three deterministic model networks (Koch networks, uniform recursive trees, and pseudo-fractal scale-free networks), further validating the effectiveness and efficiency of the algorithm. ### Experimental Results The experimental results show that for medium-sized graphs, the running time of the approximation algorithm is significantly lower than that of the exact algorithm. Especially when dealing with large-scale networks, the approximation algorithm can effectively compute all the diagonal elements, while the exact algorithm cannot run due to memory limitations and high time costs. ### Conclusion By proposing an efficient approximation algorithm, this paper solves the high complexity problem of computing the diagonal elements of the pseudo-inverse of the graph Laplacian matrix, providing a new tool for the analysis of large-scale graphs.

Diagonal of Pseudoinverse of Graph Laplacian: Fast Estimation and Exact Results

Scalable Algorithms for Laplacian Pseudo-inverse Computation

Effective and efficient approximations of the generalized inverse of the graph Laplacian matrix with an application to current-flow betweenness centrality

On Fast Computation of Directed Graph Laplacian Pseudo-Inverse

Fast Algorithms for Indices of Nested Split Graphs Approximating Real Complex Networks

The High Relative Accuracy of Computations with Laplacian Matrices

The Laplacian Polynomial and Kirchhoff Index of Graphs Based on R-graphs.

Computing Effective Resistances on Large Graphs Based on Approximate Inverse of Cholesky Factor

Distributed Estimation of Algebraic Connectivity of Directed Networks.

Kirchhoff Index As a Measure of Edge Centrality in Weighted Networks: Nearly Linear Time Algorithms

Random Walk Laplacian and Network Centrality Measures

Efficient Approximation of Kemeny's Constant for Large Graphs.

Sublinear Time Spectral Density Estimation

Isira: Integrated Shift-Invert Residual Arnoldi Method for Graph Laplacian Matrices from Big Data

Generalized Pseudospectral Shattering and Inverse-Free Matrix Pencil Diagonalization

Resistance Eccentricity in Graphs: Distribution, Computation and Optimization

Fast Computation of Kemeny's Constant for Directed Graphs

Learning Graphs from Linear Measurements: Fundamental Trade-offs and Applications

Properties and applications of Laplacian spectra for Koch networks

The Laplacian Polynomial and Kirchhoff Index of Graphs Derived from Regular Graphs

New results about the Inverse Eigenvalue Problem of a Graph