OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
EXAMPLE
MATHEMATICA
s[n_] := Total[Ceiling[Sqrt[(p = FactorInteger[n][[;; , 1]])]]^2 - p]; m = 30; seq = Table[0, {m}]; c = 0; k = 1; While[c < m, n = s[k^2]; If[PrimeQ[n] && (i = PrimePi[n]) <= m && seq[[i]] == 0, c++; seq[[i]] = k^2]; k++]; seq (* Amiram Eldar, Dec 11 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Nov 21 2002
EXTENSIONS
Name edited by Amiram Eldar, Dec 11 2019
STATUS
approved