A051007 Continued fraction for prime constant A051006.

0, 2, 2, 2, 3, 12, 131, 1, 7, 1, 2, 1, 3, 3, 1, 2, 5, 39, 2, 1, 169, 2, 2, 2, 1, 1, 2, 5, 1, 2, 1, 199, 24, 7, 7, 1, 163, 1, 3, 2, 1, 2, 14, 1, 3, 5, 1, 1, 1, 1, 3, 3, 2, 1, 279, 1, 4, 1, 1, 4, 1, 92, 1, 16, 2, 3, 1, 5, 2, 25, 9, 1, 1, 2, 8, 1, 2, 5, 1, 1, 3, 1, 1, 2, 1, 7, 1, 1, 1, 5, 1, 2, 2, 5, 6, 1

Offset: 0

Author

Eric W. Weisstein

Keywords

• nonn
• cofr

Examples

			0.414682509851111660248109622... = 0 + 1/(2 + 1/(2 + 1/(2 + 1/(3 + ...)))). - _Harry J. Smith_, Jun 15 2009
		

Links

• Harry J. Smith, Table of n, a(n) for n = 0..20000
• Eric Weisstein's World of Mathematics, Prime Constant.
• G. Xiao, Contfrac
• Index entries for continued fractions for constants

Crossrefs

Cf. A010051, A051006. Increasing partial quotients are in A102878.

Programs

Mathematica

ContinuedFraction[ FromDigits[{{Table[ If[ PrimeQ[n], 1, 0], {n, 370}]}, 0}, 2], 95] (* Robert G. Wilson v, Jan 15 2005 *)

PARI

{ default(realprecision, 20080); x=0; m=67000; for (n=1, m, if (isprime(n), a=1, a=0); x=2*x+a; ); x=contfrac(x/2^m); for (n=0, 20000, write("b051007.txt", n, " ", x[n+1])); } \\ Harry J. Smith, Jun 15 2009