A065426 Incrementally larger terms in the continued fraction (A065645) for the twin prime constant (A005597).

0, 1, 16, 18, 21, 405, 1199, 2301, 19965

Offset

1,3

Links

Table of n, a(n) for n=1..9.

Gerhard Niklasch, Some number-theoretical constants: 1000-digit values

Mathematica

(* tpc copied from Niklasch reference *)
cof = ContinuedFraction[tpc, 969]; a = -1; k = 1; Do[ While[ cof[[k]] <= a, k++ ]; a = cof[[k]]; Print[a], {n, 1, 9} ]

Prog

(PARI) \\ Increasing lprec to 30000 gives no further term beyond 19965.
a065246(lprec) = {localprec(lprec); my (m=-1, T=prodeulerrat(1-1/(p-1)^2, 1, 3), c=contfrac(T)); for (k=1, #c, if (c[k]>m, print(c[k], ", "); m=c[k]))};
a065246(1000) \\ Hugo Pfoertner, Aug 01 2023

Crossrefs

Cf. A005597 and A065645.

Sequence in context: A161874 A095953 A084800 * A308607 A097746 A081258

Adjacent sequences: A065423 A065424 A065425 * A065427 A065428 A065429

Keyword

nonn,hard,more

Author

Robert G. Wilson v, Nov 15 2001

Status

approved