A000037 := proc(n) n +floor(1/2+sqrt(n)) ; end: A092419 :=proc(n) local i, k; k := A000037(n) ; for i from 1 do if numtheory[legendre](k, ithprime(i)) = -1 then RETURN(ithprime(i)); fi; od; end: rec := -1 ; for n from 1 do a := A092419(n) ; if a > rec then printf("%d, \n", a) ; rec := a; fi; od: # R. J. Mathar, Jul 13 2009