OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := If[OddQ[e], 2*(p^(e + 1) - 1)/(p^2 - 1), (p^(e + 2) + p^e - 2)/(p^2 - 1)]; s[1] = 1; s[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[5000], Divisible[#, s[#]] &] (* Amiram Eldar, Sep 03 2020 *)
PROG
(PARI) for(n=1, 10000, if(n%sumdiv(n, d, d/core(d))==0, print1(n, ", ")))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 08 2002
STATUS
approved