|
|
A108841
|
|
Palindromic primes in which all internal digits are 4.
|
|
8
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(10) has 1255 digits, and thus is too big for a b-file. - Robert Israel, Oct 25 2018
|
|
LINKS
|
|
|
MAPLE
|
select(isprime, [seq(seq(x*(10^d+1)+40*(10^(d-1)-1)/9, x=1..9, 2), d=2..66)]); # Robert Israel, Oct 25 2018
|
|
MATHEMATICA
|
Sort[Select[FromDigits/@Flatten[Table[Join[{n}, PadRight[{}, i, 4], {n}], {n, {1, 3, 7, 9}}, {i, 150}], 1], PrimeQ]] (* Harvey P. Dale, Nov 26 2014 *)
|
|
PROG
|
(PARI) n10np1(n, d) = { local(x, y, k); for(x=1, n, for(k=1, 8, y=10^(x+1)*k+floor(10^x*d/9)*10+k; if(isprime(y), print1(y", ")) ) ) }
(Python)
from sympy import isprime
from itertools import count, islice
def agen(): yield from (t for i in count(1) for f in "1379" if isprime(t:=int(f + "4"*i + f)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|