|
|
A132971
|
|
a(2*n) = a(n), a(4*n+1) = -a(n), a(4*n+3) = 0, with a(0) = 1.
|
|
8
|
|
|
1, -1, -1, 0, -1, 1, 0, 0, -1, 1, 1, 0, 0, 0, 0, 0, -1, 1, 1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 1, 0, 1, -1, 0, 0, 1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 1, 0, 1, -1, 0, 0, 1, -1, -1, 0, 0, 0, 0, 0, 1, -1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
If binary(n) has adjacent 1 bits then a(n) = 0 else a(n) = (-1)^A000120(n).
|
|
LINKS
|
|
|
FORMULA
|
A024490(n) = number of solutions to 2^n <= k < 2^(n+1) and a(k) = 1.
A005252(n) = number of solutions to 2^n <= k < 2^(n+1) and a(k) = -1.
A027935(n-1) = number of solutions to 2^n <= k < 2^(n+1) and a(k) = 0.
G.f. A(x) satisfies A(x) = A(x^2) - x * A(x^4).
G.f. B(x) of A000621 satisfies B(x) = x * A(x^2) / A(x).
|
|
EXAMPLE
|
G.f. = 1 - x - x^2 - x^4 + x^5 - x^8 + x^9 + x^10 - x^16 + x^17 + x^18 + ...
|
|
MATHEMATICA
|
m = 100; A[_] = 1;
Do[A[x_] = A[x^2] - x A[x^4] + O[x]^m // Normal, {m}];
|
|
PROG
|
(PARI) {a(n) = if( n<1, n==0, if( n%2, if( n%4 > 1, 0, -a((n-1)/4) ), a(n/2) ) )};
(PARI) {a(n) = my(A, m); if( n<0, 0, m = 1; A = 1 + O(x); while( m<=n, m *= 2; A = subst(A, x, x^2) - x * subst(A, x, x^4) ); polcoeff(A, n)) };
(Python)
from sympy import mobius, prime, log
import math
def A(n): return n - 2**int(math.floor(log(n, 2)))
def b(n): return n + 1 if n<2 else prime(1 + (len(bin(n)[2:]) - bin(n)[2:].count("1"))) * b(A(n))
|
|
CROSSREFS
|
Cf. A085357 (gives the absolute values: -1 -> 1), A286576 (when reduced modulo 3: -1 -> 2).
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|