OFFSET
1,1
COMMENTS
It can be shown that this sequence is monotonic.
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
J. C. Lagarias, E. M. Rains and N. J. A. Sloane, The EKG sequence, Exper. Math. 11 (2002), 437-446.
FORMULA
a(n) = A064664(p(n)).
MATHEMATICA
ekg[s_] := Block[{m = s[[-1]], k = 3}, While[MemberQ[s, k] || GCD[m, k] == 1, k++]; Append[s, k]];
A064413 = Nest[ekg, {1, 2}, 1000];
Position[A064413, _?PrimeQ] // Flatten (* Jean-François Alcover, Nov 03 2018, after Robert G. Wilson v in 064413 *)
PROG
(Haskell)
import Data.List (elemIndex)
import Data.Maybe (fromJust)
a064955 n = a064955_list !! (n-1)
a064955_list =
map ((+ 1) . fromJust . (`elemIndex` a064413_list)) a000040_list
-- Reinhard Zumkeller, Sep 17 2001
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Oct 30 2001
STATUS
approved