Balls into bins problem

Short description: Balanced or random resource allocation

The balls into bins (or balanced allocations) problem is a classic problem in probability theory that has many applications in computer science. The problem involves m balls and n boxes (or "bins"). Each time, a single ball is placed into one of the bins. After all balls are in the bins, we look at the number of balls in each bin; we call this number the load on the bin. The problem can be modelled using a Multinomial distribution, and may involve asking a question such as: What is the expected number of bins with a ball in them?

Obviously, it is possible to make the load as small as m/n by putting each ball into the least loaded bin. The interesting case is when the bin is selected at random, or at least partially at random. A powerful balls-into-bins paradigm is the "power of two random choices" where each ball chooses two (or more) random bins and is placed in the lesser-loaded bin. This paradigm has found wide practical applications in shared-memory emulations, efficient hashing schemes, randomized load balancing of tasks on servers, and routing of packets within parallel networks and data centers.

Random allocation

When the bin for each ball is selected at random, independent of other choices, the maximum load might be as large as [math]\displaystyle{ m }[/math]. However, it is possible to calculate a tighter bound that holds with high probability. A "high probability" is a probability [math]\displaystyle{ 1-o(1) }[/math], i.e. the probability tends to [math]\displaystyle{ 1 }[/math] when [math]\displaystyle{ n }[/math] grows to infinity.

For the case [math]\displaystyle{ m=n }[/math], with probability [math]\displaystyle{ 1 - o(1) }[/math] the maximum load is:

[math]\displaystyle{ \frac{\log n}{\log \log n}\cdot(1+o(1)) }[/math].

Gonnet gave a tight bound for the expected value of the maximum load, which for [math]\displaystyle{ m=n }[/math] is [math]\displaystyle{ \Gamma^{-1}(n) - \frac{3}{2} + o(1) }[/math], where [math]\displaystyle{ \Gamma^{-1} }[/math] is the inverse of the gamma function, and it is known that [math]\displaystyle{ \Gamma^{-1}(n) = \frac{\log n}{\log\log n} \cdot (1 + o(1)) }[/math].

The maximum load can also be calculated for [math]\displaystyle{ m \ne n }[/math], and for example, for [math]\displaystyle{ m \gt n \log n }[/math] it is [math]\displaystyle{ \frac{m}{n}+\Theta\left(\sqrt{\frac{m \log n}{n}}\right) }[/math], and for [math]\displaystyle{ m \lt n/\log n }[/math] it is [math]\displaystyle{ \Theta\left(\frac{\log n}{\log(n/m)}\right) }[/math], with high probability.

Exact probabilities for small [math]\displaystyle{ m = n }[/math] can be computed as [math]\displaystyle{ a(n)/n^n }[/math] for [math]\displaystyle{ a(n) }[/math] defined in OEIS A208250.

Partially random allocation

Instead of just selecting a random bin for each ball, it is possible to select two or more bins for each ball and then put the ball in the least loaded bin. This is a compromise between a deterministic allocation, in which all bins are checked and the least loaded bin is selected, and a totally random allocation, in which a single bin is selected without checking other bins. This paradigm often called the "power of two random choices" has been studied in a number of settings below.

In the simplest case, if one allocates [math]\displaystyle{ m }[/math] balls into [math]\displaystyle{ n }[/math] bins (with [math]\displaystyle{ m=n }[/math]) sequentially one by one, and for each ball one chooses [math]\displaystyle{ d \ge 2 }[/math] random bins at each step and then allocates the ball into the least loaded of the selected bins (ties broken arbitrarily), then with high probability the maximum load is:

[math]\displaystyle{ \frac{\log \log n}{\log d}+\Theta(1) }[/math]

which is almost exponentially less than with totally random allocation.

This result can be generalized to the case [math]\displaystyle{ m \ge n }[/math] (with [math]\displaystyle{ d \ge 2 }[/math]), when with high probability the maximum load is:

[math]\displaystyle{ \frac{m}{n}+\frac{\log \log n}{\log d}+\Theta(1) }[/math]

which is tight up to an additive constant. (All the bounds hold with probability at least [math]\displaystyle{ 1 - 1/n^c }[/math] for any constant [math]\displaystyle{ c\gt 0 }[/math].) Note that for [math]\displaystyle{ m \gt n \log n }[/math], the random allocation process gives only the maximum load of [math]\displaystyle{ \frac{m}{n}+O\left(\log \log n\right) }[/math] with high probability, so the improvement between these two processes is especially visible for large values of [math]\displaystyle{ m }[/math].

Other key variants of the paradigm are "parallel balls-into-bins" where [math]\displaystyle{ n }[/math] balls choose [math]\displaystyle{ d }[/math] random bins in parallel, "weighted balls-into-bins" where balls have non-unit weights, and "balls-into-bins with deletions" where balls can be added as well as deleted.

Infinite stream of balls

Instead of just putting m balls, it is possible to consider an infinite process in which, at each time step, a single ball is added and a single ball is taken, such that the number of balls remains constant. For m=n, after a sufficiently long time, with high probability the maximum load is similar to the finite version, both with random allocation and with partially random allocation.

Repeated balls-into-bins

In a repeated variant of the process, [math]\displaystyle{ m }[/math] balls are initially distributed in [math]\displaystyle{ n }[/math] bins in an arbitrary way and then, in every subsequent step of a discrete-time process, one ball is chosen from each non-empty bin and re-assigned to one of the [math]\displaystyle{ n }[/math] bins uniformly at random. When [math]\displaystyle{ m=n }[/math], it has been shown that with high probability the process converges to a configuration with maximum load [math]\displaystyle{ \mathcal O(\log(n)) }[/math] after [math]\displaystyle{ \mathcal O(n) }[/math] steps.

Applications

Online Load Balancing: consider a set of n identical computers. There are n users who need computing services. The users are not coordinated - each users comes on his own and selects which computer to use. Each user would of course like to select the least loaded computer, but this requires to check the load on each computer, which might take a long time. Another option is to select a computer at random; this leads, with high probability, to a maximum load of [math]\displaystyle{ \frac{\log n}{\log \log n} \cdot (1 + o(1)) }[/math]. A possible compromise is that the user will check only two computers, and use the lesser loaded of the two. This leads, with high probability, to a much smaller maximum load of [math]\displaystyle{ \log_2 \log n + \Theta(1) }[/math].

Hashing: consider a hash table in which all keys mapped to the same location are stored in a linked list. The efficiency of accessing a key depends on the length of its list. If we use a single hash function which selects locations with uniform probability, with high probability the longest chain has [math]\displaystyle{ O\left(\frac{\log n}{\log \log n}\right) }[/math] keys. A possible improvement is to use two hash functions, and put each new key in the shorter of the two lists. In this case, with high probability the longest chain has only [math]\displaystyle{ O(\log \log n) }[/math] elements.

Fair cake-cutting: consider the problem of creating a partially proportional division of a heterogeneous resource among [math]\displaystyle{ n }[/math] people, such that each person receives a part of the resource which that person values as at least [math]\displaystyle{ 1/an }[/math] of the total, where [math]\displaystyle{ a }[/math] is some sufficiently large constant. The Edmonds–Pruhs protocol is a randomized algorithm whose analysis make use of balls-into-bins arguments.