cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 14 results. Next

A380582 a(n) = [x^n] G(x)^n, where G(x) = Product_{k >= 1} ((1 + x^k)/(1 - x^k))^(k^2) is the g.f. of A206622.

Original entry on oeis.org

1, 2, 24, 236, 2432, 25752, 277152, 3019088, 33186816, 367378814, 4089875024, 45741207228, 513537853952, 5784253405192, 65332622356032, 739706089046736, 8392732289277952, 95401363286044260, 1086232605119042424, 12386037358495697292, 141422619808922418432, 1616691574828234720352
Offset: 0

Views

Author

Peter Bala, Jan 27 2025

Keywords

Comments

Given an integer sequence {f(n) : n >= 0} with f(0) = 1, there is a unique power series F(x) with rational coefficients, where F(0) = 1, such that f(n) = [x^n] F(x)^n. F(x) is given by F(x) = series_reversion(x/E(x)), where E(x) = exp(Sum_{n >= 1} f(n)*x^n/n). Furthermore, if the series E(x) has integer coefficients then the series F(x) also has integer coefficients and the sequence {f(n)} satisfies the Gauss congruences: f(n*p^r) == f(n*p^(r-1)) (mod p^r) for all primes p and positive integers n and r (by Stanley, Ch. 5, Ex. 5.2(a), p. 72 and the Lagrange inversion formula).
Thus the present sequence satisfies the Gauss congruences. In fact, stronger congruences appear to hold for this sequence.
We conjecture that a(p) == 1 (mod p^3) for all primes p >= 5 (checked up to p = 61).
More generally, we conjecture that the supercongruence a(n*p^r) == a(n*p^(r-1)) (mod p^(3*r)) holds for all primes p >= 5 and positive integers n and r. Some examples are given below.
Let A, B be integers and let C be a positive integer. Define u(n) = [x^(C*n)] Product_{k >= 1} ((1 + x^k)^A * (1 - x^k)^B)^(k^2). The present sequence is the case A = 1, B = -1 and C = 1. We conjecture that the above supercongruences also hold for the sequence {u(n)} for all primes p >= 7.

Examples

			Examples of supercongruences:
a(7) - a(1) = 3019088 - 2 = 2*(3^3)*(7^3)*163 == 0 (mod 7^3)
a(13) - a(1) = 5784253405192 - 2 = 2*5*(13^4)*20252279 == 0 (mod 13^4)
a(2*11) - a(2) = 18501616629347623668448 - 24 = (2^3)*(11^3)*17*1951*4243*9817*1257719 == 0 (mod 11^3)
a(5^2) - a(5) = 1884578634304981694792832319004 - 256504 = (2^2)*(5^6)*193381* 155926684363405438573 == 0 (mod 5^6)
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    G(x) := series(exp(add( (1/4)*(sigma[3](2*k) - sigma[3](k))*x^k/k, k = 1..23 )),x,24):
    seq(coeftayl(G(x)^n, x = 0, n), n = 0..23);

Formula

a(n) = [x^n] exp( n*Sum_{k >= 1} (sigma_3(2*k) - sigma_3(k))/4 * x^k/k ).

A156616 G.f.: Product_{n>0} ((1+x^n)/(1-x^n))^n.

Original entry on oeis.org

1, 2, 6, 16, 38, 88, 196, 420, 878, 1794, 3584, 7032, 13572, 25792, 48352, 89512, 163774, 296444, 531234, 943072, 1659560, 2896376, 5015700, 8622108, 14718652, 24960138, 42062200, 70458160, 117349856, 194381704, 320295312, 525123604
Offset: 0

Views

Author

R. J. Mathar, Feb 11 2009

Keywords

Comments

