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-10 of 17 results. Next

A053529 a(n) = n! * number of partitions of n.

Original entry on oeis.org

1, 1, 4, 18, 120, 840, 7920, 75600, 887040, 10886400, 152409600, 2235340800, 36883123200, 628929100800, 11769069312000, 230150688768000, 4833164464128000, 105639166144512000, 2464913876705280000, 59606099200327680000, 1525429559126753280000, 40464026199993876480000
Offset: 0

Views

Author

N. J. A. Sloane, Jan 16 2000

Keywords

Comments

Commuting permutations: number of ordered pairs (g, h) in Sym(n) such that gh = hg.
Equivalently sum of the order of all normalizers of all cyclic subgroups of Sym(n). - Olivier Gérard, Apr 04 2012
From Gus Wiseman, Jan 16 2019: (Start)
Also the number of Young tableaux with distinct entries from 1 to n, where a Young tableau is an array obtained by replacing the dots in the Ferrers diagram of an integer partition of n with positive integers. For example, the a(3) = 18 tableaux are:
123 213 132 312 231 321
.
12 21 13 31 23 32
3 3 2 2 1 1
.
1 2 1 3 2 3
2 1 3 1 3 2
3 3 2 2 1 1
(End)

References

  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.12, solution.

Crossrefs

Column k=2 of A362827.
Sequences counting pairs of functions from an n-set to itself: A053529, A181162, A239749-A239785, A239836-A239841.

Programs

  • Magma
    a:= func< n | NumberOfPartitions(n)*Factorial(n) >; [ a(n) : n in [0..25]]; // Vincenzo Librandi, Jan 17 2019
    
  • Maple
    seq(count(Permutation(n))*count(Partition(n)),n=1..20); # Zerinvary Lajos, Oct 16 2006
    with(combinat): A053529 := proc(n): n! * numbpart(n) end: seq(A053529(n), n=0..20); # Johannes W. Meijer, Jul 28 2016
  • Mathematica
    Table[PartitionsP[n] n!, {n, 0, 20}] (* T. D. Noe, Jun 19 2012 *)
  • PARI
    N=66; x='x+O('x^N); Vec(serlaplace(exp(sum(k=1, N, x^k/(1-x^k)/k)))) \\ Joerg Arndt, Apr 16 2010
    
  • PARI
    N=66; x='x+O('x^N); Vec(serlaplace(sum(n=0, N, x^n/prod(k=1,n,1-x^k)))) \\ Joerg Arndt, Jan 29 2011
    
  • PARI
    a(n) = n!*numbpart(n); \\ Michel Marcus, Jul 28 2016
    
  • Python
    from math import factorial
    from sympy import npartitions
    def A053529(n): return factorial(n)*npartitions(n) # Chai Wah Wu, Jul 10 2023

Formula

E.g.f: Sum_{n>=0} x^n/(Product_{k=1..n} 1-x^k) = exp(Sum_{n>=1} (x^n/n)/(1-x^n)). - Joerg Arndt, Jan 29 2011
a(n) = Sum{k=1..n} (((n-1)!/(n-k)!)*sigma(k)*a(n-k)), n > 0, and a(0)=1. See A274760. - Johannes W. Meijer, Jul 28 2016
a(n) ~ sqrt(Pi/6)*exp(sqrt(2/3)*Pi*sqrt(n))*n^n/(2*exp(n)*sqrt(n)). - Ilya Gutkovskiy, Jul 28 2016

A319646 Number of non-isomorphic weight-n chains of distinct multisets whose dual is also a chain of distinct multisets.

Original entry on oeis.org

1, 1, 1, 4, 4, 9, 17, 28, 41, 75, 122, 192, 314, 484, 771, 1216, 1861, 2848, 4395, 6610, 10037
Offset: 0

Views

Author

Gus Wiseman, Sep 25 2018

Keywords

Comments

The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
From Gus Wiseman, Jan 17 2019: (Start)
Also the number of plane partitions of n with no repeated rows or columns. For example, the a(6) = 17 plane partitions are:
6 51 42 321
.
5 4 41 31 32 31 22 221 211
1 2 1 2 1 11 2 1 11
.
3 21 21 111
2 2 11 11
1 1 1 1
(End)

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 9 chains:
1: {{1}}
2: {{1,1}}
3: {{1,1,1}}
   {{1,2,2}}
   {{1},{1,1}}
   {{2},{1,2}}
4: {{1,1,1,1}}
   {{1,2,2,2}}
   {{1},{1,1,1}}
   {{2},{1,2,2}}
