editing
approved
editing
approved
a(n) ~ 4^n / (3*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 07 2021
approved
editing
proposed
approved
editing
proposed
a[n_] := a[n] = If[n < 4, n!, ((2*(12*n^3 - 91*n^2 + 213*n - 149))*a[n-1] - (3*(21*n^3 - 162*n^2 + 392*n - 291))*a[n-2] + (2*(33*n^3 - 257*n^2 + 633*n - 484))*a[n-3] - (4*(2*n - 7))*(3*n^2 - 13*n + 13)*a[n-4]) / ((n - 1)*(3*n^2 - 19*n + 29))]; Array[a, 30, 0] (* Jean-François Alcover, Nov 06 2017, after Alois P. Heinz *)
approved
editing
editing
approved
Alois P. Heinz, <a href="/A279563/b279563.txt">Table of n, a(n) for n = 0..1664</a>
a:= proc(n) option remember; `if`(n<4, n!,
((2*(12*n^3-91*n^2+213*n-149))*a(n-1)
-(3*(21*n^3-162*n^2+392*n-291))*a(n-2)
+(2*(33*n^3-257*n^2+633*n-484))*a(n-3)
-(4*(2*n-7))*(3*n^2-13*n+13)*a(n-4))
/ ((n-1)*(3*n^2-19*n+29)))
end:
seq(a(n), n=0..30); # Alois P. Heinz, Feb 22 2017
nonn,more,new
a(10)-a(26) from Alois P. Heinz, Feb 22 2017
1, 1, 2, 6, 22, 85, 328, 1253, 4754, 17994, 68158, 258808, 985906, 3768466, 14451386, 55585014, 214377618, 828795169, 3211030684, 12464308997, 48465092366, 188733879657, 735977084412, 2873525548315, 11231884145434, 43947466923095, 172115939825516
1, 1, 2, 6, 22, 85, 328, 1253, 4754, 17994, 68158
approved
editing