%I #31 Nov 19 2024 02:49:26
%S 1,10,300,11440,485100,21841260,1022041020,49128552000,2408829328620,
%T 119918393838100,6042249840712800,307438844121252480,
%U 15770112362658517500,814459593645444166560,42308586942403276440000,2208850973597860123741440,115825519836558228435979500
%N a(n) = binomial(3*n,3*n/2)*binomial(2*n,n)*binomial(5*n/2,n/2)/binomial(n,n/2).
%C Right-hand side of the binomial sum identity Sum_{k = 0..n} (-1)^(n+k)* binomial(3*n+k,3*n-k)*binomial(2*k,k)*binomial(2*n-k,n) = binomial(3*n,3*n/2)*binomial(2*n,n)*binomial(5*n/2,n/2)/binomial(n,n/2).
%C We also have Sum_{k = 0..3*n} (-1)^k*binomial(3*n+k,3*n-k)* binomial(2*k,k)*binomial(2*n-k,n) = binomial(3*n,3*n/2)*binomial(2*n,n)* binomial(5*n/2,n/2) /binomial(n,n/2).
%C Compare with Sum_{k = 0..n} (-1)^(n+k)*binomial(2*n+k,2*n-k)* binomial(2*k,k)*binomial(2*n-k,n) = binomial(2*n,n)^2 = A002894(n). See also A275653, A275654 and A275655.
%F a(n) = (3*n)!*(5*n/2)!*(n/2)!/((3*n/2)!^2*n!^3).
%F Recurrence: a(n) = 5*(3*n - 1)*(3*n - 5)*(5*n - 2)*(5*n - 4)*(5*n - 6)*(5*n - 8)/(n^2*(n - 1)^2*(3*n - 2)*(3*n - 4)) * a(n-2).
%F a(n) = [x^n] G(x)^(5*n), where G(x) = 1 + 2*x + 12*x^2 + 184*x^3 + 3811*x^4 + 92796*x^5 + 2497358*x^6 + ... appears to have integer coefficients.
%F exp(Sum_{n >= 1} a(n)*x^n/n) = F(x)^5, where F(x) = 1 + 2*x + 32*x^2 + 824*x^3 + 26291*x^4 + 947506*x^5 + 36934522*x^6 + ... appears to have integer coefficients.
%F a(n) ~ sqrt(5/3)*5^(5*n/2)/(2*Pi*n). - _Ilya Gutkovskiy_, Aug 07 2016
%F From _Peter Bala_, Mar 22 2022: (Start)
%F For n >= 1, a(n) = (5/3)*binomial(m*n,2*n)*binomial(m*n/2,2*n)*binomial(2*n,n)^2/ binomial(m*n/2,n)^2 at m = -1. See A352651 for the case m = 1.
%F a(n) = Sum_{k = 0..n} binomial(2*n-k-1,n-k)*binomial(3*n,k)^2.
%F a(n) = [x^n] (1 - x)^(2*n) * P(3*n,(1 + x)/(1 - x)), where P(n,x) denotes the n-th Legendre polynomial. Cf. A245086.
%F a(p) == a(1) (mod p^3) for prime p >= 5.
%F Conjecture: the supercongruences a(n*p^k) == a(n*p^(k-1)) (mod p^(3*k)) hold for all primes p >= 5 and all positive integers n and k. (End)
%F Row 1 of A365025. - _Peter Bala_, Aug 18 2023
%p seq(simplify(factorial(3*n)*factorial(n/2)*factorial(5*n/2)/(factorial(n)^3*factorial(3*n/2)^2)), n = 0 .. 20);
%t Table[Binomial[3 n, 3 n/2] Binomial[2 n, n] Binomial[5 n/2, n/2] / Binomial[n, n/2], {n, 0, 16}] (* _Michael De Vlieger_, Aug 07 2016 *)
%o (PARI) a(n) = sum(k = 0, n, binomial(2*n-k-1,n-k)*binomial(3*n,k)^2); \\ _Michel Marcus_, Apr 21 2022
%o (Python)
%o from math import factorial
%o from sympy import factorial2
%o def A275652(n): return int(factorial(3*n)*factorial2(5*n)*factorial2(n)//factorial2(3*n)**2//factorial(n)**3) # _Chai Wah Wu_, Aug 08 2023
%Y Cf. A002894, A245086, A275653, A275654, A275655, A276098, A276100, A276101, A276102, A352651, A352652, A365025.
%K nonn,easy
%O 0,2
%A _Peter Bala_, Aug 04 2016