Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #10 Oct 26 2019 17:53:57
%S 20,53,341,536,2312,2477,3380,3665,3686,4547,5009,5105,6458,6488,6731,
%T 6845,7499,7508,7562,7835,8411,8831,9032,9386,9764,9839,11027,11885,
%U 14990,19589,20498,21080,22844,24821,25220,26393,27593,29864,29921
%N Numbers n such that 1n1, 3n3, 7n7 and 9n9 are all primes.
%e 2312 is a term because 123121, 323123, 723127 and 923129 are all primes.
%t Select[ Range[ 30000 ], PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "1", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "3", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "7", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "9", {1, -1} ] ] ] & ]
%t enclose[n_]:=Table[FromDigits[Join[{i},IntegerDigits[n],{i}]],{i, {1,3,7,9}}]; Select[Range[30000],AllTrue[enclose[#],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Mar 21 2015 *)
%Y Cf. A032682, A059693, A059694.
%K base,nonn
%O 1,1
%A _Harvey P. Dale_, Feb 05 2001
%E More terms from _Patrick De Geest_, Feb 07 2001