OFFSET
0,2
COMMENTS
Obtained by formally "differentiating the Dedekind eta-function half a time".
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Don Zagier, Vassiliev invariants and a strange identity related to the Dedekind eta-function, Topology, vol.40, pp.945-960 (2001). See Eq. (20).
FORMULA
a(n) = b(24*n + 1) where b() is multiplicative and b(2^e) = b(3^e) = 0^e, b(p^e) = (1 + (-1)^e)/2 * p^(e/2) if p == 1, 11 (mod 12), b(p^e) = (1 + (-1)^e)/2 * (-p)^(e/2) if p == 5, 7 (mod 12).
G.f.: Sum_{k>0} Kronecker(12, k) * k * x^((k^2 - 1) / 24).
a(n) = sqrt(24*n + 1) * A010815(n).
EXAMPLE
G.f. = 1 - 5*x - 7*x^2 + 11*x^5 + 13*x^7 - 17*x^12 - 19*x^15 + 23*x^22 + ...
G.f. = q - 5*q^25 - 7*q^49 + 11*q^121 + 13*q^169 - 17*q^289 - 19*q^361 + ...
MATHEMATICA
a[ n_] := With[ {m = Sqrt[24 n + 1]}, If[ IntegerQ @ m, m KroneckerSymbol[ 12, m], 0]]; (* Michael Somos, Oct 15 2015 *)
a[ n_] := If[ n < 1, Boole[n == 0], Times @@ (If[ # < 5, 0, (1 - Mod[#2, 2]) (# KroneckerSymbol[ 12, #])^(#2/2)] & @@@ FactorInteger[ 24 n + 1])]; (* Michael Somos, Oct 15 2015 *)
s = QPochhammer[q] + O[q]^100; A010815 = CoefficientList[s, q]; nn = Range[0, Length[A010815]-1]; A134756 = Sqrt[24*nn+1]*A010815 (* Jean-François Alcover, Dec 01 2015 *)
PROG
(PARI) {a(n) = if( issquare( 24*n+1, &n), n * kronecker( 12, n), 0)};
(PARI) {a(n) = my(A, p, e); if( n<1, n==0, A = factor(24*n+1); prod(k = 1, matsize(A)[1], [p, e] = A[k, ]; if( (p<5) || (e%2), 0, (kronecker( 12, p) * p)^(e\2))))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Nov 08 2007
STATUS
approved