OFFSET
0,2
COMMENTS
EXAMPLE
Triangle A166884 begins:
1;
1, 1;
3, 2, 1;
15, 9, 3, 1;
114, 62, 18, 4, 1;
1159, 593, 157, 30, 5, 1;
14838, 7266, 1812, 316, 45, 6, 1;
229401, 108720, 25989, 4271, 555, 63, 7, 1; ...
of which the row sums form this sequence.
PROG
(PARI) {a(n, k)=local(F=x, M, N, P, m=max(n, k), A166884); M=matrix(m+2, m+2, r, c, F=x; for(i=1, r+c-2, F=subst(F, x, x+x^2+x^3+x*O(x^(m+2)))); polcoeff(F, c)); N=matrix(m+1, m+1, r, c, M[r, c]); P=matrix(m+1, m+1, r, c, M[r+1, c]); A166884=P~*(N~)^-1; sum(k=0, n, A166884[n+1, k+1])}
for(n=0, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Sep 13 2013
STATUS
approved