|
|
A019589
|
|
Number of nondecreasing sequences that are differences of two permutations of 1,2,...,n.
|
|
7
|
|
|
1, 1, 2, 5, 16, 59, 246, 1105, 5270, 26231, 135036, 713898, 3857113, 21220020, 118547774, 671074583
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Also, number of nondecreasing sequences that are sums of two permutations of order n. If nondecreasing requirement is dropped, the sequence becomes A175176. - Max Alekseyev, Jun 19 2023
Number of classes of permutations arrays giving the same partition by the following transformation (equivalent in this case to X-rays but more general): on the matrix representation of a permutation of order n, the sum (i.e., here, number of ones) in the i-th antidiagonal is the number of copies of i in the partition.
This gives an injection of permutations of n into partitions with parts at most 2n-1. The first or the last antidiagonal can be omitted, reducing the size of parts to 2n-2 without changing the number of classes.
This sequence is called Lambda_{n,1} in Louck's paper and appears explicitly on p. 758. Terms up to 10 were computed by Myron Stein (arXiv).
This is similar to the number of Schur functions studied by Di Francesco et al. (A007747) related to the powers of the Vandermonde determinant. Also number of classes of straight (monotonic) crossing bi-permutations. (End)
Also number of monomials in expansion of permanent of an n X n Hankel matrix [t(i+j)] in terms of its entries (cf. A019448). - Vaclav Kotesovec, Mar 29 2019
|
|
REFERENCES
|
Olivier Gérard and Karol Penson, Set partitions, multiset permutations and bi-permutations, in preparation.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MAPLE
|
with(LinearAlgebra): f:=n->nops([coeffs(Permanent(Matrix(n, (i, j) -> a[i+j])))]): [seq(f(n), n=1..10)]; # Vaclav Kotesovec, Mar 29 2019
|
|
MATHEMATICA
|
a[n_] := Table[b[i+j], {i, n}, {j, n}] // Permanent // Expand // Length;
|
|
PROG
|
(Python)
import itertools
def a019589(n):
s = set()
for p in itertools.permutations(range(n)):
s.add(tuple(sorted([k - p[k] for k in range(n)])))
return len(s)
print([a019589(n) for n in range(10)])
(PARI) a(n) = my(l=List(), v=[1..n]); for(i=0, n!-1, listput(l, vecsort(v-numtoperm(n, i)))); listsort(l, 1); #l
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more,nice
|
|
AUTHOR
|
Alex Postnikov (apost(AT)math.mit.edu)
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|