Las Vegas algorithm

Example

<span id="mwJw"></span><span class="kt" id="mwKA">int</span><span class="w" id="mwKQ"> </span><span class="nf" id="mwKg">getRandomInteger</span><span class="p" id="mwKw">(</span><span class="kt" id="mwLA">int</span><span class="w" id="mwLQ"> </span><span class="n" id="mwLg">n</span><span class="p" id="mwLw">)</span><span class="w" id="mwMA"> </span><span class="p" id="mwMQ">{</span> <span class="w" id="mwMg"> </span><span class="n" id="mwMw">Random</span><span class="w" id="mwNA"> </span><span class="n" id="mwNQ">rand</span><span class="w" id="mwNg"> </span><span class="o" id="mwNw">=</span><span class="w" id="mwOA"> </span><span class="k" id="mwOQ">new</span><span class="w" id="mwOg"> </span><span class="n" id="mwOw">Random</span><span class="p" id="mwPA">();</span> <span class="w" id="mwPQ"> </span><span class="k" id="mwPg">return</span><span class="w" id="mwPw"> </span><span class="n" id="mwQA">rand</span><span class="p" id="mwQQ">.</span><span class="na" id="mwQg">nextInt</span><span class="p" id="mwQw">(</span><span class="n" id="mwRA">n</span><span class="p" id="mwRQ">);</span> <span class="p" id="mwRg">}</span> <span class="c1" id="mwRw">// Las Vegas algorithm, assuming a is array of length n.</span> <span class="c1" id="mwSA">// Goal: return a correct index k such that a[k] == 1.</span> <span class="kt" id="mwSQ">int</span><span class="w" id="mwSg"> </span><span class="nf" id="mwSw">lasVegasAlgorithm</span><span class="p" id="mwTA">(</span><span class="kt" id="mwTQ">int</span><span class="o" id="mwTg">[]</span><span class="w" id="mwTw"> </span><span class="n" id="mwUA">a</span><span class="p" id="mwUQ">)</span><span class="w" id="mwUg"> </span><span class="p" id="mwUw">{</span> <span class="w" id="mwVA"> </span><span class="kt" id="mwVQ">int</span><span class="w" id="mwVg"> </span><span class="n" id="mwVw">n</span><span class="w" id="mwWA"> </span><span class="o" id="mwWQ">=</span><span class="w" id="mwWg"> </span><span class="n" id="mwWw">a</span><span class="p" id="mwXA">.</span><span class="na" id="mwXQ">length</span><span class="p" id="mwXg">;</span> <span class="w" id="mwXw"> </span><span class="k" id="mwYA">while</span><span class="w" id="mwYQ"> </span><span class="p" id="mwYg">(</span><span class="kc" id="mwYw">true</span><span class="p" id="mwZA">)</span><span class="w" id="mwZQ"> </span><span class="p" id="mwZg">{</span> <span class="w" id="mwZw"> </span><span class="kt" id="mwaA">int</span><span class="w" id="mwaQ"> </span><span class="n" id="mwag">k</span><span class="w" id="mwaw"> </span><span class="o" id="mwbA">=</span><span class="w" id="mwbQ"> </span><span class="n" id="mwbg">getRandomInteger</span><span class="p" id="mwbw">(</span><span class="n" id="mwcA">n</span><span class="p" id="mwcQ">);</span> <span class="w" id="mwcg"> </span><span class="k" id="mwcw">if</span><span class="w" id="mwdA"> </span><span class="p" id="mwdQ">(</span><span class="n" id="mwdg">a</span><span class="o" id="mwdw">[</span><span class="n" id="mweA">k</span><span class="o" id="mweQ">]</span><span class="w" id="mweg"> </span><span class="o" id="mwew">==</span><span class="w" id="mwfA"> </span><span class="mi" id="mwfQ">1</span><span class="p" id="mwfg">)</span><span class="w" id="mwfw"> </span><span class="p" id="mwgA">{</span> <span class="w" id="mwgQ"> </span><span class="k" id="mwgg">return</span><span class="w" id="mwgw"> </span><span class="n" id="mwhA">k</span><span class="p" id="mwhQ">;</span> <span class="w" id="mwhg"> </span><span class="p" id="mwhw">}</span> <span class="w" id="mwiA"> </span><span class="p" id="mwiQ">}</span> <span class="p" id="mwig">}</span>

As mentioned above, Las Vegas algorithms always return correct results. The code above illustrates this property. The goal is to find an index in the array a that contains the value 1. A variable k is generated randomly; after k is generated, k is used to index the array a. If this index contains the value 1, then k is returned; otherwise, the algorithm repeats this process until it finds 1. Although this Las Vegas algorithm is guaranteed to find the correct answer, it does not have a fixed runtime; due to the randomization, it is possible for arbitrarily much time to elapse before the algorithm terminates.

