login
A113093
Main diagonal of square table A113092; also, a(n) equals the n-th term in column 0 of the matrix n-th power of triangle A113095.
1
1, 1, 13, 693, 146821, 124626530, 426524622399, 5893207147435867, 328422072384464274577, 73719657441008064407836359, 66567306698774377126527799872190
OFFSET
0,3
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((M^n)[n+1, 1])}
CROSSREFS
Sequence in context: A012746 A211094 A290166 * A195576 A195555 A108223
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 14 2005
STATUS
approved