cp's OEIS Frontend

A180185 Triangle read by rows: T(n,k) is the number of permutations of [n] having no 3-sequences and having k successions (0 <= k <= floor(n/2)); a succession of a permutation p is a position i such that p(i +1) - p(i) = 1.

%I A180185 #13 May 23 2020 14:30:57
%S A180185 1,1,1,1,3,2,11,9,1,53,44,9,309,265,66,3,2119,1854,530,44,16687,14833,
%T A180185 4635,530,11,148329,133496,44499,6180,265,1468457,1334961,467236,
%U A180185 74165,4635,53,16019531,14684570,5339844,934472,74165,1854,190899411
%N A180185 Triangle read by rows: T(n,k) is the number of permutations of [n] having no 3-sequences and having k successions (0 <= k <= floor(n/2)); a succession of a permutation p is a position i such that p(i +1) - p(i) = 1.
%C A180185 Row n has 1+floor(n/2) entries.
%C A180185 Sum of entries in row n is A002628(n).
%F A180185 T(n,k) = binomial(n-k,k)*(d(n-k) + d(n-k-1)), where d(j) = A000166(j) are the derangement numbers.
%F A180185 T(n,0) = d(n) + d(n-1) = A000255(n-1).
%F A180185 T(n,1) = d(n).
%F A180185 Sum_{k>=0} k*T(n,k) = A002629(n+1).
%e A180185 T(6,3)=3 because we have 125634, 341256, and 563412.
%e A180185 Triangle starts:
%e A180185      1;
%e A180185      1;
%e A180185      1,    1;
%e A180185      3,    2;
%e A180185     11,    9,    1;
%e A180185     53,   44,    9;
%e A180185    309,  265,   66,    3;
%e A180185   2119, 1854,  530,   44;
%p A180185 d[0] := 1: for n to 51 do d[n] := n*d[n-1]+(-1)^n end do: a := proc (n, k) if n = 0 and k = 0 then 1 elif k <= (1/2)*n then binomial(n-k, k)*d[n+1-k]/(n-k) else 0 end if end proc: for n from 0 to 12 do seq(a(n, k), k = 0 .. (1/2)*n) end do; # yields sequence in triangular form
%t A180185 d[0] = 1; d[n_] := d[n] = n d[n - 1] + (-1)^n;
%t A180185 T[n_, k_] := If[n == 0 && k == 0, 1, If[k <= n/2, Binomial[n - k, k] d[n + 1 - k]/(n - k), 0]];
%t A180185 Table[T[n, k], {n, 0, 20}, {k, 0, Quotient[n, 2]}] // Flatten (* _Jean-François Alcover_, May 23 2020 *)
%o A180185 (PARI) d(n) = if(n<2, !n , round(n!/exp(1)));
%o A180185 for(n=0, 20, for(k=0, (n\2), print1(binomial(n - k, k)*(d(n - k) + d(n - k - 1)),", ");); print();) \\ _Indranil Ghosh_, Apr 12 2017
%Y A180185 Cf. A000166, A002628, A002629, A000255.
%K A180185 nonn,tabf
%O A180185 0,5
%A A180185 _Emeric Deutsch_, Sep 06 2010