OFFSET
1,2
COMMENTS
The squarefree part of x, core(x), is the smallest integer such that x*core(x) is a square.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000
MATHEMATICA
core[n_]:=Block[{f=FactorInteger[n], p, e}, {p, e} = Transpose[f]; Times@@ (p^Mod[e, 2])]; Select[Range[10^4], core[#] == Plus @@ IntegerDigits[#] &] (* Giovanni Resta, Apr 21 2017 *)
PROG
(PARI) list(lim)=my(v=List()); forfactored(n=1, lim\1, if(core(n)==sumdigits(n[1]), listput(v, n[1]))); Vec(v) \\ Charles R Greathouse IV, Nov 05 2017
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Benoit Cloitre, May 09 2002
EXTENSIONS
Data corrected by Giovanni Resta, Apr 21 2017
STATUS
approved