A066108 - OEIS

%I #21 Jun 05 2021 06:29:16

%S 1,6,30,276,3130,46914,823550,16781384,387421227,10000100110,

%T 285311670622,8916103456860,302875106592266,11112006930971730,

%U 437893890381622140,18446744078004584720,827240261886336764194,39346408075494930884190,1978419655660313589123998

%N Sum n^d over all divisors of n.

%C This is neither the Moebius transform nor the inverse Moebius transform of n^n, although it is close to them.

%H Harry J. Smith, <a href="/A066108/b066108.txt">Table of n, a(n) for n = 1..100</a>

%F a(n) = Sum_{d|n} n^d.

%F a(n) ~ n^n. - _Vaclav Kotesovec_, Jun 05 2021

%e n = 12: a(12) = A066106(12) = 8916103456860 = 8916100448256+2985984+20736+1728+144+12.

%e For comparison: M-transform of n^n at 12 = 8916100401348 = 8916100448256-46656-256+0+4+0 = A062793(12);

%e Inverse M-transform of n^n at 12 = 8916100495200 = 8916100448256+46656+256+27+4+1 = A062796(12).

%t Table[Sum[n^d, {d, Divisors@ n}], {n, 19}] (* _Michael De Vlieger_, Dec 20 2015 *)

%o (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

%o (PARI) a(n)=sumdiv(n,d, n^d );  /* _Joerg Arndt_, Oct 07 2012 */

%Y Cf. A062793, A062796.

%K nonn

%O 1,2

%A _Labos Elemer_, Dec 05 2001