TwoPi
← Back
Rows of π's digits with one block highlighted, next to a bar chart showing the expected 10% limiting frequency for each digit if π is normal

The numbers that contain everything

You’ve probably heard the claim that somewhere in the digits of π\pi is your phone number. Your birthday too, and, if you encode text as digits, every book ever written.

That claim is often presented as a fact about π\pi, but it isn’t one. We don’t actually know whether π\pi has the property needed to make it true.

The property is called normality. Roughly speaking, a number is normal in base 10 if its digits behave the way we’d expect random digits to behave in the long run: each digit occurs with frequency 1/101/10, each two-digit string with frequency 1/1001/100, each three-digit string with frequency 1/10001/1000, and so on.

That last part matters. If every finite string occurs with the expected frequency, then every finite string occurs infinitely many times. So a normal number really would contain your phone number, your birthday, and any finite piece of text you care to encode.

What makes normal numbers interesting is the mismatch between what we know in general and what we can prove in particular. Almost every real number is normal. But for familiar constants such as π\pi, ee, and 2\sqrt{2}, normality is still unproved.

Making “normal” precise

Saying that the digits of a number “look random” is vague, so let’s make it precise.

Fix a base b2b \ge 2, and write a number x[0,1)x \in [0,1) in base bb as

x=0.d1d2d3,x = 0.d_1d_2d_3\dots,

where each digit did_i belongs to {0,1,,b1}\{0,1,\dots,b-1\}.

The first thing we might ask is whether each individual digit occurs equally often.

Simply normal

A number xx is simply normal in base bb if every digit occurs with limiting frequency 1/b1/b. In other words, for every digit aa,

limn1n#{in:di=a}=1b.\lim_{n\to\infty} \frac{1}{n} \#\{\,i\le n : d_i=a\,\} = \frac{1}{b}.

In decimal, that means the proportion of each digit 0,1,,90,1,\dots,9 tends to 10%10\%.

But this isn’t enough. A sequence could contain the right number of zeroes and ones while arranging them in a very non-random way. Knowing how often individual digits occur tells us nothing about how often blocks such as 0000, 0101, or 314159314159 occur.

That leads to the stronger definition.

Normal in base b

A number is normal in base bb if every block of digits occurs with the frequency we would expect from independent, uniformly chosen digits.

For a block s=s1s2sks=s_1s_2\dots s_k of length kk,

limn1n#{in:didi+1di+k1=s}=1bk.\lim_{n\to\infty} \frac{1}{n} \#\{\,i\le n : d_id_{i+1}\dots d_{i+k-1}=s\,\} = \frac{1}{b^k}.

So in decimal, each digit has frequency 1/101/10, each two-digit block has frequency 1/1001/100, each three-digit block has frequency 1/10001/1000, and so on.

This is the property behind the claim from the introduction. If a number is normal, then every finite string of digits occurs not merely once, but infinitely many times.

Multiplying xx by bb shifts its base-bb digits one place to the left, and taking the result mod 11 discards whatever fell off the front – so bnxmod1b^n x \bmod 1 has digit sequence dn+1,dn+2,d_{n+1}, d_{n+2}, \dots. Equidistribution of that sequence on [0,1)[0,1) is exactly the statement that every digit block appears with its expected frequency. In short, xx is normal in base bb exactly when (bnxmod1)n0(b^n x \bmod 1)_{n \ge 0} is equidistributed on [0,1)[0,1)1.

Finally, normality depends on the base.

Absolutely normal

A number is absolutely normal if it is normal in every integer base b2b\ge2.

Thus absolutely normal     \implies normal in base bb     \implies simply normal in base bb, and neither implication can be reversed in general.

There is one important relation between bases. If bb and cc are multiplicatively dependent – meaning

bm=cnb^m=c^n

for some positive integers mm and nn – then normality in one base is equivalent to normality in the other2. So, for example, being normal in base 22 is equivalent to being normal in base 44 or base 88.

For multiplicatively independent bases, no such implication holds: there are numbers normal in one base but not in the other3.

A small technical point: we defined normality for numbers in [0,1)[0,1), while π=3.14159\pi = 3.14159\dots clearly is not.

