login
Continued fraction for binary expansion of Liouville's number interpreted in base 2 (A092874).
6

%I #52 Jun 19 2021 09:09:06

%S 0,1,3,3,1,2,1,4095,3,1,3,3,1,4722366482869645213695,4,3,1,3,4095,1,2,

%T 1,3,3,1

%N Continued fraction for binary expansion of Liouville's number interpreted in base 2 (A092874).

%C The continued fraction of the number obtained by reading A012245 as binary fraction.

%C Except for the first term, the only values that occur in this sequence are 1, 2, 3, 4 and values 2^((m-1)*m!) - 1 for m > 2. The probability of occurrence P(a(n) = k) are given by:

%C P(a(n) = 1) = 1/3,

%C P(a(n) = 2) = 1/12,

%C P(a(n) = 3) = 1/3,

%C P(a(n) = 4) = 1/12 and

%C P(a(n) = 2^((m-1)*m!)-1) = 1/(3*2^(m-1)) for m > 2.

%C The next term is roughly 3.12174855*10^144 (see b-file for precise value).

%H A.H.M. Smeets, <a href="/A317413/b317413.txt">Table of n, a(n) for n = 0..48</a>

%F a(n) = 1 if and only if n in A317538.

%F a(n) = 2 if and only if n in {24*m - 19 | m > 0} union {24*m - 4 | m > 0}.

%F a(n) = 3 if and only if n in A317539.

%F a(n) = 4 if and only if n in {12*m + A014710(m-1) - 2*(A014710(m-1) mod 2) | m > 0}

%F a(n) = 2^((m-1)*m!)-1 if and only if n in {3*2^(m-2)*(1+k*4) - 1 | k >= 0} union {3*2^m-2)*(3+k*4) | k >= 0} for m > 2.

%e 0.76562505... = 0 +1/( 1+ 1/(3+1/(3+1/(1+1/(2+...))))) - _R. J. Mathar_, Jun 19 2021

%p with(numtheory): cfrac(add(1/2^factorial(n),n=1..7),24,'quotients'); # _Muniru A Asiru_, Aug 11 2018

%t ContinuedFraction[ FromDigits[ RealDigits[ Sum[1/10^n!, {n, 8}], 10, 10000], 2], 60] (* _Robert G. Wilson v_, Aug 09 2018 *)

%o (Python)

%o n,f,i,p,q,base = 1,1,0,0,1,2

%o while i < 100000:

%o ....i,p,q = i+1,p*base,q*base

%o ....if i == f:

%o ........p,n = p+1,n+1

%o ........f = f*n

%o n,a,j = 0,0,0

%o while p%q > 0:

%o ....a,f,p,q = a+1,p//q,q,p%q

%o ....print(a-1,f)

%Y Cf. A012245, A014710, A317538, A317539.

%Y Cf. A058304 (in base 10), A317414 (in base 3).

%K nonn,base,cofr

%O 0,3

%A _A.H.M. Smeets_, Jul 27 2018