OFFSET
3,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 3..300
FORMULA
E.g.f.: (G(x)^2-H(x))/2 where G(x) = Sum {k>=1} (exp(x^k/k!)-1) and H(x) = Sum {k>=1} (exp(x^k/k!)-1)^2. - Vladeta Jovovic, Sep 18 2007
MAPLE
with(numtheory): with(combinat):
a:= n-> add(add(add(multinomial(n, i$j, d$((n-i*j)/d))/j!/((n-i*j)/d)!,
d=select(x->x<i, divisors(n-i*j))), j=1..n/i), i=2..n-1):
seq(a(n), n=0..30); # Alois P. Heinz, Feb 01 2014
MATHEMATICA
max = 25; G[x_] = Sum[Exp[x^k/k!]-1, {k, 1, max}]; H[x_] = Sum[(Exp[x^k/k!]-1)^2, {k, 1, max}]; Drop[CoefficientList[(G[x]^2-H[x])/2 + O[x]^max, x]*Range[0, max-1]!, 3] (* Jean-François Alcover, Jul 01 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Nov 02 2003
STATUS
approved