Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A095149
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A095149 Triangle read by rows: Aitken's array (A011971) but with a leading diagonal before it given by the Bell numbers (A000110), 1,1,2,5,15,52,... +0
5
1, 1, 1, 2, 1, 2, 5, 2, 3, 5, 15, 5, 7, 10, 15, 52, 15, 20, 27, 37, 52, 203, 52, 67, 87, 114, 151, 203, 877, 203, 255, 322, 409, 523, 674, 877, 4140, 877, 1080, 1335, 1657, 2066, 2589, 3263, 4140, 21147, 4140, 5017, 6097, 7432, 9089, 11155, 13744, 17007, 21147 (list; table; graph; listen)
OFFSET

0,4

COMMENT

Or, prefix Aitken's array (A011971) with a leading diagonal of zeros and take the differences of each row to get the new triangle.

With offset 1, triangle read by rows: T(n,k) is the number of partitions of the set {1,2,...,n} in which k is the largest entry in the block containing 1 (1<=k<=n). For example - Emeric Deutsch (deutsch(AT)duke.poly.edu), Oct 29 2006

Row term sums = the Bell numbers starting with A000110(1): 1, 2, 5, 15...

FORMULA

With offset 1, T(n,1)=T(n,n)=T(n+1,2)=B(n-1)=A000110(n-1) (the Bell numbers). T(n,k)=T(n,k-1)+T(n-1,k-1) for n>=k>=3. T(n,n-1)=B(n-1)-B(n-2)=A005493(n-3) for n>=3 (B(q) are the Bell numbers A000110). T(n,k)=A011971(n-2,k-2) for n>=k>=2. In other words, deleting first row and first column we obtain Aitken's array A011971 (also called Bell or Pierce triangle; offset in A011971 is 0). - Emeric Deutsch (deutsch(AT)duke.poly.edu), Oct 29 2006

T(n,1)=B(n-1); T(n,2)=B(n-2) for n>=2; T(n,k)=Sum(binom(k-2,i)*B(n-2-i), i=0..k-2) for 3<=k<=n, where B(q) are the Bell numbers (A000110). Generating polynomial of row n is P[n](t)=Q[n](t,1), where Q[n](t,s)=t^n*Q[n-1](1,s)+s*dQ[n-1](t,s)/ds +(s-1) Q[n-1](t,s) ; Q[1](t,s)=ts. - Emeric Deutsch (deutsch(AT)duke.poly.edu), Oct 29 2006

EXAMPLE

Triangle starts:

1;

1,1;

2,1,2;

5,2,3,5;

15,5,7,10,15;

52,15,20,27,37,52;

MAPLE

A011971 := proc(n, k) option remember ; if k = 0 then if n=0 then 1; else A011971(n-1, n-1) ; fi ; else A011971(n, k-1)+A011971(n-1, k-1) ; fi ; end: A000110 := proc(n) option remember; if n<=1 then 1 ; else add( binomial(n-1, i)*A000110(n-1-i), i=0..n-1) ; fi ; end: A095149 := proc(n, k) option remember ; if k = 0 then A000110(n) ; else A011971(n-1, k-1) ; fi ; end: for n from 0 to 11 do for k from 0 to n do printf("%d, ", A095149(n, k)) ; od ; od ; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 05 2007

with(combinat): T:=proc(n, k) if k=1 then bell(n-1) elif k=2 and n>=2 then bell(n-2) elif k<=n then add(binomial(k-2, i)*bell(n-2-i), i=0..k-2) else 0 fi end: matrix(8, 8, T): for n from 1 to 11 do seq(T(n, k), k=1..n) od; # yields sequence in triangular form Q[1]:=t*s: for n from 2 to 11 do Q[n]:=expand(t^n*subs(t=1, Q[n-1])+s*diff(Q[n-1], s)-Q[n-1]+s*Q[n-1]) od: for n from 1 to 11 do P[n]:=sort(subs(s=1, Q[n])) od: for n from 1 to 11 do seq(coeff(P[n], t, k), k=1..n) od; # yields sequence in triangular form - Emeric Deutsch (deutsch(AT)duke.poly.edu), Oct 29 2006

CROSSREFS

Cf. A000110, A005493, A011971.

Sequence in context: A153206 A144155 A109631 this_sequence A064192 A124218 A025165

Adjacent sequences: A095146 A095147 A095148 this_sequence A095150 A095151 A095152

KEYWORD

nonn,tabl,easy,nice

AUTHOR

Gary W. Adamson (qntmpkt(AT)yahoo.com), May 30 2004

EXTENSIONS

Corrected and extended by R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 05 2007

Entry revised by N. J. A. Sloane (njas(AT)research.att.com), Jun 01 2005, Jun 16 2007

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 | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 27 22:38 EST 2009. Contains 167602 sequences.


AT&T Labs Research