OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..200
FORMULA
MAPLE
isA002113 := proc(n) local dgs, i ; dgs := convert(n, base, 10) ; for i from 1 to nops(dgs)/2 do if op(i, dgs) <> op(-i, dgs) then RETURN(false): fi; od: RETURN(true) ; end: ndigs := proc(n) max(1, ilog10(n)+1) ; end: A031346 := proc(n) local p, nred; p := 0 ; nred := n ; while ndigs(nred) > 1 do nred := mul(d, d=convert(nred, base, 10) ) ; p := p+1 ; od: p ; end: isA159613 := proc(n) RETURN(isprime(n) and isA002113(n) and A031346(n) = 1) ; end: for n from 1 to 100000 do p := ithprime(n) ; if isA159613(p) then printf("%d, ", p) ; fi; od: # R. J. Mathar, Apr 21 2009
MATHEMATICA
palpQ[n_]:=Module[{idn=IntegerDigits[n]}, idn==Reverse[idn] && Length[ NestWhileList[Times@@IntegerDigits[#]&, n, #>9&]]==2]; Select[Prime[ Range[82000]], palpQ] (* Harvey P. Dale, Dec 16 2011 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Lekraj Beedassy, Apr 17 2009
EXTENSIONS
Extended by R. J. Mathar, Apr 21 2009
STATUS
approved