OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
MAPLE
revdigs := proc(n)
local L, nL, j;
L:= convert(n, base, 10);
nL:= nops(L);
add(L[i]*10^(nL-i), i=1..nL);
end:
select(isprime, [3, seq(seq(n*10^(d+1)+3*10^d + revdigs(n), n=10^(d-1) .. 10^d-1), d= 1..4)]); # Robert Israel, Nov 11 2015
MATHEMATICA
ppmd3Q[n_]:=Module[{idn=IntegerDigits[n], len}, len=Length[idn]; OddQ[len] && idn==Reverse[idn]&&idn[[(len+1)/2]]==3]; Select[Prime[ Range[ 120000]], ppmd3Q] (* Harvey P. Dale, Feb 12 2015 *)
PROG
(Magma) [ p: p in PrimesUpTo(200000000) | IsOdd(d) and D[(d+1) div 2] eq 3 and D eq Reverse(D) where d is #D where D is Intseq(p) ]; // Vincenzo Librandi, Apr 12 2011
(Python)
from gmpy2 import is_prime
A082439_list = [3]
for i in range(1, 10**6):
s = str(i)
n = int(s+'3'+s[::-1])
if is_prime(n):
A082439_list.append(n) # Chai Wah Wu, Nov 11 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Lekraj Beedassy, Apr 25 2003
EXTENSIONS
Added a(4), a(13), a(14), a(28) by Vincenzo Librandi, Apr 12 2011
STATUS
approved