The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A203421 Reciprocal of Vandermonde determinant of (1,1/2,...,1/n). 10
1, 1, -2, -18, 1152, 720000, -5598720000, -658683809280000, 1381360067999170560000, 59463021447701323327733760000, -59463021447701323327733760000000000000, -1542317635347398938581016812202229760000000000000 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Each term divides its successor, as in A000169.
LINKS
FORMULA
G.f.: G(0)/(2*x) -1/x, where G(k)= 1 + 1/(1 - x/(x + (2*k+1)/((2*k+1)^(2*k+1))/(1 + 1/(1 - x/(x - (2*k+2)/((2*k+2)^(2*k+2))/G(k+1)))))); (continued fraction). - Sergei N. Gladkovskii, Jun 03 2013
a(n) = (-1)^floor(n/2) * hyperfactorial(n)/n! = A057077(n) * A002109(n)/n!. - Paul J. Harvey, Feb 08 2014
a(n) = Product_{i=2..n} (-i)^(i-1). - Kevin Ryde, Apr 17 2022
abs(a(n)) ~ A * n^(n*(n-1)/2 - 5/12) / (sqrt(2*Pi) * exp(n^2/4 - n)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Nov 20 2023
a(n) = (-1)^binomial(n,2) * (n!)^n / BarnesG(n+2). - G. C. Greubel, Dec 07 2023
MATHEMATICA
(* First program *)
f[j_] := 1/j; z = 12;
v[n_] := Product[Product[f[k] - f[j], {j, 1, k - 1}], {k, 2, n}]
Table[v[n], {n, 1, z}]
1/% (* A203421 *)
Table[v[n]/v[n + 1], {n, 1, z}] (* A000169 signed *)
(* Additional programs *)
Table[(-1)^Floor[n/2]*Product[(k + 1)^k, {k, 0, n-1}], {n, 1, 10}] (* Vaclav Kotesovec, Oct 18 2015 *)
Table[(-1)^Binomial[n, 2]*(n!)^n/BarnesG[n+2], {n, 20}] (* G. C. Greubel, Dec 07 2023 *)
PROG
(PARI) a(n) = prod(i=2, n, (-i)^(i-1)); \\ Kevin Ryde, Apr 17 2022
(Magma)
BarnesG:= func< n | (&*[Factorial(k): k in [0..n-2]]) >;
A203421:= func< n | (-1)^Binomial(n, 2)*(Factorial(n))^n/BarnesG(n+2) >;
[A203421(n): n in [1..20]]; // G. C. Greubel, Dec 07 2023
(SageMath)
def BarnesG(n): return product(factorial(k) for k in range(n-1))
def A203421(n): return (-1)^binomial(n, 2)*(gamma(n+1))^n/BarnesG(n+2)
[A203421(n) for n in range(1, 21)] # G. C. Greubel, Dec 07 2023
CROSSREFS
Sequence in context: A268051 A268074 A202544 * A139111 A090766 A069651
KEYWORD
sign,easy
AUTHOR
Clark Kimberling, Jan 02 2012
EXTENSIONS
a(0)=1 prepended by Alois P. Heinz, Apr 13 2024
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified May 19 14:26 EDT 2024. Contains 372697 sequences. (Running on oeis4.)