login
A286576
a(n) = A132971(n) mod 3.
2
1, 2, 2, 0, 2, 1, 0, 0, 2, 1, 1, 0, 0, 0, 0, 0, 2, 1, 1, 0, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 1, 0, 1, 2, 0, 0, 1, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 1, 0, 1, 2, 0, 0, 1, 2, 2, 0, 0, 0, 0, 0, 1, 2, 2, 0, 2, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
OFFSET
0,2
FORMULA
a(n) = A010872(A132971(n)) = A132971(n) mod 3.
PROG
(Scheme) (define (A286576 n) (modulo (A132971 n) 3))
(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))
def a(n): return mobius(b(n))%3 # Indranil Ghosh, May 30 2017
CROSSREFS
Sequence in context: A193863 A363566 A273496 * A322523 A333210 A285193
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, May 30 2017
STATUS
approved