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.

A019589 Number of nondecreasing sequences that are differences of two permutations of 1,2,...,n.

Original entry on oeis.org

1, 1, 2, 5, 16, 59, 246, 1105, 5270, 26231, 135036, 713898, 3857113, 21220020, 118547774, 671074583
Offset: 0

Views

Author

Alex Postnikov (apost(AT)math.mit.edu)

Keywords

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
From Olivier Gérard, Sep 18 2007: (Start)
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.

Crossrefs

Programs

  • 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;
    Array[a, 10, 0] (* Jean-François Alcover, May 29 2019, after Vaclav Kotesovec *)
  • 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
  • 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)])
    # Bert Dobbelaere, Jan 19 2019
    

Formula

a(n) <= A007747(n) <= A362967(n). - Max Alekseyev, Jun 19 2023

Extensions

More terms from Olivier Gérard, Sep 18 2007
Two more terms from Vladeta Jovovic, Oct 04 2007
a(0)=1 prepended by Alois P. Heinz, Jul 24 2017
a(13)-a(14) from Bert Dobbelaere, Jan 19 2019
a(15) from Max Alekseyev, Jun 28 2023