login
Integers divisible by their number of digits in binary.
2

%I #26 May 19 2021 16:04:38

%S 1,2,6,8,12,20,25,30,36,42,48,54,60,70,77,84,91,98,105,112,119,126,

%T 128,136,144,152,160,168,176,184,192,200,208,216,224,232,240,248,261,

%U 270,279,288,297,306,315,324,333,342,351,360,369,378,387,396,405

%N Integers divisible by their number of digits in binary.

%H Daniel Arribas, <a href="/A177869/b177869.txt">Table of n, a(n) for n = 1..1000</a>

%e 105 is 1101001 in base 2 (length of 7); 105 / 7 is 15.

%t Select[Range[410], IntegerQ[#/Length[IntegerDigits[#, 2]]] &] (* _Alonso del Arte_, Dec 13 2010 *)

%o (Python)

%o import math

%o for n in range(1, 1000):

%o if not n % int(math.log(n, 2) + 1): print(n) # Garcia

%o (Haskell)

%o base_weight b g n | n == 0 = 0 | otherwise = (base_weight b g (n `div` b)) + (g $ n `mod` b)

%o interesting b g = filter f [1..] where f n = n `mod` (base_weight b g n) == 0

%o bin_interesting g = interesting 2 g

%o weights l n | (n >=0) && ((length l) > fromInteger n) = l !! fromInteger n | otherwise = 0

%o cnst = weights [1, 1]

%o let sequence = bin_interesting cnst -- Victor S. Miller, Oct 17 2011

%o (PARI) for(d=1, 9, forstep(n=(2^(d-1)+d-1)\d*d, 2^d-1, d, print1(n", "))) \\ _Charles R Greathouse IV_, Oct 17 2011

%Y Base 10 equivalent is A098952.

%K nonn,base

%O 1,2

%A _Grant Garcia_, Dec 13 2010