%I #10 Sep 28 2023 02:03:45
%S 1,-6,-486,839808,42515280000,-80335512599040000,
%T -6890065294166289123840000,31601087581187838970614157148160000,
%U 8925080517850366815864624583251321642024960000
%N Reciprocal of Vandermonde determinant of (1/3,1/6,...,1/(3n)).
%C Each term divides its successor, as in A203429.
%H G. C. Greubel, <a href="/A203428/b203428.txt">Table of n, a(n) for n = 1..33</a>
%F a(n) = (-3)^binomial(n,2) * (Gamma(n+1))^(n-1) / BarnesG(n+1). - _G. C. Greubel_, Sep 28 2023
%t (* First program *)
%t f[j_]:= 1/(3*j); z = 16;
%t v[n_]:= Product[Product[f[k] - f[j], {j,k-1}], {k,2,n}]
%t 1/Table[v[n], {n,z}] (* A203428 *)
%t Table[v[n]/(3*v[n+1]), {n,z}] (* A203429 *)
%t (* Second program *)
%t Table[(-3)^Binomial[n,2]*(Gamma[n+1])^(n-1)/BarnesG[n+1], {n,20}] (* _G. C. Greubel_, Sep 28 2023 *)
%o (Magma)
%o Barnes:= func< n | (&*[Factorial(j): j in [1..n-1]]) >;
%o A203428:= func< n | (-3)^Binomial(n,2)*(Factorial(n))^n/Barnes(n+1) >;
%o [A203428(n): n in [1..25]]; // _G. C. Greubel_, Sep 28 2023
%o (SageMath)
%o def barnes(n): return product(factorial(j) for j in range(n))
%o def A203428(n): return (-3)^binomial(n,2)*(factorial(n))^n/barnes(n+1)
%o [A203428(n) for n in range(1,21)] # _G. C. Greubel_, Sep 28 2023
%Y Cf. A203421, A203424, A203429.
%K sign
%O 1,2
%A _Clark Kimberling_, Jan 02 2012
|