OFFSET
1,2
COMMENTS
Multiplicative because it is the Dirichlet convolution of A000578 = n^3 and A101455 = [1 0 -1 0 1 0 -1 ...], which are both multiplicative. - Christian G. Bower, May 17 2005
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).
FORMULA
G.f.: Sum_{n>=1} n^3*x^n/(1+x^(2*n)). - Vladeta Jovovic, Oct 16 2002
From Amiram Eldar, Nov 04 2023: (Start)
Sum_{k=1..n} a(k) ~ c * n^4 / 4, where c = A175572. (End)
a(n) = Sum_{d|n} (n/d)^3*sin(d*Pi/2). - Ridouane Oudra, Sep 26 2024
MATHEMATICA
max = 40; s = Sum[n^3*x^(n-1)/(1+x^(2*n)), {n, 1, max}] + O[x]^max; CoefficientList[s, x] (* Jean-François Alcover, Dec 02 2015, after Vladeta Jovovic *)
s[n_] := If[OddQ[n], (-1)^((n-1)/2), 0]; (* A101455 *)
f[p_, e_] := (p^(3*e+3) - s[p]^(e+1))/(p^3 - s[p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 04 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, d^3*(((n/d) % 4)==1)) - sumdiv(n, d, d^3*(((n/d) % 4)==3)); \\ Michel Marcus, Feb 16 2015
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
N. J. A. Sloane, Dec 23 1999
EXTENSIONS
Offset changed from 0 to 1 by R. J. Mathar, Jul 15 2010
STATUS
approved