s := proc(z) 1/6*z^2*hypergeom([1, 1, 3/2], [2, 5/2], z^2); end proc;
f := proc(n, m) coeftayl(exp(-n*s(z)), z = 0, m); end proc;
g := proc(n) local v, m; v := 1; for m to n do v := lcm(v, denom(coeftayl(exp(-n*s(z)), z = 0, m))); end do; v; end proc;
cc := proc(n) local vv, m; if 0 = n mod 2 then vv := seq(g(n)*f(n, 2*m), m = 0 .. 1/2*n); end if; if 1 = n mod 2 then vv := seq(g(n)*f(n, 2*m), m = 0 .. 1/2*n - 1/2); end if; vv; end proc;
for i from 0 to 10 do cc(i); end do;