OFFSET
1,1
COMMENTS
Originally incorrectly named "primes which are squares mod 86", which is sequence A106891. - M. F. Hasler, Jan 15 2016
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[200]], JacobiSymbol[#, 86]==1&]
PROG
(Magma) [p: p in PrimesUpTo(653) | KroneckerSymbol(p, 86) eq 1]; // Vincenzo Librandi, Sep 11 2012
(PARI) select(p->kronecker(p, 86)==1&&isprime(p), [1..1000]) \\ This is to provide a generic characteristic function ("is_A191051") as 1st arg of select(), there are other ways to produce the sequence more efficiently. - M. F. Hasler, Jan 15 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
T. D. Noe, May 25 2011
EXTENSIONS
Definition corrected by M. F. Hasler, Jan 15 2016
STATUS
approved