OFFSET
1,2
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..10000
EXAMPLE
a(2) = 3, since the smallest number greater than all the previous terms which gives a prime when added to 2^2 is 3.
MAPLE
R:= 1: t:= 1:
for n from 2 to 100 do
t:= nextprime(t+n^2)-n^2;
R:= R, t
od:
R; # Robert Israel, Apr 11 2023
MATHEMATICA
a[n_] := a[n] = Module[{k = a[n - 1] + 1}, While[! PrimeQ[n^2 + k], k++]; k]; a[0] = 0; Array[a, 100] (* Amiram Eldar, Apr 12 2023 *)
PROG
(PARI) seq(n)={my(a=vector(n), p=0); for(n=1, #a, p++; while(!isprime(p+n^2), p++); a[n]=p); a} \\ Andrew Howroyd, Apr 11 2023
(Python)
from sympy import nextprime
from itertools import count, islice
def agen(): # generator of terms
an = 1
for n in count(2):
yield an
an = nextprime(an + n**2) - n**2
print(list(islice(agen(), 62))) # Michael S. Branicky, Apr 16 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Angad Singh, Apr 11 2023
STATUS
approved