|
%I #10 Nov 21 2023 08:41:00
%S 3,868,137683728,33559369419469824,30814027871378630714538393600,
%T 228763757339598304909903916639768291573760000,
%U 26160851477909352213716196682766580796112417641624416288768000000
%N v(n+1)/v(n), where v=A203687.
%C See A093883 for a discussion and guide to related sequences.
%F From _Vaclav Kotesovec_, Nov 21 2023: (Start)
%F a(n) ~ (n+1)!^(2*n).
%F a(n) ~ (2*Pi)^n * n^(2*n^2 + 3*n) / exp(2*n^2 - 13/6). (End)
%t f[j_] := j!; z = 8;
%t u[n_] := Product[f[j]^2 - f[j] f[k] + f[k]^2,
%t {j, 1, k - 1}]
%t v[n_] := Product[u[n], {k, 2, n}]
%t Table[v[n], {n, 1, z}] (* A203687 *)
%t Table[v[n + 1]/v[n], {n, 1, z}] (* A203688 *)
%t Table[Product[k!^2 - k!*(n+1)! + (n+1)!^2, {k, 1, n}], {n, 1, 10}] (* _Vaclav Kotesovec_, Nov 21 2023 *)
%K nonn
%O 1,1
%A _Clark Kimberling_, Jan 04 2012
|
