(PARI) \\ Here C(n) gives columns of
A361950 as vector of e.g.f.'s.
S(M)={matrix(#M, #M, i, j, sum(k=0, i-j, 2^((j-1)*k)*M[i-j+1, k+1])/(j-1)! )}
C(n, m=n)={my(M=matrix(n+1, n+1), c=vector(m+1), A=O(x*x^n)); M[1, 1]=1; c[1]=1+A; for(h=1, m, M=S(M); c[h+1]=sum(i=0, n, vecsum(M[i+1, ])*x^i, A)); c}
T(n)={my(c=C(n), b=vector(n+1, h, c[h]/c[max(h-1, 1)])); Mat(vector(n+1, h, Col(serlaplace(b[h]-if(h>1, b[h-1])), -n-1)))}
{ my(A=T(7)); for(n=1, #A, print(A[n, 1..n])) }