|
|
A108539
|
|
a(n) = that prime p such that p/prime(n) is nearest to phi, the golden ratio: phi=(1+sqrt(5))/2.
|
|
10
|
|
|
3, 5, 7, 11, 17, 23, 29, 31, 37, 47, 53, 59, 67, 71, 79, 83, 97, 97, 109, 113, 113, 127, 137, 149, 157, 163, 167, 173, 179, 181, 211, 211, 223, 223, 241, 241, 257, 263, 271, 281, 293, 293, 311, 313, 317, 317, 337, 359, 367, 373, 379, 389, 389, 409, 419, 421, 433
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(n)*prime(n) is a term of A108540 iff abs(phi-a(n)/prime(n))<1.
|
|
LINKS
|
|
|
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
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|