OFFSET
1,3
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..202
Eric Weisstein's World of Mathematics, Jacobi Theta Functions
FORMULA
a(n) = A370040(n,0) for n >= 1.
G.f. A(x) = Sum_{n>=1} a(n)*x^n satisfies the following formulas.
(1) A(x) = (1 - theta_4(x))/2 / Product_{n>=1} (1 - x^(2*n))^3.
(2) A(x) = P(x)/Q(x) where P(x) = Sum_{n>=1} (-1)^(n-1) * x^(n^2) and Q(x) = Sum_{n>=0} (-1)^n * (2*n+1) * x^(n*(n+1)).
(3) A(x) = F(x,0) where F(x,y) is the g.f. of triangle A370040 (see comment).
EXAMPLE
G.f.: A(x) = x + 3*x^3 - x^4 + 9*x^5 - 3*x^6 + 22*x^7 - 9*x^8 + 52*x^9 - 22*x^10 + 111*x^11 - 51*x^12 + 230*x^13 - 108*x^14 + 451*x^15 - 222*x^16 + ...
which equals A(x) = P(x) / Q(x)
where
P(x) = x - x^4 + x^9 - x^16 + x^25 - x^36 + x^49 + ...
Q(x) = 1 - 3*x^2 + 5*x^6 - 7*x^12 + 9*x^20 - 11*x^30 + 13*x^42 + ...
PROG
(PARI) {a(n) = my(P = sum(m=1, sqrtint(n+1), (-1)^(m-1) * x^(m^2) +x*O(x^n)),
Q = sum(m=0, sqrtint(n+1), (-1)^m * (2*m+1) * x^(m*(m+1)) +x*O(x^n)));
polcoeff(P/Q, n)}
for(n=1, 50, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Feb 10 2024
STATUS
approved