Definition

This section provides the conditions that characterize an algorithm's being of Las Vegas type.

An algorithm A is a Las Vegas algorithm for problem class X, if

1. whenever for a given problem instance x∈X it returns a solution s, s is guaranteed to be a valid solution of x
2. on each given instance x, the run-time of A is a random variable RT_A,x

There are three notions of completeness for Las Vegas algorithms:

• complete Las Vegas algorithms can be guaranteed to solve each solvable problem within run-time tmax, where tmax is an instance-dependent constant.

Let P(RT_A,x ≤ t) denote the probability that A finds a solution for a soluble instance x in time within t, then A is complete exactly if for each x there exists

some tmax such that P(RT_A,x ≤ t_max) = 1.

• approximately complete Las Vegas algorithms solve each problem with a probability converging to 1 as the run-time approaches infinity. Thus, A is approximately complete, if for each instance x, lim_t→∞ P(RT_A,x ≤ t) = 1.
• essentially incomplete Las Vegas algorithms are Las Vegas algorithms that are not approximately complete.

Approximate completeness is primarily of theoretical interest, as the time limits for finding solutions are usually too large to be of practical use.

Application scenarios

Las Vegas algorithms have different criteria for the evaluation based on the problem setting. These criteria are divided into three categories with different time limits since Las Vegas algorithms do not have set time complexity. Here are some possible application scenarios:

• Type 1: There are no time limits, which means the algorithm runs until it finds the solution.
• Type 2: There is a time limit t_max for finding the outcome.
• Type 3: The utility of a solution is determined by the time required to find the solution.

(Type 1 and Type 2 are special cases of Type 3.)

For Type 1 where there is no time limit, the average run-time can represent the run-time behavior. This is not the same case for Type 2.

Here, P(RT ≤ t_max), which is the probability of finding a solution within time, describes its run-time behavior.

In case of Type 3, its run-time behavior can only be represented by the run-time distribution function rtd: R → [0,1] defined as rtd(t) = P(RT ≤ t) or its approximation.

The run-time distribution (RTD) is the distinctive way to describe the run-time behavior of a Las Vegas algorithm.

With this data, we can easily get other criteria such as the mean run-time, standard deviation, median, percentiles, or success probabilities P(RT ≤ t) for arbitrary time-limits t.

Applications

Analogy

Las Vegas algorithms arise frequently in search problems. For example, one looking for some information online might search related websites for the desired information. The time complexity thus ranges from getting "lucky" and finding the content immediately, to being "unlucky" and spending large amounts of time. Once the right website is found, then there is no possibility of error.

Randomized quicksort

<span id="mw3A"></span><span class="k" id="mw3Q">def</span><span class="w" id="mw3g"> </span><span class="nf" id="mw3w">randomized_quicksort</span><span class="p" id="mw4A">(</span><span class="n" id="mw4Q">a</span><span class="p" id="mw4g">:</span> <span class="nb" id="mw4w">list</span><span class="p" id="mw5A">[</span><span class="nb" id="mw5Q">int</span><span class="p" id="mw5g">])</span> <span class="o" id="mw5w">-></span> <span class="kc" id="mw6A">None</span><span class="p" id="mw6Q">:</span> <span class="k" id="mw6g">if</span> <span class="nb" id="mw6w">len</span><span class="p" id="mw7A">(</span><span class="n" id="mw7Q">a</span><span class="p" id="mw7g">)</span> <span class="o" id="mw7w">==</span> <span class="mi" id="mw8A">1</span><span class="p" id="mw8Q">:</span> <span class="k" id="mw8g">return</span> <span class="n" id="mw8w">A</span> <span class="c1" id="mw9A"># A is sorted.</span> <span class="k" id="mw9Q">else</span><span class="p" id="mw9g">:</span> <span class="n" id="mw9w">i</span><span class="p" id="mw-A">:</span> <span class="nb" id="mw-Q">int</span> <span class="o" id="mw-g">=</span> <span class="n" id="mw-w">random</span><span class="o" id="mw_A">.</span><span class="n" id="mw_Q">randrange</span><span class="p" id="mw_g">(</span><span class="mi" id="mw_w">1</span><span class="p" id="mwAQA">,</span> <span class="nb" id="mwAQE">len</span><span class="p" id="mwAQI">(</span><span class="n" id="mwAQM">a</span><span class="p" id="mwAQQ">))</span> <span class="c1" id="mwAQU"># Will take a random number in the range [1, len(a))</span> <span class="n" id="mwAQY">x</span><span class="p" id="mwAQc">:</span> <span class="nb" id="mwAQg">int</span> <span class="o" id="mwAQk">=</span> <span class="n" id="mwAQo">a</span><span class="p" id="mwAQs">[</span><span class="n" id="mwAQw">i</span><span class="p" id="mwAQ0">]</span> <span class="c1" id="mwAQ4"># The pivot element</span> <span class="c1" id="mwAQ8"># Partition a into elements < x, x, and > x, as shown in the figure above.</span> <span class="c1" id="mwARA"># Execute Quicksort on a[1 : i- 1] and A[i + 1 : n].</span> <span class="c1" id="mwARE"># Combine the responses in order to obtain a sorted array.</span>

