login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A113094
Diagonal of square table A113092; also, a(n) equals the n-th term in column 0 of the matrix (n+1)-th power of triangle A113095.
1
1, 2, 27, 1504, 330745, 289031301, 1011348629263, 14213347986246578, 802722082112213275116, 182118530044524172384716760, 165892108866362877173717099499469
OFFSET
1,2
PROG
(PARI) {a(n, q=4)=local(M=matrix(n+1, n+1)); for(r=1, n+1, for(c=1, r, M[r, c]=if(r==c, 1, if(c>1, (M^q)[r-1, c-1])+(M^q)[r-1, c]))); return(if(n<1, 0, (M^n)[n, 1]))}
CROSSREFS
Sequence in context: A067075 A015217 A320417 * A327128 A320669 A299116
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 14 2005
STATUS
approved