0,2

(PARI) {a(n, q=5)=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]))}

Cf. A113103, A113106.

Adjacent sequences: A113102 A113103 A113104 this_sequence A113106 A113107 A113108

Sequence in context: A045984 A077657 A090335 this_sequence A132567 A083459 A034173

nonn

Paul D. Hanna (pauldhanna(AT)juno.com), Oct 14 2005