OFFSET
1,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..10000
EXAMPLE
Table showing n, A033203(n), a(n), a(n)^2 + 2:
1 2 0 2
2 3 1 3
3 11 3 11
4 17 7 51 = 17*3
5 19 6 38 = 19*2
6 41 11 123 = 41*3
7 43 16 258 = 43*3*2
8 59 23 531 = 59*3^2
9 67 20 402 = 67*3*2
10 73 12 146 = 73*2
...
PROG
(PARI) lista(n) = { my(L=List(), p=0); while(#L<n, p=nextprime(p+1); my(r=p%8); if(r<>5&&r<>7, my(k=0); while(vecmax(factor(k^2 + 2)[, 1]) <> p, k++); listput(L, k) )); Vec(L) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Dec 22 2024
STATUS
approved
