OFFSET
1,2
COMMENTS
a(22) > 1.3*10^5. All numbers up to this bound were sieved using newpgen and sr1sieve. Remaining numbers were checked for primality using Jean Penné's LLR application (BLS (N-1/N+1) test).
EXAMPLE
1 is a term because 1234567891 is prime.
2 is not a term because 12345678901 is composite (it is divisible by 857).
MATHEMATICA
Select[Range@ 1400, PrimeQ[123456789*10^# + 1] &] (* Michael De Vlieger, Jan 04 2019 *)
PROG
(Python)
from sympy.ntheory.primetest import isprime
for n in range(1, 1000):
if isprime(123456789*10**n+1):
print(n, end=', ') # Stefano Spezia, Jan 05 2019
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Matthias Baur, Jan 01 2019
STATUS
approved