A simple example is randomized quicksort, where the pivot is chosen randomly, and divides the elements into three partitions: elements less than pivot, elements equal to pivot, and elements greater than pivot. QuickSort always generates the solution, which in this case the sorted array. Unfortunately, the time complexity is not immediately obvious, as the runtime depends on which element is picked as a pivot.

• The worst case is Θ(n²) when the pivot is the smallest or the largest element.

${\displaystyle T(n)=T(0)+T(n-1)+\mathrm{\Theta }(n)}$

${\displaystyle T(n)=\mathrm{\Theta }(1)+T(n-1)+\mathrm{\Theta }(n)}$

${\displaystyle T(n)=T(n-1)+\mathrm{\Theta }(n)}$

${\displaystyle T(n)=\mathrm{\Theta }(n^2)}$

• However, through randomization, where the pivot is randomly picked and is exactly a middle value each time, the QuickSort can be done in Θ(nlogn).

${\displaystyle T(n)\le 2\ast T(n/2)+\mathrm{\Theta }(n)}$

${\displaystyle T(n)=\mathrm{\Theta }(n\log (n))}$

The runtime of quicksort depends heavily on how well the pivot is selected. If a value of pivot is either too big or small, the partition will be unbalanced, resulting in a poor runtime efficiency. However, if the value of pivot is near the middle of the array, then the split will be reasonably well balanced, yielding a faster runtime. Since the pivot is randomly picked, the running time will be good most of the time and bad occasionally.

In the case of average case, it is hard to determine since the analysis does not depend on the input distribution but on the random choices that the algorithm makes. The average of quicksort is computed over all possible random choices that the algorithm might make when choosing the pivot.

Although the worst-case runtime is Θ(n²), the average-case runtime is Θ(nlogn). It turns out that the worst-case does not happen often. For large values of n, the runtime is Θ(nlogn) with a high probability.

Note that the probability that the pivot is the middle value element each time is one out of n numbers, which is very rare. However, it is still the same runtime when the split is 10%-90% instead of a 50%–50% because the depth of the recursion tree will still be O(logn) with O(n) times taken each level of recursion.

Randomized greedy algorithm for the Eight queens problem

The eight queens problem is usually solved with a backtracking algorithm. However, a Las Vegas algorithm can be applied; in fact, it is more efficient than backtracking.

Place 8 queens on a chessboard so that no one attacks another. Remember that a queen attacks other pieces on the same row, column and diagonals.

Assume that k rows, 0 ≤ k ≤ 8, are successfully occupied by queens.

If k = 8, then stop with success. Otherwise, proceed to occupy row k + 1.

Calculate all positions on this row not attacked by existing queens. If there are none, then fail. Otherwise, pick one at random, increment k and repeat.

Note that the algorithm simply fails if a queen cannot be placed. But the process can be repeated and every time will generate different arrangement.

Complexity class

The complexity class of decision problems that have Las Vegas algorithms with expected polynomial runtime is ZPP.

It turns out that

${\displaystyle \mathsf{\text{ZPP}}=\mathsf{\text{RP}}\cap \mathsf{\text{co-RP}}}$

which is intimately connected with the way Las Vegas algorithms are sometimes constructed. Namely the class RP consists of all decision problems for which a randomized polynomial-time algorithm exists that always answers correctly when the correct answer is "no", but is allowed to be wrong with a certain probability bounded away from one when the answer is "yes". When such an algorithm exists for both a problem and its complement (with the answers "yes" and "no" swapped), the two algorithms can be run simultaneously and repeatedly: run each for a constant number of steps, taking turns, until one of them returns a definitive answer. This is the standard way to construct a Las Vegas algorithm that runs in expected polynomial time. Note that in general there is no worst case upper bound on the run time of a Las Vegas algorithm.

Optimal Las Vegas algorithm

In order to make a Las Vegas algorithm optimal, the expected run time should be minimized. This can be done by:

