OFFSET
0,6
COMMENTS
a(0) = 0 (an empty sum).
LINKS
FORMULA
EXAMPLE
PROG
(Python)
from sympy import factorint, prime, log
import math
def wt(n): return bin(n).count("1")
def a037445(n):
f=factorint(n)
return 2**sum([wt(f[i]) for i in f])
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 a286575(n): return a037445(b(n))
def a(n): return int(math.floor(log(a286575(n), 2))) # Indranil Ghosh, May 30 2017
(Python)
# uses RLT function from A278159
def A286574(n): return len(bin(RLT(n, lambda m: 2**(bin(m).count('1')))))-3 # Chai Wah Wu, Feb 04 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, May 28 2017
STATUS
approved