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

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

Original entry on oeis.org

1, 1, 3, 8, 15, 34, 67, 133, 255, 486, 901, 1649, 2984, 5312, 9373, 16342, 28221, 48283, 81928, 137858, 230278, 381919, 629156, 1029933, 1675856, 2711288, 4362575, 6983196, 11122327, 17630798, 27820283, 43706461, 68375137, 106534093, 165340844, 255643289
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 07 2015

Keywords

Crossrefs

Programs

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

Formula

a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4 / (2592*Zeta(3)) - Pi^2 * n^(1/3) / (12*(3*Zeta(3))^(1/3)) + 3^(4/3)/2 * Zeta(3)^(1/3) * n^(2/3)) / (2^(1/6) * 3^(1/3) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117.
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k*(1 + x^k)/(k*(1 - x^k)^2)). - Ilya Gutkovskiy, Jun 07 2018

A219555 Number of bipartite partitions of (i,j) with i+j = n into distinct pairs.

Original entry on oeis.org

1, 2, 4, 10, 19, 38, 73, 134, 242, 430, 749, 1282, 2171, 3622, 5979, 9770, 15802, 25334, 40288, 63560, 99554, 154884, 239397, 367800, 561846, 853584, 1290107, 1940304, 2904447, 4328184, 6422164, 9489940, 13967783, 20480534, 29920277, 43557272, 63194864
Offset: 0

Views

Author

Alois P. Heinz, Nov 22 2012

Keywords

