|
%I #32 Sep 06 2021 21:34:24
%S 3,5,6,7,10,11,12,13,14,17,20,22,23,24,25,26,28,29,31,34,37,40,41,43,
%T 44,46,47,48,50,52,53,55,56,58,61,62,67,68,71,73,74,77,80,82,83,86,88,
%U 91,92,94,96,97,100,101,104,106,107,110,112,113,115,116,121
%N Numbers whose binary reversal is prime.
%C Base-2 analog of A095179.
%C If k is a term, then 2*k is a term too. - _Michel Marcus_, Apr 19 2020
%H Michel Marcus, <a href="/A204232/b204232.txt">Table of n, a(n) for n = 1..1000</a>
%e 3, 5 and 7 are in the sequence because their binary reversal, equal to themselves, is prime.
%e a(3)=6 is in the sequence, because 6=110[2] (written in base 2), whose reversal 011[2]=3 is prime.
%e a(5)=11 is in the sequence, because 11=1011[2] (written in base 2), whose reversal 1101[2]=13 is prime.
%t Select[Range[170], PrimeQ[FromDigits[Reverse[IntegerDigits[#, 2]], 2]] &] (* _Alonso del Arte_, Jan 13 2012 *)
%o (PARI) for(n=1,1e2,isprime((t=binary(n))*vector(#t,i,1<<i)~\2) & print1(n","))
%o (Python)
%o from sympy import isprime
%o def ok(n): return isprime(int(bin(n)[2:][::-1], 2))
%o print(list(filter(ok, range(1, 122)))) # _Michael S. Branicky_, Sep 06 2021
%Y Cf. A030101, A095179.
%Y Positions of 2's in A227864.
%K nonn,base,easy
%O 1,1
%A _M. F. Hasler_, Jan 13 2012
|
