login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

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 56th year, we are closing in on 350,000 sequences, and we’ve crossed 9,700 citations (which often say “discovered thanks to the OEIS”).

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A191685 Eighth diagonal a(n) = s(n,n-7) of the unsigned Stirling numbers of the first kind with n>7 2
5040, 109584, 1172700, 8409500, 45995730, 206070150, 790943153, 2681453775, 8207628000, 23057159840, 60202693980, 147560703732, 342252511900, 756111184500, 1599718388730, 3256091103430, 6400590336096, 12191224980000, 22563937825000, 40681506808800 (list; graph; refs; listen; history; text; internal format)
OFFSET

8,1

COMMENTS

The Maple programs under I generate the sequence. The Maple program under II generates explicit formulas for a(n+1) = s(n+1,n+1-c) with c>=1 and n>=c.

REFERENCES

K. Seidel, Variation der Binomialkoeffizienten, Bild

der Wissenschaft, 6 (1980), pp. 127-128.

LINKS

T. D. Noe, Table of n, a(n) for n = 8..1000

FORMULA

a(n+1) = A130534(T(n,n-7)) = s(n+1,n+1-7)

a(n+1) = binomial(n+1,8)*(80*n+114*n^2-23*n^3-75*n^4-9*n^5+9*n^6)/144

EXAMPLE

c=1; a(n+1) = binomial(n+1,2)

c=2; a(n+1) = binomial(n+1,3)*(2+3*n)/4

c=3; a(n+1) = binomial(n+1,4)*(n+n^2)/2

c=4; a(n+1) = binomial(n+1,5)*(-8-10*n+15*n^2 +15*n^3)/48

c=5; a(n+1) = binomial(n+1,6)*(-6*n-7*n^2+2*n^3+ 3*n^4)/16

c=6; a(n+1) = binomial(n+1,7)*(96+140*n-224*n^2-315*n^3+63*n^5)/576

c=7; a(n+1) = binomial(n+1,8)*(80*n+114*n^2-23*n^3-75*n^4-9*n^5+9*n^6)/144

c=8; a(n+1) = binomial(n+1,9)*(-1152-1936*n+2820*n^2+

  5320*n^3+735*n^4-1575*n^5-315*n^6+135*n^7)/3840

c=9; a(n+1) = binomial(n+1,10)*(-1008*n-1676*n^2 +100*n^3+1295*n^4+392*n^5-210*n^6-60*n^7 +15*n^8)/768

MAPLE

I: programs generate the sequence:

with(combinat): c:=7; a:= proc(n) a(n):=abs(stirling1(n, n-c)); end: seq(a(n), n=c+1..28);

for n from 7 to 27 do a(n+1) := binomial(n+1, 8)*(80*n+ 114*n^2- 23*n^3- 75*n^4- 9*n^5+ 9*n^6)/144 end do: seq(a(n), n=8..28);

II: program generates explicit formulas for a(n+1) =  s(n+1, n+1-c):

k[1, 0]:=1: v:=1:

for c from 2 to 10 do

  c1:=c-1: c2:=c-2: p0:=0:

  for j from 0 to c2 do p0:=p0+k[c1, j]*m^j: end do:

  f:=expand(2*c*(m+1)*p0/v):

  p1:=0: p2:=0:

  for j from 0 to c1 do

    p1:=p1+k[c, j]*(m+1)^j:

    p2:=p2+k[c, j]*m^j:

  end do:

  g:=collect((m+2)*p1-(m-c1)*p2-f, m):

  kh[0]:=rem(g, m, m): Mk:=[kh[0]]: Mv:=[k[c, 0]]:

  for j from 1 to c1 do

    kh[j]:=coeff(g, m^j):

    Mk:=[op(Mk), kh[j]]: Mv:=[k[c, j], op(Mv)]:

  end do:

  sol:=solve(Mk, Mv):

  v:=1:

  for j from 1 to c do

    k[c, c-j]:=eval(k[c, c-j], sol[1, j]):

    nen[j]:=denom(k[c, c-j]):

    v:=ilcm(v, nen[j]):

  end do:

  for j from 0 to c1 do k[c, j]:=k[c, j]*v:

    printf("%8d", k[c, j]): end do:

  p3:=0:

  for j from 0 to c1 do p3:=p3+k[c, j]*n^j: end do:

  s[n+1, n+1-c]:=binomial(n+1, c+1)*(c+1)*p3/(2^c*k[c, c1]):

end do:

for c from 2 to 10 do print("%a\n", s[n+1, n+1-c]):

end do:

CROSSREFS

Cf. A130534, A000012 (c=0; 1st diagonal), A000217 (c=1; 2nd diagonal), A000914 (c=2; 3rd diagonal), A001303 (c=3; 4th diagonal), A000915 (c=4; 5th diagonal), A053567 (c=5; 6th diagonal), A112002 (c=6; 7th diagonal), A191685 (c=7; 8th diagonal).

Sequence in context: A061122 A029575 A179965 * A213877 A135456 A254080

Adjacent sequences:  A191682 A191683 A191684 * A191686 A191687 A191688

KEYWORD

nonn,easy

AUTHOR

Paul Weisenhorn, Jun 11 2011

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 6 01:46 EST 2021. Contains 349558 sequences. (Running on oeis4.)