OFFSET
0,3
MAPLE
a := n -> local j; add(binomial(n, j)*euler(j, 0)*(2*n)^j, j = 0..n):
seq(a(n), n = 0..16);
MATHEMATICA
a[0] := 1; a[n_] := Sum[Binomial[n, j] EulerE[j, 0] (2n)^j, {j, 0, n}];
Table[a[n], {n, 0, 16}]
PROG
(Python)
from math import comb
from sympy import euler
def A377664(n): return sum(comb(n, j)*euler(j, 0)*(n<<1)**j for j in range(n+1)) # Chai Wah Wu, Nov 14 2024
CROSSREFS
KEYWORD
sign
AUTHOR
Peter Luschny, Nov 14 2024
STATUS
approved