OFFSET
1,2
COMMENTS
LINKS
R. J. Mathar, Table of n, a(n) for n = 1..404
MAPLE
n := 1 :
for k from 1 to 10000 do
numtheory[sigma](k) ;
numtheory[sigma](%) ;
numtheory[sigma](%) ;
if type(%, 'odd') then
printf("%d %d\n", n, k) ;
n := n+1 ;
end if;
end do: # R. J. Mathar, Oct 21 2014
MATHEMATICA
sssoQ[n_]:=OddQ[Nest[DivisorSigma[1, #]&, n, 3]]; Select[Range[800], sssoQ] (* Harvey P. Dale, Jul 06 2019 *)
PROG
(PARI) for(n=1, 999, sigma(sigma(sigma(n)))%2 && print1(n", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Dec 28 2013
STATUS
approved