OFFSET
1,1
PROG
(PARI) select(p->issquare(Mod(p, 31))&&isprime(p), [1..1000]) \\ This is to provide a generic characteristic function as 1st arg of select(), there are other ways to produce the sequence more efficiently. - M. F. Hasler, Jan 15 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jan 15 2016
STATUS
approved