a(n,m)=sum(product(S2(4;j,1)^e(n,m,q,j),j=1..n),q=1..p(n,m)) if n>=m>=1, else 0. Here p(n,m)=A008284(n,m), the number of m parts partitions of n and e(n,m,q,j) is the exponent of j in the q-th m part partition of n. S2(4;j,1)= A007559(j) = A035469(j,1) = (3*j-2)!!!.