|
| |
|
|
A013957
|
|
sigma_9(n), the sum of the 9th powers of the divisors of n.
|
|
5
| |
|
|
1, 513, 19684, 262657, 1953126, 10097892, 40353608, 134480385, 387440173, 1001953638, 2357947692, 5170140388, 10604499374, 20701400904, 38445332184, 68853957121, 118587876498, 198756808749, 322687697780, 513002215782
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,2
|
|
|
COMMENTS
| If the canonical factorization of n into prime powers is the product of p^e(p) then sigma_k(n) = Product_p ((p^((e(p)+1)*k))-1)/(p^k-1).
Sum_{d|n} 1/d^k is equal to sigma_k(n)/n^k. So sequences A017665-A017712 also give the numerators and denominators of sigma_k(n)/n^k for k = 1..24. The power sums sigma_k(n) are in sequences A000203 (k=1), A001157-A001160 (k=2,3,4,5), A013954-A013972 for k = 6,7,...,24. - comment from Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 05 2001.
|
|
|
FORMULA
| G.f. sum(k>=1, k^9*x^k/(1-x^k)). - Benoit Cloitre (benoit7848c(AT)orange.fr), Apr 21 2003
|
|
|
MATHEMATICA
| lst={}; Do[AppendTo[lst, DivisorSigma[9, n]], {n, 5!}]; lst [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Mar 11 2009]
|
|
|
PROG
| (PARI) a(n)=if(n<1, 0, sigma(n, 9))
(Other) sage: [sigma(n, 9)for n in xrange(1, 21)] # [From Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Jun 04 2009]
|
|
|
CROSSREFS
| Sequence in context: A076338 A111344 A017681 * A036087 A007487 A023878
Adjacent sequences: A013954 A013955 A013956 * A013958 A013959 A013960
|
|
|
KEYWORD
| nonn,mult
|
|
|
AUTHOR
| N. J. A. Sloane (njas(AT)research.att.com).
|
| |
|
|
