Estimating Random-Walk Probabilities in Directed Graphs
Abstract
We study discounted random walks in directed graphs. In each step, the walk either terminates with a constant probability , or proceeds to a random out-neighbor. Our goal is to estimate the probability that a discounted random walk starting from terminates at . This probability is also known as the Personalized PageRank (PPR) score, which measures the relevance of to , for instance, when and are web pages on the Internet. We aim to estimate within a constant relative error with constant probability. A variety of algorithms have been developed for several problem variants, such as single-pair, single-source, single-target, and single-node estimation, under both worst-case and average-case settings, and for different combinations of allowed graph queries. However, in many important cases, there remain polynomial gaps between known upper and lower bounds. In this paper, we establish tight upper and lower bounds (up to logarithmic factors of ) for all problem variants and query combinations, closing all existing gaps in both the worst-case and average-case settings. Below we give some examples for the worst-case settings. As an upper-bound example, the classic power method estimates if it is above a threshold in time but can be as small as . For contrast, we propose algorithms that deterministically estimate arbitrarily small in time. As a lower-bound example, we improve the lower bound for the single-pair problem from to , which is optimal (up to logarithmic factors) since a simple Monte Carlo estimate takes time.
Cite
@article{arxiv.2504.16481,
title = {Estimating Random-Walk Probabilities in Directed Graphs},
author = {Christian Bertram and Mads Vestergaard Jensen and Mikkel Thorup and Hanzhi Wang and Shuyi Yan},
journal= {arXiv preprint arXiv:2504.16481},
year = {2026}
}
Comments
v4: Providing an $O(m\log n)$ upper bound for estimating $\pi(s,t)$ regardless of how small $\pi(s,t)$ is (i.e., addressing the case where the relative error threshold $\delta = 0$). The previous upper bound was $O(m\log{1/\delta})$