OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
k = 96: divisors(96) = {1,2,3,4,6,8,12,16,24,32,48,96}, 12 divisors; divisors(sigma(96)) = {1,2,3,4,6,7,9,12,14,18,21,28,36,42,63,84,126,252}, 18 divisors; 12 < 18, so 96 is a term.
MAPLE
filter:= proc(n) uses numtheory; tau(n) < tau(sigma(n)) end proc:
select(filter, [$1..100]); # Robert Israel, Aug 03 2020
MATHEMATICA
Do[s=DivisorSigma[0, DivisorSigma[1, n]]; s0=DivisorSigma[0, n]; If[Greater[s0, s], Print[n]], {n, 1, 1000}]
Select[Range[100], DivisorSigma[0, #]<DivisorSigma[0, DivisorSigma[1, #]]&] (* Harvey P. Dale, Sep 22 2019 *)
PROG
(PARI) isok(k) = {my(f = factor(k)); numdiv(f) < numdiv(sigma(f)); } \\ Amiram Eldar, Mar 07 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Aug 13 2002
STATUS
approved
