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)).

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

Original entry on oeis.org

1, 5, 200, 5125, 177500, 3952500, 150715625, 3185187500, 112844843750, 2783033593750, 86330708203125, 2019237027343750, 72195817812500000, 1591910699609375000, 50158322275878906250, 1322261581989501953125, 39183430287559814453125, 946961406814801025390625
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): A000041 (b=1), A271235 (b=2), A271236 (b=3), A303135 (b=4), this sequence (b=5).

Programs

  • Mathematica
    CoefficientList[Series[1/QPochhammer[25*x]^(1/5), {x, 0, 20}], x] (* Vaclav Kotesovec, Apr 19 2018 *)
    CoefficientList[Series[Product[(1-(25x)^n)^(-1/5),{n,20}],{x,0,20}],x] (* Harvey P. Dale, Nov 04 2021 *)

Formula

a(n) ~ exp(Pi*sqrt(2*n/15)) * 5^(2*n - 3/10) / (2^(7/5) * 3^(3/10) * n^(4/5)). - Vaclav Kotesovec, Apr 19 2018

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

Original entry on oeis.org

1, -4, -88, -992, -19360, -97152, -4296448, 4539392, -568015360, -127621120, -39357927424, 2424998313984, -38804685471744, 799759166930944, 4879962868940800, 41563181340426240, 585185165832486912, 55834295603426754560, -75535223925056208896
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): A010815 (b=1), A298411 (b=2), A303152 (b=3), this sequence (b=4), A303154 (b=5).

Programs

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

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

Original entry on oeis.org

1, -4, -24, -1248, 1632, -267136, -669440, -56925184, 597165568, -19934894080, 61831327744, -3209599664128, 47593545383936, -840449808072704, 8113679782510592, -350055154021040128, 5703847053344768000, -57129722970675609600, 704939718429511778304
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): A081362 (b=1), A298993 (b=2), A303130 (b=3), this sequence (b=4), A303132 (b=5).

Programs

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

Formula

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