login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A113082
Main diagonal of square table A113081; also, a(n) equals the n-th term in column 0 of the matrix n-th power of triangle A113084.
0
1, 1, 10, 331, 33476, 10204145, 9378590446, 26026690264407, 218132378185337416, 5518274388618175447069, 421034872020570533423509010, 96809747319527667989371938562883
OFFSET
0,3
PROG
(PARI) {a(n, q=3)=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: A257132 A218996 A352798 * A288684 A046747 A338799
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 14 2005
STATUS
approved