A128612 - OEIS

%I #20 Oct 11 2019 03:02:36

%S 1,0,1,0,2,1,1,5,5,1,1,14,30,14,1,0,28,155,147,29,1,0,56,605,1208,586,

%T 64,1,1,127,2133,7819,7819,2133,127,1,1,262,7288,44074,78190,44074,

%U 7288,262,1,0,496,23947,227623,655039,655315,227569,23893,517,1,0,992,76305,1102068,4868556,7862124,4869558,1101420,76332,1044,1

%N Triangle T(n,k) read by rows: number of permutations in [n] with exactly k ascents that have an even number of inversions.

%H Alois P. Heinz, <a href="/A128612/b128612.txt">Rows n = 1..141, flattened</a>

%H Jason Fulman, Gene B. Kim, Sangchul Lee, T. Kyle Petersen, <a href="https://arxiv.org/abs/1910.04258">On the joint distribution of descents and signs of permutations</a>, arXiv:1910.04258 [math.CO], 2019.

%H S. Tanimoto, <a href="http://arXiv.org/abs/math.CO/0602263">A new approach to signed Eulerian numbers</a>, arXiv:math/0602263 [math.CO], 2006.

%F T(n,k) = (1/2) * [A008292(n,n-k)+A049061(n,n-k)], n>=1, 0<=k<n. - _R. J. Mathar_, Nov 01 2007

%e Triangle starts:

%e 1;

%e 0,1;

%e 0,2,1;

%e 1,5,5,1;

%e 1,14,30,14,1;

%e 0,28,155,147,29,1;

%e 0,56,605,1208,586,64,1;

%e 1,127,2133,7819,7819,2133,127,1;

%p A008292 := proc(n,k) local j; add( (-1)^j*(k-j)^n*binomial(n+1,j),j=0..k) ; end: A049061 := proc(n,k) if k <= 0 or n <=0 or k > n then 0; elif n = 1 then 1 ; elif n mod 2 = 0 then A049061(n-1,k)-A049061(n-1,k-1) ; else k*A049061(n-1,k)+(n-k+1)*A049061(n-1,k-1) ; fi ; end: A128612 := proc(n,k) (A008292(n,n-k)+A049061(n,n-k))/2 ; end: for n from 1 to 11 do for k from 0 to n-1 do printf("%d,",A128612(n,k)) ; od: od: # _R. J. Mathar_, Nov 01 2007

%p # second Maple program:

%p b:= proc(u, o, i) option remember; expand(`if`(u+o=0, 1-i,

%p        add(b(u+j-1, o-j, irem(i+u+j-1, 2)), j=1..o)*x+

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

%p     end:

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

%p seq(T(n), n=1..14);  # _Alois P. Heinz_, May 02 2017

%t b[u_, o_, i_] := b[u, o, i] = Expand[If[u + o == 0, 1 - i, Sum[b[u + j - 1, o - j, Mod[i + u + j - 1, 2]], {j, 1, o}]*x + Sum[b[u - j, o + j - 1, Mod[i + u - j, 2]], {j, 1, u}]]];

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

%t Table[T[n], {n, 1, 14}] // Flatten (* _Jean-François Alcover_, Jul 25 2017, after _Alois P. Heinz_ *)

%Y Cf. A008292, A049061, A128613.

%K nonn,tabl

%O 1,5

%A _Ralf Stephan_, May 08 2007

%E More terms from _R. J. Mathar_, Nov 01 2007