We call a number “holy” if it contains no 666 in its decimal expansion, and “unholy” otherwise. For instance, 12366621 and 666609 are unholy, while 7777 and 66166266366 are holy.

Question: Is the set {2n | n∈N,2n is holy} infinite?Of course, tons of similar questions can be asked by changing the number 666, the base 2, and the base for extension (we asked for decimal, so the default was 10). I do not feel that I am the first one asking this, and I appreciate it if someone gives me references if applicable.

But my thought is the following:

Conjecture: No.I will share my reasoning at the end of the post, but let us first see some facts:

Smallest unholy instances:

2157=1826877046663628647754606040895353774569915678722192=6277101735386680763835789423207666416102355444464034512896Then, we witnessed a cluster of unholy powers: 2218,2220,2222,2224,2226, and then kept holy for a while, until we hit the unholy 2243.

Largest holy instances: I did not throw in a lot of CPU time to pursue holy numbers, nor did I try hard enough to optimize my programs, but among the 3715 holy powers of 2, the largest of them are

225357,225896,226051,226667,229784.I tested up to around 2110000, but that is all I got. It probably will be reasonable for an average computer to test up to say 2106 or 2107, but I will be surprised to see a new holy number.

Statistics: For an integer n, let H(n) be the number of holy powers of 2 up to 2n.`n | H(n) || n | H(n) || n | H(n) 1000 | 875 || 11000 | 3567 || 21000 | 3700 2000 | 1560 || 12000 | 3602 || 22000 | 3703 3000 | 2059 || 13000 | 3621 || 23000 | 3705 4000 | 2442 || 14000 | 3645 || 24000 | 3707 5000 | 2747 || 15000 | 3655 || 25000 | 3709 6000 | 2984 || 16000 | 3670 || 26000 | 3712 7000 | 3171 || 17000 | 3682 || 27000 | 3714 8000 | 3332 || 18000 | 3689 || 28000 | 3714 9000 | 3440 || 19000 | 3693 || 29000 | 3714 10000 | 3514 || 20000 | 3695 || 30000 | 3715`

The heuristics of the conjecture:This is definitely not close to a proof at all, and I still hope if rigorous arguments exists:

The idea is that we want to estimate, for an integer n, the probability P(n) that 2n is holy, and then compute ∑∞n=1P(n).

We know that 2n has O(nln2) decimal digits, so there are O(nln2) groups of three. For each group there is a 1−10−3 chance to be not 666, so very roughly

P(n)=(1−10−3)nln2≈e−10−3 nln2.And note that

∞∑n=1P(n)≈∫∞n=0e−10−3 xln2dx<∞.The red "estimation line" in the figure above follows from this integral.

Of course, one can easily argue the properness of the heuristic above:

- The distribution of the digits close to the left are not uniform; they are affected by the growth of logarithmic functions.
- The distribution of the digits close to the right are not uniform; they are affected by the pattern of 2^n \pmod{10^k}.
- P(n) and P(n+i) are not independent, partially because of the awful choice of the number 666: 6\cdots 6 \times 2^2 = 26\cdots 64.
Any thoughts are appreciated.

**Answer**

This is not an answer, but only a few comments.

Your conjecture is indeed a known open problem, and in fact we expect that there is some N>0 such that 666 occurs among the digits of the decimal expansion of 2^n, for every n>N ! Of course, one can replace k=666 by any positive integer k (more accurately, any sequence of digits).

See also

this blog post (for the case k=0),

Problem 24 in *Unsolved Problems in Number Theory* by Richard Guy,

here,

or those OEIS sequences : A035064,

A030706,

A007377,

A035062.

These are closely related questions :

(1),

(2),

(3).

Here is a relevant open problem by Fürstenberg, which apparently gives an approach to the above conjecture (even though I'm not sure to understand why), according to the page 189 of *Ten Lectures on the Interface between Analytic Number Theory and Harmonic Analysis* (reference given in John Cook's post linked above).

See also problems 10.25 and 10.27 in Yann Bugeaud's *Distribution Modulo One and Diophantine Approximation*.

According to this paper, Fürstenberg's conjecture is just "far out of reach" !

Similarly, in §6.6, p. 138 of Bugeaud's book, it is mentioned that "Erdős asked how often the ternary expansion of 2^n omits the digit 2. It is likely that there are only finitely many such integers n but the problem is still beyond reach"^{[1]} !

For the sake of completeness, here is an excerpt of Fürstenberg's paper explaining the relation between the two problems (Conjectures 2 and 2' in his paper). Reference: *Intersections of Cantor sets and transversality of semigroups* in "Problems in Analysis. A Symposium in Honor of Salomon Bochner".

To sum up, we *expect* any given sequence of digits to occur in the decimal expansion of all but finitely many powers of 2, but this is currently an unsolved problem which is out of reach!

^{[1]} One of the best result known so far is given as exercise 6.4 in Bugeaud's book. Namely, if N(X) denotes the number of positive integers m \leq X such that 2^m has no "2" in its ternary expansion, then

N(X) \leq 1.62 X^{\log(2) / \log (3)}. The (intractable) conjecture is that N(X) is bounded.

**Attribution***Source : Link , Question Author : Hw Chu , Answer Author : Watson*