login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A089505
Triangle of signed numbers used for the computation of the column sequences of triangle A089504.
7
1, -1, 4, 1, -24, 50, -1, 114, -950, 1350, 31, -15504, 400520, -1897200, 2052855, -9269, 19612560, -1431859000, 17333030000, -56265334125, 49236404224, 342953, -3011508588, 594221485000, -16634292228125, 123422029355625, -302409994743808, 222337901418633, -9945637
OFFSET
1,3
COMMENTS
A089504(n+m,m)= sum(a(m,p)*((p+2)*(p+1)*p)^n,p=1..m)/D(m) with D(m) := A089506(m); m=1,2,..., n>=0.
FORMULA
a(n, m)= D(n)*((-1)^(n-m))*(((m+2)*(m+1)*m)^(n-1))/(product(fallfac(m+2, 3)-fallfac(r+2, 3), r=1..m-1)*product(fallfac(r+2, 3)-fallfac(m+2, 3), r=m+1..n)), with D(n) := A089506(n) and fallfac(n, m) := A008279(n, m) (falling factorials), 1<=m<=n else 0. (Replace in the denominator the first product by 1 if m=1 and the second one by 1 if m=n.)
a(n, m)= A089506(n)*((-1)^(n-m))*(fallfac(m+2, 3)^(n-1))*(3*m^2+6*m+2)/((n-m)!*(m-1)!*product(fallfac(m+r+2, 2)-r*m, r=1..n)), n>=m>=1.
EXAMPLE
[1]; [ -1,4]; [1,-24,50]; [ -1,114,-950,1350]; ...
a(3,2)= -24 = 27*(-1)*((4*3*2)^2)/((4*3*2-3*2*1)*(5*4*3-4*3*2)).
A089504(2+3,3) = A089513(2) = 6156 = (1*(3*2*1)^2 - 24*(4*3*2)^2 + 50*(5*4*3)^2)/27.
MATHEMATICA
b[n_, m_] := (-1)^(n - m)*FactorialPower[m + 2, 3]^(n - 1)/(Product[ FactorialPower[m + 2, 3] - FactorialPower[r + 2, 3], {r, 1, m - 1}] * Product[ FactorialPower[r + 2, 3] - FactorialPower[m + 2, 3], {r, m + 1, n}]); den[n_] := LCM @@ Table[ Denominator[b[n, m]], {m, 1, n}]; a[n_, m_] := den[n]*b[n, m]; Table[a[n, m], {n, 1, 10}, {m, 1, n}] // Flatten (* Jean-François Alcover, Sep 02 2016 *)
CROSSREFS
Companion denominator sequence is A089506.
Sequence in context: A285061 A285066 A046860 * A300083 A062328 A378199
KEYWORD
sign,easy,tabl
AUTHOR
Wolfdieter Lang, Dec 01 2003
STATUS
approved