OFFSET
1,1
COMMENTS
a(n)*prime(n) is a term of A108540 iff abs(phi-a(n)/prime(n))<1.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Golden Ratio
MATHEMATICA
a[n_] := Module[{x = GoldenRatio * Prime[n]}, p1 = NextPrime[x, -1]; p2 = NextPrime[p1]; If[x - p1 < p2 - x, p1, p2]]; Array[a, 57] (* Amiram Eldar, Nov 28 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jun 09 2005
EXTENSIONS
Corrected by T. D. Noe, Oct 25 2006
STATUS
approved