login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A362112
a(0)=1; thereafter a(n) = 2*A110501(n+1) - A005439(n).
1
1, 1, 4, 26, 254, 3538, 67014, 1660866, 52230550, 2033261906, 96018823814, 5409008246626, 358368831222006, 27589872391918194, 2442595357421865574, 246430234111929035906, 28106918525950072081622, 3598669462582938225587602, 513978991104098010878849094
OFFSET
0,3
PROG
(Python)
from math import comb
from sympy import bernoulli
def A362112(n): return ((4<<(m:=n+1<<1))-4)*abs(bernoulli(m))-abs(sum(comb(n, k)*(2-(2<<n+k+1))*bernoulli(n+k+1) for k in range(n+1))) # Chai Wah Wu, Apr 14 2023
CROSSREFS
Sequence in context: A160886 A192546 A213438 * A056786 A350164 A006056
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Apr 14 2023
EXTENSIONS
More terms from Chai Wah Wu, Apr 14 2023
STATUS
approved