(PARI) a(n)=(-1/2)^(n-2)*sum(k=0, n, binomial(n, k)*(1-2^(n+k+1))*bernfrac(n+k+1))
(PARI) {a(n)=local(CF=1+x*O(x^n)); if(n<1, return(0), for(k=1, n, CF=1/(1-((n-k)\2+1)*((n-k)\2+2)/2*x*CF)); return(Vec(CF)[n]))} (Hanna)
(PARI) {a(n)=polcoeff( x*sum(m=0, n, m!*(m+1)!*(x/2)^m / prod(k=1, m, 1 + k*(k+1)*x/2 +x*O(x^n)) ), n)} \\
Paul D. Hanna, Feb 03 2013
(SageMath) # Algorithm of L. Seidel (1877)
# n -> [a(1), ..., a(n)] for n >= 1.
D = [0]*(n+2); D[1] = 1
R = []; z = 1/2; b = False
for i in(0..2*n-1) :
h = i//2 + 1
if b :
for k in range(h-1, 0, -1) : D[k] += D[k+1]
z *= 2
else :
for k in range(1, h+1, 1) : D[k] += D[k-1]
b = not b
if not b : R.append(D[1]/z)
return R
(Python)
from math import comb
from sympy import bernoulli
def
A000366(n): return (-1 if n&1 else 1)*sum(comb(n, k)*(1-(1<<n+k+1))*bernoulli(n+k+1) for k in range(n+1))>>n-2 if n>1 else 1 #
Chai Wah Wu, Apr 14 2023