English

Estimating Random-Walk Probabilities in Directed Graphs

Data Structures and Algorithms 2026-05-19 v4

Abstract

We study discounted random walks in directed graphs. In each step, the walk either terminates with a constant probability α\alpha, or proceeds to a random out-neighbor. Our goal is to estimate the probability π(s,t)\pi(s, t) that a discounted random walk starting from ss terminates at tt. This probability is also known as the Personalized PageRank (PPR) score, which measures the relevance of tt to ss, for instance, when ss and tt are web pages on the Internet. We aim to estimate π(s,t)\pi(s, t) 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 nn) 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 π(s,t)\pi(s,t) if it is above a threshold δ\delta in time O(mlog(1/δ))O(m\log(1/\delta)) but π(s,t)\pi(s,t) can be as small as 1/nΘ(n)1/n^{\Theta(n)}. For contrast, we propose algorithms that deterministically estimate arbitrarily small π(s,t)\pi(s,t) in O(mlogn)O(m\log n) time. As a lower-bound example, we improve the lower bound for the single-pair problem from Ω(min{n,1/δ})\Omega(\min\{n,1/\delta\}) to Ω(min{m,1/δ})\Omega(\min\{m,1/\delta\}), which is optimal (up to logarithmic factors) since a simple Monte Carlo estimate takes O(1/δ)O(1/\delta) time.

Keywords

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})$

R2 v1 2026-06-28T23:08:11.343Z