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.

A260883 Number of m-shape ordered set partitions, square array read by ascending antidiagonals, A(m, n) for m, n >= 0.

This page as a plain text file.
%I A260883 #39 Mar 16 2020 08:51:50
%S A260883 1,1,1,1,1,3,1,1,3,9,1,1,7,13,35,1,1,21,121,75,161,1,1,71,1849,3907,
%T A260883 541,913,1,1,253,35641,426405,202741,4683,6103,1,1,925,762763,
%U A260883 65782211,203374081,15430207,47293,47319,1,1,3433,17190265,11872636325,323213457781,173959321557
%N A260883 Number of m-shape ordered set partitions, square array read by ascending antidiagonals, A(m, n) for m, n >= 0.
%C A260883 A set partition of m-shape is a partition of a set with cardinality m*n for some n >= 0 such that the sizes of the blocks are m times the parts of the integer partitions of n. It is ordered if the positions of the blocks are taken into account.
%C A260883 If m = 0, all possible sizes are zero. Thus the number of ordered set partitions of 0-shape is the number of ordered partitions of n (partition numbers A101880).
%C A260883 If m = 1, the set is {1, 2, ..., n} and the set of all possible sizes are the integer partitions of n. Thus the number of ordered set partitions of 1-shape is a Fubini number (sequence A000670).
%C A260883 If m = 2, the set is {1, 2, ..., 2n} and the number of ordered set partitions of 2-shape is also the number of 2-packed words of degree n (sequence A094088).
%F A260883 From _Petros Hadjicostas_, Aug 02 2019: (Start)
%F A260883 Conjecture: For n >= 0, let P be the set of all possible lists (a_1, ..., a_n) of nonnegative integers such that a_1*1 + a_2*2 + ... + a_n*n = n. Consider terms of the form multinomial(n*m, m*[1,..., 1, 2,..., 2,..., n,..., n]) * multinomial(a_1 + ... + a_n, [a_1,..., a_n]), where in the list [1,..., 1, 2,..., 2,..., n,..., n] the number 1 occurs a_1 times, 2 occurs a_2 times, ..., and n occurs a_n times. (Here a_n = 0 or 1.) Summing these terms over P we get A(m, n) provided m >= 1. (End)
%e A260883 [ n ] [0  1   2      3         4            5                  6]
%e A260883 [ m ] -----------------------------------------------------------
%e A260883 [ 0 ] [1, 1,  3,     9,       35,          161,              913]  A101880
%e A260883 [ 1 ] [1, 1,  3,    13,       75,          541,             4683]  A000670
%e A260883 [ 2 ] [1, 1,  7,   121,     3907,       202741,         15430207]  A094088
%e A260883 [ 3 ] [1, 1, 21,  1849,   426405,    203374081,     173959321557]  A243664
%e A260883 [ 4 ] [1, 1, 71, 35641, 65782211, 323213457781, 3482943541940351]  A243665
%e A260883         A244174
%e A260883 For example the number of ordered set partitions of {1,2,...,9} with sizes in [9], [6,3] and [3,3,3] is 1, 168 and 1680 respectively. Thus A(3,3) = 1849.
%e A260883 Formatted as a triangle:
%e A260883 [1]
%e A260883 [1, 1]
%e A260883 [1, 1, 3]
%e A260883 [1, 1, 3, 9]
%e A260883 [1, 1, 7, 13, 35]
%e A260883 [1, 1, 21, 121, 75, 161]
%e A260883 [1, 1, 71, 1849, 3907, 541, 913]
%e A260883 [1, 1, 253, 35641, 426405, 202741, 4683, 6103]
%o A260883 (Sage)
%o A260883 def A260883(m, n):
%o A260883     shapes = ([x*m for x in p] for p in Partitions(n))
%o A260883     return sum(factorial(len(s))*SetPartitions(sum(s), s).cardinality() for s in shapes)
%o A260883 for m in (0..4): print([A260883(m, n) for n in (0..6)])
%Y A260883 Without order: A260876.
%Y A260883 Cf. A000670, A094088, A101880, A243664, A243665, A243666, A244174.
%K A260883 nonn,tabl
%O A260883 1,6
%A A260883 _Peter Luschny_, Aug 02 2015