login
Distributes the number of permutations in the alternating group; cf. A060351.
2

%I #20 Feb 14 2021 13:18:26

%S 1,1,1,0,1,1,1,0,1,1,3,1,1,3,1,1,1,2,5,2,5,8,5,2,2,5,8,5,2,5,2,1,1,2,

%T 8,4,9,18,12,6,8,19,31,17,12,21,9,3,2,10,19,14,18,30,21,6,4,14,17,10,

%U 6,6,3,0

%N Distributes the number of permutations in the alternating group; cf. A060351.

%C The symmetric group distribution of permutation descents is summarized in Table A008292; for example 1 57 302 302 57 1 sums the following A060351 values:

%C 1.......5.......10.......10.......5.......1

%C .......14.......35.......35.......14.......

%C .......19.......26.......40.......19.......

%C .......14.......40.......19.......14.......

%C ........5.......61.......26.......5.......

%C ................26.......61..............

%C ................19.......40..............

%C ................40.......26..............

%C ................35.......35..............

%C ................10.......10..............

%H Alois P. Heinz, <a href="/A124921/b124921.txt">Rows n = 0..14, flattened</a>

%e The distribution is based on the frequency of descents; for example, when permuting four symbols the 12 patterns are ddd ddu dud udu dud duu udd udu dud udu uud and uuu yielding the frequency distribution 1 1 3 1 1 3 1 1.

%e Triangle T(n,k) begins:

%e 1;

%e 1;

%e 1, 0;

%e 1, 1, 1, 0;

%e 1, 1, 3, 1, 1, 3, 1, 1;

%e 1, 2, 5, 2, 5, 8, 5, 2, 2, 5, 8, 5, 2, 5, 2, 1;

%e ...

%p b:= proc(u, o, t, h) option remember; expand(`if`(u+o=0, h,

%p add(b(u-j, o+j-1, t+1, irem(h+u-j, 2))*x^floor(2^(t-1)), j=1..u)+

%p add(b(u+j-1, o-j, t+1, irem(h+u+j-1, 2)), j=1..o)))

%p end:

%p T:= n-> (p-> seq(coeff(p, x, i), i=0..ceil(2^(n-1)-1)))(b(n, 0$2, 1)):

%p seq(T(n), n=0..7); # _Alois P. Heinz_, Sep 09 2020

%t b[u_, o_, t_, h_] := b[u, o, t, h] = Expand[If[u+o == 0, h,

%t Sum[b[u-j, o+j-1, t+1, Mod[h+u-j, 2]]*x^Floor[2^(t-1)], {j, 1, u}]+

%t Sum[b[u+j-1, o-j, t+1, Mod[h+u+j-1, 2]], {j, 1, o}]]];

%t T[n_] := With[{p = b[n, 0, 0, 1]}, Table[Coefficient[p, x, i],

%t {i, 0, Ceiling[2^(n-1)-1]}]];

%t T /@ Range[0, 7] // Flatten (* _Jean-François Alcover_, Feb 14 2021, after _Alois P. Heinz_ *)

%Y Cf. A001710 (row sums), A008292, A060351, A011782 (row lengths).

%K tabf,look,nonn

%O 0,11

%A _Alford Arnold_, Nov 16 2006

%E More terms from _Alois P. Heinz_, Sep 09 2020