A064955 - OEIS

%I #16 Nov 03 2018 15:57:59

%S 2,5,10,14,20,28,33,37,43,57,61,67,74,81,89,100,107,115,128,134,138,

%T 151,160,167,182,189,197,203,207,216,236,253,259,264,279,287,297,305,

%U 314,328,336,344,363,371,377,381,401,420,430,438,444,458,462,474,487,501,510,517,530,540,549,557,581,587,599,606,629,639,655,664,670,681,699,707,724,730,736,756,766,781,798,802,814,819,833,848,857,874,882,889,898,911,927,942,953,961,971,997,1004,1029,1041,1059,1072,1080,1087,1099,1118,1135,1142,1150,1156,1175,1181,1190,1203,1223,1232,1242,1249,1258,1266,1287,1298,1306,1324,1350,1357,1378,1391,1398,1410,1425,1433,1442,1456,1470,1478,1504,1516,1542,1546,1564,1568,1578,1586,1610,1626,1638,1646,1652,1680,1686,1693,1702,1734,1739,1760

%N Position of n-th prime in A064413.

%C It can be shown that this sequence is monotonic.

%C A073734(a(n)) = A000040(n) for n > 1. [_Reinhard Zumkeller_, Sep 17 2001]

%H Reinhard Zumkeller, <a href="/A064955/b064955.txt">Table of n, a(n) for n = 1..10000</a>

%H J. C. Lagarias, E. M. Rains and N. J. A. Sloane, <a href="http://arXiv.org/abs/math.NT/0204011">The EKG sequence</a>, Exper. Math. 11 (2002), 437-446.

%H <a href="/index/Ed#EKG">Index entries for sequences related to EKG sequence</a>

%F a(n) = A064664(p(n)).

%t ekg[s_] := Block[{m = s[[-1]], k = 3}, While[MemberQ[s, k] || GCD[m, k] == 1, k++]; Append[s, k]];

%t A064413 = Nest[ekg, {1, 2}, 1000];

%t Position[A064413, _?PrimeQ] // Flatten (* _Jean-François Alcover_, Nov 03 2018, after _Robert G. Wilson v_ in 064413 *)

%o (Haskell)

%o import Data.List (elemIndex)

%o import Data.Maybe (fromJust)

%o a064955 n = a064955_list !! (n-1)

%o a064955_list =

%o    map ((+ 1) . fromJust . (`elemIndex` a064413_list)) a000040_list

%o -- _Reinhard Zumkeller_, Sep 17 2001

%Y Cf. A064413, A064425.

%K nonn

%O 1,1

%A _N. J. A. Sloane_, Oct 30 2001