OFFSET
0,3
COMMENTS
G.f. A(x) satisfies x=(1+4*A(x))A(A(x)).
FORMULA
a(n)=T(n,1), T(n,m)=1/2*(sum(k=1..n-m, 4^k*T(n-m,k)*binomial(k+m-1,m-1)*(-1)^(k))-sum(k=m+1..n-1, T(n,k)*T(k,m))), n>m, T(n,n)=1. [Vladimir Kruchinin, May 04 2012]
PROG
(PARI) a(n)=local(A); if(n<1, 0, A=x; for(k=1, n, A=Pol(A+serreverse(A+x*O(x^k))/(1+4*x))/2); polcoeff(A, n))
(Maxima) T(n, m):=if n=m then 1 else 1/2*(sum(4^k*T(n-m, k)*binomial(k+m-1, m-1)*(-1)^(k), k, 1, n-m)-sum(T(n, k)*T(k, m), k, m+1, n-1)); makelist(T(n, 1), n, 1, 10); [Vladimir Kruchinin, May 04 2012]
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Oct 04 2003
STATUS
approved