This does not matter. The integer part contributes only finitely many digits, and finitely many digits cannot change a limiting frequency. Asking whether π\pi is normal is therefore the same as asking whether 0.141592650.14159265\dots is normal.

More generally, when talking about normality, we can ignore the integer part and look only at the digits after the radix point.

Borel’s theorem: almost all numbers are normal

In 1909, Émile Borel4 proved that almost every real number is absolutely normal – normal in every integer base at once.

Here, “almost every” has a precise meaning. The exceptions form a set of Lebesgue measure zero. Equivalently, if xx is chosen uniformly at random from [0,1)[0,1), then xx is absolutely normal with probability 11.

What follows is a modern probabilistic proof, not Borel’s original 1909 argument, but it establishes the same theorem. The basic idea is simple: treat the digits of a randomly chosen number as independent random variables, apply the strong law of large numbers, and then use the fact that there are only countably many bases and finite digit strings to check.

The digits of a random number behave like independent rolls. Fix a base bb and choose xx uniformly from [0,1)[0,1). Write

x=0.d1d2d3x = 0.d_1d_2d_3\dots

in base bb.

Suppose we ask for the first three decimal digits to be 33, 11, 44. The numbers with that prefix are exactly those in

[0.314, 0.315),[0.314,\ 0.315),

an interval of length 0.001=1030.001 = 10^{-3}. Since xx is uniform on [0,1)[0,1), the probability of landing in that interval – that is, P(d1d2d3=314)P(d_1d_2d_3 = 314) – is exactly 10310^{-3}, the same number you’d get from three independent, uniform digit draws: (1/10)3(1/10)^3.

Nothing about 314314 was special. Prescribing any kk digits a1aka_1 \dots a_k pins xx to an interval of length bkb^{-k}, so

P(d1=a1,,dk=ak)=bk=(1b)k=P(d1=a1)P(dk=ak).P(d_1=a_1,\dots,d_k=a_k) = b^{-k} = \left(\frac{1}{b}\right)^k = P(d_1=a_1)\cdots P(d_k=a_k).

The joint probability factors exactly as it should for independent variables – knowing the earlier digits tells you nothing about the next one. The same counting argument works for digits at any kk positions, not just a prefix. Each fixed digit cuts the surviving set of xx’s by a factor of 1/b1/b, regardless of what the unfixed digits in between are free to do. So fixing di1,,dikd_{i_1}, \dots, d_{i_k} restricts xx to a union of intervals of total length bkb^{-k}, giving the same factorization. Apart from the usual ambiguity at numbers with two base-bb expansions (like 0.4999=0.50000.4999\dots = 0.5000\dots in base 10 – a countable, hence measure-zero, set of exceptions), the digits

d1,d2,d3,d_1,d_2,d_3,\dots

are independent and uniformly distributed over {0,1,,b1}\{0,1,\dots,b-1\}.

Now fix one block. Let ss be some string of kk digits. We want to know how often ss appears in the expansion of xx.

There’s a small complication, though: consecutive occurrences can overlap. Take k=3k=3 and s=314s = 314: the block starting at position 11 uses d1d2d3d_1d_2d_3, while the block starting at position 22 uses d2d3d4d_2d_3d_4 – they share two digits, so the corresponding events aren’t independent, and the law of large numbers doesn’t directly apply.

Split the starting positions into kk groups by residue modulo kk. For k=3k=3, that’s 1,4,7,10,1,4,7,10,\dots, then 2,5,8,11,2,5,8,11,\dots, then 3,6,9,12,3,6,9,12,\dots. Within the first group, the blocks d1d2d3, d4d5d6, d7d8d9,d_1d_2d_3,\ d_4d_5d_6,\ d_7d_8d_9,\dots use completely disjoint digits, so they’re independent – each equal to 314314 with probability 1/10001/1000. The strong law of large numbers says the fraction of these blocks equal to 314314 converges to 1/10001/1000 almost surely, and the same holds for the other two groups. Since the three groups split the starting positions evenly, their combined frequency is still 1/10001/1000.

The same argument applies for a block of any length kk: split the starting positions into kk groups by residue modulo kk, apply the strong law within each group – where the blocks are disjoint and hence independent, each equal to ss with probability 1/bk1/b^k – and recombine. That gives

