OFFSET
-1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = -1..10000
FORMULA
Expansion of (eta(q) / eta(q^13))^2 in powers of q.
Euler transform of period 13 sequence [ -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 0, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^3 + v^3 - 13*u*v - 4*u*v * (u+v) - (u*v)^2.
G.f. A(x) satisfies 0 = f(A(x), A(x^3)) where f(u, v) = (u^2 - u*v + v^2)^2 - u*v * (13 + 6*u + u^2) * (13 + 6*v + v^2).
G.f. is a period 1 Fourier series which satisfies f(-1 / (13 t)) = 13 g(t) where q = exp(2 Pi i t) and g() is the g.f. of A121597.
G.f.: x^(-1) * (Product_{k>0} (1 - x^k) / (1 - x^(13*k)))^2.
a(-1) = 1, a(n) = -(2/(n+1))*Sum_{k=1..n+1} A284587(k)*a(n-k) for n > -1. - Seiichi Manyama, Mar 29 2017
EXAMPLE
G.f. = 1/q - 2 - q + 2*q^2 + q^3 + 2*q^4 - 2*q^5 - 2*q^7 - 2*q^8 + q^9 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ (1/q) (QPochhammer[ q] / QPochhammer[ q^13])^2, {q, 0, n}]; (* Michael Somos, Jan 17 2015 *)
a[ n_] := With[{A = q Product[ (1 - q^k)^KroneckerSymbol[ 13, k], {k, n + 1}]}, SeriesCoefficient[ 1/A - 3 - A, {q, 0, n}]]; (* Michael Somos, Jan 17 2015 *)
PROG
(PARI) {a(n) = my(A); if( n<-1, 0, n++; A = x * O(x^n); polcoeff( (eta(x + A) / eta(x^13 + A))^2, n))};
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Sep 11 2007
STATUS
approved