OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..388
FORMULA
a(n) = Sum_{p|n} p^n. - Wesley Ivan Hurt, Jun 14 2021
EXAMPLE
a(6) = 793 because the distinct prime divisors of 6 are 2 and 3, and 2^6 + 3^6 = 793.
MAPLE
f:= proc(n) local p;
add(p^n, p = numtheory:-factorset(n))
end proc:
map(f, [$1..30]); # Robert Israel, Feb 20 2024
MATHEMATICA
Prepend[Array[Plus@@First[Transpose[FactorInteger[#]^#]]&, 100, 2], 0]
Join[{0}, Table[Total[FactorInteger[n][[All, 1]]^n], {n, 2, 25}]] (* Harvey P. Dale, Jan 23 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Nov 22 2011
STATUS
approved