limn1n#{in:didi+1di+k1=s}=1bk\lim_{n\to\infty} \frac{1}{n} \#\{\,i\le n : d_id_{i+1}\dots d_{i+k-1}=s\,\} = \frac{1}{b^k}

almost surely, where nn counts starting positions rather than digits used (counting only blocks that fit entirely within the first nn digits, via ink+1i \le n-k+1, gives the same limit). For any fixed base bb, block length kk, and block ss, the set of numbers for which the required frequency fails has measure zero.

There are only countably many conditions. Normality asks for this statement to hold for every finite block, and absolute normality asks for it to hold in every integer base.

But there are only countably many choices involved: countably many bases b2b\ge2, countably many lengths k1k\ge1, and finitely many strings of each fixed length in each base. Altogether, there are still only countably many conditions.

Each condition has a measure-zero exceptional set. A countable union of measure-zero sets still has measure zero.

Outside that union, every finite block has the expected limiting frequency in every base. In other words, almost every real number is absolutely normal.

That is Borel’s theorem.

The catch: we can’t prove it for the numbers we care about

Borel’s theorem tells us that almost every real number is normal, and in fact absolutely normal. But take the numbers that turn up naturally in mathematics – π\pi, ee, 2\sqrt{2} – and the situation is very different.

We do not know whether any of them is normal in base 10. In base 2, we cannot even prove the weaker statement that their zeroes and ones occur equally often5.

Explicit normal numbers are known – if we are allowed to build a number specifically for the purpose. The simplest is Champernowne’s constant6,

C10=0.12345678910111213141516,C_{10}=0.12345678910111213141516\dots,

formed by writing the positive integers one after another. Champernowne proved in 1933 that this number is normal in base 10. Explicit absolutely normal numbers are also known7.

And it is easy to write down numbers that are not normal. Every rational number has an eventually periodic expansion in any integer base; a periodic expansion can only ever realize a bounded number of distinct length-kk blocks, while there are bkb^k possible ones, so most blocks of that length never occur at all, which rules out normality. The usual Liouville constant,

L=n=110n!=0.110001000000000000000001,L=\sum_{n=1}^{\infty}10^{-n!} =0.110001000000000000000001\dots,

is another example in base 10. Its 11’s become increasingly sparse – the proportion of 11’s among the first nn digits tends to 00, not 1/101/10 – so even the individual digits fail to occur with equal limiting frequency.

The standard Liouville constant above is not normal, but being a Liouville number does not itself prevent normality: there are Liouville numbers that are normal, as well as ones that are normal in no base8. It is the particular digit structure of the constant above, not its membership in the Liouville numbers as a class, that makes it non-normal.

For π\pi, computation gives us plenty of evidence but no proof. In 2016, Peter Trueb analyzed the first 22,459,157,718,361 decimal digits of π\pi and, separately, the first 18,651,926,753,033 hexadecimal digits9, counting every block of length one, two, and three in each expansion. The observed frequencies were consistent with what normality predicts. But any finite computation can only examine a finite prefix. Normality is a statement about what happens as the number of digits tends to infinity, so no amount of checking alone can settle it.

We can compute the digits of π\pi and ee to enormous precision, prove that both are transcendental, and write down series, products, integrals, and rapidly converging algorithms for calculating them. None of those descriptions has yielded the digit-frequency estimates needed to prove normality.

The equidistribution formulation from earlier makes the difficulty concrete. To prove that π\pi is normal in base bb, we would need to show that

π,  bπ,  b2π,  b3π,(mod1)\pi,\; b\pi,\; b^2\pi,\; b^3\pi,\dots \pmod 1

is equidistributed in [0,1)[0,1). This isn’t a matter of computing more digits. It’s about understanding the long-term behavior of the fractional parts of bnπb^n\pi for arbitrarily large nn. The familiar representations of π\pi – geometric and analytic – have not yielded a proof of the required equidistribution.

Why “almost every” tells us so little about one number

There is no contradiction between saying that almost every real number is normal and admitting that we cannot prove normality for π\pi.

A set of measure zero does not have to be empty, or even small in any ordinary sense. It can be infinite and highly structured. Measure tells us how large a set is inside the continuum of real numbers; it does not tell us whether a particular number belongs to it.

