login
Numbers k which are emirps in more bases 1 < b < k than any previous number.
2

%I #17 Jan 25 2020 09:07:55

%S 2,5,7,11,13,17,23,29,31,41,43,53,67,73,79,97,113,127,157,193,223,269,

%T 277,311,379,421,431,467,487,563,613,647,743,907,937,977,1093,1193,

%U 1249,1259,1373,1483,1543,1637,1667,1933,2239,2393,2477,2521,2857,2957,3083

%N Numbers k which are emirps in more bases 1 < b < k than any previous number.

%C The corresponding numbers of bases are 0, 1, 3, 6, 8, 9, 12, 13, 17, 21, 24, ... (see the link for more values).

%H Amiram Eldar, <a href="/A331486/b331486.txt">Table of n, a(n) for n = 1..175</a>

%H Amiram Eldar, <a href="/A331486/a331486.txt">Table of n, a(n), number of bases for n = 1..175</a>

%e 2 is not emirp in any base.

%e 5 is emirp in one base, 3: 5 is 12 in base 3, and 21 in base 3 is 7 which is also a prime.

%e 7 is emirp in 3 bases, 3, 4, and 5.

%t emirpQ[n_, b_] := n != (rev = FromDigits[Reverse @ IntegerDigits[n, b], b]) && And @@ PrimeQ[{n, rev}];

%t emirpCount[n_] := Length @ Select[Range[2, n - 1], emirpQ[n, #] &];

%t seq = {}; emax = -1; Do[e1 = emirpCount[n]; If[e1 > emax, emax = e1; AppendTo[seq, n]], {n, 2, 3000}]; seq

%Y Cf. A006567, A080790, A087911, A107129, A228768.

%K nonn,base

%O 1,1

%A _Amiram Eldar_, Jan 23 2020