OFFSET
0,4
COMMENTS
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..200
FORMULA
The sequence can be obtained by shifting down triangle A084938 one row, then inserting a "1" as the first term, obtaining a new triangle Q. Then take the Lim_{n->inf.} Q^n, obtaining a shifted left column vector [1, 1, 1, 2, 6, 23, 105, 550, 3236,...].
PROG
(PARI)
EigenSeq(nn, tf)={my(v=vector(nn+1)); v[1] = 1; for(n=2, #v, v[n] = sum(k=0, n-2, tf(n-2, k)*v[k+1]); ); v}
T084938(n)={my(v=vector(n+1)); v[1] = [1]; for(n=1, n, v[1+n]=vector(1+n, k, k--; if(k>0, sum(j=0, n-k, j!*v[n-j][k])))); v}
seq(nn) = { my(T=T084938(nn-1)); EigenSeq(nn, (n, k)->T[1+n][1+k]) } \\ Andrew Howroyd, Sep 23 2025
CROSSREFS
KEYWORD
eigen,nonn
AUTHOR
Gary W. Adamson, Sep 20 2009
EXTENSIONS
a(0)=1 prepended and a(11) onwards from Andrew Howroyd, Sep 23 2025
STATUS
approved
