OFFSET
1,1
PROG
(PARI) listar(nn) = {my(va = vector(nn), rec = 0); va[1] = 5; for (n=2, nn, p = factor(va[n-1]^2-4)[1, 1]; if (p > rec, print1(p, ", "); rec = p); va[n] = va[n-1] + p - 1; ); } \\ Michel Marcus, Dec 14 2018
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Vladimir Shevelev, May 15 2010
EXTENSIONS
Sequence extended by R. J. Mathar, Jun 30 2010
More terms from Michel Marcus, Dec 14 2018
STATUS
approved