cp's OEIS Frontend

A201636 Triangle read by rows, n>=0, k>=0, T(0,0) = 1, T(n,k) = Sum_{j=0..k} (C(n+k,k-j)*(-1)^(k-j)*2^(n-j)*Sum_{i=0..j} (C(n+j,i)*|S(n+j-i,j-i)|)), S = Stirling number of first kind.

%I A201636 #17 Jun 29 2019 02:17:02
%S A201636 1,0,1,0,4,3,0,24,40,15,0,192,520,420,105,0,1920,7392,9520,5040,945,0,
%T A201636 23040,116928,211456,176400,69300,10395,0,322560,2055168,4858560,
%U A201636 5642560,3465000,1081080,135135,0,5160960,39896064,117722880,177580480,150870720,73153080,18918900,2027025
%N A201636 Triangle read by rows, n>=0, k>=0, T(0,0) = 1, T(n,k) = Sum_{j=0..k} (C(n+k,k-j)*(-1)^(k-j)*2^(n-j)*Sum_{i=0..j} (C(n+j,i)*|S(n+j-i,j-i)|)),  S = Stirling number of first kind.
%C A201636 This triangle was inspired by a formula of Vladimir Kruchinin given in A001662.
%F A201636 T(n,1) = A002866(n) for n>0.
%F A201636 T(n,n) = A001147(n).
%F A201636 Sum((-1)^(n-k)*T(n,k)) = A001662(n+1).
%e A201636 [n\k 0,    1,    2,    3,    4,   5]
%e A201636 [0]  1,
%e A201636 [1]  0,    1,
%e A201636 [2]  0,    4,    3,
%e A201636 [3]  0,   24,   40,   15,
%e A201636 [4]  0,  192,  520,  420,  105,
%e A201636 [5]  0, 1920, 7392, 9520, 5040, 945,
%p A201636 A201636 := proc(n,k) if n=0 and k=0 then 1 else
%p A201636 add(binomial(n+k,k-j)*(-1)^(n+k-j)*2^(n-j)*
%p A201636 add(binomial(n+j,i)*stirling1(n+j-i,j-i),i=0..j),j=0..k) fi end:
%p A201636 for n from 0 to 8 do print(seq(A201636(n,k),k=0..n)) od;
%t A201636 T[0, 0] = 1; T[n_, k_] := T[n, k] = Sum[Binomial[n+k, k-j]*(-1)^(n+k-j)* 2^(n-j)*Sum[Binomial[n+j, i]*StirlingS1[n+j-i, j-i], {i, 0, j}], {j, 0, k}];
%t A201636 Table[T[n, k], {n, 0, 8}, {k, 0, n}] (* _Jean-François Alcover_, Jun 29 2019 *)
%o A201636 (Sage)
%o A201636 def A201636(n,k) :
%o A201636     if n==0 and k==0: return 1
%o A201636     return add(binomial(n+k,k-j)*(-1)^(k-j)*2^(n-j)*add(binomial(n+j,i)* stirling_number1(n+j-i,j-i) for i in (0..j)) for j in (0..k))
%Y A201636 Cf. A001662, A001147, A002866.
%K A201636 nonn,tabl
%O A201636 0,5
%A A201636 _Peter Luschny_, Nov 13 2012