Algorithmically random sequence

Richard von Mises Richard von Mises formalized the notion of a test for randomness in order to define a random sequence as one that passed all tests for randomness. He defined a "collective" (kollektiv) to be an infinite binary string x_{1:\infty} defined such that • There exists a limit \lim_n \frac 1n \sum_{i=1}^n x_i = p \in (0, 1) . • For any "admissible" rule, such that it picks out an infinite subsequence (x_{m_i})_i from the string, we still have \lim_n \frac 1n \sum_{i=1}^n x_{m_i} = p . He called this principle "impossibility of a gambling system". To pick out a subsequence, first pick a binary function \phi, such that given any binary string x_{1:k}, it outputs either 0 or 1. If it outputs 1, then we add x_{k+1} to the subsequence, else we continue. In this definition, some admissible rules might abstain forever on some sequences, and thus fail to pick out an infinite subsequence. We only consider those that do pick an infinite subsequence. Stated in another way, each infinite binary string is a coin-flip game, and an admissible rule is a way for a gambler to decide when to place bets. A collective is a coin-flip game where there is no way for one gambler to do better than another over the long run. That is, there is no gambling system that works for the game. The definition generalizes from binary alphabet to countable alphabet: • The frequency of each letter converges to a limit greater than zero. • For any "admissible" rule, such that it picks out an infinite subsequence (x_{m_i})_i from the string, the frequency of each letter in the subsequence still converges to the same limit. Usually the admissible rules are defined to be rules computable by a Turing machine, and we require p=1/2. With this, we have the Mises–Wald–Church random sequences. This is not a restriction, since given a sequence with p=1/2, we can construct random sequences with any other computable p \in (0, 1). (Here, "Church" refers to Alonzo Church, whose 1940 paper proposed using Turing-computable rules.) {{block indent| Theorem (Abraham Wald, 1936, 1937) If there are only countably many admissible rules, then almost any sequence is a collective. Proof sketch: Use measure-theoretic probability. Fix one admissible rule. Sample a random sequence from Bernoulli space. With probability 1 (use martingales), the subsequence picked by the admissible rule still has \lim_n \frac 1n \sum_{i=1}^n x_{m_i} = p . Now add all the countably many rules. With probability 1, each subsequence picked by each rule still has \lim_n \frac 1n \sum_{i=1}^n x_{m_i} = p . }} However, this definition was found not to be strong enough. Intuitively, the long-time average of a random sequence should oscillate on both sides of p, like how a random walk should cross the origin infinitely many times. However, Jean Ville showed that, even with countably many rules, there exists a binary sequence that tends towards p fraction of ones, but, for every finite prefix, the fraction of ones is less than p. {{block indent| '''Ville's Construction''' (Jean Ville, 1939) There exists a collective with countably many admissible rules such that, for all n, \frac 1n \sum_{k=1}^n x_k \leq p. }} Per Martin-Löf The Ville construction suggests that the Mises–Wald–Church sense of randomness is not good enough, because some random sequences do not satisfy some laws of randomness. For example, the Ville construction does not satisfy one of the laws of the iterated logarithm: \limsup_{n \to \infty} \frac{-\sum_{k=1}^n (x_k - 1/2)}{\sqrt{2n \log\log n}} \neq 1 Naively, one can fix this by requiring a sequence to satisfy all possible laws of randomness, where a "law of randomness" is a property that is satisfied by all sequences with probability 1. However, for each infinite sequence y_{1:\infty} \in 2^\N, we have a law of randomness that x_{1:\infty} \neq y_{1:\infty} , leading to the conclusion that there are no random sequences. (Per Martin-Löf, 1966) defined "Martin-Löf randomness" by only allowing laws of randomness that are Turing-computable. In other words, a sequence is random iff it passes all Turing-computable tests of randomness. The thesis that the definition of Martin-Löf randomness "correctly" captures the intuitive notion of randomness has been called the Martin-Löf–Chaitin Thesis; it is somewhat similar to the Church–Turing thesis. == Three equivalent definitions ==

