OFFSET
1,1
FORMULA
a(n) = sqrt(A065408(n)).
MATHEMATICA
bb={}; Do[If[PrimeQ[Apply[Plus, IntegerDigits[n^2]]], bb=Append[bb, n]], {n, 500}]; bb
PROG
(PARI) isok(n) = isprime(sumdigits(n^2)); \\ Michel Marcus, Jan 09 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, May 20 2005
STATUS
approved