OFFSET
1,1
COMMENTS
The first digit of each term must be 2, 6, or 8. - Harvey P. Dale, Jun 13 2020
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..40
FORMULA
a(x) = 10^(x+1)*n+floor(10^x*3/9)*10+n+1. Output if a(x) is prime.
MATHEMATICA
Select[Sort[Flatten[Table[10FromDigits[PadRight[{i}, n, 3]]+i+1, {i, {2, 6, 8}}, {n, 2, 25}]]], PrimeQ] (* Harvey P. Dale, Jun 13 2020 *)
PROG
(PARI) n10np1(n) = { local(x, y, k); for(x=1, n, for(k=1, 8, y=10^(x+1)*k+floor(10^x*3/9)*10+k+1; if(isprime(y), print1(y", ")) ) ) }
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Cino Hilliard, Jul 11 2005
EXTENSIONS
Corrected by T. D. Noe, Nov 15 2006
STATUS
approved