OFFSET
1,1
COMMENTS
How can one prove that the sequence is infinite?
Probabilistic arguments imply that the number of terms not exceeding x is not less than 5/9*log(10)/log(6)*x^(log(6)/log(10))/log(x) = 0.7139...*x^0.778.../log(x).
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[Prime@ Range@ 310, Total@ Drop[First /@ Partition[DigitCount@ #, 2, 2], {3}] == 1 &] (* Michael De Vlieger, Sep 21 2015 *)
PROG
(PARI) nbd(vd, d) = #select(i->(i == d), vd);
lista(nn) = {forprime(p=2, nn, vd = digits(p); if (nbd(vd, 1) + nbd(vd, 3) + nbd(vd, 7) + nbd(vd, 9) == 1, print1(p, ", ")); ); } \\ Michel Marcus, Sep 22 2015
(PARI) list(lim)=my(v=List([3])); forprime(p=7, lim, if(#setintersect(Set(digits(p\10)), [1, 3, 7, 9])==0, listput(v, p))); Vec(v) \\ Charles R Greathouse IV, Sep 22 2015
(Perl) use ntheory ":all"; say join ", ", grep { tr/1379// == 1 } @{primes(3000)}; # Dana Jacobsen, Oct 13 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Vladimir Shevelev, Sep 21 2015
STATUS
approved