A291697 -id:A291697 - cp's OEIS frontend

A270919 Coefficient of x^n in Product_{k>=1} ((1 + x^k) / (1 - x^k))^n.

1, 2, 12, 80, 552, 3912, 28224, 206208, 1520784, 11297546, 84413912, 633713808, 4776117216, 36115518376, 273868321536, 2081866609920, 15859616674336, 121046064563376, 925411686479820, 7085465166635440, 54323193841192752, 416993869451825424, 3204447137019290944

Offset: 0

Vaclav Kotesovec, Mar 25 2016

nonn

From Peter Bala, Apr 18 2023: (Start)
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and all positive integers n and k.
Conjecture: the supercongruence a(p) == 2*p + 2 (mod p^2) holds for all primes p. Cf. A291697. (End)

Seiichi Manyama, Table of n, a(n) for n = 0..1118 (terms 0..500 from Vaclav Kotesovec)

Main diagonal of A288515.

Cf. A008485, A008705, A156616, A270913, A270923, A270924, A291697, A309986, A362408.

Table[SeriesCoefficient[Product[((1+x^k)/(1-x^k))^n, {k, 1, n}], {x, 0, n}], {n, 0, 25}]
Table[SeriesCoefficient[(QPochhammer[-1, x]/QPochhammer[x, x])^n, {x, 0, n}]/2^n, {n, 0, 25}]
(* Calculation of constants {d,c}: *) eq = FindRoot[{2*s*QPochhammer[r*s] == QPochhammer[-1, r*s], (Log[1 - r*s] + QPolyGamma[0, 1, r*s])/Log[r*s] + r*((Derivative[0, 1][QPochhammer][-1, r*s] - 2*s*Derivative[0, 1][QPochhammer][r*s, r*s]) / (2*QPochhammer[r*s])) == 1}, {r, 1/8}, {s, 2}, WorkingPrecision -> 1000]; {N[1/r /. eq, 120], val = Sqrt[(1 - r*s)*Log[r*s]^2*(QPochhammer[r*s] / (Pi*(-r*s*(-1 + r*s) * Log[r*s]*(4*(2*ArcTanh[1 - 2*r*s] + QPolyGamma[0, 1, r*s])* Derivative[0, 1][QPochhammer][r*s, r*s] + r*Log[r*s]*(Derivative[0, 2][QPochhammer][-1, r*s] - 2*s*Derivative[0, 2][QPochhammer][r*s, r*s])) + 2*QPochhammer[r*s] * (4*r*s*ArcTanh[1 - 2*r*s] + 2*(-1 + (-1 + r*s)*ArcTanh[1 - 2*r*s])*Log[1 - r*s] - (-1 + r*s)*(-2 + Log[r*s] - 2*Log[1 - r*s])*QPolyGamma[0, 1, r*s] + (-1 + r*s) * QPolyGamma[0, 1, r*s]^2 + (-1 + r*s)*(QPolyGamma[1, 1, r*s] - 2*r*s*Log[r*s]* Derivative[0, 0, 1][QPolyGamma][0, 1, r*s])))))] /. eq; N[Chop[val], -Floor[Log[10, Abs[Im[val]]]] - 3]} (* Vaclav Kotesovec, Oct 03 2023 *)

a(n) ~ c * d^n / sqrt(n), where d = 7.862983395705905261519347909953827161057584... and c = 0.299856802806668079413694689903953367699319...

a(n) = [x^n] 1/theta_4(x)^n, where theta_4() is the Jacobi theta function. - Ilya Gutkovskiy, Nov 03 2017

A289522 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=0} ((1 + x^(2j+1))/(1 - x^(2j+1)))^k.

Original entry on oeis.org

1, 1, 0, 1, 2, 0, 1, 4, 2, 0, 1, 6, 8, 4, 0, 1, 8, 18, 16, 6, 0, 1, 10, 32, 44, 32, 8, 0, 1, 12, 50, 96, 102, 56, 12, 0, 1, 14, 72, 180, 256, 216, 96, 16, 0, 1, 16, 98, 304, 550, 624, 428, 160, 22, 0, 1, 18, 128, 476, 1056, 1512, 1408, 816, 256, 30, 0, 1, 20, 162, 704, 1862, 3240, 3820, 3008, 1494, 404, 40, 0

Offset: 0

Ilya Gutkovskiy, Jul 07 2017

			Square array begins:
1,  1,   1,    1,    1,     1,  ...
0,  2,   4,    6,    8,    10,  ...
0,  2,   8,   18,   32,    50,  ...
0,  4,  16,   44,   96,   180,  ...
0,  6,  32,  102,  256,   550,  ...
0,  8,  56,  216,  624,  1512,  ...

Columns k=0-6 give: A000007, A080054, A007096, A261647, A014969, A261648, A014970.
Rows n=0-3 give: A000012, A005843, A001105, A217873.
Main diagonal gives A291697.

Table[Function[k, SeriesCoefficient[Product[((1 + x^(2 i + 1))/(1 - x^(2 i + 1)))^k, {i, 0, n}], {x, 0, n}]][j - n], {j, 0, 11}, {n, 0, j}] // Flatten
Table[Function[k, SeriesCoefficient[(QPochhammer[-x, x^2]/QPochhammer[x, x^2])^k, {x, 0, n}]][j - n], {j, 0, 11}, {n, 0, j}] // Flatten

