|
|
A132228
|
|
a(n) = Sum_{k=0..n} C(n-1,k)^2*a(k)*a(n-k-1) for n>0 with a(0)=1.
|
|
2
|
|
|
1, 1, 2, 8, 52, 504, 6808, 122304, 2820048, 81183200, 2853990496, 120321094656, 5991955466560, 347996920977664, 23312947041336960, 1784445116557881344, 154767015393810489600, 15098457734490931766784
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Let A(x) = Sum_{n>=0} a(n) * x^n / n!^2. Then A(x)^2 = A'(x) + x * A''(x). - Michael Somos, May 20 2018
|
|
LINKS
|
|
|
PROG
|
(PARI) a(n)=if(n==0, 1, sum(k=0, n-1, a(k)*a(n-k-1)*binomial(n-1, k)^2 ))
(PARI) {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( k=0, n, A = 1 + intformal( intformal(A^2) / x)); n!^2 * polcoeff(A, n))}; /* Michael Somos, May 20 2018 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|