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

A047967 Number of partitions of n with some part repeated.

Original entry on oeis.org

0, 0, 1, 1, 3, 4, 7, 10, 16, 22, 32, 44, 62, 83, 113, 149, 199, 259, 339, 436, 563, 716, 913, 1151, 1453, 1816, 2271, 2818, 3496, 4309, 5308, 6502, 7959, 9695, 11798, 14298, 17309, 20877, 25151, 30203, 36225, 43323, 51748, 61651, 73359, 87086, 103254, 122164
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Also number of partitions of n with at least one even part. - Vladeta Jovovic, Sep 10 2003. Example: a(5)=4 because we have [4,1], [3,2], [2,2,1] and [2,1,1,1] ([5], [3,1,1] and [1,1,1,1,1] do not qualify). - Emeric Deutsch, Mar 30 2006
Also number of partitions of n (where it is assumed that the least part is 0) such that at least one difference is at least two. Example: a(5)=4 because we have [5,0], [4,1,0], [3,2,0] and [3,1,1,0] ([2,2,1,0], [2,1,1,1,0] and [1,1,1,1,1,0] do not qualify). - Emeric Deutsch, Mar 30 2006
The Heinz numbers of these partitions (with some part repeated) are given by A013929. Equivalent to Vladeta Jovovic's comment, a(n) is also the number of integer partitions whose product of parts is even. The Heinz numbers of these latter partitions are given by A324929. - Gus Wiseman, Mar 23 2019

Examples

			a(5) = 4 because we have [3,1,1], [2,2,1], [2,1,1,1] and [1,1,1,1,1] ([5], [4,1] and [3,2] do not qualify).

Links

Crossrefs

Cf. A038348, A261982.
Column k=1 of A320264.
Cf. A324847, A324929, A324966, A324967.

Programs

  • Maple
    g:=sum(x^(2*k)*product(1+x^j,j=k+1..70)/product(1-x^j,j=1..k),k=1..40): gser:=series(g,x=0,50): seq(coeff(gser,x,n),n=0..44); # Emeric Deutsch, Mar 30 2006
  • Mathematica
    Table[PartitionsP[n]-PartitionsQ[n],{n,0,50}] (* Harvey P. Dale, Jan 17 2019 *)
  • PARI
    x='x+O('x^66); concat([0,0], Vec(1/eta(x)-eta(x^2)/eta(x))) \\ Joerg Arndt, Jun 21 2011

Formula

a(n) = A000041(n) - A000009(n).
G.f.: Sum_{k>=1} x^(2*k)*(Product_{j>=k+1} (1+x^j)) / Product_{j=1..k} (1-x^j) = Sum_{k>=1} x^(2*k)/(Product_{j=1..2*k} (1-x^j)*Product_{j>=k} (1-x^(2*j+1))). - Emeric Deutsch, Mar 30 2006
G.f.: 1/P(x) - P(x^2)/P(x) where P(x) = Product_{k>=1} (1-x^k). - Joerg Arndt, Jun 21 2011
a(n) = p(n-2)+p(n-4)-p(n-10)-p(n-14)+...+(-1^(j-1))*p(n-j*(3*j-1)) + (-1^(j-1))*p(n-j*(3*j+1))+..., where p(n) = A000041(n). - Gregory L. Simay, Aug 28 2023

A006152 Exponential generating function x*exp(x/(1-x)).

Original entry on oeis.org

1, 2, 9, 52, 365, 3006, 28357, 301064, 3549177, 45965530, 648352001, 9888877692, 162112109029, 2841669616982, 53025262866045, 1049180850990736, 21937381717388657, 483239096122434354, 11184035897992673017, 271287473871771163460, 6881656485607798743261
Offset: 1

Views

Author

Simon Plouffe

Keywords

Comments

a(n) is the number of labeled rooted trees with every non-root vertex of degree 1 or 2. - Geoffrey Critzer, May 21 2012.
Total number of unit length lists in all sets of lists, cf. A000262. - Alois P. Heinz, May 10 2016

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Links