Martin-Löf's original definition of a random sequence was in terms of constructive null covers; he defined a sequence to be random if it is not contained in any such cover. Gregory Chaitin, Leonid Levin and Claus Peter Schnorr proved a characterization in terms of algorithmic complexity: a sequence is random if there is a uniform bound on the compressibility of its initial segments. Schnorr gave a third equivalent definition in terms of martingales. Li and Vitanyi's book An Introduction to Kolmogorov Complexity and Its Applications is the standard introduction to these ideas. • Algorithmic complexity (Chaitin 1969, Schnorr 1973, Levin 1973): Algorithmic complexity (also known as (prefix-free) Kolmogorov complexity or program-size complexity) can be thought of as a lower bound on the algorithmic compressibility of a finite sequence (of characters or binary digits). It assigns to each such sequence w a natural number K(w) that, intuitively, measures the minimum length of a computer program (written in some fixed programming language) that takes no input and will output w when run. Given a natural number c and a sequence w, we say that w is '''c-incompressible''' if K(w) \geq |w| - c . :An infinite sequence S is Martin-Löf random if and only if there is a constant c such that all of S finite prefixes are c-incompressible. More succinctly, K(w) \geq |w| - O(1) . • Constructive null covers (Martin-Löf 1966): This is Martin-Löf's original definition. For a finite binary string w we let Cw denote the cylinder generated by w. This is the set of all infinite sequences beginning with w, which is a basic open set in Cantor space. The product measure μ(Cw) of the cylinder generated by w is defined to be 2−|w|. Every open subset of Cantor space is the union of a countable sequence of disjoint basic open sets, and the measure of an open set is the sum of the measures of any such sequence. An effective open set is an open set that is the union of the sequence of basic open sets determined by a recursively enumerable sequence of binary strings. A constructive null cover or effective measure 0 set is a recursively enumerable sequence U_i of effective open sets such that U_{i+1} \subseteq U_i and \mu (U_i) \leq 2^{-i} for each natural number i. Every effective null cover determines a G_\delta set of measure 0, namely the intersection of the sets U_i. :A sequence is defined to be Martin-Löf random if it is not contained in any G_\delta set determined by a constructive null cover. • Constructive martingales (Schnorr 1971): A martingale is a function d:\{0,1\}^*\to[0,\infty) such that, for all finite strings w, d(w) = (d(w^\smallfrown 0) + d(w^\smallfrown 1))/2, where a^\smallfrown b is the concatenation of the strings a and b. This is called the "fairness condition": if a martingale is viewed as a betting strategy, then the above condition requires that the bettor plays against fair odds. A martingale d is said to succeed on a sequence S if \limsup_{n\to\infty} d(S \upharpoonright n) = \infty, where S \upharpoonright n is the first n bits of S. A martingale d is constructive (also known as weakly computable, lower semi-computable) if there exists a computable function \widehat{d}:\{0,1\}^*\times\N\to{\mathbb{Q}} such that, for all finite binary strings w • \widehat{d}(w,t) \leq \widehat{d}(w,t+1) for all positive integers t, • \lim_{t\to\infty} \widehat{d}(w,t) = d(w). :A sequence is Martin-Löf random if and only if no constructive martingale succeeds on it. == Interpretations of the definitions ==

