OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..1100
EXAMPLE
a(3) = 3 because 3! divides sigma(3!) = 12.
MAPLE
N:= 500: # to get a(1) .. a(N)
k:= 1; skf:= 1;
for n from 1 to N do
nf:= n!;
while skf mod nf <> 0 do
k:= k+1;
skf:= numtheory:-sigma(k!);
od:
A[n]:= k;
od:
seq(A[i], i=1..N); # Robert Israel, Aug 09 2016
MATHEMATICA
Table[k = 1; While[! Divisible[DivisorSigma[1, k!], n!], k++]; k, {n, 58}] (* Michael De Vlieger, Aug 08 2016 *)
PROG
(PARI) a(n) = {my(k = 1); while(sigma(k!) % n! != 0, k++); k; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Jul 29 2016
STATUS
approved