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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A080350 Costé prime expansion of 1/Pi. 1
5, 2, 7, 5, 3, 5, 11, 3, 7, 5, 2, 7, 79, 67, 19, 29, 59, 11, 5, 11, 3, 47, 31, 607, 419, 47, 19, 43, 37, 7, 3, 29, 7, 7, 5, 5, 157, 53, 89, 17, 5, 37, 7, 11, 3, 17, 5, 3, 3, 29, 127, 19, 11, 11, 5, 3, 13, 41, 13, 11, 3, 11, 11, 5, 11, 3, 3, 5, 2, 5, 2, 5, 3, 7, 3, 17, 5, 11, 3, 11, 17, 5, 3, 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

Table of n, a(n) for n=0..83.

A. Costé, Sur un système fibré lié à la suite des nombres premiers, Exper. Math., 11 (2002), 383-405.

MAPLE

Digits := 500: 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(1/Pi);

CROSSREFS

Cf. A079385, A079386, A079366, A079367, A079368, A080348, A080349.

Sequence in context: A084340 A065925 A250720 * A204899 A253545 A195343

Adjacent sequences:  A080347 A080348 A080349 * A080351 A080352 A080353

KEYWORD

nonn,easy

AUTHOR

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 3 04:30 EDT 2020. Contains 336197 sequences. (Running on oeis4.)