%I #12 Oct 17 2017 14:33:03
%S 1,4,34,410,6170,109506,2227882,50987378,1293386410,35933906946,
%T 1083069266634,35146231883122,1220316225699786,45099576850647106,
%U 1766204747811044266,73012853049788936050,3175271515525778312810
%N E.g.f. exp(sum_{d|M} (exp(d*x)-1)/d), M=10.
%H Vincenzo Librandi, <a href="/A141007/b141007.txt">Table of n, a(n) for n = 0..200</a>
%t u[0, j_]: = 1; u[k_, j_]: = u[k, j] = Sum[Binomial[k-1, i-1]Plus@@(u[k-i, j]#^(i-1)&/@Divisors[j]), {i, k}]; Table[u[n, 10], {n, 0, 30}] (* _Vincenzo Librandi_, Dec 12 2012 *)
%Y Column k=10 of A162663.
%K nonn
%O 0,2
%A _R. J. Mathar_, Jul 11 2008