OFFSET
1,2
COMMENTS
a(n) = sqrt(A171671(n)).
LINKS
MATHEMATICA
nn=5*10^5; list=Table[n + Total[IntegerDigits[n]], {n, nn}]; Select[Sqrt[Complement[Range[nn], list]], IntegerQ[#] &] (* Jayanta Basu, May 06 2013 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Zak Seidov, Dec 15 2009
EXTENSIONS
Changed the word "safe" in this entry to "self". - N. J. A. Sloane, Feb 26 2017
STATUS
approved