cp's OEIS Frontend

A138784 a(n) = n*(n-1)*(n-2)(27*n^3 - 174*n^2 + 511*n - 552)*n!/48.

%I A138784 #18 Sep 08 2022 08:45:33
%S A138784 0,0,108,5232,154200,3747600,82908000,1756661760,36715472640,
%T A138784 770648256000,16426162368000,358029745920000,8015066233574400,
%U A138784 184788598259865600,4394609711331840000,107899664364011520000
%N A138784 a(n) = n*(n-1)*(n-2)(27*n^3 - 174*n^2 + 511*n - 552)*n!/48.
%C A138784 a(n) = Sum [f(L)^2 Sum h(u)^2*h(v)^2h(w)^2], where L is a partition of n, f(L) is the number of standard Young tableaux of shape L, h(z) is the hook length of the box z in L (i.e., in the Ferrers diagram of L), the inner summation is over all unordered triples of distinct boxes u, v and w in L and the outer summation is over all partitions of n. Example: a(3)=108 because for the partitions L=(3), (2,1), (1,1,1) of n=3 the values of f(L) are 1, 2, 1, respectively, the hook length multi-sets are {3,2,1}, {3,1,1},{3,2,1}, respectively, Sum h(u)^2*h(v)^2*h(w)^2 = 36, 9, 36, respectively and now a(n) 1^2*36 + 2^2*9 + 1^2*36 = 108.
%C A138784 In Proposition 6.12 of the Han paper the number 600 should be replaced by 552. - _Emeric Deutsch_, Dec 07 2015
%H A138784 Guo-Niu Han, <a href="http://arxiv.org/abs/0804.1849">An explicit expansion formula for the powers of the Euler product in terms of partition hook lengths</a>, arXiv:0804.1849v3 [math.CO] 9 May 2008 (p. 29).
%F A138784 D-finite 9*(n-3) *(4179399*n -20778487) *a(n) +(-37614591*n^3 -15080676*n^2 +317158871*n +3150537632) *a(n-1) +2*(44621643*n +200235296) *(-1+n)^2*a(n-2)=0. - _R. J. Mathar_, Dec 10 2015
%p A138784 A138784 := proc(n)
%p A138784     n*(n-1)*(n-2)*(27*n^3-174*n^2+511*n-552)*factorial(n)/48 ;
%p A138784 end proc:
%p A138784 seq(A138784(n),n=1..30) ;
%t A138784 Table[n (n - 1) (n - 2) (27 n^3 - 174 n^2 + 511 n - 552) n!/48, {n, 25}] (* _Vincenzo Librandi_, Dec 08 2015 *)
%o A138784 (PARI) a(n) = n*(n-1)*(n-2)*(27*n^3-174*n^2+511*n-552)*n!/48; \\ _Altug Alkan_, Dec 07 2015
%o A138784 (Magma) [n*(n-1)*(n-2)*(27*n^3-174*n^2+511*n-552)*Factorial(n)/48: n in [1..25]]; // _Vincenzo Librandi_, Dec 08 2015
%Y A138784 Cf. A138782, A138783.
%K A138784 nonn,easy
%O A138784 1,3
%A A138784 _Emeric Deutsch_, May 15 2008