Abstract:This paper addresses a fundamental problem in random variate generation: given access to a random source that emits a stream of independent fair bits, what is the most accurate and entropy-efficient algorithm for sampling from a discrete probability distribution $(p_1, \dots, p_n)$, where the probabilities of the output distribution $(\hat{p}_1, \dots, \hat{p}_n)$ of the sampling algorithm must be specified using at most $k$ bits of precision? We present a theoretical framework for formulating this problem and provide new techniques for finding sampling algorithms that are optimal both statistically (in the sense of sampling accuracy) and information-theoretically (in the sense of entropy consumption). We leverage these results to build a system that, for a broad family of measures of statistical accuracy, delivers a sampling algorithm whose expected entropy usage is minimal among those that induce the same distribution (i.e., is "entropy-optimal") and whose output distribution $(\hat{p}_1, \dots, \hat{p}_n)$ is a closest approximation to the target distribution $(p_1, \dots, p_n)$ among all entropy-optimal sampling algorithms that operate within the specified $k$-bit precision. This optimal approximate sampler is also a closer approximation than any (possibly entropy-suboptimal) sampler that consumes a bounded amount of entropy with the specified precision, a class which includes floating-point implementations of inversion sampling and related methods found in many software libraries. We evaluate the accuracy, entropy consumption, precision requirements, and wall-clock runtime of our optimal approximate sampling algorithms on a broad set of distributions, demonstrating the ways that they are superior to existing approximate samplers and establishing that they often consume significantly fewer resources than are needed by exact samplers.

Approximating Perpetuities

Fast Perfect Simulation of Vervaat Perpetutities

Perfect simulation of Vervaat perpetuities

Approximating the Permanent with Belief Propagation

Approximation Algorithms for Perishable Inventory Systems

Approximating partial differential equations without boundary conditions

Approximation methods for piecewise deterministic Markov processes and their costs

A Tight Analysis of Bethe Approximation for Permanent

Approximating the Bethe partition function

Pathwise and distributional approximations of semi-Markov processes

A $(2+\varepsilon)$-approximation algorithm for preemptive weighted flow time on a single machine

Approximations in Probabilistic Programs

Optimal Approximate Sampling from Discrete Probability Distributions

Greedy Algorithms for Optimal Distribution Approximation

Approximate optimality and the risk/reward tradeoff given repeated gambles

Some Functionals and Approximation Operators Associated with a Family of Discrete Probability Distributions

Density in Approximation Theory

On perpetuities with gamma-like tails

Optimal distributions for randomized unbiased estimators with an infinite horizon and an adaptive algorithm

Approximation of Hysteresis Functional

On discrete approximations of stable distributions