OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Jean-Marie De Koninck and Rafael Jakimczuk, Summing the largest prime factor over integer sequences, Revista de la Unión Matemática Argentina, Vol. 67, No. 1 (2024), pp. 27-35.
Michael De Vlieger, Log log scatterplot of pi(a(n)), n = 2..2^20, highlighting in red pi(a(n)) such that b(n) is a prime power, where b(n) = A001694(n). The remaining terms are such that b(n) is in A286708. The topmost red line corresponds with b(n) in A001248 (prime squares), the row of red dots at the bottom corresponds with b(n) in A000079 (powers of 2), the topmost blue line corresponds with b(n) > 16 in A069262 (4*p^2, with prime p).
FORMULA
MATHEMATICA
s[n_] := Module[{f = FactorInteger[n]}, If[n == 1 || AllTrue[f[[;; , 2]], # > 1 &], f[[-1, 1]], Nothing]]; Array[s, 4000]
PROG
(PARI) lista(kmax) = {my(f); print1(1, ", "); for(k = 2, kmax, f = factor(k); if(vecmin(f[, 2]) > 1, print1(f[#f~, 1], ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Mar 03 2024
STATUS
approved