%I M5382 N2338 #64 Aug 05 2023 13:08:47
%S 1,120,7308,303660,11098780,389449060,13642629000,486591585480,
%T 17856935296200,678103775949600,26726282654771700,1094862336960892500,
%U 46641683693715610500,2066075391660447667500,95122549872697437090000
%N Second-order reciprocal Stirling number (Fekete) a(n) = [[2n+4, n]]. The number of n-orbit permutations of a (2n+4)-set with at least 2 elements in each orbit. Also known as associated Stirling numbers of the first kind (e.g., Comtet).
%D L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 256.
%D C. Jordan, Calculus of Finite Differences. Budapest, 1939, p. 152.
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H A. E. Fekete, <a href="http://www.jstor.org/stable/2974533">Apropos two notes on notation</a>, Amer. Math. Monthly, 101 (1994), 771-778.
%H H. W. Gould, Harris Kwong, and Jocelyn Quaintance, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL18/Kwong/kwong9.html">On Certain Sums of Stirling Numbers with Binomial Coefficients</a>, J. Integer Sequences, 18 (2015), #15.9.6.
%H C. Jordan, <a href="https://www.jstage.jst.go.jp/article/tmj1911/37/0/37_0_254/_pdf">On Stirling's Numbers</a>, Tohoku Math. J., 37 (1933), 254-278.
%F a(n) = [[2n+4, n]] = Sum_{i=0..n} (-1)^i*binomial(2n+4, 2n+4-i)*[2n+4-i, n-i] where [n, k] is the unsigned Stirling number of the first kind. - Barbara Haas Margolius (margolius(AT)math.csuohio.edu), Dec 14 2000
%F Recurrence: 30*(n-1)*(116*n+75)*a(n) + (-6960*n^3-49760*n^2-112691*n-80787)*a(n-1) + (n+1)*(2*n+1)*(20*n+21)*a(n-2) = 0. - _R. J. Mathar_, Jul 18 2015
%F For n>0, a(n) = (1113 + 1447*n + 600*n^2 + 80*n^3)*(2*n+4)!/(1215*2^(n+3)*(n-1)!). - _Vaclav Kotesovec_, Jan 17 2016
%F Recurrence (for n>1): (n-1)*(80*n^3 + 360*n^2 + 487*n + 186)*a(n) = (n+2)*(2*n+3)*(80*n^3 + 600*n^2 + 1447*n + 1113)*a(n-1). - _Vaclav Kotesovec_, Jan 18 2016
%p with(combinat):s1 := (n,k)->sum((-1)^i*binomial(n,i)*abs(stirling1(n-i,k-i)),i=0..n); for j from 1 to 20 do s1(2*j+4,j); od; # Barbara Haas Margolius (margolius(AT)math.csuohio.edu), Dec 14 2000
%t Prepend[Table[Sum[(-1)^i Binomial[2 n + 4, 2 n + 4 - i] Abs@ StirlingS1[2 n + 4 - i, n - i], {i, 0, n}], {n, 14}] , 1] (* _Michael De Vlieger_, Jan 04 2016 *)
%o (PARI) a(n) = if (!n, 1, sum(i=0, n, (-1)^i*binomial(2*n+4, 2*n+4-i)*abs(stirling(2*n+4-i, n-i, 1)))); \\ _Michel Marcus_, Jan 04 2016
%o (Magma) [1] cat [(1113+1447*n+600*n^2+80*n^3)*Factorial(2*n+4)/(1215*2^(n+ 3)*Factorial(n-1)): n in [1..15]]; // _Vincenzo Librandi_, Jan 18 2016
%Y Cf. A000907, A000483, A001784.
%K nonn
%O 0,2
%A _N. J. A. Sloane_
%E More terms from Barbara Haas Margolius (margolius(AT)math.csuohio.edu), Dec 14 2000
%E Offset changed to 0 by _Michel Marcus_, Jan 04 2016
|