OFFSET
1,2
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..200
EXAMPLE
a(6) = 12 because the divisors of 12 are 1, 2, 3, 4, 6, 12 and 1^6 + 2^6 + 3^6 + 4^6 + 6^6 + 12^6 = 3037530 = 6*506255.
MAPLE
with(numtheory):
a:= proc(n) local k;
for k while irem (add (d &^n mod n, d=divisors(k)), n)<>0
do od; k
end:
seq (a(n), n=1..63);
MATHEMATICA
snk[n_]:=Module[{k=1}, While[!Divisible[DivisorSigma[n, k], n], k++]; k]; Array[ snk, 70] (* Harvey P. Dale, Jun 07 2016 *)
PROG
(PARI) a(n)=my(k); while(sigma(k++, n)%n, ); k \\ Charles R Greathouse IV, Nov 03 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Nov 03 2011
STATUS
approved