A307522 Expansion of Product_{k>=1} ((1 + x)^k - x^k)/((1 + x)^k + x^k).

1, -2, 2, -2, 4, -10, 22, -42, 72, -116, 188, -332, 662, -1432, 3148, -6736, 13784, -26894, 50254, -90782, 160856, -285230, 518170, -983710, 1964800, -4090002, 8705322, -18582722, 39219572, -81148034, 163946630, -323136562, 622125982, -1173528562, 2179230066

Offset: 0

Author

Seiichi Manyama, Apr 12 2019

Keywords

• sign
• easy

Crossrefs

Cf. A002448, A103198, A307520, A307521, A318570.

Programs

Maple

a(n) := 2*(-1)^n*add( binomial(n-1, n-k^2), k = 1..floor(sqrt(n))):
print(1, seq(a(n), n = 1..40)); # Peter Bala, Dec 31 2024

Mathematica

m = 34; CoefficientList[Series[Product[((1 + x)^k - x^k)/((1 + x)^k + x^k), {k, 1, m}], {x, 0, m}], x] (* Amiram Eldar, May 14 2021 *)

PARI

my(N=66, x='x+O('x^N)); Vec(prod(k=1, N, ((1+x)^k-x^k)/((1+x)^k+x^k)))

Formula

G.f.: theta_4(x/(1 + x)), where theta_4() is the Jacobi theta function.

From Peter Bala, Dec 31 2024: (Start)

For n >= 1, a(n) = 2 * (-1)^n * Sum_{k = 1..floor(sqrt(n))} binomial(n-1, n-k^2).

For n >= 1, |a(n)| = 2 * A103198(n). (End)