Crazy pattern in the simple continued fraction for $\sum_{k=1}^\infty \frac{1}{(2^k)!}$

The continued fraction of this series exhibits a truly crazy pattern and I found no reference for it so far. We have:

$$\sum_{k=1}^\infty \frac{1}{(2^k)!}=0.5416914682540160487415778421$$

But the continued fraction is just beautiful:

[1, 1, 5, 2, 69, 1, 1, 5, 1, 1, 12869, 2, 5, 1, 1, 69, 2, 5, 1, 1, 601080389, 2, 5, 2, 69, 1, 1, 5, 2, 12869, 1, 1, 5, 1, 1, 69, 2, 5, 1, 1, 1832624140942590533, 2, 5, 2, 69, 1, 1, 5, 1, 1, 12869, 2, 5, 1, 1, 69, 2, 5, 2, 601080389, 1, 1, 5, 2, 69, 1, 1, 5, 2, 12869, 1, 1, 5, 1, 1, 69, 2, 5, 1, 1, 23951146041928082866135587776380551749, 2, 5, 2, 69, 1, 1, 5, 1, 1, 12869, 2, 5, 1, 1, 69, 2, 5, 1, 1, 601080389, 2, 5, 2, 69, 1, 1, 5, 2, 12869, 1, 1, 5, 1, 1, 69, 2, 5, 2, 1832624140942590533, 1, 1, 5, 2, 69, 1, 1, 5, 1, 1, 12869, 2, 5, 1, 1, 69, 2, 5, 2, 601080389, 1, 1, 5, 2, 69, 1, 1, 5, 2,...]

All of these large numbers are not just random - they have a simple closed form:

