login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A059259 Triangle read by rows giving coefficient T(i,j) of x^i y^j in 1/(1-x-x*y-y^2) = 1/((1+y)(1-x-y)) for (i,j) = (0,0), (1,0), (0,1), (2,0), (1,1), (0,2), ... 14
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(n-1)+k(k+1)a(n-2), 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/(1-x^2) and subsequent rows are successive partial sums given by 1/((1-x)^n)(1-x^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/(1-x-x*y-y^2).

As a square array read by antidiagonals, this is T(n, k) = sum{i=0..n, (-1)^(n-i)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(n-1,k) + T(n-2,k-1) + T(n-2,k-2), 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(n-1,k) + T(n-1,k-1) 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/(1-x-x*y-y^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(n-1, k-1)-sum(prec(n, k+i-1) for i in (2..n-k+1))

    return [(-1)^(n-k+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

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

License Agreements, Terms of Use, Privacy Policy .

Last modified June 27 11:25 EDT 2017. Contains 288788 sequences.