OFFSET
1,2
EXAMPLE
sum_{k=1..13}(k^k) = 312086923782437; the digits of 312086923782437 sum to 65, which is divisible by 13, so 13 is in the sequence.
MATHEMATICA
s = 0; Do[s += n^n; k = Plus @@ IntegerDigits[s]; If[Mod[k, n] == 0, Print[n]], {n, 1, 10^5}]
Module[{nn=2000, pwr}, pwr=Accumulate[Table[n^n, {n, nn}]]; Select[Thread[ {Range[nn], pwr}], Divisible[Total[IntegerDigits[#[[2]]]], #[[1]]]&]][[All, 1]] (* The program generates the first 13 terms of the sequence. To generate more, increase the nn constant, but the program may take a long time to run and the computer may run out of memory. *) (* Harvey P. Dale, Dec 19 2021 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Ryan Propper, Aug 06 2005
EXTENSIONS
More terms from Ryan Propper, Nov 13 2005
STATUS
approved