%I #22 Jan 17 2019 09:21:39
%S 1,3,16,111,931,9066,99925,1224579,16466524,240481659,3783550243,
%T 63705870894,1141642250521,21673475232051,434151377912680,
%U 9144596750676663,201921333704487283,4661560653703827474,112247618536314408013,2813147392487637246483
%N Expansion of e.g.f. exp(Sum_{d|M} (exp(d*x)-1)/d), M=4.
%H Vincenzo Librandi, <a href="/A141003/b141003.txt">Table of n, a(n) for n = 0..200</a>
%H J. Pasukonis, S. Ramgoolam, <a href="https://arxiv.org/abs/1010.1683">From counting to construction for BPS states in N=4SYM</a>, arXiv:1010.1683 [hep-th], 2010, (E.4).
%H J. Pasukonis, S. Ramgoolam, <a href="https://doi.org/10.1007/JHEP02(2011)078">From counting to construction for BPS states in N=4SYM</a>, J. High En. Phys. 2011 (2) (2011), (E.4).
%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,4],{n,0,12}] (* _Wouter Meeussen_, Dec 06 2008 *)
%Y Column 4 of A162663. u[n,j] generates the j-th column of A162663. - _Wouter Meeussen_, Dec 06 2008; _Danny Rorabaugh_, Oct 29 2015
%K nonn
%O 0,2
%A _R. J. Mathar_, Jul 11 2008