OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 2000 terms from Zak Seidov)
FORMULA
a(n) = Sum_{i=1..n+1} A139033(i). - Alois P. Heinz, Jan 27 2023
MAPLE
s:= proc(n) option remember; `if`(n<1, 0, b(n)^2+s(n-1)) end:
b:= proc(n) option remember; local k, m;
k:= s(n-1); for m from 1+b(n-1)
while not isprime(k+m^2) do od; m
end: b(1):=1:
a:= n-> s(n+1):
seq(a(n), n=1..34); # Alois P. Heinz, Jan 27 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Apr 07 2008
EXTENSIONS
Terms corrected by Alois P. Heinz, Jan 27 2023
STATUS
approved