The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A079387 Costé prime expansion of sqrt(3) - 1. 3
 2, 3, 3, 7, 5, 7, 3, 37, 7, 3, 149, 19, 41, 17, 7, 3, 11, 2, 11, 17, 23, 19, 11, 5, 3, 5, 3, 3, 5, 2, 5, 2, 23, 7, 13, 13, 19, 37, 7, 41, 29, 11, 2, 11, 3, 3, 7, 7, 3, 23, 7, 19, 11, 11, 17, 11, 7, 5, 7, 5, 5, 3, 5, 2, 5, 7, 19, 31, 19, 17, 7, 5, 11, 3, 3, 3, 103, 853, 211, 23, 19, 17, 11, 7, 5 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,1 COMMENTS For x in (0,1], define P(x) = min{p: p prime, 1/x < p}, Phi(x) = P(x)x - 1. Costé prime expansion of x(0) is sequence a(0), a(1), ... given by x(n) = Phi(x(n-1)) (n>0), a(n) = P(x(n)) (n >= 0). LINKS G. C. Greubel, Table of n, a(n) for n = 0..2000 A. Costé, Sur un système fibré lié à la suite des nombres premiers, Exper. Math., 11 (2002), 383-405. MAPLE Digits := 200: P := proc(x) local y; y := ceil(evalf(1/x)); if isprime(y) then y else nextprime(y); fi; end; F := proc(x) local y, i, t1; y := x; t1 := []; for i from 1 to 100 do p := P(y); t1 := [op(t1), p]; y := p*y-1; od; t1; end; F(sqrt(3)-1); MATHEMATICA \$MaxExtraPrecision = 500; P[x_] := Module[{y}, y = Ceiling[1/x]; If[PrimeQ[y], y, NextPrime[y]]]; F[x_] := Module[{y, i, t1}, y = x; t1 = {}; For[i = 1, i <= 100, i++, AppendTo[t1, p = P[y]]; y = p*y - 1]; t1]; F[Sqrt[3] - 1] (* G. C. Greubel, Jan 20 2019 *) CROSSREFS Cf. A079388, A079389, A079366, A079367, A079368. Sequence in context: A097359 A208525 A209574 * A141061 A256447 A076557 Adjacent sequences:  A079384 A079385 A079386 * A079388 A079389 A079390 KEYWORD nonn,easy AUTHOR N. J. A. Sloane, Feb 16 2003 EXTENSIONS More terms from Mark Hudson (mrmarkhudson(AT)hotmail.com), Feb 17 2003 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 17 09:44 EST 2021. Contains 340214 sequences. (Running on oeis4.)