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”).

A113083
Diagonal of square table A113081; also, a(n) equals the n-th term in column 0 of the matrix (n+1)-th power of triangle A113084.
0
1, 2, 21, 724, 75695, 23694838, 22239639177, 62747494950248, 532868670719193651, 13624738004791751175370, 1048678107774203901392276461, 242892250870416811233766661498812
OFFSET
1,2
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(if(n<1, 0, (M^n)[n, 1]))}
CROSSREFS
Sequence in context: A015209 A342267 A376689 * A174964 A210830 A139164
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 14 2005
STATUS
approved