OFFSET
1,2
COMMENTS
This is neither the Moebius transform nor the inverse Moebius transform of n^n, although it is close to them.
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..100
FORMULA
a(n) = Sum_{d|n} n^d.
a(n) ~ n^n. - Vaclav Kotesovec, Jun 05 2021
EXAMPLE
MATHEMATICA
Table[Sum[n^d, {d, Divisors@ n}], {n, 19}] (* Michael De Vlieger, Dec 20 2015 *)
PROG
(PARI) { for (n=1, 100, d=divisors(n); a=sum(i=1, length(d), n^d[i]); write("b066108.txt", n, " ", a) ) } \\ Harry J. Smith, Nov 15 2009
(PARI) a(n)=sumdiv(n, d, n^d ); /* Joerg Arndt, Oct 07 2012 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Dec 05 2001
STATUS
approved