(PARI)x='x+O('x^66);Vec(1/2(2*x)*(1/(1-4*x)^(3/2)-1))\\约尔格·阿恩特2013年5月1日
(哈斯克尔)
a033876 n=总和$zipWith(!!)zss[0..n]其中
zss=取(n+1)$g(取(n+1)(1:[0,0..]),其中
g us=(取(n+1)$g’us):g(0:init us)
g’vs=最后$take(2*n+3)$
map snd$iterate h(0,vs++reverse vs)
h(p,ws)=(1-p,删除p$zipWith(+)([0]++ws)(ws++[0]))
(岩浆)[(2*n+3)*二项式(2*n+1,n):n in[0..40]]//韦斯利·伊万·赫特2017年11月30日