OFFSET
1,1
COMMENTS
a(n-1) is a term in row a(n) of A046071. - Reinhard Zumkeller, May 10 2015
LINKS
T. D. Noe, Table of n, a(n) for n=1..1000
MAPLE
f:= proc(t) local i, p;
p:= t;
do
p:= nextprime(p);
if numtheory:-jacobi(t, p) = 1 then return p fi
od
end proc:
A[1]:= 2:
for n from 2 to 100 do A[n]:= f(A[n-1]) od:
seq(A[i], i=1..100); # Robert Israel, Sep 04 2016
MATHEMATICA
f[n_] := Block[{k = PrimePi[n] + 1}, While[ JacobiSymbol[n, Prime[k]] == -1, k++ ]; Prime[k]]; NestList[f, 2, 54] (* Robert G. Wilson v, Mar 16 2004 *)
PROG
(Haskell)
a034794 n = a034794_list !! (n-1)
a034794_list = 2 : f 2 (tail a000040_list) where
f x (p:ps) = if elem x $ a046071_row p then p : f p ps else f x ps
-- Reinhard Zumkeller, May 10 2015
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
EXTENSIONS
Mathematica updated by Jean-François Alcover, Jul 04 2013
Name corrected by Robert Israel, Sep 04 2016
STATUS
approved