OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
EXAMPLE
168 = sigma(60) = sigma(78) = sigma(92) = sigma(123) = sigma(143) = sigma(167).
MAPLE
N:= 10^4: # for terms <= N
V:= Vector(N):
for n from 1 to N-1 do
s:= numtheory:-sigma(n);
if s <= N then V[s]:= V[s]+1 fi
od:
select(t -> V[t]=6, [$1..N]); # Robert Israel, Nov 21 2019
MATHEMATICA
a = Table[ 0, {10000} ]; Do[ s = DivisorSigma[ 1, n ]; If[ s < 10001, a[ [ s ] ]++ ], {n, 1, 10000} ]; Select[ Range[ 10000 ], a[ [ # ] ] == 6 & ]
PROG
(PARI) is(k) = invsigmaNum(k) == 6 \\ Amiram Eldar, Nov 18 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Apr 18 2001
STATUS
approved