1. The Las Vegas algorithm A(x) runs repeatedly for some number t₁ steps. If A(x) stops during the run time then A(x) is done; otherwise, repeat the process from the beginning for another t₂ steps, and so on.
2. Designing a strategy that is optimal among all strategies for A(x), given the full information about the distribution of T_A(x).

The existence of the optimal strategy might be a fascinating theoretical observation. However, it is not practical in real life because it is not easy to find the information of distribution of T_A(x). Furthermore, there is no point of running the experiment repeatedly to obtain the information about the distribution since most of the time, the answer is needed only once for any x.

Relation to Monte Carlo algorithms

Las Vegas algorithms can be contrasted with Monte Carlo algorithms, in which the resources used are bounded but the answer may be incorrect with a certain (typically small) probability. A Las Vegas algorithm can be converted into a Monte Carlo algorithm by running it for set time and generating a random answer when it fails to terminate. By an application of Markov's inequality, we can set the bound on the probability that the Las Vegas algorithm would go over the fixed limit.

Here is a table comparing Las Vegas and Monte Carlo algorithms:

If a deterministic way to test for correctness is available, then it is possible to turn a Monte Carlo algorithm into a Las Vegas algorithm. However, it is hard to convert a Monte Carlo algorithm to a Las Vegas algorithm without a way to test the algorithm. On the other hand, changing a Las Vegas algorithm to a Monte Carlo algorithm is easy. This can be done by running a Las Vegas algorithm for a specific period of time given by confidence parameter. If the algorithm finds the solution within the time, then it is success and if not then output can simply be "sorry".

This is an example of Las Vegas and Monte Carlo algorithms for comparison:

Assume that there is an array with the length of even n. Half of the entries in the array are 0's and the remaining half are 1's. The goal here is to find an index that contains a 1.

