OFFSET
1,3
FORMULA
a(n, k) = 2^max{0, n-1} * sum_{d|k} (2^d - 1) * mu(k/d)
EXAMPLE
The binary expansion of 7/24 = 0.010(01)... has 3 preperiodic bits (to the right of the binary point) followed by 2 periodic (i.e., repeating) bits, while 1/2 = 0.1(0)... has one bit of each type. The preperiodic and periodic parts are both chosen to be as short as possible.
a(2, 2) = |{1/12 = 0.00(01)..., 5/12 = 0.01(10)..., 7/12 = 0.10(01)..., 11/12 = 0.11(10)...}| = 4
Table begins:
1 2 6 12
1 2 6 12
2 4 12 24
4 8 24 48
MATHEMATICA
Table[2^ Max[0, n-1](Plus@@((2^Divisors[k]-1)MoebiusMu[k/Divisors[k]])), {n, 0, 1 0}, {k, 1, 10}]
CROSSREFS
KEYWORD
AUTHOR
Brad Chalfan (brad(AT)chalfan.net), May 28, 2006
STATUS
approved