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-5 of 5 results.

A271235 G.f. equals the square root of P(4*x), where P(x) is the g.f. of the partition numbers (A000041).

Original entry on oeis.org

1, 2, 14, 68, 406, 1820, 10892, 48008, 266214, 1248044, 6454116, 29642424, 156638076, 707729176, 3551518936, 16671232784, 81685862790, 375557689292, 1843995831412, 8437648295384, 40779718859796, 188104838512840, 891508943457064, 4091507664092016, 19457793452994012, 88760334081132280, 415942096027738728, 1905990594266105648, 8875964207106121784, 40416438507461834160
Offset: 0

Views

Author

Paul D. Hanna, Apr 02 2016

Keywords

Comments

More formulas and information can be derived from entry A000041.
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1/2, g(n) = 4^n. - Seiichi Manyama, Apr 20 2018

Examples

			G.f.: A(x) = 1 + 2*x + 14*x^2 + 68*x^3 + 406*x^4 + 1820*x^5 + 10892*x^6 + 48008*x^7 + 266214*x^8 + 1248044*x^9 + 6454116*x^10 +...
where A(x)^2 = P(4*x).
RELATED SERIES.
P(x) = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 7*x^5 + 11*x^6 + 15*x^7 + 22*x^8 + 30*x^9 + 42*x^10 + 56*x^11 + 77*x^12 + 101*x^13 + 135*x^14 +...+ A000041(n)*x^n +...
1/A(x)^6 = 1 - 12*x + 320*x^3 - 28672*x^6 + 9437184*x^10 - 11811160064*x^15  + 57174604644352*x^21 +...+ (-1)^n*(2*n+1)*(4*x)^(n*(n+1)/2) +...
		

Crossrefs

Expansion of Product_{n>=1} (1 - ((b^2)*x)^n)^(-1/b): A000041 (b=1), this sequence (b=2), A271236 (b=3), A303135 (b=4), A303136 (b=5).

Programs

  • Mathematica
    nmax = 30; CoefficientList[Series[Product[1/Sqrt[1 - (4*x)^k], {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Apr 02 2016 *)
  • PARI
    {a(n) = polcoeff( prod(k=1,n, 1/sqrt(1 - (4*x)^k +x*O(x^n))),n)}
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sum( k=1, sqrtint(n+1), (4*x)^(k^2) / prod(j=1, k, 1 - (4*x)^j, 1 + x*O(x^n))^2, 1)^(1/2), n))};
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    N=99; x='x+O('x^N); Vec(prod(k=1, N, 1/(1-(4*x)^k)^(1/2))) \\ Altug Alkan, Apr 20 2018

Formula

G.f.: Product_{n>=1} 1/sqrt(1 - (4*x)^n).
Sum_{k=0..n} a(k) * a(n-k) = 4^n * A000041(n), for n>=0, where A000041(n) equals the number of partitions of n.
a(n) ~ 4^(n-1) * exp(sqrt(n/3)*Pi) / (3^(3/8) * n^(7/8)). - Vaclav Kotesovec, Apr 02 2016

A271236 G.f.: Product_{k>=1} 1/(1 - (9*x)^k)^(1/3).

Original entry on oeis.org

1, 3, 45, 450, 5805, 52326, 705591, 6190425, 77219325, 751178610, 8522919063, 80502824835, 975122402985, 8949951461925, 100088881882830, 1003346683458480, 10828622925516312, 104307212166072165, 1152197107898173875, 11048535008792967825, 119509353627934830327
Offset: 0

Views

Author

Vaclav Kotesovec, Apr 02 2016

Keywords

Comments

In general, for h>=1, if g.f. = Product_{k>=1} 1/(1-(h^2*x)^k)^(1/h), then a(n) ~ h^(2*n) * exp(Pi*sqrt(2*n/(3*h))) / (2^(5*h+3) * 3^(h+1) * h^(h+1) * n^(3*h+1))^(1/(4*h)).
This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1/3, g(n) = 9^n. - Seiichi Manyama, Apr 20 2018

Crossrefs

Expansion of Product_{n>=1} (1 - ((b^2)*x)^n)^(-1/b): A000041 (b=1), A271235 (b=2), this sequence (b=3), A303135 (b=4), A303136 (b=5).

