%I #8 Dec 07 2023 02:20:13
%S 1,4,1120,79833600,3173289799680000,123650071173117090201600000,
%T 7337799401269093351612002462597120000000,
%U 951792703318385182295191545713146608287219712000000000000
%N a(n) = Product_{i=1..n, j=1..n} (i + 3*j).
%F a(n) ~ A^(1/3) * 2^(16*n*(n+1)/3 + 13/18) * n^(n^2 - 19/36) / (Pi^(1/3) * Gamma(1/3)^(1/3) * 3^(n*(3*n+4)/2 + 11/36) * exp(3*n^2/2 + 1/36)), where A = A074962 is the Glaisher-Kinkelin constant.
%t Table[Product[i + 3*j, {i, 1, n}, {j, 1, n}], {n, 0, 10}]
%Y Cf. A074962, A079478, A324402, A367957, A367958.
%K nonn
%O 0,2
%A _Vaclav Kotesovec_, Dec 06 2023