cp's OEIS Frontend

A191936 Triangle read by rows of Legendre-Stirling numbers of the first kind.

%I A191936 #29 Jun 07 2021 04:38:54
%S A191936 1,1,0,1,-2,0,1,-8,12,0,1,-20,108,-144,0,1,-40,508,-2304,2880,0,1,-70,
%T A191936 1708,-17544,72000,-86400,0,1,-112,4648,-89280,808848,-3110400,
%U A191936 3628800,0,1,-168,10920,-349568,5808528,-48405888,177811200,-203212800,0
%N A191936 Triangle read by rows of Legendre-Stirling numbers of the first kind.
%C A191936 Apparently this is the mirror of triangle A129467. - _Omar E. Pol_, Jan 10 2012
%H A191936 G. C. Greubel, <a href="/A191936/b191936.txt">Rows n = 1..50 of the triangle, flattened</a>
%H A191936 G. E. Andrews, W. Gawronski and L. L. Littlejohn, <a href="http://www.math.psu.edu/andrews/pdf/283.pdf">The Legendre-Stirling Numbers</a>
%H A191936 G. E. Andrews et al., <a href="https://doi.org/10.1016/j.disc.2011.02.028">The Legendre-Stirling numbers</a>, Discrete Math., 311 (2011), 1255-1272.
%H A191936 J. Pan, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL16/Pan/pan21.html">Convolution Properties of the Generalized Stirling Numbers and the Jacobi-Stirling Numbers of the First Kind</a>, Journal of Integer Sequences, 16 (2013), #13.9.2.
%F A191936 T(n, k) = ps(n-1, n-k), where ps(n, k) = ps(n-1, k-1) - n*(n-1)*ps(n-1, k), ps(n, 0) = 0, and ps(n, n) = 1. - _G. C. Greubel_, Jun 07 2021
%e A191936 Triangle begins:
%e A191936   1;
%e A191936   1,    0;
%e A191936   1,   -2,    0;
%e A191936   1,   -8,   12,      0;
%e A191936   1,  -20,  108,   -144,      0;
%e A191936   1,  -40,  508,  -2304,   2880,        0;
%e A191936   1,  -70, 1708, -17544,  72000,   -86400,       0;
%e A191936   1, -112, 4648, -89280, 808848, -3110400, 3628800, 0;
%e A191936   ...
%t A191936 ps[n_, k_]:= ps[n, k]= If[k==n, 1, If[k==0, 0, ps[n-1, k-1] - n*(n-1)*ps[n-1, k]]];
%t A191936 Table[ps[n-1, n-k], {n,12}, {k,n}]//Flatten (* _G. C. Greubel_, Jun 07 2021 *)
%o A191936 (Sage)
%o A191936 @CachedFunction
%o A191936 def ps(n, k):
%o A191936     if (k==n): return 1
%o A191936     elif (k==0): return 0
%o A191936     else: return ps(n-1, k-1) - n*(n-1)*ps(n-1, k)
%o A191936 flatten([[ps(n-1, n-k) for k in (1..n)] for n in (1..12)]) # _G. C. Greubel_, Jun 07 2021
%Y A191936 Cf. A191935.
%K A191936 sign,tabl
%O A191936 1,5
%A A191936 _N. J. A. Sloane_, Jun 19 2011
%E A191936 More terms from _Omar E. Pol_, Jan 10 2012