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.

Previous Showing 11-20 of 46 results. Next

A321722 Number of non-normal magic squares whose entries are nonnegative integers summing to n.

Original entry on oeis.org

1, 1, 1, 1, 10, 21, 97, 657, 5618, 48918, 494530, 5383553, 65112565, 840566081, 11834555867, 176621056393, 2838064404989, 48060623405313
Offset: 0

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

A non-normal magic square is a square matrix with row sums, column sums, and both diagonals all equal to d, for some d|n.

Examples

			The a(4) = 10 magic squares:
  [4]
.
  [1 1]
  [1 1]
.
  [1 0 0 0][1 0 0 0][0 1 0 0][0 1 0 0][0 0 1 0][0 0 1 0][0 0 0 1][0 0 0 1]
  [0 0 1 0][0 0 0 1][0 0 1 0][0 0 0 1][1 0 0 0][0 1 0 0][1 0 0 0][0 1 0 0]
  [0 0 0 1][0 1 0 0][1 0 0 0][0 0 1 0][0 1 0 0][0 0 0 1][0 0 1 0][1 0 0 0]
  [0 1 0 0][0 0 1 0][0 0 0 1][1 0 0 0][0 0 0 1][1 0 0 0][0 1 0 0][0 0 1 0]
		

Crossrefs

