OFFSET
0,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,1).
FORMULA
Euler transform of length 10 sequence [0, 1, -1, 0, 1, 0, 0, 1, 0, -1].
Moebius transform is length 8 sequence [1, -1, 0, -1, 0, 0, 0, 1].
a(n) is multiplicative with a(2) = 0, a(4) = -1, a(2^e) = 0 if e > 2, a(p^e) = 1 if p > 2.
G.f.: x / (1 - x^2) - x^4 / (1 - x^8).
G.f.: (x + x^3 - x^4 + x^5 + x^7) / (1 - x^8).
G.f.: x * (1 - x^2) * (1 - x^5) * (1 - x^8) / ((1 - x^3) * (1 - x^10)).
G.f.: x / (1 - x) - x^2 / (1 - x^2) - x^4 / (1 - x^4) + x^8 / (1 - x^8).
G.f. A(x) = Sum_{k>0} F(x^k) where P(x) = x - x^2 - x^4 + x^8 = x * (1 - x^3) * (1 - x - x^4).
a(n) = a(-n) = a(n+8) for all n in Z.
a(n) = -(-1)^n * A259044(n).
a(2*n + 1) = 1. a(4*n + 2) = 0. a(8*n) = 0. a(8*n + 4) = -1.
Dirichlet g.f.: zeta(s)*(2^s-1)^2*(2^s+1)/8^s. - Amiram Eldar, Dec 29 2022
EXAMPLE
G.f. = x + x^3 - x^4 + x^5 + x^7 + x^9 + x^11 - x^12 + x^13 + x^15 + ...
MATHEMATICA
a[ n_] := Mod[n, 2] - Boole[ Mod[n, 8] == 4];
a[ n_] := {0, 1, 0, 1, -1, 1, 0, 1}[[Mod[n, 8] + 1]];
a[ n_] := SeriesCoefficient[ x / (1 - x^2) - x^4 / (1 - x^8), {x, 0, Abs@n}];
PadRight[{}, 120, {0, 1, 0, 1, -1, 1, 0, 1}] (* Harvey P. Dale, May 11 2019 *)
PROG
(PARI) {a(n) = n%2 - (n%8==4)};
(PARI) {a(n) = [0, 1, 0, 1, -1, 1, 0, 1][n%8 + 1]};
(PARI) {a(n) = 1 - (n%2==0) - (n%4==0) + (n%8==0)};
(PARI) {a(n) = polcoeff( x / (1 - x^2) - x^4 / (1 - x^8) + x * O(x^abs(n)), abs(n))};
(Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(x*(1 +x^2-x^3+x^4+x^6)/(1-x^8))); // G. C. Greubel, Jul 29 2018
CROSSREFS
KEYWORD
sign,mult,easy
AUTHOR
Michael Somos, Dec 29 2016
STATUS
approved