OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..20000
MAPLE
q:= m-> is(m mod ilog2(m)=1):
select(q, [$2..361])[]; # Alois P. Heinz, Jul 03 2026
PROG
(Python)
from oeis_sequences.OEISsequences import bisection
def A112251(n):
def f(x):
c, l, a = n+2+x, x.bit_length(), 2
for i in range(1, l-1):
c += (a-2)//i
a <<= 1
c -= (a-2)//i
if l>1:
c -= (x-1)//(l-1)-(a-2)//(l-1)
return c
return bisection(f, n, n) # Chai Wah Wu, Jul 02 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Aug 30 2005
EXTENSIONS
Name clarified by Sean A. Irvine, Jul 02 2026
STATUS
approved