OFFSET
0,2
COMMENTS
Diagonal of rational function 1/(1 - (x^4 + y^4 + z^4 - t^4 + 4*x*y*z*t)). - Gheorghe Coserea, Aug 04 2018
FORMULA
Let H(m, n, x) = m^n*hypergeom([(k-n)/m for k=0..m-1], [1 for k=0..m-2], x) then a(n) = H(4, n, -1).
From Peter Bala, Sep 26 2025: (Start)
a(n) = Sum_{k = 0..floor(n/4)} (-1)^k * 4^(n-4*k) * n!/((n-4*k)!*k!^4).
E.g.f.: exp(4*x) * Sum_{k >= 0} (-1)^k * x^(4*k)/k!^4 = 1 + 4*x + 16*x^2/2! + 64*x^3/3! + 232*x^4/4! + 544*x^5/5! - ....
n^3*a(n) = 4*(2*n - 1)*(2*n^2 - 2*n + 1)*a(n-1) - 16*(n - 1)*(6*n^2 - 12*n + 7)*a(n-2) + 128*(n - 1)*(n - 2)*(2*n - 3)*a(n-3) - 512*(n - 1)*(n - 2)*(n - 3)*a(n-4) with a(0) = 1, a(1) = 4, a(2) = 16 and a(3) = 64. (End)
MAPLE
T := (m, n, x) -> m^n*hypergeom([seq((k-n)/m, k=0..m-1)], [seq(1, k=0..m-2)], x):
lprint(seq(simplify(T(4, n, -1)), n=0..39));
MATHEMATICA
Table[4^n * HypergeometricPFQ[{-n/4, (1-n)/4, (2-n)/4, (3-n)/4}, {1, 1, 1}, -1], {n, 0, 20}] (* Vaclav Kotesovec, Nov 02 2017 *)
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Peter Luschny, Nov 02 2017
STATUS
approved