Examples

			a(2) = 4: [(2,0)], [(1,1)], [(1,0),(0,1)], [(0,2)].
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, add(
           b(n-i*j, min(n-i*j, i-1))*binomial(i+1, j), j=0..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..42);  # Alois P. Heinz, Sep 19 2019
  • Mathematica
    nmax=50; CoefficientList[Series[Product[(1+x^k)^(k+1),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Mar 07 2015 *)

Formula

a(n) = Sum_{i+j=n} [x^i*y^j] 1/2 * Product_{k,m>=0} (1+x^k*y^m).
G.f.: Product_{k>=1} (1+x^k)^(k+1). - Vaclav Kotesovec, Mar 07 2015
a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4 / (1296*Zeta(3)) + Pi^2 * n^(1/3) / (2^(5/3) * 3^(4/3) * Zeta(3)^(1/3)) + (3/2)^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(5/4) * 3^(1/3) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117. - Vaclav Kotesovec, Mar 07 2015
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k*(2 - x^k)/(k*(1 - x^k)^2)). - Ilya Gutkovskiy, Aug 11 2018

A052812 A simple grammar: power set of pairs of sequences.

Original entry on oeis.org

1, 0, 1, 2, 3, 6, 9, 16, 24, 42, 63, 102, 157, 244, 373, 570, 858, 1290, 1930, 2858, 4228, 6208, 9084, 13216, 19175, 27666, 39804, 57020, 81412, 115820, 164264, 232178, 327220, 459796, 644232, 900214, 1254554, 1743896, 2418071, 3344896, 4616026
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Number of partitions of n objects of two colors into distinct parts, where each part must contain at least one of each color. - Franklin T. Adams-Watters, Dec 28 2006

Crossrefs

Programs

  • Maple
    spec := [S,{B=Sequence(Z,1 <= card),C=Prod(B,B),S= PowerSet(C)},unlabeled]: seq(combstruct[count](spec,size=n), n=0..20);
  • Mathematica
    nmax=50; CoefficientList[Series[Product[(1+x^k)^(k-1),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Mar 07 2015 *)

Formula

G.f.: exp(Sum((-1)^(j[1]+1)*(x^j[1])^2/(x^j[1]-1)^2/j[1], j[1]=1 .. infinity))
G.f.: Product_{k>=1} (1+x^k)^(k-1). - Vladeta Jovovic, Sep 17 2002
Weigh transform of b(n) = n-1. - Franklin T. Adams-Watters, Dec 28 2006
a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4/(1296*Zeta(3)) - Pi^2 * n^(1/3) / (3^(4/3) * 2^(5/3) * Zeta(3)^(1/3)) + (3/2)^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(1/4) * 3^(1/3) * n^(2/3) * sqrt(Pi)), where Zeta(3) = A002117. - Vaclav Kotesovec, Mar 07 2015

Extensions

More terms from Vladeta Jovovic, Sep 17 2002

A255836 G.f.: Product_{k>=1} (1+x^k)^(3*k+1).

Original entry on oeis.org

1, 4, 13, 42, 117, 310, 785, 1896, 4433, 10062, 22248, 48080, 101821, 211682, 432795, 871520, 1730491, 3391894, 6568996, 12580316, 23841774, 44742634, 83193865, 153347110, 280336704, 508499474, 915540681, 1636805438, 2906642396, 5128530946, 8993376689
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 07 2015

Keywords

Crossrefs

Programs

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

Formula

a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4 / (3888*Zeta(3)) + Pi^2 * n^(1/3) / (6^(5/3) * Zeta(3)^(1/3)) + 3^(5/3)/2^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(17/12) * 3^(1/6) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117.

A120844 Number of multi-trace BPS operators for the quiver gauge theory of the orbifold C^2/Z_2.

Original entry on oeis.org

1, 3, 11, 32, 90, 231, 576, 1363, 3141, 7003, 15261, 32468, 67788, 138892, 280103, 556302, 1089991, 2108332, 4030649, 7620671, 14261450, 26431346, 48544170, 88393064, 159654022, 286149924, 509137464, 899603036, 1579014769
Offset: 0

Views

Author

Amihay Hanany (hanany(AT)mit.edu), Aug 25 2006

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> 2*n+1): seq(a(n), n=0..50); # Vaclav Kotesovec, Mar 06 2015 after Alois P. Heinz
    # alternative program
    with(numtheory):
    series(exp(add((2*sigma[2](k) + sigma[1](k))*x^k/k, k = 1..30)), x, 31):
    seq(coeftayl(%, x = 0, n), n = 0..30); # Peter Bala, Jan 16 2025
  • Mathematica
    nmax=50; CoefficientList[Series[Product[1/(1-x^k)^(2*k+1),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Feb 27 2015 *)

Formula

G.f.: exp( Sum_{n>0} (3*x^n - x^(2*n)) / (n*(1-x^n)^2) ).
a(n) ~ Zeta(3)^(7/18) * exp(1/6 - Pi^4/(864*Zeta(3)) + Pi^2 * n^(1/3)/(3 * 2^(5/3) * Zeta(3)^(1/3)) + 3 * (Zeta(3)/2)^(1/3) * n^(2/3)) / (A^2 * 2^(2/9) * 3^(1/2) * Pi * n^(8/9)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... . - Vaclav Kotesovec, Mar 07 2015
From Peter Bala, Jan 16 2025: (Start)
G.f.: 1/Product_{k >= 1} (1 - x^k)^(2*k+1).
G.f.: exp(Sum_{k >= 1} (2*sigma_2(k) + sigma_1(k))*x^k/k) = 1 + 3*x + 11*x^2 + 32*x^3 + 90*x^4 + 231*x^5 + .... (End)

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

Original entry on oeis.org

1, 5, 23, 86, 295, 926, 2748, 7732, 20891, 54401, 137355, 337249, 808043, 1893402, 4348634, 9805669, 21741925, 47463473, 102133056, 216841459, 454648373, 942113618, 1930779697, 3915946921, 7864385266, 15647363323, 30858285440, 60345383394, 117065924679
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 07 2015

Keywords

Comments

In general, if g.f. = Product_{k>=1} 1/(1-x^k)^(m*k+c), m > 0, then a(n) ~ (m*Zeta(3))^(m/36 + c/6 + 1/6) * exp(m/12 - c^2 * Pi^4 / (432*m*Zeta(3)) + c * Pi^2 * n^(1/3) / (3 * 2^(4/3) * (m*Zeta(3))^(1/3)) + 3 * (m*Zeta(3))^(1/3) * n^(2/3) / 2^(2/3)) / (A^m * 2^(c/3 + 1/3 - m/36) * 3^(1/2) * Pi^((c+1)/2) * n^(m/36 + c/6 + 2/3)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Mar 08 2015

Crossrefs

Cf. A000219 (k), A005380 (k+1), A052847 (k-1), A120844 (2k+1), A253289 (2k-1), A255802 (2k+3), A255271 (3k+1).

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> 3*n+2): seq(a(n), n=0..50); # after Alois P. Heinz
    with(numtheory):
    series(exp(add((3*sigma[2](k) + 2*sigma[1](k))*x^k/k, k = 1..30)), x, 31):
    seq(coeftayl(%, x = 0, n), n = 0..30); # Peter Bala, Jan 16 2025
  • Mathematica
    nmax=50; CoefficientList[Series[Product[1/(1-x^k)^(3*k+2),{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ Zeta(3)^(7/12) * 3^(1/12) * exp(1/4 - Pi^4/(324*Zeta(3)) + Pi^2 * n^(1/3) / (3^(4/3) * (2*Zeta(3))^(1/3)) + 3^(4/3) * Zeta(3)^(1/3) * n^(2/3) / 2^(2/3)) / (A^3 * 2^(11/12) * Pi^(3/2) * n^(13/12)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... .
G.f.: exp(Sum_{k >= 1} (3*sigma_2(k) + 2*sigma_1(k))*x^k/k) = 1 + 5*x + 23*x^2 + 86*x^3 + 295*x^4 + .... - Peter Bala, Jan 16 2025

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

Original entry on oeis.org

1, 5, 22, 79, 259, 777, 2201, 5911, 15239, 37865, 91224, 213741, 488759, 1093173, 2396934, 5160756, 10928181, 22787949, 46848176, 95046026, 190466354, 377295743, 739319876, 1433974869, 2754597217, 5243308562, 9894376295, 18517966608, 34386781020, 63378252332
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 07 2015

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr(n-> 2*n+3): seq(a(n), n=0..50); # after Alois P. Heinz
    with(numtheory):
    series(exp(add((2*sigma[2](k) + 3*sigma[1](k))*x^k/k, k = 1..30)), x, 31):
    seq(coeftayl(%, x = 0, n), n = 0..30); # Peter Bala, Jan 16 2025
  • Mathematica
    nmax=50; CoefficientList[Series[Product[1/(1-x^k)^(2*k+3),{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ Zeta(3)^(13/18) * exp(1/6 - Pi^4/(96*Zeta(3)) + Pi^2 * n^(1/3) / (2^(5/3) * Zeta(3)^(1/3)) + 3 * (Zeta(3)/2)^(1/3) * n^(2/3)) / (A^2 * 2^(5/9) * 3^(1/2) * Pi^2 * n^(11/9)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant and Zeta(3) = A002117 = 1.202056903... .
G.f.: exp(Sum_{k >= 1} (2*sigma_2(k) + 3*sigma_1(k))*x^k/k) = 1 + 5*x + 22*x^2 + 29*x^3 + 777*x^4 + .... - Peter Bala, Jan 16 2025

A255837 G.f.: Product_{k>=1} (1+x^k)^(3*k+2).

Original entry on oeis.org

1, 5, 18, 61, 182, 506, 1338, 3369, 8172, 19197, 43833, 97636, 212748, 454461, 953505, 1968095, 4001627, 8024295, 15885484, 31074351, 60111277, 115071431, 218126868, 409662895, 762679151, 1408172844, 2579599582, 4690277001, 8467363674, 15182486586
Offset: 0

Views

Author

Vaclav Kotesovec, Mar 07 2015

Keywords

Comments

In general, if g.f. = Product_{k>=1} (1+x^k)^(m*k+c), m > 0, then a(n) ~ (m*Zeta(3))^(1/6) * exp(-c^2 * Pi^4 / (1296*m*Zeta(3)) + (c * Pi^2 * n^(1/3)) / (2^(5/3) * 3^(4/3) * (m*Zeta(3))^(1/3)) + 3^(4/3) * (m*Zeta(3))^(1/3) * n^(2/3) / 2^(4/3)) / (2^(m/12 + c/2 + 2/3) * 3^(1/3) * sqrt(Pi) * n^(2/3)). - Vaclav Kotesovec, Mar 08 2015

Crossrefs

Cf. A026007 (k), A219555 (k+1), A052812 (k-1), A255834 (2*k+1), A255835 (2*k-1), A255836 (3*k+1).
Cf. A255803.

Programs

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

Formula

a(n) ~ Zeta(3)^(1/6) * exp(-Pi^4/(972*Zeta(3)) + Pi^2 * n^(1/3) / (2^(2/3) * 3^(5/3) * Zeta(3)^(1/3)) + 3^(5/3)/2^(4/3) * Zeta(3)^(1/3) * n^(2/3)) / (2^(23/12) * 3^(1/6) * sqrt(Pi) * n^(2/3)), where Zeta(3) = A002117.

A363600 Number of partitions of n into distinct parts where there are k^2+1 kinds of part k.

Original entry on oeis.org

1, 2, 6, 20, 52, 140, 356, 880, 2123, 5016, 11610, 26400, 59130, 130476, 284216, 611592, 1301344, 2740194, 5713930, 11806144, 24184908, 49142504, 99091244, 198360536, 394342884, 778818658, 1528531702, 2982017956, 5784365082, 11158728448, 21413292868
Offset: 0

Views

Author

Seiichi Manyama, Jun 10 2023

Keywords

Crossrefs

Programs

  • PARI
    my(N=40, x='x+O('x^N)); Vec(prod(k=1, N, (1+x^k)^(k^2+1)))

Formula

G.f.: Product_{k>=1} (1+x^k)^(k^2+1).
a(0) = 1; a(n) = (-1/n) * Sum_{k=1..n} ( Sum_{d|k} (-1)^(k/d) * d * (d^2+1) ) * a(n-k).

Extensions

Name suggested by Joerg Arndt, Jun 11 2023
Showing 1-9 of 9 results.