5: {{1,1,1,1,1}}
   {{1,1,2,2,2}}
   {{1,2,2,2,2}}
   {{1},{1,1,1,1}}
   {{2},{1,1,2,2}}
   {{2},{1,2,2,2}}
   {{1,1},{1,1,1}}
   {{1,2},{1,2,2}}
   {{2,2},{1,2,2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnplane[n_]:=Union[Map[Reverse@*primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Sum[Length[Select[ptnplane[Times@@Prime/@y],And[UnsameQ@@#,UnsameQ@@Transpose[PadRight[#]],And@@GreaterEqual@@@#,And@@(GreaterEqual@@@Transpose[PadRight[#]])]&]],{y,IntegerPartitions[n]}],{n,10}] (* Gus Wiseman, Jan 18 2019 *)

Extensions

a(11)-a(17) from Gus Wiseman, Jan 18 2019
a(18)-a(21) from Robert Price, Jun 21 2021

A003293 Number of planar partitions of n decreasing across rows.

Original entry on oeis.org

1, 1, 2, 4, 7, 12, 21, 34, 56, 90, 143, 223, 348, 532, 811, 1224, 1834, 2725, 4031, 5914, 8638, 12540, 18116, 26035, 37262, 53070, 75292, 106377, 149738, 209980, 293473, 408734, 567484, 785409, 1083817, 1491247, 2046233, 2800125, 3821959, 5203515
Offset: 0

Views

Author

Keywords

Comments

Also number of planar partitions monotonically decreasing down antidiagonals (i.e., with b(n,k) <= b(n-1,k+1)). Transpose (to get planar partitions decreasing down columns), then take the conjugate of each row. - Franklin T. Adams-Watters, May 15 2006
Also number of partitions into one kind of 1's and 2's, two kinds of 3's and 4's, three kinds of 5's and 6's, etc. - Joerg Arndt, May 01 2013
Also count of semistandard Young tableaux with sum of entries equal to n (row sums of A228125). - Wouter Meeussen, Aug 11 2013

Examples

			From _Gus Wiseman_, Jan 17 2019: (Start)
The a(6) = 21 plane partitions with strictly decreasing columns (the count is the same as for strictly decreasing rows):
  6   51   42   411   33   321   3111   222   2211   21111   111111
.
  5   4   41   31   32   311   22   221   2111
  1   2   1    2    1    1     11   1     1
.
  3
  2
  1
(End)
		

References

  • D. M. Bressoud, Proofs and Confirmations, Camb. Univ. Press, 1999; p. 133.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

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-> `if`(modp(n,2)=0,n,n+1)/2): seq(a(n), n=0..45);  # Alois P. Heinz, Sep 08 2008
  • Mathematica
    CoefficientList[Series[Product[(1-x^k)^(-Ceiling[k/2]), {k, 1, 40}], {x, 0, 40}], x][[1 ;; 40]] (* Jean-François Alcover, Apr 18 2011, after Michael Somos *)
    nmax=50; CoefficientList[Series[Product[1/(1-x^k)^((2*k+1-(-1)^k)/4),{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Feb 28 2015 *)
    nmax = 50; CoefficientList[Series[Product[1/((1-x^(2*k-1))*(1-x^(2*k)))^k, {k, 1, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 02 2015 *)
  • PARI
    {a(n)=if(n<0, 0, polcoeff( prod(k=1, n, (1-x^k+x*O(x^n))^-ceil(k/2)), n))} /* Michael Somos, Sep 19 2006 */

Formula

G.f.: Product_(1 - x^k)^{-c(k)}, c(k) = 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ....
Euler transform of A110654. - Michael Somos, Sep 19 2006
a(n) ~ 2^(-3/4) * (3*Pi*Zeta(3))^(-1/2) * (n/Zeta(3))^(-49/72) * exp(3/2*Zeta(3) * (n/Zeta(3))^(2/3) + Pi^2*(n/Zeta(3))^(1/3)/24 - Pi^4/(3456*Zeta(3)) + Zeta'(-1)/2) [Basil Gordon and Lorne Houten, 1969]. - Vaclav Kotesovec, Feb 28 2015

Extensions

More terms from James Sellers, Feb 06 2000
Additional comments from Michael Somos, May 19 2000

A323429 Number of rectangular plane partitions of n.

Original entry on oeis.org

1, 1, 3, 5, 10, 14, 26, 35, 58, 81, 124, 169, 257, 345, 501, 684, 968, 1304, 1830, 2452, 3387, 4541, 6188, 8257, 11193, 14865, 19968, 26481, 35341, 46674, 62007, 81611, 107860, 141602, 186292, 243800, 319610, 416984, 544601, 708690, 922472, 1197018, 1553442
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Comments

Number of ways to fill a (not necessarily square) matrix with the parts of an integer partition of n so that the rows and columns are weakly decreasing.

Examples

			The a(5) = 14 matrices:
  [5] [4 1] [3 2] [3 1 1] [2 2 1] [2 1 1 1] [1 1 1 1 1]
.
  [4] [3] [2 1]
  [1] [2] [1 1]
.
  [3] [2]
  [1] [2]
  [1] [1]
.
  [2]
  [1]
  [1]
  [1]
.
  [1]
  [1]
  [1]
  [1]
  [1]
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Length[Select[Union[Sort/@Tuples[IntegerPartitions[#,{k}]&/@ptn]],And@@OrderedQ/@Transpose[#]&]],{ptn,IntegerPartitions[n]},{k,Min[ptn]}],{n,30}]

A323438 Number of ways to fill a Young diagram with the prime indices of n such that all rows and columns are weakly increasing.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 5, 1, 3, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 4, 1, 7, 2, 2, 2, 8, 1, 2, 2, 7, 1, 4, 1, 4, 4, 2, 1, 12, 2, 3, 2, 4, 1, 5, 2, 7, 2, 2, 1, 10, 1, 2, 4, 11, 2, 4, 1, 4, 2, 4, 1, 13, 1, 2, 3, 4, 2, 4, 1, 12, 5, 2, 1, 10, 2
Offset: 1

Views

Author

Gus Wiseman, Jan 16 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The a(96) = 19 tableaux:
  111112
.
  111   1111   1112   11111   11112
  112   12     11     2       1
.
  11   111   111   112   1111   1112
  11   11    12    11    1      1
  12   2     1     1     2      1
.
  11   11   111   112
  11   12   1     1
  1    1    1     1
  2    1    2     1
.
  11   12
  1    1
  1    1
  1    1
  2    1
.
  1
  1
  1
  1
  1
  2
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnplane[n_]:=Union[Map[primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Length[Select[ptnplane[y],And[And@@LessEqual@@@#,And@@(LessEqual@@@Transpose[PadRight[#]/.(0->Infinity)])]&]],{y,100}]

Formula

Sum_{A056239(n) = k} a(k) = A323450(n).

A323430 Number of rectangular plane partitions of n with strictly decreasing rows and columns.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 7, 9, 12, 16, 22, 27, 36, 44, 57, 72, 89, 110, 139, 170, 210, 261, 318, 390, 478, 581, 705, 860, 1036, 1252, 1511, 1816, 2178, 2618, 3127, 3743, 4471, 5330, 6347, 7564, 8984, 10674, 12669, 15016, 17780, 21050, 24868, 29371, 34655, 40836, 48080
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Comments

Number of ways to fill a (not necessarily square) matrix with the parts of an integer partition of n so that the rows and columns are strictly decreasing.

Examples

			The a(8) = 12 matrices:
  [8] [7 1] [6 2] [5 3] [5 2 1] [4 3 1]
.
  [7] [6] [5] [3 2]
  [1] [2] [3] [2 1]
.
  [5] [4]
  [2] [3]
  [1] [1]
The a(10) = 22 matrices:
  [10] [9 1] [8 2] [7 3] [7 2 1] [6 4] [6 3 1] [5 4 1] [5 3 2] [4 3 2 1]
.
  [9] [8] [7] [6] [5 2] [4 2] [4 3]
  [1] [2] [3] [4] [2 1] [3 1] [2 1]
.
  [7] [6] [5] [5]
  [2] [3] [4] [3]
  [1] [1] [1] [2]
.
  [4]
  [3]
  [2]
  [1]
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Length[Select[Union[Tuples[Select[IntegerPartitions[#,{k}],UnsameQ@@#&]&/@ptn]],And@@(OrderedQ[#,Greater]&/@Transpose[#])&]],{ptn,IntegerPartitions[n]},{k,Min[ptn]}],{n,30}]

A323432 Number of semistandard rectangular plane partitions of n.

Original entry on oeis.org

1, 1, 2, 4, 6, 9, 15, 20, 30, 42, 59, 79, 112, 146, 199, 264, 350, 455, 603, 774, 1010, 1297, 1668, 2124, 2724, 3441, 4372, 5513, 6955, 8718, 10960, 13670, 17091, 21264, 26454, 32786, 40667, 50215, 62048, 76435, 94126
Offset: 0

Views

Author

Gus Wiseman, Jan 16 2019

Keywords

Comments

Number of ways to fill a (not necessarily square) matrix with the parts of an integer partition of n so that the rows are weakly decreasing and the columns are strictly decreasing.

Examples

			The a(6) = 15 matrices:
  [6] [51] [42] [411] [33] [321] [3111] [222] [2211] [21111] [111111]
.
  [5] [4] [22]
  [1] [2] [11]
.
  [3]
  [2]
  [1]
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Length[Select[Union[Tuples[IntegerPartitions[#,{k}]&/@ptn]],And@@(OrderedQ[#,Greater]&/@Transpose[#])&]],{ptn,IntegerPartitions[n]},{k,Min[ptn]}],{n,30}]

A323437 Number of semistandard Young tableaux whose entries are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 4, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 4, 1, 2, 2, 4, 1, 3, 1, 2, 2, 2, 2, 4, 1, 2, 1, 2, 1, 5, 2, 2, 2
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Comments

Number of ways to fill a Young diagram with the prime indices of n such that all rows are weakly increasing and all columns are strictly increasing.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Is this a duplicate of A339887? - R. J. Mathar, Feb 03 2021

Examples

			The a(60) = 5 tableaux:
  1123
.
  11   112   113
  23   3     2
.
  11
  2
  3
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnplane[n_]:=Union[Map[primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Length[Select[ptnplane[y],And[And@@Less@@@#,And@@(LessEqual@@@Transpose[PadRight[#]/.(0->Infinity)])]&]],{y,100}]

Formula

Sum_{A056239(n) = k} a(k) = A003293(n).

A323439 Number of ways to fill a Young diagram with the prime indices of n such that all rows and columns are strictly increasing.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 0, 1, 2, 2, 0, 1, 1, 1, 0, 2, 2, 1, 0, 0, 2, 0, 0, 1, 4, 1, 0, 2, 2, 2, 0, 1, 2, 2, 0, 1, 4, 1, 0, 0, 2, 1, 0, 0, 1, 2, 0, 1, 0, 2, 0, 2, 2, 1, 0, 1, 2, 0, 0, 2, 4, 1, 0, 2, 4, 1, 0, 1, 2, 1, 0, 2, 4, 1, 0, 0, 2, 1, 0, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jan 16 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			The a(630) = 8 tableaux:
  123   124   1234
  24    23    2
.
  12   12   123   124
  23   24   2     2
  4    3    4     3
.
  12
  2
  3
  4
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnplane[n_]:=Union[Map[primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Length[Select[ptnplane[y],And[And@@Less@@@#,And@@(Less@@@DeleteCases[Transpose[PadRight[#]],0,{2}]),And@@(LessEqual@@@Transpose[PadRight[#]/.(0->Infinity)])]&]],{y,100}]

Formula

Sum_{A056239(n) = k} a(k) = A323451(n).

A323431 Number of strict rectangular plane partitions of n.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 7, 9, 11, 15, 21, 25, 33, 41, 53, 65, 81, 97, 121, 143, 173, 215, 255, 305, 367, 441, 527, 637, 751, 899, 1067, 1269, 1491, 1775, 2071, 2439, 2875, 3357, 3911, 4577, 5309, 6177, 7171, 8305, 9609, 11151
Offset: 0

Views

Author

Gus Wiseman, Jan 16 2019

Keywords

Comments

Number of ways to fill a (not necessarily square) matrix with the parts of a strict integer partition of n so that the rows and columns are strictly decreasing.

Examples

			The a(10) = 21 matrices:
  [10] [9 1] [8 2] [7 3] [7 2 1] [6 4] [6 3 1] [5 4 1] [5 3 2] [4 3 2 1]
.
  [9] [8] [7] [6] [4 2] [4 3]
  [1] [2] [3] [4] [3 1] [2 1]
.
  [7] [6] [5] [5]
  [2] [3] [4] [3]
  [1] [1] [1] [2]
.
  [4]
  [3]
  [2]
  [1]
		

Crossrefs

Programs

  • Mathematica
    Table[Sum[Length[Select[Union[Sort/@Tuples[IntegerPartitions[#,{k}]&/@ptn]],UnsameQ@@Join@@#&&And@@OrderedQ/@Transpose[#]&]],{ptn,IntegerPartitions[n]},{k,Min[ptn]}],{n,30}]
Showing 1-10 of 17 results. Next