peven[i_]:=(Sum[(-1)^j*Binomial[n-i-1, j]/(n-i-1)!*(n-i+1-j)^(n/2)*(n-i-j)^(n/2-1), {j, 0, n-i-1}]);
poddblack[i_]:=(Sum[(-1)^j*Binomial[n-i-1, j]/(n-i-1)!*(n-i+1-j)^((n+1)/2)*(n-i-j)^((n-3)/2), {j, 0, n-i-1}]);
poddwhite[i_]:=(Sum[(-1)^j*Binomial[n-i-1, j]/(n-i-1)!*(n-i+1-j)^((n-1)/2)*(n-i-j)^((n-1)/2), {j, 0, n-i-1}]);
Table[If[n==1, 1, Sum[If[EvenQ[n], peven[i]*peven[n-i], poddblack[i]*poddwhite[n-i]], {i, 1, n-1}]], {n, 1, 50}]
(* Alternative formula with Stirling numbers of the second kind: *)
Table[If[n==1, 1, Sum[Sum[Binomial[Floor[(n+1)/2], j] * StirlingS2[j+Floor[n/2], n-i], {j, 0, Floor[(n+1)/2]}] * Sum[Binomial[Floor[n/2], j] * StirlingS2[j+Floor[(n+1)/2], i], {j, 0, Floor[n/2]}], {i, 1, n-1}]], {n, 1, 50}] (*
Vaclav Kotesovec, Mar 23 2011 *)