OFFSET
1,1
COMMENTS
abs(phi*a(n) - A108542(n)) < 1, where phi = golden ratio = (1+sqrt(5))/2.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_] := Module[{x = GoldenRatio * p}, p1 = NextPrime[x, -1]; p2 = NextPrime[p1]; q = If[x - p1 < p2 - x, p1, p2]; If[Abs[q - x] < 1, q, 0]]; seq = {}; p=1; Do[p = NextPrime[p]; If[f[p] > 0, AppendTo[seq, p]], {200}]; seq (* Amiram Eldar, Nov 28 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jun 09 2005
STATUS
approved