OFFSET
1,1
COMMENTS
See A090711 for a similar sequence whose definition works "in the opposite direction". - M. F. Hasler, Jan 03 2014
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
A000040(10)=29 in base 11 is 2*11^1+9*11^0=31 prime, therefore 29 is a term.
MAPLE
filter:= proc(n) local L;
if not isprime(n) then return false fi;
L:= convert(n, base, 10);
isprime(add(L[i]*11^(i-1), i=1..nops(L)))
end proc:
select(filter, [2, seq(i, i=3..10000, 2)]); # Robert Israel, Jan 28 2018
MATHEMATICA
Select[Prime@ Range@ 250, PrimeQ@ FromDigits[IntegerDigits@ #, 11] &] (* Michael De Vlieger, Aug 29 2015 *)
PROG
(PARI) is(p, b=11)={my(d=digits(p)); isprime(vector(#d, i, b^(#d-i))*d~)&&isprime(p)} \\ M. F. Hasler, Jan 03 2014
(Magma) [n:n in PrimesUpTo(1600)| IsPrime(Seqint(Intseq(n), 11))]; // Marius A. Burtea, Jun 30 2019
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Reinhard Zumkeller, Feb 13 2004
EXTENSIONS
Corrected by Zak Seidov, Feb 25 2004
STATUS
approved