OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
core[n_] := Times @@ Power @@@ ({#[[1]], Mod[ #[[2]], 2]} & /@ FactorInteger[n]); Select[Range[10^5], DivisorSigma[0, #] == DivisorSigma[1, core[#]] &] (* Amiram Eldar, Jul 11 2019 after Zak Seidov at A007913 *)
PROG
(PARI) isok(n) = sigma(core(n)) == numdiv(n); \\ Michel Marcus, Aug 09 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 28 2002
STATUS
approved