Generating function for a sum over strict plane partitions weighted with 2 powered to their number of connected components.
The inverse Euler transform is apparently 2, 3, 6, 6, 10, 9, 14, 12, 18, 15, 22, 18, 26, 21, ..., A016825 interlaced with A008585. - R. J. Mathar, Apr 23 2009
In general, for m >= 1, if g.f. = Product_{k>=1} ((1+x^k)/(1-x^k))^(m*k), then a(n) ~ exp(m/12 + 3/2 * (7*m*Zeta(3)/2)^(1/3) * n^(2/3)) * m^(1/6 + m/36) * (7*Zeta(3))^(1/6 + m/36) / (A^m * 2^(2/3 + m/9) * sqrt(3*Pi) * n^(2/3 + m/36)), where Zeta(3) = A002117 and A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Aug 17 2015
In general, for m >= 0, if g.f. = Product_{k>=1} ((1+x^k)/(1-x^k))^(k^m), then a(n) ~ ((2^(m+2)-1) * Gamma(m+2) * Zeta(m+2) / (2^(2*m+3) * n))^((1-2*Zeta(-m))/(2*m+4)) * exp((m+2)/(m+1) * ((2^(m+2)-1) * n^(m+1) * Gamma(m+2) * Zeta(m+2) / 2^(m+1))^(1/(m+2)) + Zeta'(-m)) / sqrt((m+2)*Pi*n). - Vaclav Kotesovec, Aug 19 2015

Crossrefs

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[((1+x^k)/(1-x^k))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 17 2015 *)
  • PARI
    {a(n)=polcoeff(exp(sum(m=1,n,(sigma(2*m,2)-sigma(m,2))/2*x^m/m)+x*O(x^n)),n)} \\ Paul D. Hanna, May 01 2010

Formula

Convolve A000219 with A026007.
O.g.f.: exp( Sum_{n>=1} (sigma_2(2n) - sigma_2(n))/2 *x^n/n ), where sigma_2(n) is the sum of squares of divisors of n (A001157). - Paul D. Hanna, May 01 2010
a(n) ~ exp(1/12 + 3 * 2^(-4/3) * (7*Zeta(3))^(1/3) * n^(2/3)) * (7*Zeta(3))^(7/36) / (A * 2^(7/9) * sqrt(3*Pi) * n^(25/36)), where Zeta(3) = A002117 and A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Aug 17 2015
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A076577(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 30 2017
G.f.: A(x) = exp( 2*Sum_{n >= 0} x^(2*n+1)/((2*n+1)*(1 - x^(2*n+1))^2) ). Cf. A000122 and A302237. - Peter Bala, Dec 23 2021

A206623 G.f.: Product_{n>0} ( (1+x^n)/(1-x^n) )^(n^3).

Original entry on oeis.org

1, 2, 18, 88, 398, 1768, 7508, 30644, 121310, 467234, 1756080, 6457168, 23274788, 82381584, 286760344, 982874120, 3320800590, 11070619228, 36446345198, 118581503192, 381552358872, 1214868568728, 3829841265428, 11959828895612, 37013411304892, 113570015855642
Offset: 0

Views

Author

Paul D. Hanna, Feb 12 2012

Keywords

Comments

Convolution of A023872 and A248882. - Vaclav Kotesovec, Aug 19 2015

Examples

			G.f.: A(x) = 1 + 2*x + 18*x^2 + 88*x^3 + 398*x^4 + 1768*x^5 + 7508*x^6 +...
where A(x) = (1+x)/(1-x) * (1+x^2)^8/(1-x^2)^8 * (1+x^3)^27/(1-x^3)^27 *...
Also, A(x) = Euler transform of [2,15,54,120,250,405,686,960,1458,...]:
A(x) = 1/((1-x)^2*(1-x^2)^15*(1-x^3)^54*(1-x^4)^120*(1-x^5)^250*(1-x^6)^405*...).
		

Crossrefs

Cf. A156616, A206622, A206624, A001159 (sigma_4).

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[((1+x^k)/(1-x^k))^(k^3), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 19 2015 *)
  • PARI
    {a(n)=polcoeff(prod(m=1,n+1,((1+x^m)/(1-x^m+x*O(x^n)))^(m^3)),n)}
    
  • PARI
    {a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 4)-sigma(m, 4))/8*x^m/m)+x*O(x^n)), n)}
    
  • PARI
    {a(n)=local(InvEulerGF=x*(2+15*x+46*x^2+60*x^3+46*x^4+15*x^5+2*x^6)/(1-x^2+x*O(x^n))^4);polcoeff(1/prod(k=1,n,(1-x^k+x*O(x^n))^polcoeff(InvEulerGF,k)),n)}
    for(n=0,30,print1(a(n),", "))

Formula

