|
%I #23 Nov 21 2018 02:48:50
%S 1,1,7,64,843,13876,276792,6438797,170938483,5091463423,167965714273,
%T 6074571662270,238837895468954,10138497426332796,461941179848628434,
%U 22478593443737857695,1163160397700757351363,63760710281671647692688,3690276585886363643056992
%N Number of multisets of nonempty words with a total of n letters over n-ary alphabet.
%H Alois P. Heinz, <a href="/A252654/b252654.txt">Table of n, a(n) for n = 0..350</a>
%F a(n) = [x^n] Product_{j>=1} 1/(1-x^j)^(n^j).
%F a(n) = n-th term of the Euler transform of the powers of n.
%F a(n) ~ n^(n-3/4) * exp(2*sqrt(n) - 1/2) / (2*sqrt(Pi)). - _Vaclav Kotesovec_, Mar 14 2015
%F a(n) = [x^n] exp(n*Sum_{k>=1} x^k/(k*(1 - n*x^k))). - _Ilya Gutkovskiy_, Nov 20 2018
%e a(2) = 7: {aa}, {ab}, {ba}, {bb}, {a,a}, {a,b}, {b,b}.
%p with(numtheory):
%p A:= proc(n, k) option remember; `if`(n=0, 1, add(add(
%p d*k^d, d=divisors(j)) *A(n-j, k), j=1..n)/n)
%p end:
%p a:= n-> A(n$2):
%p seq(a(n), n=0..25);
%t A[n_, k_] := A[n, k] = If[n == 0, 1, Sum[Sum[d*k^d, {d, Divisors[j]}]*A[n - j, k], {j, 1, n}]/n];
%t a[n_] := A[n, n];
%t Table[a[n], {n, 0, 25}] (* _Jean-François Alcover_, Mar 21 2017, translated from Maple *)
%Y Main diagonal of A144074.
%Y Cf. A292805, A292873.
%K nonn
%O 0,3
%A _Alois P. Heinz_, Dec 19 2014
%E New name from comment by _Alois P. Heinz_, Sep 21 2018
|
