0,2
(PARI) {a(n)=local(M=matrix(n+1, n+1, r, c, if(r>=c, binomial((r-1)*r*(r+1)/3!-(c-1)*c*(c+1)/3!, r-c))), N=matrix(n+1, n+1, r, c, if(r>=c, binomial((r-1)*r*(r+1)/3!-(c-1)*c*(c+1)/3!+1, r-c)))); sum(k=0, n, (N*M^-1)[n+1, k+1])}
Cf. A126465, A126466, A126467, A126468.
Sequence in context: A185052 A009457 A175172 * A054859 A076725 A059917
Adjacent sequences: A126466 A126467 A126468 * A126470 A126471 A126472
nonn
Paul D. Hanna (pauldhanna(AT)juno.com), Dec 27 2006
