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

A281425 a(n) = [q^n] (1 - q)^n / Product_{j=1..n} (1 - q^j).

Original entry on oeis.org

1, 0, 1, -1, 2, -4, 9, -21, 49, -112, 249, -539, 1143, -2396, 5013, -10550, 22420, -48086, 103703, -223806, 481388, -1029507, 2187944, -4625058, 9742223, -20490753, 43111808, -90840465, 191773014, -405523635, 858378825, -1817304609, 3845492204, -8129023694, 17162802918, -36191083386
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 05 2017

Keywords

Comments

a(n) is n-th term of the Euler transform of -n + 1, 1, 1, 1, ...
Inverse zero-based binomial transform of A000041. The version for strict partitions is A380412, or A293467 up to sign. - Gus Wiseman, Feb 06 2025

Crossrefs

Programs

  • Maple
    b:= proc(n, k) option remember; `if`(k=0,
          combinat[numbpart](n), b(n, k-1)-b(n-1, k-1))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..35);  # Alois P. Heinz, Dec 21 2024
  • Mathematica
    Table[SeriesCoefficient[(1 - q)^n / Product[(1 - q^j), {j, 1, n}], {q, 0, n}], {n, 0, 35}]
    Table[SeriesCoefficient[(1 - q)^n QPochhammer[q^(1 + n), q]/QPochhammer[q, q], {q, 0, n}], {n, 0, 35}]
    Table[SeriesCoefficient[1/QFactorial[n, q], {q, 0, n}], {n, 0, 35}]
    Table[Differences[PartitionsP[Range[0, n]], n], {n, 0, 35}] // Flatten
    Table[Sum[(-1)^j*Binomial[n, j]*PartitionsP[n-j], {j, 0, n}], {n, 0, 30}] (* Vaclav Kotesovec, Oct 06 2017 *)

Formula

a(n) = [q^n] 1/((1 + q)*(1 + q + q^2)*...*(1 + q + ... + q^(n-1))).
a(n) = Sum_{j=0..n} (-1)^j * binomial(n, j) * A000041(n-j). - Vaclav Kotesovec, Oct 06 2017
a(n) ~ (-1)^n * 2^(n - 3/2) * exp(Pi*sqrt(n/12) + Pi^2/96) / (sqrt(3)*n). - Vaclav Kotesovec, May 07 2018

A292463 Number of partitions of n with n kinds of 1.

Original entry on oeis.org

1, 1, 4, 14, 51, 188, 702, 2644, 10026, 38223, 146359, 562456, 2168134, 8379539, 32459199, 125984039, 489837300, 1907490728, 7438346255, 29042470132, 113522618066, 444199913556, 1739735079466, 6819657196928, 26753893533257, 105034060120469, 412637434996367
Offset: 0

Views

Author

Alois P. Heinz, Sep 16 2017

Keywords

Examples

			a(2) = 4: 2, 1a1a, 1a1b, 1b1b.
		

Crossrefs

Main diagonal of A292508.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0 or i<2,
          binomial(k+n-1, n), add(b(n-i*j, i-1, k), j=0..n/i))
        end:
    a:= n-> b(n$3):
    seq(a(n), n=0..30);
    # second Maple program:
    b:= proc(n, k) option remember; `if`(n=0, 1, add(
          (numtheory[sigma](j)+k-1)*b(n-j, k), j=1..n)/n)
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..30);
    # third Maple program:
    b:= proc(n, k) option remember; `if`(n=0, 1, `if`(k=1,
          combinat[numbpart](n), b(n-1, k) +b(n, k-1)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..30);
  • Mathematica
    Table[SeriesCoefficient[1/(1-x)^(n-1) * Product[1/(1-x^k), {k,1,n}], {x,0,n}], {n,0,30}] (* Vaclav Kotesovec, Sep 19 2017 *)

Formula

a(n) = [x^n] 1/(1-x)^n * 1/Product_{j=2..n} (1-x^j).
a(n) is n-th term of the Euler transform of n,1,1,1,... .
a(n) ~ c * 4^n / sqrt(n), where c = QPochhammer[-1, 1/2] / (8*sqrt(Pi) * QPochhammer[1/4, 1/4]) = 0.48841139329043831428669851139824427133317... - Vaclav Kotesovec, Sep 19 2017
Equivalently, c = 1/(4*sqrt(Pi)*QPochhammer(1/2)). - Vaclav Kotesovec, Mar 17 2024
Showing 1-2 of 2 results.