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 12 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

A323433 Number of ways to split an integer partition of n into consecutive subsequences of equal length.

Original entry on oeis.org

1, 1, 3, 5, 10, 14, 25, 34, 54, 74, 109, 146, 211, 276, 381, 501, 675, 871, 1156, 1477, 1926, 2447, 3142, 3957, 5038, 6291, 7918, 9839, 12277, 15148, 18773, 23027, 28333, 34587, 42284, 51357, 62466, 75503, 91344, 109971, 132421, 158755, 190365, 227354, 271511
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Examples

			The a(5) = 14 split partitions:
  [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

  • Maple
    b:= proc(n, i, t) option remember; `if`(n=0 or i=1, numtheory
          [tau](t+n), b(n, i-1, t)+b(n-i, min(n-i, i), t+1))
        end:
    a:= n-> `if`(n=0, 1, b(n$2, 0)):
    seq(a(n), n=0..50);  # Alois P. Heinz, Jan 15 2019
  • Mathematica
    Table[Sum[Length[Divisors[Length[ptn]]],{ptn,IntegerPartitions[n]}],{n,30}]
    (* Second program: *)
    b[n_, i_, t_] := b[n, i, t] = If[n == 0 || i == 1,
         DivisorSigma[0, t+n], b[n, i-1, t] + b[n-i, Min[n-i, i], t+1]];
    a[n_] := If[n == 0, 1, b[n, n, 0]];
    a /@ Range[0, 50] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)
  • PARI
    my(N=66, x='x+O('x^N)); Vec(1+sum(k=1, N, numdiv(k)*x^k/prod(j=1, k, 1-x^j))) \\ Seiichi Manyama, Jan 21 2022
    
  • PARI
    my(N=66, x='x+O('x^N)); Vec(1+sum(i=1, N, sum(j=1, N\i, x^(i*j)/prod(k=1, i*j, 1-x^k)))) \\ Seiichi Manyama, Jan 21 2022

Formula

a(n) = Sum_y A000005(k), where the sum is over all integer partitions of n and k is the number of parts.
From Seiichi Manyama, Jan 21 2022: (Start)
G.f.: 1 + Sum_{k>=1} A000005(k) * x^k/Product_{j=1..k} (1-x^j).
G.f.: 1 + Sum_{i>=1} Sum_{j>=1} x^(i*j)/Product_{k=1..i*j} (1-x^k). (End)
a(n) = Sum_{i=1..n} Sum_{j=1..n} A008284(n,i*j). - Ridouane Oudra, Apr 13 2023

A323300 Number of ways to fill a matrix with the parts of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 4, 1, 2, 2, 4, 1, 6, 1, 4, 4, 3, 1, 6, 1, 6, 4, 4, 1, 12, 2, 4, 2, 6, 1, 12, 1, 2, 4, 4, 4, 18, 1, 4, 4, 12, 1, 12, 1, 6, 6, 4, 1, 10, 2, 6, 4, 6, 1, 12, 4, 12, 4, 4, 1, 36, 1, 4, 6, 4, 4, 12, 1, 6, 4, 12, 1, 20, 1, 4, 6, 6, 4, 12, 1, 10, 3, 4
Offset: 1

Views

Author

Gus Wiseman, Jan 12 2019

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The a(24) = 12 matrices whose entries are (2,1,1,1):
  [1 1 1 2] [1 1 2 1] [1 2 1 1] [2 1 1 1]
.
  [1 1] [1 1] [1 2] [2 1]
  [1 2] [2 1] [1 1] [1 1]
.
  [1] [1] [1] [2]
  [1] [1] [2] [1]
  [1] [2] [1] [1]
  [2] [1] [1] [1]
		

Crossrefs

