OFFSET
0,2
COMMENTS
Radius of convergence of g.f. A(x) is r = 2^4/3^6 = 16/729.
EXAMPLE
G.f.: A(x) = 1 + 3*x + 72*x^2 + 2208*x^3 + 75531*x^4 + 2748957*x^5 + 104125542*x^6 + 4055630148*x^7 + 161248468944*x^8 + 6513248563281*x^9 + ...
such that
A(x)^2 + A(x) - 1 = 1 + 9*x + 225*x^2 + 7056*x^3 + 245025*x^4 + 9018009*x^5 + 344622096*x^6 + 13521038400*x^7 + 540917591841*x^8 + 21966328580625*x^9 + ... + binomial(3*n,n)^2 * x^n + ...
PROG
(PARI) {S(n) = sum(m=0, n, binomial(3*m, m)^2 * x^m ) +x*O(x^n)}
{A(n) = (sqrt(4*S(n) + 5) - 1)/2 }
{a(n) = polcoeff( A(n), n)}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 07 2018
STATUS
approved