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