|
%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
|
