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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A135886 Column 1 of triangle Q = A135885; also equals column 0 of Q^2. 3
1, 4, 20, 126, 980, 9186, 101492, 1296934, 18868652, 308478492, 5605768476, 112198139500, 2454071216496, 58267971181456, 1493114371576942, 41084194594171729, 1208473333806735096, 37849717704435895370 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

LINKS

Table of n, a(n) for n=0..17.

EXAMPLE

Triangle Q = A135885 begins:

1;

2, 1;

6, 4, 1;

25, 20, 6, 1;

138, 126, 42, 8, 1;

970, 980, 351, 72, 10, 1;

8390, 9186, 3470, 748, 110, 12, 1; ...

where column k of Q equals column 0 of Q^(k+1) such that

column 0 of Q equals column 0 of P=A135880 shift left and Q=P^2.

PROG

(PARI) {a(n)=local(P=Mat(1), R, PShR); if(n==0, 1, for(i=0, n+1, PShR=matrix(#P, #P, r, c, if(r>=c, if(r==c, 1, if(c==1, 0, P[r-1, c-1])))); R=P*PShR; R=matrix(#P+1, #P+1, r, c, if(r>=c, if(r<#P+1, R[r, c], if(c==1, (P^2)[ #P, 1], (P^(2*c-1))[r-c+1, 1])))); P=matrix(#R, #R, r, c, if(r>=c, if(r<#R, P[r, c], (R^c)[r-c+1, 1])))); (P^2)[n+2, 2])}

CROSSREFS

Cf. A135885; other columns: A135881, A135887.

Sequence in context: A151341 A285868 A266490 * A007550 A080795 A126674

Adjacent sequences:  A135883 A135884 A135885 * A135887 A135888 A135889

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Dec 15 2007

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 May 12 10:54 EDT 2021. Contains 343821 sequences. (Running on oeis4.)