OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
N:= 100: # for a(1)..a(N)
A:=Vector(N): A[1]:= 3: n:= 1:
for d from 3 while n < N do
for m from 1 to d-1 while n < N do
if igcd(m, d)=1 then
n:= n+1;
for k from ceil((A[n-1]+1 - m)/d) do
q:= d*k+m;
if isprime(q) then A[n]:= q; break fi
od
fi
od od:
convert(A, list);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 22 2021
STATUS
approved