

A059259


Triangle read by rows giving coefficient T(i,j) of x^i y^j in 1/(1xx*yy^2) = 1/((1+y)(1xy)) for (i,j) = (0,0), (1,0), (0,1), (2,0), (1,1), (0,2), ...


13



1, 1, 0, 1, 1, 1, 1, 2, 2, 0, 1, 3, 4, 2, 1, 1, 4, 7, 6, 3, 0, 1, 5, 11, 13, 9, 3, 1, 1, 6, 16, 24, 22, 12, 4, 0, 1, 7, 22, 40, 46, 34, 16, 4, 1, 1, 8, 29, 62, 86, 80, 50, 20, 5, 0, 1, 9, 37, 91, 148, 166, 130, 70, 25, 5, 1, 1, 10, 46, 128, 239, 314, 296, 200
(list;
table;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

0,8


COMMENTS

This sequence provides the general solution to the recurrence a(n)=a(n1)+k(k+1)a(n2), a(0)=a(1)=1. The solution is (1,1,k^2+k+1,2k^2+2k+1, .....) whose coefficients can be read from the rows of the triangle. The row sums of the triangle are given by the case k=1. These are the Jacobsthal numbers, A001045. Viewed as a square array, its first row is (1,0,1,0,1,...) with E.g.f. cosh(x), G.f. 1/(1x^2) and subsequent rows are successive partial sums given by 1/((1x)^n)(1x^2)).  Paul Barry, Mar 17 2003
Conjecture: every second column of this triangle is identical to a column in the square array A071921. For example, column 4 of A059259 (1, 3, 9, 22, 46, ...) appears to be the same as column 3 of A071921; column 6 of A059259 (1, 4, 16, 50, 130, 296, ...) appears to be the same as column 4 of A071921; and in general column 2k of A059259 appears to be the same as column k+1 of A071921. Furthermore, since A225010 is a transposition of A071921 (ignoring the latter's top row and two leftmost columns), there appears to be a correspondence between column 2k of A059259 and row k of A225010.  Mathew Englander, May 17 2014


LINKS

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened


FORMULA

G.f.: 1/(1xx*yy^2).
As a square array read by antidiagonals, this is T(n, k) = sum{i=0..n, (1)^(ni)C(i+k, k)}.  Paul Barry, Jul 01 2003
T(2*n,n) = A026641(n).  Philippe Deléham, Mar 08 2007
T(n,k) = T(n1,k) + T(n2,k1) + T(n2,k2), T(0,0) = T(1,0) = T(2,0) = T(2,1) = T(2,2)=1, T(1,1)=0, T(n,k)=0 if k<0 or if k>n.  Philippe Deléham, Nov 24 2013
T(n,0) = 1, T(n,n) = (1+(1)^n)/2, and T(n,k) = T(n1,k) + T(n1,k1) for 0 < k < n.  Mathew Englander, May 24 2014


EXAMPLE

Triangle begins:
1;
1,0;
1,1,1;
1,2,2,0;
1,3,4,2,1;
1,4,7,6,3,0;
1,5,11,13,9,3,1;
1,6,16,24,22,12,4,0;
1,7,22,40,46,34,16,4,1;
1,8,29,62,86,80,50,20,5,0;
1,9,37,91,148,166,130,70,25,5,1;
1,10,46,128,239,314,296,200,95,30,6,0;
...


MAPLE

read transforms; 1/(1xx*yy^2); SERIES2(%, x, y, 12); SERIES2TOLIST(%, x, y, 12);


MATHEMATICA

T[n_, 0] := 1; T[n_, n_] := (1 + (1)^n)/2; T[n_, k_] := T[n, k] = T[n  1, k] + T[n  1, k  1]; Table[T[n, k], {n, 0, 10} , {k, 0, n}]//Flatten (* G. C. Greubel, Jan 03 2017 *)


PROG

(Sage)
def A059259_row(n):
@cached_function
def prec(n, k):
if k==n: return (1)^n
if k==0: return 0
return prec(n1, k1)sum(prec(n, k+i1) for i in (2..nk+1))
return [(1)^(nk+1)*prec(n+1, k) for k in (1..n)]
for n in (1..12): print A059259_row(n) # Peter Luschny, Mar 16 2016


CROSSREFS

See A059260 for an explicit formula.
Diagonals of this triangle are given by A006498.
Similar to the triangles A035317, A080242, A108561, A112555.
Sequence in context: A167637 A109754 A220074 * A124394 A086460 A136431
Adjacent sequences: A059256 A059257 A059258 * A059260 A059261 A059262


KEYWORD

nonn,tabl


AUTHOR

N. J. A. Sloane, Jan 23 2001


STATUS

approved



