OFFSET
1,1
COMMENTS
Originally incorrectly named "Primes that are squares mod 91", which is sequence A106989, and a subsequence of this one, because these primes have (p|7) = (p|13) = 1. The present sequence includes also primes with (p|7) = (p|13) = -1. - David Broadhurst and M. F. Hasler, Jan 15 2016
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[200]], JacobiSymbol[#, 91]==1&]
PROG
(Magma) [p: p in PrimesUpTo(641) | JacobiSymbol(p, 91) eq 1]; // Vincenzo Librandi, Sep 10 2012
(PARI) is(p)=kronecker(p, 91)==1&&isprime(p) \\ M. F. Hasler, Jan 15 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
T. D. Noe, May 25 2011
EXTENSIONS
Definition corrected, following a suggestion from David Broadhurst, by M. F. Hasler, Jan 15 2016
STATUS
approved