OFFSET
1,1
COMMENTS
Originally incorrectly named "Primes which are squares mod 58", which is sequence A038901. - M. F. Hasler, Jan 15 2016
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MAPLE
select(t -> isprime(t) and numtheory:-jacobi(t, 58)=1, [seq(i, i=3..1000, 2)]); # Robert Israel, Jan 15 2016
MATHEMATICA
Select[Prime[Range[200]], JacobiSymbol[#, 58]==1&]
PROG
(Magma) [p: p in PrimesUpTo(647) | KroneckerSymbol(p, 58) eq 1]; // Vincenzo Librandi, Sep 11 2012
(PARI) select(p->kronecker(p, 58)==1&&isprime(p), [1..1000]) \\ This is to provide a generic characteristic function ("is_A191037") 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 (following an observation by David Broadhurst) by M. F. Hasler, Jan 15 2016
STATUS
approved