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

A103261 Number of partitions of 2n into parts with 10 types c^1 c^2...C^10 of each part. The even parts appear with multiplicity 1 for each type . The odd parts occur with multiplicity 2 for each part.

Original entry on oeis.org

1, 20, 200, 1360, 7200, 32024, 125280, 443680, 1450240, 4435940, 12827888, 35346800, 93377920, 237675640, 585229760, 1398704736, 3253934080, 7386124520, 16392493800, 35634450320, 75992326592, 159199081600, 328027789600
Offset: 0

Views

Author

Noureddine Chair, Feb 16 2005

Keywords

Comments

This is also Sequence(A080054)^(10) or sequence(A007096)^(5).
In general, if j > 0 and g.f. = Product_{k>=0} ((1 + x^(2*k+1))/(1 - x^(2*k+1)))^j, then a(n) ~ exp(Pi*sqrt(j*n/2)) * j^(1/4) / (2^(j/2 + 7/4) * n^(3/4)). - Vaclav Kotesovec, Aug 28 2015

Examples

			a(2)=200 because we have 10 types of 4, 45 ways of writing 4 in terms of ten of 2's only or ten of 11's only and 100 ways of writing 2's combined with 11's so the total number of ways of writing 4 is 200.

Crossrefs

Cf. A080054 (j=1), A007096 (j=2), A261647 (j=3), A014969 (j=4), A261648 (j=5), A014970 (j=6), A014972 (j=8).

Programs

  • Maple
    series(product(((1+x^k)*(1-x^(2*k)))^(10)/((1-x^k)*(1+x^(2*k)))^(10),k=1..100),x=0,100);
  • Mathematica
    nmax=60; CoefficientList[Series[Product[((1+x^(2*k+1))/(1-x^(2*k+1)))^10,{k,0,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Aug 28 2015 *)

Formula

G.f.:(theta_4(0, x^2)/theta_4(0, x))^10= (theta_3(0, x)/theta_4(0, x))^5.
a(n) ~ exp(Pi*sqrt(5*n)) * 5^(1/4) / (64 * sqrt(2) * n^(3/4)). - Vaclav Kotesovec, Aug 28 2015

Extensions

Example corrected by Vaclav Kotesovec, Sep 01 2015

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

Original entry on oeis.org

1, 6, 18, 44, 102, 216, 428, 816, 1494, 2650, 4584, 7740, 12804, 20808, 33264, 52400, 81462, 125100, 189966, 285516, 425016, 627040, 917436, 1331856, 1919332, 2746926, 3905784, 5519352, 7754064, 10833192, 15055216, 20817600, 28647414, 39241336, 53517060
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Links

Crossrefs

Cf. A015128, A156616.
Cf. A080054, A007096, A014969, A261648, A014970, A014972, A103261.
Cf. A261610, A261649, A261651.
Cf. A261611, A261650, A261652.

Programs

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

Formula

a(n) ~ exp(Pi*sqrt(3*n/2)) * 3^(1/4) / (8 * 2^(1/4) * n^(3/4)).

A261652 Expansion of Product_{k>=0} ((1+x^(4*k+1))/(1-x^(4*k+1)))^3.

Original entry on oeis.org

1, 6, 18, 38, 66, 108, 182, 306, 486, 728, 1068, 1578, 2318, 3312, 4614, 6388, 8862, 12192, 16488, 22038, 29400, 39156, 51702, 67554, 87810, 113982, 147384, 189200, 241446, 307356, 390408, 493662, 621006, 778712, 974628, 1216284, 1511756, 1872840, 2315538
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 28 2015

Keywords

Comments

In general, if j > 0, a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} ((1 + x^(a*k+b))/(1 - x^(a*k+b)))^j, then a(n) ~ Gamma(b/a)^j * 2^(j/2 - 3/2 - 2*b*j/a) * a^(-j/4 - 1/4 + b*j/(2*a)) * exp(Pi*sqrt(j*n/a)) * j^(1/4 - j/4 + b*j/(2*a)) * Pi^(b*j/a - j) * n^(j/4 - 3/4 - b*j/(2*a)).

Crossrefs

Cf. A015128 (a=1, b=1, j=1), A156616.
Cf. A080054 (a=2, b=1, j=1), A007096 (a=2, b=1, j=2), A261647 (a=2, b=1, j=3), A014969 (a=2, b=1, j=4), A261648 (a=2, b=1, j=5), A014970 (a=2, b=1, j=6), A014972 (a=2, b=1, j=8), A103261 (a=2, b=1, j=10).
Cf. A261610 (a=3, b=1, j=1), A261649 (a=3, b=1, j=2), A261651 (a=3, b=1, j=3).
Cf. A261611 (a=4, b=1, j=1), A261650 (a=4, b=1, j=2), A261652 (a=4, b=1, j=3).

Programs

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

Formula

a(n) ~ exp(Pi*sqrt(3*n)/2) * 2^(1/4) * Gamma(1/4)^3 / (8 * 3^(1/8) * Pi^(9/4) * n^(3/8)).

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.

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

Views

Author

Ilya Gutkovskiy, Jul 07 2017

Keywords

Examples

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

Crossrefs

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.

Programs

  • Mathematica
    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

Formula

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.
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.