%N A positive integer k is included if the value of (the reversal of k's representation in binary) divides k.
%C By "reversal" of k's representation in binary, it is meant: write k in binary, reverse the order of its digits, and read the result as a binary value.
%C It seems (verified for the first 120000 entries) that a(n) = A057890(n+1). - _R. J. Mathar_, Jun 18 2009
%C Indeed, this is A057890 (palindromes with optional trailing zeros) without the initial term. In other bases this does not have to be so, as illustrated by A071687. - _Ivan Neretin_, Sep 04 2015
%e 24 in binary is 11000. The reversal of this is 11 (ignoring leading 0's), which is 3 in decimal. Since 3 divides 24, then 24 is included in this sequence.
%p A030101 := proc(n) local bdgs ; bdgs := convert(n,base,2) ; add( op(-i,bdgs)*2^(i-1),i=1..nops(bdgs)) ; end: isA161604 := proc(n) if ( n mod A030101(n) ) = 0 then true ; else false; fi; end: for n from 1 to 600 do if isA161604(n) then printf("%d,",n) ; fi; od: # _R. J. Mathar_, Jun 18 2009
%t Select[Range@ 153, Divisible[#, FromDigits[Reverse@ IntegerDigits[#, 2], 2]] &] (* _Michael De Vlieger_, Sep 04 2015 *)
%Y Cf. A030101, A057890, A071687.
%A _Leroy Quet_, Jun 14 2009
%E Extended by _R. J. Mathar_, Jun 18 2009