The Kolmogorov complexity characterization conveys the intuition that a random sequence is incompressible: no prefix can be produced by a program much shorter than the prefix. The null cover characterization conveys the intuition that a random real number should not have any property that is "uncommon". Each measure 0 set can be thought of as an uncommon property. It is not possible for a sequence to lie in no measure 0 sets, because each one-point set has measure 0. Martin-Löf's idea was to limit the definition to measure 0 sets that are effectively describable; the definition of an effective null cover determines a countable collection of effectively describable measure 0 sets and defines a sequence to be random if it does not lie in any of these particular measure 0 sets. Since the union of a countable collection of measure 0 sets has measure 0, this definition immediately leads to the theorem that there is a measure 1 set of random sequences. Note that if we identify the Cantor space of binary sequences with the interval [0,1] of real numbers, the measure on Cantor space agrees with Lebesgue measure. An effective measure 0 set can be interpreted as a Turing machine that is able to tell, given an infinite binary string, whether the string looks random at levels of statistical significance. The set is the intersection of shrinking sets U_1 \supset U_2 \supset U_3 \supset \cdots, and since each set U_n is specified by an enumerable sequence of prefixes, given any infinite binary string, if it is in U_n, then the Turing machine can decide in finite time that the string does fall inside U_n. Therefore, it can "reject the hypothesis that the string is random at significance level 2^{-n}". If the Turing machine can reject the hypothesis at all significance levels, then the string is not random. A random string is one that, for each Turing-computable test of randomness, manages to remain forever un-rejected at some significance level. The martingale characterization conveys the intuition that no effective procedure should be able to make money betting against a random sequence. A martingale d is a betting strategy. d reads a finite string w and bets money on the next bit. It bets some fraction of its money that the next bit will be 0, and then remainder of its money that the next bit will be 1. d doubles the money it placed on the bit that actually occurred, and it loses the rest. d(w) is the amount of money it has after seeing the string w. Since the bet placed after seeing the string w can be calculated from the values d(w), d(w0), and d(w1), calculating the amount of money it has is equivalent to calculating the bet. The martingale characterization says that no betting strategy implementable by any computer (even in the weak sense of constructive strategies, which are not necessarily computable) can make money betting on a random sequence. ==Properties and examples of Martin-Löf random sequences==

Universality There is a universal constructive martingale d. This martingale is universal in the sense that, given any constructive martingale d, if d succeeds on a sequence, then d succeeds on that sequence as well. Thus, d succeeds on every sequence in RANDc (but, since d is constructive, it succeeds on no sequence in RAND). (Schnorr 1971) There is a constructive null cover of RANDc. This means that all effective tests for randomness (that is, constructive null covers) are, in a sense, subsumed by this universal test for randomness, since any sequence that passes this single test for randomness will pass all tests for randomness. (Martin-Löf 1966) Intuitively, this universal test for randomness says "If the sequence has increasingly long prefixes that can be increasingly well-compressed on this universal Turing machine", then it is not random." -- see next section. Construction sketch: Enumerate the effective null covers as ((U_{m, n})_n)_m. The enumeration is also effective (enumerated by a modified universal Turing machine). Now we have a universal effective null cover by diagonalization: (\cup_{n}U_{n, n+k+1})_k. Passing randomness tests If a sequence fails an algorithmic randomness test, then it is algorithmically compressible. Conversely, if it is algorithmically compressible, then it fails an algorithmic randomness test. Construction sketch: Suppose the sequence fails a randomness test, then it can be compressed by lexicographically enumerating all sequences that fails the test, then code for the location of the sequence in the list of all such sequences. This is called "enumerative source encoding". Conversely, if the sequence is compressible, then by the pigeonhole principle, only a vanishingly small fraction of sequences are like that, so we can define a new test for randomness by "has a compression by this universal Turing machine". Incidentally, this is the universal test for randomness. For example, consider a binary sequence sampled IID from the Bernoulli distribution. After taking a large number N of samples, we should have about M \approx pN ones. We can code for this sequence as "Generate all binary sequences with length N, and M ones. Of those, the i-th sequence in lexicographic order.". By Stirling approximation, \log_2 \binom{N}{pN} \approx N H(p) where H is the binary entropy function. Thus, the number of bits in this description is 2(1 + \epsilon) \log_2 N + (1 + \epsilon) N H(p) + O(1) The first term is for prefix-coding the numbers N and M. The second term is for prefix-coding the number i. (Use Elias omega coding.) The third term is for prefix-coding the rest of the description. When N is large, this description has just \sim H(p) N bits, and so it is compressible, with compression ratio \sim H(p) . In particular, the compression ratio is exactly one (incompressible) only when p = 1/2. (Example 14.2.8 ) Impossibility of a gambling system Consider a casino offering fair odds at a roulette table. The roulette table generates a sequence of random numbers. If this sequence is algorithmically random, then there is no lower semi-computable strategy to win, which in turn implies that there is no computable strategy to win. That is, for any gambling algorithm, the long-term log-payoff is zero (neither positive nor negative). Conversely, if this sequence is not algorithmically random, then there is a lower semi-computable strategy to win. Examples • Chaitin's halting probability Ω is an example of a random sequence. • Every random sequence is not computable. • Every random sequence is normal, satisfies the law of large numbers, and satisfies all Turing-computable properties satisfied by an IID stream of uniformly random numbers. (Theorem 14.5.2 ) Relation to the arithmetic hierarchy • RANDc (the complement of RAND) is a measure 0 subset of the set of all infinite sequences. This is implied by the fact that each constructive null cover covers a measure 0 set, there are only countably many constructive null covers, and a countable union of measure 0 sets has measure 0. This implies that RAND is a measure 1 subset of the set of all infinite sequences. • The class RAND is a \Sigma^0_2 subset of Cantor space, where \Sigma^0_2 refers to the second level of the arithmetical hierarchy. This is because a sequence S is in RAND if and only if there is some open set in the universal effective null cover that does not contain S; this property can be seen to be definable by a \Sigma^0_2 formula. • There is a random sequence which is \Delta^0_2, that is, computable relative to an oracle for the Halting problem. (Schnorr 1971) Chaitin's Ω is an example of such a sequence. • No random sequence is decidable, computably enumerable, or co-computably-enumerable. Since these correspond to the \Delta_1^0, \Sigma_1^0, and \Pi_1^0 levels of the arithmetical hierarchy, this means that \Delta_2^0 is the lowest level in the arithmetical hierarchy where random sequences can be found. • Every sequence is Turing reducible to some random sequence. (Kučera 1985/1989, Gács 1986). Thus there are random sequences of arbitrarily high Turing degree. == Relative randomness ==

