A076271 a(1) = 1, a(2) = 2, and for n > 2, a(n) = a(n-1) + gpf(a(n-1)), where gpf = greatest prime factor = A006530.
1, 2, 4, 6, 9, 12, 15, 20, 25, 30, 35, 42, 49, 56, 63, 70, 77, 88, 99, 110, 121, 132, 143, 156, 169, 182, 195, 208, 221, 238, 255, 272, 289, 306, 323, 342, 361, 380, 399, 418, 437, 460, 483, 506, 529, 552, 575, 598, 621, 644, 667, 696, 725, 754, 783, 812, 841
Offset: 1
Keywords
Links
- T. D. Noe, Table of n, a(n) for n = 1..10000
Programs
-
Haskell
a076271 n = a076271_list !! (n-1) a076271_list = iterate a070229 1 -- Reinhard Zumkeller, Nov 07 2015
-
Mathematica
NestList[#+FactorInteger[#][[-1,1]]&,1,60] (* Harvey P. Dale, May 11 2015 *)
-
PARI
print1(n=1);for(i=1,199,print1(","n+=A006530(n))) \\ M. F. Hasler, Apr 08 2015
Formula
a(n+1) = A070229(a(n)). - Reinhard Zumkeller, Nov 07 2015
Extensions
Edited by M. F. Hasler, Apr 08 2015
Comments