OFFSET
0,2
LINKS
N. J. A. Sloane, Transforms
FORMULA
G.f.: exp(Sum_{k>=1} ( Sum_{d|k} ((-1)^(k/d+1) + 1)*d^(d+1) ) * x^k/k).
a(n) ~ 2 * n^n * (1 + 2*exp(-1)/n + (exp(-1) + 10*exp(-2))/n^2). - Vaclav Kotesovec, Nov 09 2018
MAPLE
a:=series(mul(((1+x^k)/(1-x^k))^(k^k), k=1..100), x=0, 21): seq(coeff(a, x, n), n=0..20); # Paolo P. Lava, Apr 02 2019
MATHEMATICA
nmax = 20; CoefficientList[Series[Product[((1 + x^k)/(1 - x^k))^(k^k), {k, 1, nmax}], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[((-1)^(k/d + 1) + 1) d^(d + 1), {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 20}]
PROG
(PARI) seq(n)={Vec(exp(sum(k=1, n, sumdiv(k, d, ((-1)^(k/d+1) + 1)*d^(d+1) ) * x^k/k) + O(x*x^n)))} \\ Andrew Howroyd, Nov 09 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 08 2018
STATUS
approved