login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = Sum_{j=0..n, j even} binomial(n, j) * oddfactorial(j/2) * n^j, where oddfactorial(n) = (2*n)! / (2^n*n!).
2

%I #9 Jan 18 2023 09:34:20

%S 1,1,5,28,865,9626,758701,12606280,1872570113,41351249980,

%T 9925656304501,273345587759696,96567039881462305,3185756105692821688,

%U 1555524449985942662045,59790093545794928817376,38565845285812075675023361,1692346747225524397926264080,1393672439437011815394433559653

%N a(n) = Sum_{j=0..n, j even} binomial(n, j) * oddfactorial(j/2) * n^j, where oddfactorial(n) = (2*n)! / (2^n*n!).

%F Let K(n, x) = 2^(n/2)*(-1/x^2)^(-n/2)*KummerU(-n/2, 1/2, -1/(2*x^2)) denote the Kummer polynomials, defined in A359760.

%F a(n) = K(n, n) for n >= 1.

%p A359739 := n -> ifelse(n=0, 1, KummerU(-n/2, 1/2, -1/(2*n^2))*(-1/(2*n^2))^(-n/2)): seq(simplify(A359739(n)), n = 0..18);

%o (Python)

%o from math import factorial, comb

%o def oddfactorial(n: int) -> int: return factorial(2 * n) // (2**n * factorial(n))

%o def a(n: int) -> int:

%o return sum(comb(n, j) * oddfactorial(j//2) * n**j for j in range(0, n+1, 2))

%o print([a(n) for n in range(19)])

%Y Cf. A359760.

%K nonn

%O 0,3

%A _Peter Luschny_, Jan 12 2023