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!)
A107727 Matrix inverse of A107719. 3
1, -1, 1, -3, -2, 1, -21, -7, -3, 1, -219, -53, -13, -4, 1, -2973, -583, -115, -21, -5, 1, -49323, -8249, -1437, -217, -31, -6, 1, -964173, -141655, -22715, -3101, -369, -43, -7, 1, -21680571, -2853185, -430877, -55251, -5975, -581, -57, -8, 1, -551173053, -65887783, -9505707, -1168349, -119137 (list; table; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

Matrix square is A107728. Matrix cube is A107726. Column 0 is negative A107716 shift right.

LINKS

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

EXAMPLE

Triangle begins:

1;

-1,1;

-3,-2,1;

-21,-7,-3,1;

-219,-53,-13,-4,1;

-2973,-583,-115,-21,-5,1;

-49323,-8249,-1437,-217,-31,-6,1;

-964173,-141655,-22715,-3101,-369,-43,-7,1; ...

PROG

(PARI) {T(n, k)=local(L, N, M=matrix(n+1, n+1, m, j, if(m>=j, if(m==j, 1, if(m==j+1, -3*j, polcoeff(1/sum(i=0, m-j, prod(r=0, i-1, 3*r+1)*x^i)+O(x^m), m-j)))))^-1); L=sum(i=1, #M, (M^0-M)^i/i)/3; N=sum(i=0, #L, L^i/i!); return(if(n<0, 0, N[n+1, k+1]))}

CROSSREFS

Cf. A107716, A107719, A107726, A107728.

Sequence in context: A291080 A107862 A117265 * A346743 A087041 A152790

Adjacent sequences:  A107724 A107725 A107726 * A107728 A107729 A107730

KEYWORD

sign,tabl

AUTHOR

Paul D. Hanna, May 30 2005

STATUS

approved

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 May 26 07:48 EDT 2022. Contains 354077 sequences. (Running on oeis4.)