a(n) = f(3*n, 2*n), where f(m, n) = 2^(-m) * Sum_{j=0..n} (-1)^j*n!*m!*(2*m-2*j)!/(j!*(m-j)!*(n-j)!*6^(n-j)).
a(n) = ((6*n)!/(288)^n)*Sum_{j=0..2*n} b(2*n,j)*b(3*n,j)*(-6)^j/(j!*b(2*j, j)*b(6*n,2*j)), where b(x,y) = binomomial(x,y).
a(n) = (6*n)!/(288)^n * Hypergeometric1F1([-2*n], [1/2-3*n], -3/2). (End)
a(n) ~ sqrt(Pi) * 2^(n+1) * 3^(4*n + 1/2) * n^(6*n + 1/2) / exp(6*n+1). -
Vaclav Kotesovec, Oct 21 2023