As each of the equivalent definitions of a Martin-Löf random sequence is based on what is computable by some Turing machine, one can naturally ask what is computable by a Turing oracle machine. For a fixed oracle A, a sequence B which is not only random but in fact, satisfies the equivalent definitions for computability relative to A (e.g., no martingale which is constructive relative to the oracle A succeeds on B) is said to be random relative to A. Two sequences, while themselves random, may contain very similar information, and therefore neither will be random relative to the other. Any time there is a Turing reduction from one sequence to another, the second sequence cannot be random relative to the first, just as computable sequences are themselves nonrandom; in particular, this means that Chaitin's Ω is not random relative to the halting problem. An important result relating to relative randomness is van Lambalgen's theorem, which states that if C is the sequence composed from A and B by interleaving the first bit of A, the first bit of B, the second bit of A, the second bit of B, and so on, then C is algorithmically random if and only if A is algorithmically random, and B is algorithmically random relative to A. A closely related consequence is that if A and B are both random themselves, then A is random relative to B if and only if B is random relative to A. == Stronger than Martin-Löf randomness ==

Relative randomness gives us the first notion which is stronger than Martin-Löf randomness, which is randomness relative to some fixed oracle A. For any oracle, this is at least as strong, and for most oracles, it is strictly stronger, since there will be Martin-Löf random sequences which are not random relative to the oracle A. Important oracles often considered are the halting problem, \emptyset ', and the nth jump oracle, \emptyset^{(n)}, as these oracles are able to answer specific questions which naturally arise. A sequence which is random relative to the oracle \emptyset^{(n-1)} is called n-random; a sequence is 1-random, therefore, if and only if it is Martin-Löf random. A sequence which is n-random for every n is called arithmetically random. The n-random sequences sometimes arise when considering more complicated properties. For example, there are only countably many \Delta^0_2 sets, so one might think that these should be non-random. However, the halting probability Ω is \Delta^0_2 and 1-random; it is only after 2-randomness is reached that it is impossible for a random set to be \Delta^0_2. == Weaker than Martin-Löf randomness ==

Additionally, there are several notions of randomness which are weaker than Martin-Löf randomness. Some of these are weak 1-randomness, Schnorr randomness, computable randomness, partial computable randomness. Yongge Wang showed that Schnorr randomness is different from computable randomness. Additionally, Kolmogorov–Loveland randomness is known to be no stronger than Martin-Löf randomness, but it is not known whether it is actually weaker. At the opposite end of the randomness spectrum there is the notion of a K-trivial set. These sets are anti-random in that all initial segment is logarithmically compressible (i.e., K(w) \leq K(|w|) + b for each initial segment w), but they are not computable. ==See also==