Nearly optimal Bernoulli factories for linear functions

Abstract

Suppose that X1,X2,… are independent identically distributed Bernoulli random variables with mean p. A Bernoulli factory for a function f takes as input X1,X2,… and outputs a random variable that is Bernoulli with mean f(p). A fast algorithm is a function that only depends on the values of X1,…,XT, where T is a stopping time with small mean. When f(p) is a real analytic function the problem reduces to being able to draw from linear functions Cp for a constant C > 1. Also it is necessary that Cp ≤ 1 - ε for known ε > 0. Previous methods for this problem required extensive modification of the algorithm for every value of C and ε. These methods did not have explicit bounds on E[T] as a function of C and ε. This paper presents the first Bernoulli factory for f(p) = Cp with bounds on E[T] as a function of the input parameters. In fact, p ∈ [0,(1-ε)/C] E[T] ≤ 9.5Cε-1. In addition, this method is very simple to implement. Furthermore, a lower bound on the average running time of any Cp Bernoulli factory is shown. For ε ≤ 1/2, p ∈ [0,(1 - ε)/C] E[T] ≥ 0.004 C ε-1, so the new method is optimal up to a constant in the running time.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…