G.f.: exp( Sum_{n>=1} (sigma_4(2*n) - sigma_4(n))/8 * x^n/n ), where sigma_4(n) is the sum of 4th powers of divisors of n (A001159).
Inverse Euler transform has g.f.: x*(2 + 15*x + 46*x^2 + 60*x^3 + 46*x^4 + 15*x^5 + 2*x^6)/(1-x^2)^4.
a(n) ~ (93*Zeta(5))^(59/600) * exp(5/4 * (93*Zeta(5)/2)^(1/5) * n^(4/5) + Zeta'(-3)) / (2^(59/100) * sqrt(5*Pi) * n^(359/600)), where Zeta(5) = A013663, Zeta'(-3) = A259068. - Vaclav Kotesovec, Aug 19 2015

A206624 G.f.: Product_{n>0} ( (1+x^n)/(1-x^n) )^(n^4).

Original entry on oeis.org

1, 2, 34, 228, 1414, 8872, 52876, 301136, 1662614, 8929406, 46738920, 239036116, 1197187780, 5882369976, 28397283056, 134864166352, 630819797174, 2908948327780, 13236421303742, 59477002686404, 264104800719672, 1159649708139680, 5037895127964316
Offset: 0

Views

Author

Paul D. Hanna, Feb 12 2012

Keywords

Comments

Convolution of A023873 and A248883. - Vaclav Kotesovec, Aug 19 2015
In general, for m >= 0, if g.f. = Product_{k>=1} ((1+x^k)/(1-x^k))^(k^m), then a(n) ~ ((2^(m+2)-1) * Gamma(m+2) * Zeta(m+2) / (2^(2*m+3) * n))^((1-2*Zeta(-m))/(2*m+4)) * exp((m+2)/(m+1) * ((2^(m+2)-1) * n^(m+1) * Gamma(m+2) * Zeta(m+2) / 2^(m+1))^(1/(m+2)) + Zeta'(-m)) / sqrt((m+2)*Pi*n). - Vaclav Kotesovec, Aug 19 2015
If m is even and m >= 2, then can be simplified as: a(n) ~ ((2^(m+2)-1) * Gamma(m+2) * Zeta(m+2) / (2^(2*m+3) * n))^(1/(2*m+4)) * exp((m+2)/(m+1) * ((2^(m+2)-1) * n^(m+1) * Gamma(m+2) * Zeta(m+2) / 2^(m+1))^(1/(m+2)) + (-1)^(m/2) * Gamma(m+1) * Zeta(m+1) / (2^(m+1) * Pi^m)) / sqrt((m+2)*Pi*n). - Vaclav Kotesovec, Aug 19 2015

Examples

			G.f.: A(x) = 1 + 2*x + 18*x^2 + 88*x^3 + 398*x^4 + 1768*x^5 + 7508*x^6 +...
where A(x) = (1+x)/(1-x) * (1+x^2)^16/(1-x^2)^16 * (1+x^3)^81/(1-x^3)^81 *...
Also, A(x) = Euler transform of [2,31,162,496,1250,2511,4802,7936,...]:
A(x) = 1/((1-x)^2*(1-x^2)^31*(1-x^3)^162*(1-x^4)^496*(1-x^5)^1250*(1-x^6)^2511*...).
		

Crossrefs

Cf. A015128 (m=0), A156616 (m=1), A206622 (m=2), A206623 (m=3), A001160 (sigma_5).

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[Product[((1+x^k)/(1-x^k))^(k^4), {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 19 2015 *)
  • PARI
    {a(n)=polcoeff(prod(m=1,n+1,((1+x^m)/(1-x^m+x*O(x^n)))^(m^4)),n)}
    
  • PARI
    {a(n)=polcoeff(exp(sum(m=1, n, (sigma(2*m, 5)-sigma(m, 5))/16*x^m/m)+x*O(x^n)), n)}
    
  • PARI
    {a(n)=local(InvEulerGF=x*(2+31*x+152*x^2+341*x^3+460*x^4+341*x^5+152*x^6+31*x^7+2*x^8)/(1-x^2+x*O(x^n))^5); polcoeff(1/prod(k=1,n,(1-x^k+x*O(x^n))^polcoeff(InvEulerGF,k)),n)}
    for(n=0,30,print1(a(n),", "))

Formula

G.f.: exp( Sum_{n>=1} (sigma_5(2*n) - sigma_5(n))/16 * x^n/n ), where sigma_5(n) is the sum of 5th powers of divisors of n (A001160).
Inverse Euler transform has g.f.: x*(2 + 31*x + 152*x^2 + 341*x^3 + 460*x^4 + 341*x^5 + 152*x^6 + 31*x^7 + 2*x^8)/(1-x^2)^5.
a(n) ~ exp(3*2^(2/3)*Pi*n^(5/6)/5 + 3*Zeta(5)/(4*Pi^4)) / (2^(7/6) * 3^(1/2) * n^(7/12)), where Zeta(5) = A013663. - Vaclav Kotesovec, Aug 19 2015
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A096960(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 30 2017

A285988 Expansion of Product_{n>0} ((1-x^n)/(1+x^n))^(n^2) in powers of x.

Original entry on oeis.org

1, -2, -6, -4, 22, 72, 92, -48, -522, -1294, -1624, 300, 6948, 19032, 30192, 20432, -45578, -202788, -437178, -599460, -311112, 1038624, 4023532, 8423280, 11892004, 8429270, -12073032, -60747944, -139842736, -223644552, -232762256, -15050944, 636838518
Offset: 0

Views

Author

Seiichi Manyama, Apr 30 2017

Keywords

Crossrefs

Product_{n>0} ((1-x^n)/(1+x^n))^(n^m): A002448 (m=0), A285675 (m=1), this sequence (m=2), A285990 (m=3), A285991 (m=4).

Formula

a(0) = 1, a(n) = -(2/n)*Sum_{k=1..n} A007331(k)*a(n-k) for n > 0.
G.f.: exp(Sum_{k>=1} (sigma_3(k) - sigma_3(2*k))*x^k/(4*k)). - Ilya Gutkovskiy, Apr 14 2019

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

Original entry on oeis.org

1, 2, 2, 10, 18, 36, 86, 150, 326, 608, 1164, 2230, 4046, 7632, 13622, 24868, 44222, 78304, 138312, 240138, 418648, 718292, 1233494, 2097350, 3552370, 5987642, 10026088, 16745600, 27779030, 45970868, 75650248, 124100970, 202720814, 329909400, 535132036
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 08 2017

Keywords

Comments

Convolution of A294749 and A294750.
In general, if g.f. = Product_{k>=1} ((1 + x^(2*k-1))/(1 - x^(2*k-1)))^(c2*k^2 + c1*k + c0) and c2 > 0, then a(n) ~ exp(Pi*sqrt(2) * c2^(1/4) * n^(3/4) / 3 + 7*(c1+c2) * Zeta(3) * sqrt(n) / (2*sqrt(c2) * Pi^2) + (Pi*(4*c0 + 2*c1 + c2) / (8*sqrt(2) * c2^(1/4)) - 49*(c1+c2)^2 * Zeta(3)^2 / (2^(3/2) * c2^(5/4) * Pi^5)) * n^(1/4) - (7*c0 + 21*c1/4 + c2 + 7*c0*c1/c2 + 7*c1^2/(2*c2)) * Zeta(3) / (4*Pi^2) + 22411*(c1+c2)^3 * Zeta(3)^3 / (196 * c2^2 * Pi^8) - (c1+c2)/24) * A^((c1+c2)/2) * (n^((c1+c2)/96 - 5/8) / (2^(c0/2 + (11*c1 + 5*c2)/48 + 9/4) * Pi^((c1+c2)/24) * c2^((c1+c2)/96 - 1/8))), where A is the Glaisher-Kinkelin constant A074962.

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[((1+x^(2*k-1))/(1-x^(2*k-1)))^(k^2), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(sqrt(2)*Pi * n^(3/4)/3 + 7*Zeta(3) * sqrt(n) / (2*Pi^2) + (Pi / (8*sqrt(2)) - 49*Zeta(3)^2 / (2^(3/2) * Pi^5)) * n^(1/4) + 22411*Zeta(3)^3 / (196*Pi^8) - Zeta(3)/(4*Pi^2) - 1/24) * sqrt(A) / (2^(113/48) * Pi^(1/24) * n^(59/96)), where A is the Glaisher-Kinkelin constant A074962.

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

Original entry on oeis.org

1, 2, 2, 2, 10, 18, 18, 18, 50, 100, 118, 118, 206, 438, 582, 582, 806, 1606, 2344, 2506, 3122, 5322, 8202, 9498, 11130, 16844, 26110, 32272, 37018, 52274, 78018, 100098, 115986, 155026, 223190, 291674, 345132, 439518, 618734, 811790, 972846, 1204190, 1653726
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 29 2017

Keywords

Comments

Convolution of A291649 and A291655.

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[((1+x^(k^2))/(1-x^(k^2)))^(k^2), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

log(a(n)) ~ 5 * Pi^(1/5) * ((8-sqrt(2)) * Zeta(5/2))^(2/5) * n^(3/5) / (4*3^(3/5)).

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

Original entry on oeis.org

1, 2, 40, 1320, 61984, 3797560, 287566368, 25957422400, 2721948311680, 325260627848442, 43635601119149040, 6494550360714973304, 1062063969900788407680, 189301256401392643093560, 36526821128512112807216192, 7585918627122817713267856320
Offset: 0

Views

Author

Vaclav Kotesovec, May 12 2018

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 20; Table[SeriesCoefficient[Product[((1+x^k)/(1-x^k))^(n^2), {k, 1, n}], {x, 0, n}], {n, 0, nmax}]
    nmax = 20; Table[SeriesCoefficient[(QPochhammer[-1, x]/2/QPochhammer[x])^(n^2), {x, 0, n}], {n, 0, nmax}]

Formula

a(n) ~ 2^(n - 1/2) * exp(n + 1/2) * n^(n - 1/2) / sqrt(Pi).

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

Original entry on oeis.org

1, 0, 2, 6, 14, 32, 74, 166, 370, 810, 1736, 3682, 7718, 15976, 32754, 66508, 133794, 266948, 528424, 1038178, 2025456, 3925360, 7559298, 14470162, 27540598, 52130440, 98159832, 183905636, 342896254, 636384748, 1175823512, 2163221030, 3963353706, 7232529308
Offset: 0

Views

Author

Vaclav Kotesovec, Nov 08 2017

Keywords

Comments

Convolution of A027999 and A258349.

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[((1+x^k)/(1-x^k))^(k*(k-1)/2), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

a(n) ~ exp(2*Pi * n^(3/4) / 3 - 7*Zeta(3) * sqrt(n) / (2*Pi^2) - 49*Zeta(3)^2 * n^(1/4) / (4*Pi^5) - 22411 * Zeta(3)^3 / (392*Pi^8) - Zeta(3) / (8*Pi^2) - 1/24) * sqrt(A) / (2^(23/12) * Pi^(1/24) * n^(59/96)), where A is the Glaisher-Kinkelin constant A074962.

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

Original entry on oeis.org

1, 2, 8, 26, 76, 216, 590, 1554, 3988, 9988, 24464, 58794, 138866, 322808, 739658, 1672372, 3734848, 8245956, 18012114, 38952586, 83448832, 177194716, 373111970, 779430870, 1615995262, 3326484686, 6800794428, 13813260736, 27881653590, 55942340000, 111601021856
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 03 2018

Keywords

Comments

Convolution of the sequences A000294 and A028377.

Crossrefs

Programs

  • Mathematica
    nmax = 30; CoefficientList[Series[Product[((1 + x^k)/(1 - x^k))^(k (k + 1)/2), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

G.f.: Product_{k>=1} ((1 + x^k)/(1 - x^k))^A000217(k).
a(n) ~ exp(2*Pi*n^(3/4)/3 + 7*Zeta(3)*sqrt(n) / (2*Pi^2) - 49*Zeta(3)^2 * n^(1/4) / (4*Pi^5) + 22411 * Zeta(3)^3 / (392*Pi^8) - Zeta(3)/(8*Pi^2) + 1/24) * Pi^(1/24) / (sqrt(A) * 2^(25/12) * n^(61/96)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Apr 08 2018
G.f.: A(x) = exp( 2*Sum_{n >= 0} x^(2*n+1)/((2*n+1)*(1 - x^(2*n+1))^3) ). Cf. A000122 and A156616. - Peter Bala, Dec 23 2021
Showing 1-10 of 14 results. Next