OFFSET
1,1
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..370 (terms 1..210 from Robert Israel)
MAPLE
A[1]:= 3;
for n from 2 to 100 do
d:= 10^(ilog10(A[n-1])+1);
for k from 1 do
p:= A[n-1]+d*k;
if isprime(p) then
A[n]:= p;
break
fi
od
od:
seq(A[n], n=1..100); # Robert Israel, Jul 15 2014
PROG
(Python)
from sympy import isprime
from itertools import count, islice
def agen(an=3):
while True:
yield an
pow10 = 10**len(str(an))
for t in count(pow10+an, step=pow10):
if isprime(t):
an = t
break
print(list(islice(agen(), 17))) # Michael S. Branicky, Jun 23 2022
CROSSREFS
KEYWORD
base,nonn
AUTHOR
G. L. Honaker, Jr., Jan 18 2000
EXTENSIONS
Definition corrected by Robert Israel, Jul 15 2014
STATUS
approved