|
Search: id:A066108
|
|
|
| A066108 |
|
Sum n^d over all divisors of n. |
|
+0 1
|
|
| 1, 6, 30, 276, 3130, 46914, 823550, 16781384, 387421227, 10000100110, 285311670622, 8916103456860, 302875106592266, 11112006930971730, 437893890381622140, 18446744078004584720, 827240261886336764194
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
COMMENT
|
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)=Apply[Plus, n^Divisors[n]]= Sum{n^d; d runs over divisors of n}.
|
|
EXAMPLE
|
n = 12: a(12) = A066106(12) = 8916103456860 = 8916100448256+2985984+20736+1728+144+12; For comparison: M-transform of n^n at 12 = 8916100401348 = 8916100448256-46656-256+0+4+0 = A062793(12); Inverse M-transform of n^n at 12 = 8916100495200 = 8916100448256+46656+256+27+4+1 = A062796(12);
|
|
PROGRAM
|
(PARI) { for (n=1, 100, d=divisors(n); a=sum(i=1, length(d), n^d[i]); write("b066108.txt", n, " ", a) ) } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Nov 15 2009]
|
|
CROSSREFS
|
Cf. A062793, A062796.
Sequence in context: A052585 A051821 A099031 this_sequence A110374 A111876 A119634
Adjacent sequences: A066105 A066106 A066107 this_sequence A066109 A066110 A066111
|
|
KEYWORD
|
nonn,new
|
|
AUTHOR
|
Labos E. (labos(AT)ana.sote.hu), Dec 05 2001
|
|
|
Search completed in 0.002 seconds
|