(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(exp(-log(1-x)^5)))
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace((1-x)^(-log(1-x)^4)))
(PARI) a_vector(n) = my(v=vector(n+1)); v[1]=1; for(i=1, n, v[i+1]=120*sum(j=1, i, binomial(i-1, j-1)*abs(stirling(j, 5, 1))*v[i-j+1])); v;
(PARI) a(n) = sum(k=0, n\5, (5*k)!*abs(stirling(n, 5*k, 1))/k!);