OFFSET
1,1
COMMENTS
Primes with the same number of prime decimal digits as odd nonprime digits.
MAPLE
digs1or9 := proc(n) local a, c; a := 0 ; for c in convert(n, base, 10) do if c in {1, 9} then a := a+1 ; end if; end do; a; end proc: A109066 := proc(n) local a, d; p := ithprime(n) ; a := 0 ; for d in convert(p, base, 10) do if isprime(d) then a := a+1 ; end if; end do; a; end proc: for n from 1 to 400 do p := ithprime(n) ; if digs1or9(p) = A109066(n) then printf("%d, ", p) ; fi; od: # R. J. Mathar, Oct 22 2009
MATHEMATICA
ddQ[n_]:=Module[{idn=IntegerDigits[n]}, Count[idn, 1|9]==Count[idn, 2|3|5|7]]; Select[Prime[Range[250]], ddQ] (* Harvey P. Dale, Nov 29 2011 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Juri-Stepan Gerasimov, Jan 15 2009
STATUS
approved