OFFSET
1,1
COMMENTS
31 has the same number of digits 1 and 3, so it is not in the sequence.
See comment in A245626.
LINKS
Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n)
local L;
if not isprime(n) then return false fi;
L:= convert(n, base, 10);
numboccur(1, L)>numboccur(3, L)
end proc:
select(filter, [2*i+1$i=1..1000]); # Robert Israel, Aug 05 2014
PROG
(PARI) nbd(n, d) = my(v=digits(n), m=0); for(i=1, #v, if(v[i]==d, m++)); m
select(n->nbd(n, 1)>nbd(n, 3), primes(300)) \\ Colin Barker, Jul 27 2014
(Python)
import sympy
from sympy import prime
for n in range(1, 10**3):
..s = str(prime(n))
..if s.count('1') > s.count('3'):
....print(int(s), end=', ')
# Derek Orr, Jul 28 2014
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Vladimir Shevelev, Jul 27 2014
EXTENSIONS
More terms from Peter J. C. Moses, Jul 27 2014
More terms from Colin Barker, Jul 27 2014
STATUS
approved