%I #57 Feb 24 2024 11:04:49

%S 1,1,28,1,126,28,344,1,757,126,1332,28,2198,344,3528,1,4914,757,6860,

%T 126,9632,1332,12168,28,15751,2198,20440,344,24390,3528,29792,1,37296,

%U 4914,43344,757,50654,6860,61544,126,68922,9632,79508,1332,95382,12168,103824,28

%N Sum of cubes of odd divisors of n.

%C The sum of cubes of even divisors of 2*k equals 8*A001158(k), and the sum of cubes of even divisors of 2*k-1 vanishes, for k >= 1. - _Wolfdieter Lang_, Jan 07 2017

%C Sum_{k>=1} 1/a(k) diverges. - _Vaclav Kotesovec_, Sep 21 2020

%H Reinhard Zumkeller, <a href="/A051000/b051000.txt">Table of n, a(n) for n = 1..10000</a>

%H John A. Ewell, <a href="http://www.pphmj.com/abstract/2504.htm">On a relation between two divisor functions</a>, JP Journal of Algebra, Number Theory and Applications, Vol. 7, No. 2 (2007), pp. 241-243.

%H J. W. L. Glaisher, <a href="https://books.google.com/books?id=bLs9AQAAMAAJ&amp;pg=RA1-PA1">On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares</a>, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/OddDivisorFunction.html">Odd Divisor Function</a>.

%H <a href="/index/Ge#Glaisher">Index entries for sequences mentioned by Glaisher</a>.

%F Dirichlet g.f.: (1-2^(3-s))*zeta(s)*zeta(s-3). Dirichlet convolution of (-1)^n*A176415(n) and A000578. - _R. J. Mathar_, Apr 06 2011

%F a(n) = Sum_{k=1..A001227(n)} A182469(n,k)^3. - _Reinhard Zumkeller_, May 01 2012

%F G.f.: Sum_{k>=1} (2*k - 1)^3*x^(2*k-1)/(1 - x^(2*k-1)). - _Ilya Gutkovskiy_, Jan 04 2017

%F Sum_{k=1..n} a(k) ~ Pi^4 * n^4 / 720. - _Vaclav Kotesovec_, Jan 31 2019

%F Multiplicative with a(2^e) = 1 and a(p^e) = (p^(3*e+3)-1)/(p^3-1) for p > 2. - _Amiram Eldar_, Sep 14 2020

%F For k>=0, a(2^k) = 1. - _Vaclav Kotesovec_, Sep 21 2020

%F G.f.: Sum_{n >= 1} x^n*(1 + 23*x^(2*n) + 23*x^(4*n) + x^(6*n))/(1 - x^(2*n))^4. See row 4 of A060187. - _Peter Bala_, Dec 20 2021

%F a(n) = Sum_{k=0..n-1} A000203(2*n-2*k-1)*A000203(2*k+1)/A006519(n)^3 (Ewell, 2007). - _Amiram Eldar_, Feb 24 2024

%t Table[Total[Select[Divisors[n],OddQ]^3],{n,50}] (* _Harvey P. Dale_, Jun 28 2012 *)

%t f[2, e_] := 1; f[p_, e_] := (p^(3*e + 3) - 1)/(p^3 - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Sep 14 2020 *)

%o (Haskell)

%o a051000 = sum . map (^ 3) . a182469_row

%o -- _Reinhard Zumkeller_, May 01 2012

%o (PARI) a(n) = sumdiv(n, d, (d%2)*d^3); \\ _Michel Marcus_, Jan 04 2017

%o (Python)

%o from sympy import divisor_sigma

%o def A051000(n): return int(divisor_sigma(n>>(~n&n-1).bit_length(),3)) # _Chai Wah Wu_, Jul 16 2022

%Y Cf. A000203, A000593, A001227, A001158, A006519, A050999, A051001, A051002.

%K nonn,mult

%O 1,3

%A _Eric W. Weisstein_