<span id="mwAZk"></span><span class="c1" id="mwAZo">// Las Vegas algorithm</span> <span class="kt" id="mwAZs">int</span><span class="w" id="mwAZw"> </span><span class="nf" id="mwAZ0">lasVegasAlgorithm</span><span class="p" id="mwAZ4">(</span><span class="kt" id="mwAZ8">int</span><span class="o" id="mwAaA">[]</span><span class="w" id="mwAaE"> </span><span class="n" id="mwAaI">a</span><span class="p" id="mwAaM">)</span><span class="w" id="mwAaQ"> </span><span class="p" id="mwAaU">{</span> <span class="w" id="mwAaY"> </span><span class="kt" id="mwAac">int</span><span class="w" id="mwAag"> </span><span class="n" id="mwAak">n</span><span class="w" id="mwAao"> </span><span class="o" id="mwAas">=</span><span class="w" id="mwAaw"> </span><span class="n" id="mwAa0">a</span><span class="p" id="mwAa4">.</span><span class="na" id="mwAa8">length</span><span class="p" id="mwAbA">;</span> <span class="w" id="mwAbE"> </span><span class="k" id="mwAbI">while</span><span class="w" id="mwAbM"> </span><span class="p" id="mwAbQ">(</span><span class="kc" id="mwAbU">true</span><span class="p" id="mwAbY">)</span><span class="w" id="mwAbc"> </span><span class="p" id="mwAbg">{</span> <span class="w" id="mwAbk"> </span><span class="kt" id="mwAbo">int</span><span class="w" id="mwAbs"> </span><span class="n" id="mwAbw">k</span><span class="w" id="mwAb0"> </span><span class="o" id="mwAb4">=</span><span class="w" id="mwAb8"> </span><span class="n" id="mwAcA">getRandomInteger</span><span class="p" id="mwAcE">(</span><span class="n" id="mwAcI">n</span><span class="p" id="mwAcM">);</span> <span class="w" id="mwAcQ"> </span><span class="k" id="mwAcU">if</span><span class="w" id="mwAcY"> </span><span class="p" id="mwAcc">(</span><span class="n" id="mwAcg">a</span><span class="o" id="mwAck">[</span><span class="n" id="mwAco">k</span><span class="o" id="mwAcs">]</span><span class="w" id="mwAcw"> </span><span class="o" id="mwAc0">==</span><span class="w" id="mwAc4"> </span><span class="mi" id="mwAc8">1</span><span class="p" id="mwAdA">)</span><span class="w" id="mwAdE"> </span><span class="p" id="mwAdI">{</span> <span class="w" id="mwAdM"> </span><span class="k" id="mwAdQ">return</span><span class="w" id="mwAdU"> </span><span class="n" id="mwAdY">k</span><span class="p" id="mwAdc">;</span> <span class="w" id="mwAdg"> </span><span class="p" id="mwAdk">}</span> <span class="w" id="mwAdo"> </span><span class="p" id="mwAds">}</span> <span class="p" id="mwAdw">}</span> <span class="c1" id="mwAd0">// Monte Carlo algorithm</span> <span class="kt" id="mwAd4">int</span><span class="w" id="mwAd8"> </span><span class="nf" id="mwAeA">monteCarloAlgorithm</span><span class="p" id="mwAeE">(</span><span class="kt" id="mwAeI">int</span><span class="o" id="mwAeM">[]</span><span class="w" id="mwAeQ"> </span><span class="n" id="mwAeU">a</span><span class="p" id="mwAeY">)</span><span class="w" id="mwAec"> </span><span class="p" id="mwAeg">{</span> <span class="w" id="mwAek"> </span><span class="kt" id="mwAeo">int</span><span class="w" id="mwAes"> </span><span class="n" id="mwAew">n</span><span class="w" id="mwAe0"> </span><span class="o" id="mwAe4">=</span><span class="w" id="mwAe8"> </span><span class="n" id="mwAfA">a</span><span class="p" id="mwAfE">.</span><span class="na" id="mwAfI">length</span><span class="p" id="mwAfM">;</span> <span class="w" id="mwAfQ"> </span><span class="k" id="mwAfU">for</span><span class="w" id="mwAfY"> </span><span class="p" id="mwAfc">(</span><span class="kt" id="mwAfg">int</span><span class="w" id="mwAfk"> </span><span class="n" id="mwAfo">i</span><span class="w" id="mwAfs"> </span><span class="o" id="mwAfw">=</span><span class="w" id="mwAf0"> </span><span class="mi" id="mwAf4">0</span><span class="p" id="mwAf8">;</span><span class="w" id="mwAgA"> </span><span class="n" id="mwAgE">i</span><span class="w" id="mwAgI"> </span><span class="o" id="mwAgM"><</span><span class="w" id="mwAgQ"> </span><span class="mi" id="mwAgU">300</span><span class="p" id="mwAgY">;</span><span class="w" id="mwAgc"> </span><span class="o" id="mwAgg">++</span><span class="n" id="mwAgk">i</span><span class="p" id="mwAgo">)</span><span class="w" id="mwAgs"> </span><span class="p" id="mwAgw">{</span> <span class="w" id="mwAg0"> </span><span class="kt" id="mwAg4">int</span><span class="w" id="mwAg8"> </span><span class="n" id="mwAhA">k</span><span class="w" id="mwAhE"> </span><span class="o" id="mwAhI">=</span><span class="w" id="mwAhM"> </span><span class="n" id="mwAhQ">getRandomInteger</span><span class="p" id="mwAhU">(</span><span class="n" id="mwAhY">n</span><span class="p" id="mwAhc">);</span> <span class="w" id="mwAhg"> </span><span class="k" id="mwAhk">if</span><span class="w" id="mwAho"> </span><span class="p" id="mwAhs">(</span><span class="n" id="mwAhw">a</span><span class="o" id="mwAh0">[</span><span class="n" id="mwAh4">k</span><span class="o" id="mwAh8">]</span><span class="w" id="mwAiA"> </span><span class="o" id="mwAiE">==</span><span class="w" id="mwAiI"> </span><span class="mi" id="mwAiM">1</span><span class="p" id="mwAiQ">)</span><span class="w" id="mwAiU"> </span><span class="p" id="mwAiY">{</span> <span class="w" id="mwAic"> </span><span class="k" id="mwAig">return</span><span class="w" id="mwAik"> </span><span class="n" id="mwAio">k</span><span class="p" id="mwAis">;</span> <span class="w" id="mwAiw"> </span><span class="p" id="mwAi0">}</span> <span class="w" id="mwAi4"> </span><span class="p" id="mwAi8">}</span> <span class="w" id="mwAjA"> </span><span class="k" id="mwAjE">return</span><span class="w" id="mwAjI"> </span><span class="o" id="mwAjM">-</span><span class="mi" id="mwAjQ">1</span><span class="p" id="mwAjU">;</span><span class="w" id="mwAjY"> </span><span class="c1" id="mwAjc">// denotes failure</span> <span class="p" id="mwAjg">}</span>

Since Las Vegas does not end until it finds 1 in the array, it does not gamble with the correctness but run-time. On the other hand, Monte Carlo runs 300 times, which means it is impossible to know that Monte Carlo will find "1" in the array within 300 times of loops until it actually executes the code. It might find the solution or not. Therefore, unlike Las Vegas, Monte Carlo does not gamble with run-time but correctness.