OFFSET
1,1
EXAMPLE
a(1) = 5011 is a prime;
5011_2 = 1001110010011_10 is a prime;
5 + 0 + 1 + 1 = 7;
1 + 0 + 0 + 1 + 1 + 1 + 0 + 0 + 1 + 0 + 0 + 1 + 1 = 7; both the digit sums are equal.
MATHEMATICA
Select[Prime[Range[1000000]], PrimeQ[FromDigits[IntegerDigits[#, 2]]] && Plus @@ IntegerDigits[#] == Plus @@ IntegerDigits[FromDigits[IntegerDigits[#, 2]]] &]
PROG
(PARI) eva(n) = subst(Pol(n), x, 10)
is(n) = ispseudoprime(n) && ispseudoprime(eva(binary(n))) && sumdigits(n)==sumdigits(eva(binary(n))) \\ Felix Fröhlich, Jan 19 2017
CROSSREFS
KEYWORD
nonn,base
AUTHOR
K. D. Bajpai, Jan 19 2017
STATUS
approved