Coin-flipping, ball-dropping, and grass-hopping for generating random graphs from matrices of edge probabilities

TL;DR

This paper introduces efficient graph sampling techniques, grass-hopping and ball-dropping, for generating large sparse random graphs from probability matrices, improving over naive methods by reducing unnecessary computations.

Contribution

It presents novel algorithms and computational primitives for sampling from complex probability matrices, including a new connection between Kronecker graphs and Morton codes.

Findings

Grass-hopping reduces sampling complexity for sparse graphs.

New methods handle complex probability matrices efficiently.

Connections to applied math problems like coupon collector are established.

Abstract

Common models for random graphs, such as Erd\H{o}s-R\'{e}nyi and Kronecker graphs, correspond to generating random adjacency matrices where each entry is non-zero based on a large matrix of probabilities. Generating an instance of a random graph based on these models is easy, although inefficient, by flipping biased coins (i.e. sampling binomial random variables) for each possible edge. This process is inefficient because most large graph models correspond to sparse graphs where the vast majority of coin flips will result in no edges. We describe some not-entirely-well-known, but not-entirely-unknown, techniques that will enable us to sample a graph by finding only the coin flips that will produce edges. Our analogies for these procedures are ball-dropping, which is easier to implement, but may need extra work due to duplicate edges, and grass-hopping, which results in no duplicated work or extra edges. Grass-hopping does this using geometric random variables. In order to use this idea on complex probability matrices such as those in Kronecker graphs, we decompose the problem into three steps, each of which are independently useful computational primitives: (i) enumerating non-decreasing sequences, (ii) unranking multiset permutations, and (iii) decoding and encoding z-curve and Morton codes and permutations. The third step is the result of a new connection between repeated Kronecker product operations and Morton codes. Throughout, we draw connections to ideas underlying applied math and computer science including coupon collector problems.