OFFSET
0,2
LINKS
PROG
(Python)
from sympy import mobius, prime
def a053645(n): return n - 2**(n.bit_length()-1)
def b(n): return n + 1 if n<2 else prime(1 + (len(bin(n)[2:]) - bin(n)[2:].count("1"))) * b(a053645(n))
def a(n): return mobius(b(n))%3 # Indranil Ghosh, May 30 2017
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, May 30 2017
STATUS
approved