Programs

  • Mathematica
    nmax = 30; CoefficientList[Series[Product[1/(1 - (9*x)^k)^(1/3), {k, 1, nmax}], {x, 0, nmax}], x]
  • PARI
    N=99; x='x+O('x^N); Vec(prod(k=1, N, 1/(1-(9*x)^k)^(1/3))) \\ Altug Alkan, Apr 20 2018

Formula

a(n) ~ 3^(2*n - 2/3) * exp(sqrt(2*n)*Pi/3) / (2^(3/2) * n^(5/6)).

A303135 Expansion of Product_{n>=1} (1 - (16*x)^n)^(-1/4).

Original entry on oeis.org

1, 4, 104, 1760, 39520, 590720, 14285056, 205151232, 4596467200, 75375073280, 1504196046848, 23673049726976, 525315968712704, 7912159583600640, 158055039529779200, 2726833423421800448, 51889395654107463680, 840470097284214292480, 16765991910040314839040
Offset: 0

Views

Author

Seiichi Manyama, Apr 19 2018

Keywords

Comments

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1/4, g(n) = 16^n.

Crossrefs

Expansion of Product_{n>=1} (1 - ((b^2)*x)^n)^(-1/b): A000041 (b=1), A271235 (b=2), A271236 (b=3), this sequence (b=4), A303136 (b=5).

Programs

  • Mathematica
    CoefficientList[Series[1/QPochhammer[16*x]^(1/4), {x, 0, 20}], x] (* Vaclav Kotesovec, Apr 19 2018 *)

Formula

a(n) ~ exp(sqrt(n/6)*Pi) * 2^(4*n - 33/16) / (3^(5/16) * n^(13/16)). - Vaclav Kotesovec, Apr 19 2018

A303154 Expansion of Product_{n>=1} (1 - (25*x)^n)^(1/5).

Original entry on oeis.org

1, -5, -175, -3250, -100625, -1015000, -58034375, -154171875, -22257500000, -154144921875, -6824828906250, 175448177734375, -8774446542968750, 164769756689453125, 756859169189453125, 9661555852294921875, -16148589271240234375, 81663068586871337890625
Offset: 0

Views

Author

Seiichi Manyama, Apr 19 2018

Keywords

Comments

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -1/5, g(n) = 25^n.

Crossrefs

Expansion of Product_{n>=1} (1 - ((b^2)*x)^n)^(1/b): A010815 (b=1), A298411 (b=2), A303152 (b=3), A303153 (b=4), this sequence (b=5).

Programs

  • PARI
    N=66; x='x+O('x^N); Vec(prod(k=1, N, (1-(25*x)^k)^(1/5)))

A303132 Expansion of Product_{n>=1} (1 + (25*x)^n)^(-1/5).

Original entry on oeis.org

1, -5, -50, -3875, 2500, -2046250, -12409375, -1087687500, 13232343750, -907225000000, 1545669140625, -362705679687500, 6007095839843750, -224713698632812500, 2118331116210937500, -226812683210205078125, 4765872641563720703125
Offset: 0

Views

Author

Seiichi Manyama, Apr 19 2018

Keywords

Comments

This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = 1/5, g(n) = -25^n.
In general, for h>=1, if g.f. = Product_{k>=1} (1 + (h^2*x)^k)^(-1/h), then a(n) ~ (-1)^n * exp(Pi*sqrt(n/(6*h))) * h^(2*n) / (2^(7/4) * 3^(1/4) * h^(1/4) * n^(3/4)). - Vaclav Kotesovec, Apr 20 2018

Crossrefs

Expansion of Product_{n>=1} (1 + ((b^2)*x)^n)^(-1/b): A081362 (b=1), A298993 (b=2), A303130 (b=3), A303131 (b=4), this sequence (b=5).

Programs

  • Mathematica
    CoefficientList[Series[(2/QPochhammer[-1, 25*x])^(1/5), {x, 0, 20}], x] (* Vaclav Kotesovec, Apr 20 2018 *)

Formula

a(n) ~ (-1)^n * exp(Pi*sqrt(n/30)) * 5^(2*n - 1/4) / (2^(7/4) * 3^(1/4) * n^(3/4)). - Vaclav Kotesovec, Apr 20 2018
Showing 1-5 of 5 results.