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.

A328436 Number of inversion sequences of length n avoiding the consecutive patterns 000 and 001.

Original entry on oeis.org

1, 1, 2, 3, 9, 37, 190, 1181, 8564, 70914, 659810, 6811371, 77232836, 953969548, 12747856402, 183218649413, 2818050980941, 46182485773217, 803323102085452, 14781372445602234, 286838921699435184, 5854404018902152208, 125367868007259046305, 2810511319383912299122
Offset: 0

Views

Author

Juan S. Auli, Oct 17 2019

Keywords

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 000 and 001.

Examples

			The a(4)=9 length 4 inversion sequences avoiding the consecutive patterns 000 and 001 are 0100, 0110, 0120, 0101, 0121, 0102, 0122, 0103, and 0123.
		

Crossrefs

Programs

  • Maple
    # after Alois P. Heinz in A328357
    b := proc(n, x, t) option remember; `if`(n = 0, 1, add(
           `if`(t and i = x, 0, b(n - 1, i, i <= x)), i = 0 .. n - 1))
         end proc:
    a := n -> b(n, -1, 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, i <= x]], {i, 0, n - 1}]];
    a[n_] := b[n, -1, False];
    a /@ Range[0, 24] (* Jean-François Alcover, Mar 02 2020 after Alois P. Heinz in A328357 *)