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!)
A104732 Square array T[i,j]=T[i-1,j]+T[i-1,j-1], T[1,j]=j, T[i,1]=1, read by antidiagonals. 2

%I #23 Apr 26 2020 15:20:31

%S 1,2,1,3,3,1,4,5,4,1,5,7,8,5,1,6,9,12,12,6,1,7,11,16,20,17,7,1,8,13,

%T 20,28,32,23,8,1,9,15,24,36,48,49,30,9,1,10,17,28,44,64,80,72,38,10,1,

%U 11,19,32,52,80,112,129,102,47,11,1,12,21,36,60,96,144,192,201,140,57,12,1

%N Square array T[i,j]=T[i-1,j]+T[i-1,j-1], T[1,j]=j, T[i,1]=1, read by antidiagonals.

%C Original definition was "Triangle, row sums are A001924". Reading the rows of the triangle as antidiagonals of a square array allows a precise, yet simple definition and a method for computing the terms. - _M. F. Hasler_, Apr 26 2008

%C When formatted as a triangle, row sums are A001924: 1, 3, 7, 14, 26...(apply the partial sum operator twice to the Fibonacci sequence).

%H Charles R Greathouse IV, <a href="/A104732/b104732.txt">Table of n, a(n) for n = 1..10000</a>

%F The triangle is extracted from A * B; where A = [1; 2, 1; 3, 2, 1;...], B = [1; 0, 1; 0, 1, 1; 0, 0, 2, 1;...]; both infinite lower triangular matrices with the rest of the terms zeros. The sequence in "B" (1, 0, 1, 0, 1, 1, 0, 0, 2, 1...) = A026729.

%F As a square array, g.f. Sum T[i,j] x^j y^i = xy/((1-(1+x)y)*(1-x)^2). - Alec Mihailovs (alec(AT)mihailovs.com), Apr 26 2008

%e The first few rows of the triangle (= rising diagonals of the square array) are:

%e 1;

%e 2, 1;

%e 3, 3, 1;

%e 4, 5, 4, 1;

%e 5, 7, 8, 5, 1;

%e 6, 9, 12, 12, 6, 1;

%e ...

%p A104732 := proc(i,j) coeftayl(coeftayl(x*y/(1-x)^2/(1-y*(1+x)),y=0,i),x=0,j) ; end: for d from 1 to 20 do for j from d to 1 by -1 do printf("%d,",A104732(d-j+1,j)) ; od: od: # _R. J. Mathar_, May 04 2008

%t nn = 10; Map[Select[#, # > 0 &] &,Drop[CoefficientList[

%t Series[y x/(1 - x - y x + y x^3)/(1 - x), {x, 0, nn}], {x, y}],

%t 1]] // Grid (* _Geoffrey Critzer_, Mar 17 2015 *)

%o (Python)

%o def A104732_rows(n):

%o """Produces n rows of A104732 triangle"""

%o from operator import iadd

%o a,b,c = [], [1], [1]

%o for i in range(2,n+1):

%o a,b = b, [i]+list(map(iadd,a,b[:-1]))+[1]

%o c+=b

%o return c

%o # Alec Mihailovs (alec(AT)mihailovs.com), May 04 2008

%Y Cf. A001924, A026729.

%K nonn,tabl

%O 1,2

%A _Gary W. Adamson_, Mar 20 2005

%E Edited by _M. F. Hasler_, Apr 26 2008

%E More terms from _R. J. Mathar_ and Alec Mihailovs (alec(AT)mihailovs.com), May 04 2008

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 April 30 23:08 EDT 2024. Contains 372141 sequences. (Running on oeis4.)