OFFSET
1,2
PROG
(PARI) f(n)=my(k=n+(sqrtint(4*n)+1)\2); forprime(p=2, , if(!issquare(Mod(k, p)), return(p))); \\ A144294
lista(nn) = {my(v=vector(nn, n, f(n))); my(m=0, nm=0); for (n=1, nn, nm = v[n]; if (nm > m, print1(n, ", "); m = nm; ); ); } \\ Michel Marcus, Jun 25 2021
CROSSREFS
KEYWORD
nonn,more
AUTHOR
N. J. A. Sloane, Dec 03 2008
EXTENSIONS
a(9)-a(12) from R. J. Mathar, Dec 04 2008
a(13)-a(15) from Michel Marcus, Jun 25 2021
a(16)-a(18) from Michael S. Branicky, Sep 30 2024
STATUS
approved