%I #15 Sep 08 2022 08:46:23
%S 1,1,2,11,74,708,8583,127424,2239965,45514345,1049365071,27061132159,
%T 771695223819,24109698083919,818914886275467,30044684789498522,
%U 1184048086192376822,49883929845112421452,2237287911899357657492,106426388125032988691636,5352033610656721914626572
%N Expansion of Product_{k>=1} (1 + x^k)^(k^(k-1)).
%C Weigh transform of A000169.
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%F G.f.: exp(Sum_{k>=1} ( Sum_{d|k} (-1)^(k/d+1)*d^d ) * x^k/k).
%F a(n) ~ n^(n-1) * (1 + exp(-1)/n + (3*exp(-1)/2 + 2*exp(-2))/n^2). - _Vaclav Kotesovec_, Nov 09 2018
%p a:=series(mul((1+x^k)^(k^(k-1)),k=1..100),x=0,21): seq(coeff(a,x,n),n=0..20); # _Paolo P. Lava_, Apr 02 2019
%t nmax = 20; CoefficientList[Series[Product[(1 + x^k)^(k^(k - 1)), {k, 1, nmax}], {x, 0, nmax}], x]
%t a[n_] := a[n] = If[n == 0, 1, Sum[Sum[(-1)^(k/d + 1) d^d, {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 20}]
%o (PARI) seq(n)={Vec(exp(sum(k=1, n, sumdiv(k, d, (-1)^(k/d+1)*d^d ) * x^k/k) + O(x*x^n)))} \\ _Andrew Howroyd_, Nov 09 2018
%o (PARI) m=30; x='x+O('x^m); Vec(prod(k=1,m,(1+x^k)^(k^(k-1)))) \\ _G. C. Greubel_, Nov 09 2018
%o (Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1+x^k)^(k^(k-1)): k in [1..m]]) )); // _G. C. Greubel_, Nov 09 2018
%Y Cf. A000169, A023879, A261053, A283335, A321385, A321388.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, Nov 08 2018
|