login
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