With[{nn=20}, CoefficientList[Series[Log[1+Tan[Sin[x]]], {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Nov 01 2021 *)
PROG
(Maxima)
a(n):=sum(sum((((-1)^(k-m)+1)*(sum(binomial(j-1, m-1)*j!*2^(k-j-1)*stirling2(k, j)*(-1)^((m+k)/2+j), j, m, k))*((-1)^(n-k)+1)*sum((2*i-k)^n*binomial(k, i)*(-1)^((n+k)/2-i), i, 0, k/2))/(2^k*k!), k, m, n)/m*(-1)^(m-1), m, 1, n); /* Vladimir Kruchinin, May 05 2011 */