The same gap shows up with transcendence. Almost every real number is transcendental. The algebraic numbers are countable, so they have measure zero. But this fact alone tells us nothing about whether a specific number is transcendental.

Proving that ee is transcendental required Hermite’s work in 187310. Lindemann proved the transcendence of π\pi in 188211. And even today, we do not know whether a number as simple to write down as

π+e\pi + e

is transcendental12. In fact even the weaker question, whether π+e\pi + e is irrational, is open.

The theorem that almost every number is transcendental does not help with that question. It says that a randomly chosen real number will avoid the exceptional set with probability 11. π+e\pi + e is one particular point, and a measure-theoretic statement about almost all points does not tell us where that point lies.

Borel’s theorem has the same limitation. It tells us that the non-normal numbers occupy a measure-zero part of the real line. It does not give us a test that we can apply to π\pi, ee, or 2\sqrt{2}.

“Almost every number is normal” tells us a great deal about the real numbers as a whole, and almost nothing about any particular number we happen to care about.

So, does π\pi contain your phone number?

We don’t know.

If π\pi is normal, then yes: your phone number appears somewhere in its decimal expansion, and not just once, but infinitely many times. The same would be true of your birthday, this article encoded as digits, and every other finite string.

Borel’s theorem tells us that this is what happens for almost every real number. Computations tell us that the digits of π\pi behave as normality would predict, in the tests performed so far.

The claim behind the phone-number story – that π\pi contains every finite sequence of digits – is still unproved, no matter how plausible it looks.

Borel’s theorem is more than a century old. The simplest questions it raises about specific numbers are still open: is π\pi normal? Is ee normal? Is 2\sqrt{2} normal? All three are widely believed to be normal, but nobody has found a way to prove it.

Footnotes

  1. D. D. Wall, “Normal Numbers”, PhD thesis, University of California, Berkeley (1949) – the source of the normality/equidistribution equivalence.

  2. John E. Maxfield, “Normal k-tuples,” Pacific Journal of Mathematics 3, 189-196 (1953) – the dependent-base case.

  3. The independent-base existence result is due independently to J. W. S. Cassels and to Wolfgang M. Schmidt, “Über die Normalität von Zahlen zu verschiedenen Basen”, Acta Arithmetica 7, 299-309 (1962).

  4. Émile Borel, “Les probabilités dénombrables et leurs applications arithmétiques”, Rendiconti del Circolo Matematico di Palermo 27, 247-271 (1909) – the original theorem.

  5. Tanguy Rivoal, “On the Bits Counting Function of Real Numbers”, Journal of the Australian Mathematical Society 85(1), 95-111 (2008) – the abstract states explicitly that simple normality in base 2 for 2\sqrt{2}, ee, and π\pi is conjectural.

  6. D. G. Champernowne, “The Construction of Decimals Normal in the Scale of Ten”, Journal of the London Mathematical Society 8, 254-260 (1933) – the original construction and proof.

  7. Wacław Sierpiński, “Démonstration élémentaire du théorème de M. Borel sur les nombres absolument normaux et détermination effective d’un tel nombre”, Bulletin de la Société Mathématique de France 45, 125-132 (1917) – the first explicit construction, though the construction did not provide an algorithm for computing its digits. A computable reformulation was later given by Verónica Becher & Santiago Figueira, “An Example of a Computable Absolutely Normal Number”, Theoretical Computer Science 270(1-2), 947-958 (2002).

  8. Yann Bugeaud, “Nombres de Liouville et nombres normaux”, Comptes Rendus de l’Académie des Sciences 335, 117-120 (2002) – proves both that normal Liouville numbers exist and that Liouville numbers normal in no base exist.

  9. Peter Trueb, “Digit Statistics of the First 22.4 Trillion Decimal Digits of Pi” (2016) – the digit-frequency computation cited above.

  10. Charles Hermite, “Sur la fonction exponentielle,” Comptes Rendus de l’Académie des Sciences 77 (1873), 18-24, 74-79, 226-233, 285-293.

  11. Ferdinand von Lindemann, “Über die Zahl π,” Mathematische Annalen 20 (1882), 213-225.

  12. F. M. S. Lima, “Some Transcendence Results from a Harmless Irrationality Theorem” (2014), lists π+e\pi + e among the numbers for which not even an irrationality proof is known.