p1[n_] := Product[x^k + 1, {k, 1, n - 1, 2}] // Expand; p2[n_] := Product[x^k + 1, {k, 1, n/2}] // Expand; c[n_] := Coefficient[p1[n], x, n^2/8] * Coefficient[p2[n], x, n (n/2 + 1)/8]; a[n_] := c[8*n]/4; Table[a[n], {n, 1, 13}] (*
Jean-François Alcover, Jul 24 2013, after
Eric W. Weisstein *)