OFFSET
1,1
COMMENTS
Motivated by the former (incorrect) definition of A191029.
Also, primes p which have Legendre symbols (p|3) = (p|13) = 1, together with 3 and 13.
Apparently this contains the 3 plus the elements of A139494. - R. J. Mathar, May 28 2025
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..10000
MATHEMATICA
Join[{3, 13}, Select[Prime[Range[500]], JacobiSymbol[#, {3, 13}] == {1, 1} &]] (* Paolo Xausa, May 29 2025 *)
PROG
(PARI) select(p->issquare(Mod(p, 39))&&isprime(p), [1..1000])
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
M. F. Hasler, Jan 15 2016
STATUS
approved