G.f. of column k: Product_{j>=0} ((1 + x^(2*j+1))/(1 - x^(2*j+1)))^k.
G.f. of column 2k: (theta_3(x)/theta_4(x))^k, where theta_() is the Jacobi theta function.
For asymptotics of column k see comment from Vaclav Kotesovec in A261648.

A361008 G.f.: Product_{k >= 0} ((1 + x^(2k+1)) / (1 - x^(2k+1)))^k.

Original entry on oeis.org

1, 0, 0, 2, 0, 4, 2, 6, 8, 10, 20, 18, 42, 40, 78, 92, 140, 192, 258, 382, 480, 728, 902, 1334, 1698, 2404, 3148, 4292, 5742, 7608, 10304, 13430, 18192, 23592, 31720, 41144, 54766, 71188, 93762, 122156, 159420, 207820, 269380, 350726, 452434, 587520, 755446

Offset: 0

Vaclav Kotesovec, Apr 20 2023

nonn

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000

Cf. A080054, A291697, A263149, A263150.

Cf. A156616, A270919.

Table[SeriesCoefficient[Product[((1 + x^(2*k + 1))/(1 - x^(2*k + 1)))^k, {k, 0, n}], {x, 0, n}], {n, 0, 50}]

a(n) ~ sqrt(A/(3*Pi)) * (7*zeta(3))^(11/72) * exp(3*(7*zeta(3))^(1/3) * n^(2/3)/4 - Pi^2 * n^(1/3)/(8*(7*zeta(3))^(1/3)) - 1/24 - Pi^4/(1344*zeta(3))) / (2^(3/4) * n^(47/72)), where A = A074962 is the Glaisher-Kinkelin constant.

A294592 a(n) = [x^n] (theta_3(x)/theta_4(x))^n, where theta_() is the Jacobi theta function.

Original entry on oeis.org

1, 4, 32, 304, 3072, 32024, 340352, 3666016, 39878656, 437091892, 4819567552, 53401892240, 594093969408, 6631726263608, 74242911364864, 833237193123104, 9371924860764160, 105614054423502408, 1192210691317862048, 13478559927485340144, 152589996020498655232, 1729590806617202662528

Offset: 0

Ilya Gutkovskiy, Nov 03 2017

nonn

Robert Israel, Table of n, a(n) for n = 0..500
Eric Weisstein's MathWorld, Jacobi Theta Functions
Wikipedia, Theta function

Cf. A007096, A014969, A014970, A014972, A066535, A080054, A103261, A270919, A289522, A291697.

S:= series((JacobiTheta3(0,x)/JacobiTheta4(0,x))^n,x,51):
seq(coeff(S,x,n),n=0..50); # Robert Israel, Nov 03 2017

Table[SeriesCoefficient[(EllipticTheta[3, 0, x]/EllipticTheta[4, 0, x])^n, {x, 0, n}], {n, 0, 21}]
Table[SeriesCoefficient[Product[((1 + x^(2 k + 1))/(1 - x^(2 k + 1)))^(2 n), {k, 0, n}], {x, 0, n}], {n, 0, 21}]
Table[SeriesCoefficient[(QPochhammer[-x, x^2]/QPochhammer[x, x^2])^(2 n), {x, 0, n}], {n, 0, 21}]
(* Calculation of constant d: *) 1/r /. FindRoot[{s == EllipticTheta[3, 0, r*s]/EllipticTheta[4, 0, r*s], EllipticTheta[4, 0, r*s] + r*s*Derivative[0, 0, 1][EllipticTheta][4, 0, r*s] == r*Derivative[0, 0, 1][EllipticTheta][3, 0, r*s]}, {r, 1/10}, {s, 2}, WorkingPrecision -> 120] (* Vaclav Kotesovec, Jan 17 2024 *)

a(n) = [x^n] Product_{k>=0} ((1 + x^(2*k+1))/(1 - x^(2*k+1)))^(2*n).
From Vaclav Kotesovec, Nov 05 2017: (Start)
a(n) ~ c * d^n / sqrt(n), where
d = 11.61255065799699699891360038489317237925475956178123836149123386457... and
c = 0.34456510029264878768512693687607064416428117641473856418257649837... (End)

cp's OEIS Frontend

A270919 Coefficient of x^n in Product_{k>=1} ((1 + x^k) / (1 - x^k))^n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A289522 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=0} ((1 + x^(2j+1))/(1 - x^(2j+1)))^k.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

A361008 G.f.: Product_{k >= 0} ((1 + x^(2k+1)) / (1 - x^(2k+1)))^k.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A294592 a(n) = [x^n] (theta_3(x)/theta_4(x))^n, where theta_() is the Jacobi theta function.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

cp's OEIS Frontend

A270919 Coefficient of x^n in Product_{k>=1} ((1 + x^k) / (1 - x^k))^n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A289522 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=0} ((1 + x^(2*j+1))/(1 - x^(2*j+1)))^k.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

A361008 G.f.: Product_{k >= 0} ((1 + x^(2*k+1)) / (1 - x^(2*k+1)))^k.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A294592 a(n) = [x^n] (theta_3(x)/theta_4(x))^n, where theta_() is the Jacobi theta function.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A289522 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=0} ((1 + x^(2j+1))/(1 - x^(2j+1)))^k.

A361008 G.f.: Product_{k >= 0} ((1 + x^(2k+1)) / (1 - x^(2k+1)))^k.