(PARI) getmnk(m, n, k) = {if (n<k || k<0, return (0)); if (k==0, return (1)); if (! m[n, k], if (n==k, m[n, k] = 1, m[n, k] = sum(j=0, n-k, getmnk(m, n-k, j)*(2*k)^j))); m[n, k]; }
lista(nn) = {my(m=matrix(nn, nn)); for(n=1, nn, for (k=1, n, m[n, k] = getmnk(m, n, k); ); ); vector(nn, n, m[n, 1]); } \\
Michel Marcus, Jun 03 2019