OFFSET
1,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..1000
FORMULA
G.f.: Sum_{n>=1} x^n/(1 + Lucas(n)*x^n + (-1)^n*x^(2*n)) where Lucas(n) = A000204(n).
EXAMPLE
G.f.: A(x) = x + 3*x^3 - 5*x^4 + 6*x^5 - 3*x^6 + 14*x^7 - 48*x^8 + 52*x^9 +...
where A(x) = x/(1+x-x^2) + x^2/(1+3*x^2+x^4) + x^3/(1+4*x^3-x^6) + x^4/(1+7*x^4+x^8) + x^5/(1+11*x^5-x^10) + x^6/(1+18*x^6+x^12) +...+ x^n/(1 + Lucas(n)*x^n + (-1)^n*x^(2*n)) +...
Illustration of terms.
a(1) = 1; a(2) = 1*(-1/1 + 1/1) = 0; a(3) = 2*(1/1 + 1/2) = 3;
a(4) = 3*(-1/1 - 1/1 + 1/3) = -5; a(5) = 5*(1/1 + 1/5) = 6;
a(6) = 8*(-1/1 + 1/1 - 1/2 + 1/8) = -3; a(7) = 13*(1/1 + 1/13) = 14;
a(8) = 21*(-1/1 - 1/1 - 1/3 + 1/21) = -48; ...
PROG
(PARI) {a(n)=fibonacci(n) * sumdiv(n, d, -(-1)^(n/d) / fibonacci(d))}
(PARI) {Lucas(n)=fibonacci(n-1)+fibonacci(n+1)}
{a(n)=polcoeff(sum(m=1, n, x^m/(1 + Lucas(m)*x^m + (-1)^m*x^(2*m) +x*O(x^n))), n)}
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Jan 11 2012
STATUS
approved