$$A_n= \left( \begin{array}( 2^n \\ 2^{n-1} \end{array} \right) -1$$

$$A_1=1$$

$$A_2=5$$

$$A_3=69$$

$$A_4=12869$$

$$A_5=601080389$$

And so on. This sequence is not in OEIS, only the larger sequence is, which contains this one as a subsequence https://oeis.org/A014495

What is the explanation for this?

Is there a regular pattern in this continued fraction (in the positions of numbers)?

Is there generalizations for other sums of the form $\sum_{k=1}^\infty \frac{1}{(a^k)!}$?


Edit

I think a good move will be to rename the strings of small numbers:

$$a=1, 1, 5, 2,\qquad b=1, 1, 5, 1, 1,\qquad c=2,5,1,1,\qquad d=2, 5, 2$$

As a side note if we could set $1,1=2$ then all these strings will be the same.

Now we rewrite the sequence. I will denote $A_n$ by just their indices $n$:

$$[a, 3, b, 4, c, 3, c, 5, d, 3, a, 4, b, 3, c, 6, d, 3, b, 4, c, 3, d, 5, a, 3, a, 4, b, 3, c, 7, \\ d, 3, b, 4, c, 3, c, 5, d, 3, a, 4, b, 3, d, 6, a, 3, b, 4, c, 3, d, 5, a, 3, a,...]$$

$$[a3b4c3c5d3a4b3c6d3b4c3d5a3a4b3c7d3b4c3c5d3a4b3d6a3b4c3d5a3a,...]$$

Now we have new large numbers $A_n$ appear at positions $2^n$. And postitions of the same numbers are in a simple arithmetic progression with a difference $2^n$ as well.

Now we only have to figure out the pattern (if any exists) for $a,b,c,d$.

The $10~000$ terms of the continued fraction are uploaded at github here.


I also link my related question, from the iformation there we can conclude that the series above provide a greedy algorithm Egyptian fraction expansion of the number, and the number is irrational by the theorem stated in this paper.


actually your pattern is true and it is relatively easy to prove. Most of it can be found in an article from Henry Cohn in Acta Arithmetica (1996) ("Symmetry and specializability in continued fractions") where he finds similar patterns for other kind of continued fractions such as $\sum \frac{1}{10^{n!}}$. Curiously he doesn't mention your particular series although his method applies directly to it.

Let $[a_0,a_1,a_2,\dots,a_m]$ a continued fraction and as usual let $$ \frac{p_m}{q_m} = [a_0,a_1,a_2,\dots,a_m] $$ we use this lemma from the mentioned article (the proof is not difficult and only uses elementary facts about continued fractions):

[Folding Lemma] $$ \frac{p_m}{q_m} + \frac{(-1)^m}{xq_m^2} = [a_0,a_1,a_2,\dots,a_m,x,-a_m,-a_{m-1},\dots,-a_2,-a_1] $$

This involves negative integers but it can be easiy transformed in a similar expresion involving only posive numbers using the fact that for any continued fraction: $$ [\dots, a, -\beta] = [\dots, a-1, 1, \beta-1] $$ So $$ \frac{p_m}{q_m} + \frac{(-1)^m}{xq_m^2} = [a_0,a_1,a_2,\dots,a_m,x-1,1,a_m-1,a_{m-1},\dots,a_2,a_1] $$

With all this consider the write the $m$th partial sum of your series as $$ S_n = \sum_{k=1}^n \frac{1}{2^k!} = [0,a_1,a_2,\dots,a_m] = \frac{p_m}{q_m}$$

Where we can take always $m$ even (ie if $m$ is odd and $a_m>1$ then we can consider instead the continued fraction $[0,a_1,a_2,\dots,a_m-1,1]$ and so on).

Now $q_m = 2^n!$ we see it using induction on $n$: it is obvious for $n=1$ and if $S_{n-1} = P/2^{n-1}!$ then $$ S_n = \frac{P}{2^{n-1}!} + \frac{1}{2^n!} = \frac{P (2^n!/2^{n-1}!) + 1}{2^n!} $$ now any common factor of the numerator and the denominator is has to be a factor of $2^{n-1}!$ dividing also $P$, but this is impossible as both are coprime so $q_m = 2^n!$ and we are done.

Using the "positive" form of the folding lemma with $x = \binom{2^{n+1}}{2^n}$ we get:

$$ \frac{p_m}{q_m} + \frac{(-1)^m}{\binom{2^{n+1}}{2^n}(2^n!)^2} = \frac{p_m}{q_m} + \frac{1}{2^n!} = [0,a_1,a_2,\dots,a_m,\binom{2^{n+1}}{2^n}-1,1,a_m-1,a_{m-1},\dots,a_1] $$

And we get the "shape" of the continued fraction and the your $A_m$. Let's see several steps:

We start with the first term wich is $$ \frac{1}{2} = [0,2] $$ as $m$ is odd we change and use instead the continuos fraction $$ \frac{1}{2} = [0,1,1] $$ and apply the last formula getting $$ \frac{1}{2}+\frac{1}{2^2!} = [0,1,1,5,1,0,1] $$ We can dispose of the zeros using the fact that for any continued fraction: $$ [\dots, a, 0, b, \dots] = [\dots, a+b, \dots ] $$ so $$ \frac{1}{2}+\frac{1}{2^2!} = [0,1,1,5,2] $$ this time $m$ is even so we apply again the formula getting $$ \frac{1}{2}+\frac{1}{2^2!}+\frac{1}{2^3!} = [0,1,1,5,2,69,1,1,5,1,1] $$ again $m$ is even (and it will always continue even as easy to infer) so we apply again the formula getting as the next term $$ \frac{1}{2}+\frac{1}{2^2!}+\frac{1}{2^3!}+\frac{1}{2^4!} = [0,1,1,5,2,69,1,1,5,1,1,12869,1,0,1,5,1,1,69,2,5,1,1] $$ and we reduce it using the zero trick: $$ \frac{1}{2}+\frac{1}{2^2!}+\frac{1}{2^3!}+\frac{1}{2^4!} = [0,1,1,5,2,69,1,1,5,1,1,12869,2,5,1,1,69,2,5,1,1] $$

from now it is easy to see that the obtained continued fraction always has an even number of terms and we always have to remove the zero leaving a continued fraction ending again in 1,1. So the rule to obtain the continued fraction from here to an arbitrary number of termes is repeating the follewing steps: let the shape of the last continued fraction be $[0,1,1,b_1,\dots,b_k,1,1]$ then the next continued fraction will be $$[0,1,1,b_1,\dots,b_k,1,1,A_n,2,b_k,\dots,b_1,1,1 ]$$ from this you can easily derive the patterns you have found for the position of apperance of the different integers.


I have computed $10~000$ entries of the continued fraction, using Mathematica (the number itself was evaluated with $100~000$ digit precision).

The results show that the pattern is very simple for the most part.

First, we denote again:

$$A_n= \left( \begin{array}( 2^n \\ 2^{n-1} \end{array} \right) -1$$

The computed CF contains $A_n$ up to $n=13$, and the formula is numerically confirmed.

Now the positions of $A_n$ go like this ($P_n$ is the position of $A_n$ in the list of all CF entries):

$$P_2=P(5)=[3,8,13,18,23,28,\dots]$$

$$P_3=P(69)=[5,16,25,36,45,56,\dots]$$

$$P_4=[11,30,51,70,91,110,\dots]$$

$$P_5=[21,60,101,140,181,220,\dots]$$

$$P_6=[41,120,201,280,361,440,\dots]$$

But all these are just a combination of two arithmetic progressions for even and odd terms!

The first two terms are a little off, but the rest goes exactly like $P_4,P_5$, meaning for $n \geq 4$ we can write the general position for $A_n$ ($k=0,1,2,\dots$):

$$p_k(A_n)= \begin{cases} 5 \cdot 2^{n-3}(1+2 k)+1,\qquad k \text{ even} \\ 5 \cdot 2^{n-3}(1+2 k),\qquad \qquad k \text{ odd} \end{cases}$$

As special, different cases, we have:

$$p_k(5)=3+5 k,\qquad k=0,1,2,\dots$$

$$p_k(69)= \begin{cases} 2(3+5 k)-1,\qquad k \text{ even} \\ 2(3+5 k),\qquad \qquad k \text{ odd} \end{cases}$$


For $P_1=1$ and for $2$ I can't see any definite pattern so far.


Basically, we now know the explicit expression for every CF entry and all of its positions in the list, except for entries $1$ and $2$.

It's enough now to consider the positions of $2$, then we just fill the rest of the list with $1$. The positions of $2$ start like this:

[4, 12, 17, 22, 24, 29, 37, 42, 44, 52, 57, 59, 64, 69, 77, 82, 84, 92, 97, 102, 104, 109, 117, 119, 124, 132, 137, 139, 144, 149, 157, 162, 164, 172, 177, 182, 184, 189, 197, 202, 204, 212, 217, 219, 224, 229, 237, 239, 244, 252, 257, 262, 264, 269, 277, 279, 284, 292, 297, 299, 304, 309, 317, 322, 324, 332, 337, 342, 344, 349, 357, 362, 364, 372, 377, 379, 384, 389, 397, 402, 404, 412, 417, 422, 424, 429, 437, 439, 444, 452, 457, 459, 464, 469, 477, 479, 484, 492, 497, 502, 504, 509, 517, 522, 524, 532, 537, 539, 544, 549, 557, 559, 564, 572, 577, 582, 584, 589, 597,...]

So far I have found four uninterrupted patterns for $2$:

$$p_{1k}(2)=4+20k,\qquad k=0,1,2,\dots$$

$$p_{2k}(2)=17+20k,\qquad k=0,1,2,\dots$$

$$p_{3k}(2)=29+40k,\qquad k=0,1,2,\dots$$

$$p_{4k}(2)=12+40k,\qquad k=0,1,2,\dots$$

Edit

Discounting these four progressions the rest of the sequence is very close to $20k$, but some numbers are $20k+2$, while some $20k-1$ with no apparent pattern:

[22,42,59,82,102,119,139,162,182,202,219,239,262,279,299,322,342,362,379,402,422,439,459,479,502,522,539,559,582,599,619,642,662,682,699,722,742,759,779,802,822,842,859,879,902,919,939,959,982,1002,1019,1042,1062,1079,1099,1119,1142,1162,1179,1199,1222,1239,1259,1282,1302,1322,1339,1362,1382,1399,1419,1442,1462,1482,1499,1519,1542,1559,1579,1602,1622,1642,1659,1682,1702,1719,1739,1759,1782,1802,1819,1839,1862,1879,1899,1919,1942,1962,1979,2002,2022,2039,2059,2082,2102,2122,2139,2159,2182,2199,2219,2239,2262,2282,2299,2322,2342,2359,2379,2399,2422,2442,2459,2479,2502,2519,2539,2562,2582,2602,2619,2642,2662,2679,2699,2722,2742,2762,2779,2799,2822,2839,2859,2882,2902,2922,2939,2962,2982,2999,3019,3039,3062,3082,3099,3119,3142,3159,3179,3202,3222,3242,3259,3282,3302,3319,3339,3362,3382,3402,3419,3439,3462,3479,3499,3519,3542,3562,3579,3602,3622,3639,3659,3679,3702,3722,3739,3759,3782,3799,3819,3839,3862,3882,3899,3922,3942,3959,3979,4002,4022,4042,4059,4079,4102,4119,4139,4162,4182,4202,4219,4242,4262,4279,4299,4319,4342,4362,4379,4399,4422,4439,4459,4479,4502,4522,4539,4562,4582,4599,4619,4642,4662,4682,4699,4719,4742,4759,4779,4799,4822,4842,4859,4882,4902,4919,4939,4959,4982,5002,...]


The $10~000$ terms of the continued fraction are uploaded at github here. You can check my conjectures or try to obtain the full pattern.


I would also like some hints and outlines for a proof of the above conjectures.

I understand that it's quite likely that any of the patterns I found break down for some large $k$.