OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := (-1)^e*e; q[1] = True; q[n_] := Plus @@ f @@@ FactorInteger[n] == 0; Select[Range[1000], q]
PROG
(PARI) isok(n) = {my(f = factor(n)); sum(i = 1, #f~, (-1)^f[i, 2]*f[i, 2]) == 0};
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Aug 06 2022
STATUS
approved