OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
f[p_, e_] := If[OddQ[e], p, 1]; sqf[n_] := Times @@ (f @@@ FactorInteger[n]); imomega[n_] := DivisorSum[n, PrimeNu[#] &]; Select[Range[32000], sqf[#] == imomega[#] &] (* Amiram Eldar, Apr 06 2020 *)
PROG
(PARI) for(n=1, 50000, if(core(n)==sumdiv(n, d, omega(d)), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 17 2002
STATUS
approved