%I #10 Dec 28 2020 09:53:37
%S 14,300,4400,55692,657370,7488228,83752760,928406556,10254052556,
%T 113186465340,1250820198264,13852280754980,153813849202674,
%U 1712835575525140,19129590267619304,214261857777632700,2406509409480345364,27100348605141932540,305944173898725745944
%N Number of words of length 2n such that all letters of the quaternary alphabet occur at least once and are introduced in ascending order and which can be built by repeatedly inserting doublets into the initially empty word.
%H Alois P. Heinz, <a href="/A258491/b258491.txt">Table of n, a(n) for n = 4..900</a>
%F a(n) ~ 12^n / (8*sqrt(Pi)*n^(3/2)). - _Vaclav Kotesovec_, Jun 01 2015
%p A:= proc(n, k) option remember; `if`(n=0, 1, k/n*
%p add(binomial(2*n, j)*(n-j)*(k-1)^j, j=0..n-1))
%p end:
%p T:= (n, k)-> add((-1)^i*A(n, k-i)/(i!*(k-i)!), i=0..k):
%p a:= n-> T(n, 4):
%p seq(a(n), n=4..25);
%t A[n_, k_] := A[n, k] = If[n == 0, 1, (k/n) Sum[Binomial[2n, j] (n - j)* If[j == 0, 1, (k - 1)^j], {j, 0, n - 1}]];
%t T[n_, k_] := Sum[(-1)^i A[n, k - i]/(i! (k - i)!), {i, 0, k}];
%t a[n_] := T[n, 4];
%t a /@ Range[4, 25] (* _Jean-François Alcover_, Dec 28 2020, after _Alois P. Heinz_ *)
%Y Column k=4 of A256117.
%K nonn
%O 4,1
%A _Alois P. Heinz_, May 31 2015
|