|
%I #15 Jan 25 2022 12:24:14
%S 5,6,9,10,15,18,19,20,21,22,33,36,37,38,39,40,41,42,43,44,45,46,69,72,
%T 73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,
%U 141,144,145,150,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168
%N a(1) = 5; for n >1, a(n) = a(n-1) + gcd(n, a(n-1)).
%D Eric S. Rowland, A simple prime-generating recurrence, Abstracts Amer. Math. Soc., 29 (No. 1, 2008), p. 50 (Abstract 1035-11-986).
%H T. D. Noe, <a href="/A134736/b134736.txt">Table of n, a(n) for n=1..1000</a>
%H Eric S. Rowland, <a href="https://arxiv.org/abs/0710.3217">A natural prime-generating recurrence</a>, arXiv:0710.3217 [math.NT], 2007-2008.
%t a[1] = 5; a[n_] := a[n] = a[n-1] + GCD[n, a[n-1]];
%t Array[a, 66] (* _Jean-François Alcover_, Oct 01 2018 *)
%t RecurrenceTable[{a[1]==5,a[n]==a[n-1]+GCD[n,a[n-1]]},a,{n,70}] (* _Harvey P. Dale_, Nov 24 2018 *)
%o (Haskell)
%o a134736 n = a134736_list !! (n-1)
%o a134736_list =
%o 5 : zipWith (+) a134736_list (zipWith gcd a134736_list [2..])
%o -- _Reinhard Zumkeller_, Nov 15 2013
%Y See A106108 for other cross-references.
%Y Cf. A230504, A134743 (first differences), A084662, A084663.
%K nonn
%O 1,1
%A _N. J. A. Sloane_, Jan 28 2008
|
