login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A328431
Number of inversion sequences of length n avoiding the consecutive patterns 010, 021, 120, and 210.
15
1, 1, 2, 5, 15, 53, 214, 960, 4701, 24873, 141147, 853641, 5472642, 37024569, 263342224, 1962835806, 15288074104, 124120865849, 1048092680689, 9186689045482, 83435365244510, 783923558286071, 7608398620990535, 76177574145052258, 785853360840424425
OFFSET
0,3
COMMENTS
A length n inversion sequence e_1e_2...e_n is a sequence of integers such that 0 <= e_i < i. The term a(n) counts the inversion sequences of length n with no entries e_i, e_{i+1}, e_{i+2} such that e_i != e_{i+1} > e_{i+2}. This is the same as the set of inversion sequences of length n avoiding the consecutive patterns 010, 021, 120, and 210.
LINKS
EXAMPLE
The a(4)=15 length 4 inversion sequences avoiding the consecutive patterns 010, 021, 120, and 210 are 0000, 0110, 0001, 0011, 0111, 0002, 0012, 0112, 0022, 0122, 0003, 0013, 0113, 0023, and 0123.
MAPLE
b := proc(n, x, t) option remember; `if`(n = 0, 1, add(
`if`(t and i <> x, 0, b(n - 1, i, x < i)), i = 0 .. n - 1))
end proc:
a := n -> b(n, n, false):
seq(a(n), n = 0 .. 24);
MATHEMATICA
b[n_, x_, t_] := b[n, x, t] = If[n == 0, 1, Sum[If[t && i != x, 0, b[n - 1, i, x < i]], {i, 0, n - 1}]];
a[n_] := b[n, n, False];
a /@ Range[0, 24] (* Jean-François Alcover, Mar 02 2020 after Alois P. Heinz in A328357 *)
KEYWORD
nonn
AUTHOR
Juan S. Auli, Oct 16 2019
STATUS
approved