Crossrefs

Cf. A000262, A320264.

Programs

  • Mathematica
    nn = 17; a = x/(1 - x);
Range[0, nn]! CoefficientList[Series[x Exp[a], {x, 0, nn}], x]  (* Geoffrey Critzer, May 21 2012 *)
  • PARI
    a(n)=n!*polcoeff(x*exp(x/(1-x)+O(x^n)), n)

Formula

a(n) = n*A000262(n-1).
D-finite with recurrence a(n) = 2*(n-1)*a(n-1)-(n^2-5*n+5)*a(n-2)-(n-4)*(n-2)*a(n-3). - Vaclav Kotesovec, Oct 05 2012
a(n) ~ n^(n-1/4)*exp^(2*sqrt(n)-n-1/2)/sqrt(2). - Vaclav Kotesovec, Oct 05 2012
a(n) = A320264(n+1,n). - Alois P. Heinz, Oct 08 2018

Extensions

More terms from Michael Somos, Jun 07 2000

A320265 Number of proper multisets of nonempty words with a total of n letters over n-ary alphabet such that if a letter occurs in the multiset all predecessors occur at least once.

Original entry on oeis.org

1, 3, 23, 178, 1786, 20927, 282520, 4299263, 72750927, 1353700567, 27452623890, 602326265519, 14209892886819, 358576428141962, 9634718410829852, 274567642777650028, 8270000441627265937, 262464788618069324640, 8752908129221863491691, 305968679117675345995513
Offset: 2

Views

Author

Alois P. Heinz, Oct 08 2018

Keywords

Examples

			a(2) = 1: {a,a}.
a(3) = 3: {a,a,a}, {a,a,b}, {a,b,b}.
a(4) = 23: {a,a,a,a}, {a,a,aa}, {aa, aa}, {a,a,a,b}, {a,a,b,b}, {a,b,b,b}, {a,a,ab}, {a,a,ba}, {a,a,bb}, {b,b,ab}, {b,b,ba}, {b,b,aa}, {ab,ab}, {ba,ba}, {a,a,b,c}, {a,a,bc}, {a,a,cb}, {b,b,a,c}, {b,b,ac}, {b,b,ca}, {c,c,a,b}, {c,c,ab}, {c,c,ba}.

Links

Crossrefs

Row sums of A320264.
Cf. A257741, A319518.

Programs

  • Maple
    h:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(h(n-i*j, i-1, k)*binomial(k^i, j), j=0..n/i)))
    end:
g:= proc(n, k) option remember; `if`(n=0, 1, add(add(
      d*k^d, d=numtheory[divisors](j))*g(n-j, k), j=1..n)/n)
    end:
a:= n-> add(add((-1)^i*(g(n, k-i)-h(n$2, k-i))*
        binomial(k, i), i=0..k), k=1..n-1):
seq(a(n), n=2..25);
  • Mathematica
    h[n_, i_, k_] := h[n, i, k] = If[n == 0, 1, If[i < 1, 0, Sum[h[n - i*j, i - 1, k]*Binomial[k^i, j], {j, 0, n/i}]]];
g[n_, k_] := g[n, k] = If[n == 0, 1, Sum[Sum[d*k^d, {d, Divisors[j]}]*g[n - j, k], {j, 1, n}]/n];
T[n_, k_] := Sum[(-1)^i*(g[n, k-i]-h[n, n, k-i])*Binomial[k, i], {i, 0, k}];
a[n_] := Sum[T[n, k], {k, 1, n - 1}];
Table[a[n], {n, 2, 25}] (* Jean-François Alcover, Jun 01 2022, after Alois P. Heinz in A320264 *)

Formula

a(n) = Sum_{k=1..n-1} A320264(n,k).
a(n) = A257741(n) - A319518(n).
Showing 1-3 of 3 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.