Programs

  • Mathematica
    prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
    multsubs[set_,k_]:=If[k==0,{{}},Join@@Table[Prepend[#,set[[i]]]&/@multsubs[Drop[set,i-1],k-1],{i,Length[set]}]];
    Table[Length[Select[multsubs[Tuples[Range[n],2],n],And[Union[First/@#]==Range[Max@@First/@#]==Union[Last/@#],SameQ@@Join[{Tr[prs2mat[#]],Tr[Reverse[prs2mat[#]]]},Total/@prs2mat[#],Total/@Transpose[prs2mat[#]]]]&]],{n,5}]

Formula

a(p) = A007016(p) + 1 if p is prime. a(n) >= A007016(n) + 1 for n > 1. - Chai Wah Wu, Jan 15 2019

Extensions

a(7)-a(15) from Chai Wah Wu, Jan 15 2019
a(16)-a(17) from Chai Wah Wu, Jan 16 2019

A323306 Heinz numbers of integer partitions that can be arranged into a matrix with equal row-sums and equal column-sums.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 36, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 100, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193
Offset: 1

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

First differs from A137944 in lacking 120.
The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			6480 belongs to the sequence because it is the Heinz number of (3,2,2,2,2,1,1,1,1), which can be arranged in the following ways:
  [1 1 3] [1 2 2] [1 2 2] [1 3 1] [2 1 2] [2 1 2] [2 2 1] [2 2 1] [3 1 1]
  [2 2 1] [1 2 2] [3 1 1] [2 1 2] [1 3 1] [2 1 2] [1 1 3] [2 2 1] [1 2 2]
  [2 2 1] [3 1 1] [1 2 2] [2 1 2] [2 1 2] [1 3 1] [2 2 1] [1 1 3] [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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    Select[Range[100],!Select[ptnmats[#],And[SameQ@@Total/@#,SameQ@@Total/@Transpose[#]]&]=={}&]

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

A321729 Number of integer partitions of n whose Young diagram can be partitioned into vertical sections of the same sizes as the parts of the original partition.

Original entry on oeis.org

1, 1, 1, 2, 3, 4, 6, 8, 12, 16, 22, 28, 40, 51
Offset: 0

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

First differs from A046682 at a(11) = 28, A046682(11) = 29.
A vertical section is a partial Young diagram with at most one square in each row. For example, a suitable partition (shown as a coloring by positive integers) of the Young diagram of (322) is:
1 2 3
1 2
2 3
Conjecture: a(n) is the number of half-loop-graphical partitions of n. An integer partition is half-loop-graphical if it comprises the multiset of vertex-degrees of some graph with half-loops, where a half-loop is an edge with one vertex, to be distinguished from a full loop, which has two equal vertices.

Examples

			The a(1) = 1 through a(8) = 12 partitions whose Young diagram cannot be partitioned into vertical sections of the same sizes as the parts of the original partition are the same as the half-loop-graphical partitions up to n = 8:
  (1)  (11)  (21)   (22)    (221)    (222)     (322)      (332)
             (111)  (211)   (311)    (321)     (2221)     (2222)
                    (1111)  (2111)   (2211)    (3211)     (3221)
                            (11111)  (3111)    (4111)     (3311)
                                     (21111)   (22111)    (4211)
                                     (111111)  (31111)    (22211)
                                               (211111)   (32111)
                                               (1111111)  (41111)
                                                          (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
For example, the half-loop-graphs
  {{1},{1,2},{1,3},{2,3}}
  {{1},{2},{3},{1,2},{1,3}}
both have degrees y = (3,2,2), so y is counted under a(7).
		

Crossrefs

The complement is counted by A321728.
The following pertain to the conjecture.
Half-loop-graphical partitions by length are A029889 or A339843 (covering).
The version for full loops is A339656.
A027187 counts partitions of even length, ranked by A028260.
A058696 counts partitions of even numbers, ranked by A300061.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs, ranked by A340018/A340019.
A339659 is a triangle counting graphical partitions by length.

Programs

  • Mathematica
    spsu[,{}]:={{}};spsu[foo,set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@spsu[Select[foo,Complement[#,Complement[set,s]]=={}&],Complement[set,s]]]/@Cases[foo,{i,_}];
    ptnpos[y_]:=Position[Table[1,{#}]&/@y,1];
    ptnverts[y_]:=Select[Join@@Table[Subsets[ptnpos[y],{k}],{k,Reverse[Union[y]]}],UnsameQ@@First/@#&];
    Table[Length[Select[IntegerPartitions[n],Length[Select[spsu[ptnverts[#],ptnpos[#]],Function[p,Sort[Length/@p]==Sort[#]]]]>0&]],{n,8}]

Formula

a(n) is the number of integer partitions y of n such that the coefficient of m(y) in e(y) is nonzero, where m is monomial symmetric functions and e is elementary symmetric functions.
a(n) = A000041(n) - A321728(n).

A323774 Number of multiset partitions, whose parts are constant and all have the same sum, of integer partitions of n.

Original entry on oeis.org

1, 1, 3, 3, 7, 3, 12, 3, 16, 8, 14, 3, 39, 3, 16, 15, 40, 3, 50, 3, 54, 17, 20, 3, 135, 10, 22, 25, 73, 3, 129, 3, 119, 21, 26, 19, 273, 3, 28, 23, 217, 3, 203, 3, 123, 74, 32, 3, 590, 12, 106, 27, 154, 3, 370, 23, 343, 29, 38, 3, 963, 3, 40, 95, 450, 25, 467, 3
Offset: 0

Views

Author

Gus Wiseman, Jan 27 2019

Keywords

Comments

An unlabeled version of A279789.

Examples

			The a(1) = 1 through a(6) = 12 multiset partitions:
  (1)  (2)     (3)        (4)           (5)              (6)
       (11)    (111)      (22)          (11111)          (33)
       (1)(1)  (1)(1)(1)  (1111)        (1)(1)(1)(1)(1)  (222)
                          (2)(2)                         (3)(3)
                          (2)(11)                        (111111)
                          (11)(11)                       (3)(111)
                          (1)(1)(1)(1)                   (2)(2)(2)
                                                         (111)(111)
                                                         (2)(2)(11)
                                                         (2)(11)(11)
                                                         (11)(11)(11)
                                                         (1)(1)(1)(1)(1)(1)
		

Crossrefs

Cf. A001970, A006171 (constant parts), A007716, A034729, A047966 (uniform partitions), A047968, A279787, A279789 (twice-partitions version), A305551 (equal part-sums), A306017, A319056, A323766, A323775, A323776.

Programs

  • Mathematica
    Table[Length[Join@@Table[Union[Sort/@Tuples[Select[IntegerPartitions[#],SameQ@@#&]&/@ptn]],{ptn,Select[IntegerPartitions[n],SameQ@@#&]}]],{n,30}]
  • PARI
    a(n) = if (n==0, 1, sumdiv(n, d, binomial(numdiv(d) + n/d - 1, n/d))); \\ Michel Marcus, Jan 28 2019

Formula

a(0) = 1; a(n) = Sum_{d|n} binomial(tau(d) + n/d - 1, n/d), where tau = A000005.

A319071 Number of integer partitions of n whose product of parts is a perfect power and whose parts all have the same number of prime factors, counted with multiplicity.

Original entry on oeis.org

1, 0, 0, 0, 2, 0, 2, 0, 3, 2, 3, 0, 4, 1, 4, 3, 7, 1, 7, 1, 8, 6, 8, 0, 15, 5, 12, 6, 15, 4, 22, 4, 24, 12, 22, 8, 35, 7, 30, 16, 42, 9, 50, 9, 50, 30, 53, 7, 79, 22, 72, 33, 87, 21, 109, 26, 111, 55, 117, 24, 168, 40, 149, 65, 178, 59
Offset: 0

Views

Author

Gus Wiseman, Oct 10 2018

Keywords

Comments

The positions of zeros appear to be A048278.

Examples

			The a(4) = 2 through a(16) = 7 integer partitions (G = 16):
  4   33   8     9    55     66      94  77       555     G
  22  222  44    333  3322   444         5522     33333   88
           2222       22222  3333        332222   333222  664
                             222222      2222222          4444
                                                          5533
                                                          333322
                                                          22222222
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[GCD@@FactorInteger[Times@@#][[All,2]]>1,SameQ@@PrimeOmega/@#]&]],{n,30}]

A321720 Number of non-normal (0,1) semi-magic squares with sum of entries equal to n.

Original entry on oeis.org

1, 1, 2, 6, 25, 120, 726, 5040, 40410, 362881, 3630840, 39916800, 479069574, 6227020800, 87181402140, 1307674370040, 20922977418841, 355687428096000, 6402388104196400, 121645100408832000, 2432903379962038320, 51090942171778378800, 1124000886592995642000, 25852016738884976640000
Offset: 0

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

A non-normal semi-magic square is a nonnegative integer matrix with row sums and column sums all equal to d, for some d|n.

Crossrefs

Programs

  • Mathematica
    prs2mat[prs_]:=Table[Count[prs,{i,j}],{i,Union[First/@prs]},{j,Union[Last/@prs]}];
    Table[Length[Select[Subsets[Tuples[Range[n],2],{n}],And[Union[First/@#]==Union[Last/@#]==Range[Max@@First/@#],SameQ@@Total/@prs2mat[#],SameQ@@Total/@Transpose[prs2mat[#]]]&]],{n,5}]

Formula

a(p) = p! for p prime as the squares are all permutation matrices of order p and a(n) >= n! for n > 1 (see comments in A321717 and A321719). - Chai Wah Wu, Jan 13 2019
a(n) = Sum_{d|n, d<=n/d} A008300(n/d, d) for n > 0. - Andrew Howroyd, Apr 11 2020

Extensions

a(7) from Chai Wah Wu, Jan 13 2019
a(8)-a(15) from Chai Wah Wu, Jan 14 2019
a(16)-a(21) from Chai Wah Wu, Jan 16 2019
Terms a(22) and beyond from Andrew Howroyd, Apr 11 2020

A321724 Irregular triangle read by rows where T(n,k) is the number of non-isomorphic non-normal semi-magic square multiset partitions of weight n and length d = A027750(n, k).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 3, 5, 1, 1, 5, 1, 1, 3, 7, 1, 1, 1, 1, 4, 9, 12, 11, 1, 1, 1, 1, 4, 15, 1, 1, 13, 31, 1, 1, 5, 43, 22, 1, 1, 1, 1, 5, 22, 103, 30, 1, 1, 1, 1, 6, 106, 264, 42, 1, 1, 30, 383, 1, 1, 6, 56, 1, 1, 1, 1, 7, 45, 321, 2804, 1731, 77, 1
Offset: 1

Views

Author

Gus Wiseman, Nov 18 2018

Keywords

Comments

Also the number of nonnegative integer square matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, with row sums and column sums all equal to d.
A non-normal semi-magic square multiset partition of weight n is a multiset partition of weight n whose part sizes and vertex degrees are all equal to d, for some d|n.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Triangle begins:
  1
  1 1
  1 1
  1 2 1
  1 1
  1 2 3 1
  1 1
  1 3 5 1
  1 5 1
  1 3 7 1
Inequivalent representatives of the T(10,3) = 7 semi-magic squares (zeros not shown):
  [2    ] [2    ] [2    ] [2    ] [2    ] [11   ] [11   ]
  [ 2   ] [ 2   ] [ 2   ] [ 11  ] [ 11  ] [11   ] [1 1  ]
  [  2  ] [  2  ] [  11 ] [ 11  ] [ 1 1 ] [  11 ] [ 1 1 ]
  [   2 ] [   11] [  1 1] [   11] [  1 1] [  1 1] [  1 1]
  [    2] [   11] [   11] [   11] [   11] [   11] [   11]
		

Crossrefs

Formula

T(n,k) = A333733(d, n/d), where d = A027750(n, k). - Andrew Howroyd, Apr 11 2020

Extensions

a(28)-a(39) from Chai Wah Wu, Jan 16 2019
Terms a(40) and beyond from Andrew Howroyd, Apr 11 2020
Edited by Peter Munn, Mar 05 2025

A323347 Number of integer partitions of n whose parts can be arranged into a (not necessarily square) matrix with equal row-sums and equal column-sums.

Original entry on oeis.org

1, 1, 2, 2, 3, 2, 5, 2, 5, 3, 6, 2, 11, 2, 7, 7, 10, 2, 18, 2, 17, 13, 9, 2, 50, 3, 10, 24, 34, 2, 85, 2, 51, 46, 12, 9, 261, 2, 13, 80, 257, 2, 258, 2, 323, 431, 15, 2, 1533, 3, 227, 206, 1165, 2, 971, 483, 2409, 309, 18, 2
Offset: 0

Views

Author

Gus Wiseman, Jan 13 2019

Keywords

Comments

Rectangles must be of size m X k where m, k are divisors of n and mk <= n. This implies that a(p) = 2 for p prime, since the only allowable rectangles must be of size 1 X 1 corresponding to the partition (p), or 1 X p or p X 1 corresponding to the partition (1,1,...,1). Similarly, a(p^2) = 3 since the allowable rectangles must be of sizes 1 X 1 (partition (p^2)), 1 X p or p X 1 (partition (p,p,...,p)), 1 X p^2, p^2 X 1 and p X p (partition (1,1,...,1)). - Chai Wah Wu, Jan 14 2019

Examples

			The a(8) = 5 integer partitions are (8), (44), (2222), (3311), (11111111).
The a(12) = 11 integer partitions (C = 12):
  (C)
  (66)
  (444)
  (3333)
  (4422)
  (5511)
  (222222)
  (332211)
  (22221111)
  (222111111)
  (111111111111)
For example, the arrangements of (222111111) are:
  [1 1 2] [1 1 2] [1 2 1] [1 2 1] [2 1 1] [2 1 1]
  [1 2 1] [2 1 1] [1 1 2] [2 1 1] [1 1 2] [1 2 1]
  [2 1 1] [1 2 1] [2 1 1] [1 1 2] [1 2 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]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    Table[Length[Select[IntegerPartitions[n],!Select[ptnmats[Times@@Prime/@#],And[SameQ@@Total/@#,SameQ@@Total/@Transpose[#]]&]=={}&]],{n,10}]

Formula

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

Extensions

a(17)-a(53) from Chai Wah Wu, Jan 15 2019
a(54)-a(59) from Chai Wah Wu, Jan 16 2019

A320451 Number of multiset partitions of uniform integer partitions of n in which all parts have the same length.

Original entry on oeis.org

1, 1, 3, 5, 8, 7, 19, 11, 24, 26, 38, 28, 85, 46, 89, 99, 146, 110, 246, 163, 326, 305, 416, 376, 816, 591, 903, 971, 1450, 1295, 2517, 1916, 3045, 3141, 4042, 4117, 7073, 5736, 8131, 9026, 12658, 11514, 19459, 16230, 24638, 27129, 33747, 32279, 55778, 45761, 71946
Offset: 0

Views

Author

Gus Wiseman, Oct 12 2018

Keywords

Comments

An integer partitions is uniform if all parts appear with the same multiplicity.
Terms can be computed by the formula: Sum_{d|n} Sum_{i>=1} P(n/d,i) * Sum_{h|i*d} M(i*d/h, i, h, d) where P(n,k) is the number of partitions of n into k distinct parts and M(h,w,r,s) is the number of nonnegative integer h X w matrices up to row permutations with all row sums equal to r and all column sums equal to s. The cases of M(h,w,w,h) and M(n,n,k,k) are enumerated by the arrays A257462 and A257463. - Andrew Howroyd, Feb 04 2022

Examples

			The a(9) = 26 multiset partitions:
  {{9}}
  {{1,8}}
  {{2,7}}
  {{3,6}}
  {{4,5}}
  {{1,2,6}}
  {{1,3,5}}
  {{1},{8}}
  {{2,3,4}}
  {{2},{7}}
  {{3,3,3}}
  {{3},{6}}
  {{4},{5}}
  {{1},{2},{6}}
  {{1},{3},{5}}
  {{2},{3},{4}}
  {{3},{3},{3}}
  {{1,1,1,2,2,2}}
  {{1,1,1},{2,2,2}}
  {{1,1,2},{1,2,2}}
  {{1,1},{1,2},{2,2}}
  {{1,2},{1,2},{1,2}}
  {{1,1,1,1,1,1,1,1,1}}
  {{1,1,1},{1,1,1},{1,1,1}}
  {{1},{1},{1},{2},{2},{2}}
  {{1},{1},{1},{1},{1},{1},{1},{1},{1}}
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Length[Select[Join@@mps/@IntegerPartitions[n],And[SameQ@@Length/@Split[Sort[Join@@#]],SameQ@@Length/@#]&]],{n,10}]

Extensions

Terms a(11) and beyond from Andrew Howroyd, Feb 04 2022
Previous Showing 11-20 of 46 results. Next