A341445 - OEIS

%I #23 Feb 13 2021 05:01:18

%S 1,0,2,2,0,3,2,6,0,6,8,8,16,0,10,16,32,24,40,0,20,52,84,108,60,90,0,

%T 35,134,262,294,310,150,210,0,70,432,816,1008,880,816,336,448,0,126,

%U 1248,2544,3192,3208,2460,2100,784,1008,0,252

%N Triangle read by rows: T(n,k) is the number of Dyck paths of semilength n having degree of symmetry k (n >= 1, 1 <= k <= n).

%C The degree of symmetry of a Dyck path is defined as the number of steps in the first half that are mirror images of steps in the second half, with respect to the reflection along the vertical line through the midpoint of the path.

%H Sergi Elizalde, <a href="https://arxiv.org/abs/2002.12874">The degree of symmetry of lattice paths</a>, arXiv:2002.12874 [math.CO], 2020.

%H Sergi Elizalde, <a href="https://www.mat.univie.ac.at/~slc/wpapers/FPSAC2020/26.html">Measuring symmetry in lattice paths and partitions</a>, Sem. Lothar. Combin. 84B.26, 12 pp (2020).

%e For n=4 there are 6 Dyck paths with degree of symmetry equal to 2: uuuddudd, uuduuddd, uududdud, uuddudud, uduududd, ududuudd.

%e Triangle begins:

%e      1;

%e      0,    2;

%e      2,    0,    3;

%e      2,    6,    0,    6;

%e      8,    8,   16,    0,   10;

%e     16,   32,   24,   40,    0,   20;

%e     52,   84,  108,   60,   90,    0,  35;

%e    134,  262,  294,  310,  150,  210,   0,   70;

%e    432,  816, 1008,  880,  816,  336, 448,    0, 126;

%e   1248, 2544, 3192, 3208, 2460, 2100, 784, 1008,   0, 252;

%e   ...

%p b:= proc(x, y, v) option remember; expand(

%p       `if`(min(y, v, x-max(y, v))<0, 0, `if`(x=0, 1, (l-> add(add(

%p       `if`(y=v+(j-i)/2, z, 1)*b(x-1, y+i, v+j), i=l), j=l))([-1, 1]))))

%p     end:

%p g:= proc(n) option remember; add(b(n, j$2), j=0..n) end:

%p T:= (n, k)-> coeff(g(n), z, k):

%p seq(seq(T(n, k), k=1..n), n=1..10);  # _Alois P. Heinz_, Feb 12 2021

%t b[x_, y_, v_] := b[x, y, v] = Expand[If[Min[y, v, x - Max[y, v]] < 0, 0, If[x == 0, 1, Function[l, Sum[Sum[If[y == v + (j - i)/2, z, 1]*b[x - 1, y + i, v + j], {i, l}], {j, l}]][{-1, 1}]]]];

%t g[n_] := g[n] = Sum[b[n, j, j], {j, 0, n}];

%t T[n_, k_] := Coefficient[g[n], z, k];

%t Table[Table[T[n, k], {k, 1, n}], {n, 1, 10}] // Flatten (* _Jean-François Alcover_, Feb 13 2021, after _Alois P. Heinz_ *)

%Y Equivalent to A298645 with rows reversed.

%Y Row sums give A000108.

%Y Main diagonal gives A001405.

%Y Column k=1 gives A298647 (for n>2).

%Y Second subdiagonal gives 2*A191522.

%K nonn,tabl

%O 1,3

%A _Sergi Elizalde_, Feb 12 2021