OFFSET
0,2
COMMENTS
Given g.f. A(x), note that A(x)^(1/5) does not yield an integer series.
Compare to: exp( Sum_{n>=1} Lucas(n)*x^n/n ) = 1/(1-x-x^2) where Lucas(n) = Fibonacci(n-1) + Fibonacci(n+1).
FORMULA
G.f.: 1 / ( (1+x)^4 * (1 - 3*x + x^2)^3 ).
EXAMPLE
G.f.: 1 + 5*x + 25*x^2 + 100*x^3 + 380*x^4 + 1348*x^5 + 4610*x^6 +...
where
log(A(x))/5 = x + 5*x^2/2 + 10*x^3/3 + 29*x^4/4 + 73*x^5/5 + 194*x^6/6 + 505*x^7/7 + 1325*x^8/8 +...+ A069921(n-1)*x^n/n +...
PROG
(PARI) {L(n)=fibonacci(n-1)^2+fibonacci(n+1)^2}
{a(n)=polcoeff(exp(sum(m=1, n, 5*L(m)*x^m/m)+x*O(x^n)), n)}
for(n=0, 30, print1((a(n)), ", "))
(PARI) {a(n)=polcoeff(1/((1+x)^4*(1-3*x+x^2)^3+x*O(x^n)), n)}
for(n=0, 30, print1((a(n)), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Apr 05 2013
STATUS
approved