Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A060739
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A060739 (-1)^(n(n-1)/2) * Product_{k=0,...,n-1} (n+k-1)!/((k!)^2 * (n-1-k)!). +0
4
1, 1, -2, -36, 7200, 17640000, -560105280000, -239102222768640000, 1408147589778024775680000, 116620600756651855983415296000000, -137839975629646325813680872620851200000000, -2352568589682795058651211199786427114330521600000000 (list; graph; listen)
OFFSET

0,3

COMMENT

Let A_n be the matrix of size n X n defined by: A_n[i,j] = 1/(binomial coefficient i+j-2 over i-1) = 1/C(i+j-2,i-1) where 1 <= i,j <= n. The diagonals of this matrix are the reciprocals of the entries in the Pascal triangle. Then a(n) = 1/det(A_n) = det((A_n)^(-1)).

From the formula for a(n) it follows that the determinant of (A_n)^(-1) is an integer. By inspecting the values of (A_n)^(-1) for small values of n it looks like (A_n)^(-1) is actually a matrix of integers but I do not have a proof of this fact.

Let M_n be the n X n matrix with M_n(i,j)=i/(i+j); then |a(n-1)|=1/det(M_n). - Benoit Cloitre (benoit7848c(AT)orange.fr), Apr 21 2002

Also related to the multinomial coefficients (i+j)!/i!/j! : abs(a(n))=(1/detQ_n-1) where Q_n is the n X n matrix q(i,j)=i!j!/(i+j)! - Benoit Cloitre (benoit7848c(AT)orange.fr), May 30 2002

Contribution from Alexander Adamchuk (alex(AT)kolmogorov.com), Nov 14 2009: (Start)

Also a(n) = (-1)^(n(n-1)/2) * Product[ Binomial[2k,k]^2/2, {k,1,n-1} ].

It is simpler definition of a(n).

It follows from the observation that Sqrt[ Abs[ a(n+1)/a(n)/2 ] ] = {1, 3, 10, 35, 126, 462, ...} = C(2n+1, n+1) = A001700. (End)

REFERENCES

Doron Zeilberger, Reverend Charles to the aid of Major Percy and Fields-Medalist Enrico, Amer. Math. Monthly 103 (1996), 501-502.

LINKS

Harry J. Smith, Table of n, a(n) for n=0,...,43

FORMULA

If Multinomial[a, b, c] denotes the multinomial coefficient (a+b+c)! / (a! * b! * c!) (which is an integer) then : a(n) = (-1)^(n(n-1)/2) * Product k=0, ..., n-1 Multinomial[k, k, n-1-k] = (-1)^(n(n-1)/2) * product k=0, ..., n-1 (n+k-1)!/((k!)^2 * (n-1-k)!)

a(n) = (-1)^(n(n-1)/2) * Product[ Binomial[2k,k]^2/2, {k,1,n-1} ]. [From Alexander Adamchuk (alex(AT)kolmogorov.com), Nov 14 2009]

EXAMPLE

Here is the matrix A_4 for n=4: [1, 1, 1, 1; 1, 1/2, 1/3, 1/4; 1, 1/3, 1/6, 1/10; 1, 1/4, 1/10, 1/20]; a(4) = 7200 because det(A_4) = 1/7200

MAPLE

A060739 := n->(-1)^(n*(n-1)/2) * mul( (n+k-1)!/((k!)^2 * (n-1-k)!), k=0..n-1);

PROGRAM

(PARI) for(n=1, 15, print1(1/matdet(matrix(n, n, i, j, i/(j+i))), ", ")) [See Cloitre's comment]

(PARI) { for (n=0, 43, if (n<2, a=1, a=(-1)^(n\2)/matdet(matrix(n-1, n-1, i, j, i/(j+i)))); write("b060739.txt", n, " ", a); ) } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Jul 10 2009]

CROSSREFS

Cf. A005249, A067689.

Cf. A001700. [From Alexander Adamchuk (alex(AT)kolmogorov.com), Nov 14 2009]

Sequence in context: A095229 A047832 A004003 this_sequence A134366 A127234 A051459

Adjacent sequences: A060736 A060737 A060738 this_sequence A060740 A060741 A060742

KEYWORD

easy,sign,nice

AUTHOR

Noam Katz (noamkj(AT)hotmail.com), Apr 25 2001

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | The OEIS Foundation | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified March 20 09:10 EDT 2010. Contains 173642 sequences.


AT&T Labs Research