login
A094889
Odd nonunitary abundant numbers.
9
33075, 99225, 165375, 231525, 259875, 297675, 363825, 429975, 467775, 496125, 552825, 562275, 571725, 606375, 628425, 675675, 694575, 716625, 760725, 779625, 798525, 826875, 848925, 883575, 893025, 921375, 937125, 959175, 1003275
OFFSET
1,1
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1); nuabQ[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (1 + Power @@@ FactorInteger[n]) > n; Select[Range[3, 10^6, 2], nuabQ] (* Amiram Eldar, May 12 2019 *)
PROG
(PARI) { usigma(n)=local(s=1, fac, i); fac=factor(n); for(i=1, matsize(fac)[1], s=s*(1+fac[i, 1]^fac[i, 2])); return(s); } nusigma(n)=sigma(n)-usigma(n); forstep(n=1, 2^24, 2, if(nusigma(n)>n, print1(n", ")))
CROSSREFS
Cf. A064597.
Sequence in context: A118278 A118280 A062682 * A332200 A031660 A252291
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Jun 16 2004
EXTENSIONS
Corrected and extended by Jason Earls, Jun 18 2004
STATUS
approved