OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
EXAMPLE
1512 = sigma(480) = sigma(636) = sigma(736) = sigma(748) = sigma(830) = sigma(902) = sigma(1006) = sigma(1105) = sigma(1255) = sigma(1391) = sigma(1411) = sigma(1511).
MATHEMATICA
a = Table[ 0, {50000} ]; Do[ s = DivisorSigma[ 1, n ]; If[ s < 50001, a[ [ s ] ]++ ], {n, 1, 50000} ]; Select[ Range[ 50000 ], a[ [ # ] ] == 12 & ]
Take[Sort[Transpose[Select[Tally[DivisorSigma[1, Range[100000]]], #[[2]] == 12&]][[1]]], 50] (* Harvey P. Dale, Jan 18 2013 *)
PROG
(PARI) is(k) = invsigmaNum(k) == 12 \\ Amiram Eldar, Nov 18 2024, using Max Alekseyev's invphi.gp
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Robert G. Wilson v, Apr 18 2001
STATUS
approved