Positions of 1's are one and prime numbers A008578.
Positions of 2's are primes to prime powers A053810.

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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    Array[Length[ptnmats[#]]&,100]

Formula

a(n) = A008480(n) * A000005(A001222(n)).

A323351 Number of ways to fill a (not necessarily square) matrix with n zeros and ones.

Original entry on oeis.org

1, 2, 8, 16, 48, 64, 256, 256, 1024, 1536, 4096, 4096, 24576, 16384, 65536, 131072, 327680, 262144, 1572864, 1048576, 6291456, 8388608, 16777216, 16777216, 134217728, 100663296, 268435456, 536870912, 1610612736, 1073741824, 8589934592, 4294967296, 25769803776
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Examples

			The a(3) = 16 matrices:
  [000] [001] [010] [011] [100] [101] [110] [111]
.
  [0] [0] [0] [0] [1] [1] [1] [1]
  [0] [0] [1] [1] [0] [0] [1] [1]
  [0] [1] [0] [1] [0] [1] [0] [1]
		

Crossrefs

Programs

  • Mathematica
    Table[2^n*DivisorSigma[0,n],{n,10}]
  • PARI
    a(n) = if (n==0, 1, 2^n*numdiv(n)); \\ Michel Marcus, Jan 15 2019

Formula

a(n) = 2^n * A000005(n) for n > 0, a(0) = 1.
G.f.: 1 + Sum_{k>=1} 2^k*x^k/(1 - 2^k*x^k). - Ilya Gutkovskiy, May 23 2019

A323307 Number of ways to fill a matrix with the parts of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 2, 4, 2, 6, 3, 12, 18, 12, 2, 36, 4, 10, 20, 72, 2, 60, 4, 40, 60, 24, 3, 120, 80, 14, 360, 120, 4, 240, 2, 240, 42, 32, 70, 720, 6, 27, 112, 480, 2, 210, 4, 84, 420, 40, 4, 1440, 280, 280, 108, 224, 5, 1260, 224, 420, 180, 22, 2, 840, 6, 72, 1680, 2880
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

This multiset (row n of A305936) is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

Examples

			The a(22) = 24 matrices:
  [111112] [111121] [111211] [112111] [121111] [211111]
.
  [111] [111] [111] [112] [121] [211]
  [112] [121] [211] [111] [111] [111]
.
  [11] [11] [11] [11] [12] [21]
  [11] [11] [12] [21] [11] [11]
  [12] [21] [11] [11] [11] [11]
.
  [1] [1] [1] [1] [1] [2]
  [1] [1] [1] [1] [2] [1]
  [1] [1] [1] [2] [1] [1]
  [1] [1] [2] [1] [1] [1]
  [1] [2] [1] [1] [1] [1]
  [2] [1] [1] [1] [1] [1]
		

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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    Array[Length[ptnmats[Times@@Prime/@nrmptn[#]]]&,30]

Formula

a(n) = A318762(n) * A000005(A056239(n)).

A323301 Number of ways to fill a matrix with the parts of a strict integer partition of n.

Original entry on oeis.org

1, 1, 1, 5, 5, 9, 21, 25, 37, 53, 137, 153, 249, 337, 505, 845, 1085, 1497, 2061, 2785, 3661, 7589, 8849, 13329, 18033, 26017, 34225, 48773, 70805, 91977, 123765, 164761, 216373, 283205, 367913, 470889, 758793, 913825, 1264105, 1651613, 2251709, 2894793, 3927837
Offset: 0

Views

Author

Gus Wiseman, Jan 12 2019

Keywords

Examples

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

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember;
          `if`(n>i*(i+1)/2, 0, `if`(n=0, t!*numtheory[tau](t),
           b(n, i-1, t)+b(n-i, min(n-i, i-1), t+1)))
        end:
    a:= n-> `if`(n=0, 1, b(n$2, 0)):
    seq(a(n), n=0..50);  # Alois P. Heinz, Jan 15 2019
  • 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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    Table[Sum[Length[ptnmats[k]],{k,Select[Times@@Prime/@#&/@IntegerPartitions[n],SquareFreeQ]}],{n,20}]
    (* Second program: *)
    b[n_, i_, t_] := b[n, i, t] = If[n > i(i+1)/2, 0,
         If[n == 0, t!*DivisorSigma[0, t], b[n, i - 1, t] +
         b[n - i, Min[n - i, i - 1], t + 1]]];
    a[n_] := If[n == 0, 1, b[n, n, 0]];
    a /@ Range[0, 50] (* Jean-François Alcover, May 13 2021, after Alois P. Heinz *)

Formula

a(n) = Sum_{y1 + ... + yk = n, y1 > ... > yk} k! * A000005(k) for n > 0, a(0) = 1.

Extensions

a(0)=1 prepended by Alois P. Heinz, Jan 15 2019

A323349 Number of positive integer matrices with entries summing to n, with equal row-sums and equal column-sums.

Original entry on oeis.org

1, 1, 3, 3, 6, 3, 11, 3, 12, 6, 13, 3, 52, 3, 15, 30, 57, 3, 156, 3, 238, 129, 19, 3, 2221, 6, 21, 415, 3114, 3, 14921, 3, 12853, 1044, 25, 6219, 164743, 3, 27, 2220, 851476, 3, 954088, 3, 434106, 3326714, 31, 3, 24648724, 6, 22309800, 7269, 2737618, 3, 69823653
Offset: 0

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

Also the number of non-normal semi-magic rectangles summing to n with no zeros.
Matrices must be of size m X k where m, k are divisors of n and mk <= n. This implies that a(p) = 3 for p prime, since the only allowable matrices must be of size 1 X 1, 1 X p or p X 1 with only one way to fill in the entries for each matrix size. Similarly, a(p^2) = 6 with additional allowable matrices of sizes 1 X p^2, p^2 X 1 and p X p, again with only one way to fill in the entries for each size. - Chai Wah Wu, Jan 13 2019

Examples

			The a(6) = 11 matrices:
  [6] [3 3] [2 2 2] [1 1 1 1 1 1]
.
  [3] [1 2] [2 1] [1 1 1]
  [3] [2 1] [1 2] [1 1 1]
.
  [2] [1 1]
  [2] [1 1]
  [2] [1 1]
.
  [1]
  [1]
  [1]
  [1]
  [1]
  [1]
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Join@@Table[Partition[cmp,d],{cmp,Join@@Permutations/@IntegerPartitions[n]},{d,Divisors[Length[cmp]]}],And[SameQ@@Total/@#,SameQ@@Total/@Transpose[#]]&]],{n,10}]

Formula

a(p) = 3 and a(p^2) = 6 for p prime (see comment). - Chai Wah Wu, Jan 13 2019

Extensions

a(21)-a(31) from Chai Wah Wu, Jan 13 2019
a(32)-a(53) from Chai Wah Wu, Jan 14 2019
a(54) from Chai Wah Wu, Jan 16 2019

A323434 Number of ways to split a strict integer partition of n into consecutive subsequences of equal length.

Original entry on oeis.org

1, 1, 1, 3, 3, 5, 7, 9, 11, 15, 20, 24, 31, 38, 48, 59, 72, 86, 106, 125, 150, 180, 213, 250, 296, 347, 407, 477, 555, 645, 751, 869, 1003, 1161, 1334, 1534, 1763, 2018, 2306, 2637, 3002, 3418, 3886, 4409, 4994, 5659, 6390, 7214, 8135, 9160, 10300, 11580, 12990
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Examples

			The a(10) = 20 split partitions:
  [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 3]
  [1] [2] [3] [4] [2 1]
.
  [7] [6] [5] [5]
  [2] [3] [4] [3]
  [1] [1] [1] [2]
.
  [4]
  [3]
  [2]
  [1]
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n>i*(i+1)/2, 0,
          `if`(n=0, numtheory[tau](t), b(n, i-1, t)+
             b(n-i, min(n-i, i-1), t+1)))
        end:
    a:= n-> `if`(n=0, 1, b(n$2, 0)):
    seq(a(n), n=0..60);  # Alois P. Heinz, Jan 15 2019
  • Mathematica
    Table[Sum[Length[Divisors[Length[ptn]]],{ptn,Select[IntegerPartitions[n],UnsameQ@@#&]}],{n,30}]
    (* Second program: *)
    b[n_, i_, t_] := b[n, i, t] = If[n>i(i+1)/2, 0,
         If[n == 0, DivisorSigma[0, t], b[n, i-1, t] +
         b[n-i, Min[n-i, i-1], t+1]]];
    a[n_] := If[n == 0, 1, b[n, n, 0]];
    a /@ Range[0, 60] (* Jean-François Alcover, May 18 2021, after Alois P. Heinz *)

Formula

a(n) = Sum_y A000005(k), where the sum is over all strict integer partitions of n and k is the number of parts.

A323305 Number of divisors of the number of prime factors of n counted with multiplicity.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 1, 3, 2, 2, 2, 2, 1, 2, 1, 2, 2, 2, 2, 3, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 3, 1, 2, 2, 4, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 2, 2, 1, 2, 3, 2, 1, 3, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

a(1) = 1 by convention.
First differs from A036430 at a(64) = 4, A036430(64) = 3.

Crossrefs

Positions of 1's are 1 and the prime numbers A008578.
Positions of 2's are A063989.

Programs

  • Mathematica
    Array[Length@*Divisors@*PrimeOmega,100]
  • PARI
    a(n) = if (n==1, 1, numdiv(bigomega(n))); \\ Michel Marcus, Jan 13 2019

Formula

a(n) = A000005(A001222(n)).

A338864 Triangle T(n,k) defined by Sum_{k=1..n} T(n,k)*u^k*x^n/n! = Product_{j>0} ( exp(x^j/(1 - x^j)) )^u.

Original entry on oeis.org

1, 4, 1, 12, 12, 1, 72, 96, 24, 1, 240, 840, 360, 40, 1, 2880, 7200, 4920, 960, 60, 1, 10080, 70560, 65520, 19320, 2100, 84, 1, 161280, 745920, 887040, 362880, 58800, 4032, 112, 1, 1088640, 7983360, 12640320, 6652800, 1481760, 150192, 7056, 144, 1
Offset: 1

Views

Author

Seiichi Manyama, Nov 13 2020

Keywords

Comments

Also the Bell transform of A323295.

Examples

			exp(Sum_{n>0} u*d(n)*x^n) = 1 + u*x + (4*u+u^2)*x^2/2! + (12*u+12*u^2+u^3)*x^3/3! + ... .
Triangle begins:
       1;
       4,      1;
      12,     12,      1;
      72,     96,     24,      1;
     240,    840,    360,     40,     1;
    2880,   7200,   4920,    960,    60,    1;
   10080,  70560,  65520,  19320,  2100,   84,   1;
  161280, 745920, 887040, 362880, 58800, 4032, 112, 1;
  ...
		

Crossrefs

Column k=1..2 give A323295, (n!/2) * A055507(n-1).
Rows sum give A294363.

Programs

  • Mathematica
    T[n_, 0] := Boole[n == 0]; T[n_, k_] := T[n, k] = Sum[Boole[j > 0] * Binomial[n - 1, j - 1] * j! * DivisorSigma[0, j] * T[n - j, k - 1], {j, 0, n - k + 1}]; Table[T[n, k], {n, 1, 9}, {k, 1, n}] // Flatten (* Amiram Eldar, Nov 13 2020 *)
  • PARI
    {T(n, k) = my(u='u); n!*polcoef(polcoef(prod(j=1, n, exp(x^j/(1-x^j+x*O(x^n)))^u), n), k)}
    
  • PARI
    a(n) = if(n<1, 0, n!*numdiv(n));
    T(n, k) = if(k==0, 0^n, sum(j=0, n-k+1, binomial(n-1, j-1)*a(j)*T(n-j, k-1)))

Formula

E.g.f.: exp(Sum_{n>0} u*d(n)*x^n), where d(n) is the number of divisors of n.
T(n; u) = Sum_{k=1..n} T(n,k)*u^k is given by T(n; u) = u * (n-1)! * Sum_{k=1..n} k*d(k)*T(n-k; u)/(n-k)!, T(0; u) = 1.
T(n,k) = (n!/k!) * Sum_{i_1,i_2,...,i_k > 0 and i_1+i_2+...+i_k=n} Product_{j=1..k} d(i_j).
Showing 1-10 of 12 results. Next