A191037 Primes p that have Jacobi symbol (p|58) = 1.

3, 7, 11, 19, 23, 37, 43, 61, 71, 101, 103, 131, 151, 157, 163, 167, 199, 211, 223, 229, 233, 239, 241, 251, 257, 269, 281, 293, 307, 313, 317, 331, 353, 379, 383, 389, 401, 421, 431, 439, 443, 457, 461, 463, 467, 487, 491, 521, 541, 563, 593, 619, 631, 647

Offset: 1

Author

T. D. Noe, May 25 2011

Keywords

• nonn
• easy

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

Crossrefs

Cf. A191017, A191018, A191020, A191023, A191025, A191026, A191028, A191029, A191032, A191034, A191036.

Programs

Magma

[p: p in PrimesUpTo(647) | KroneckerSymbol(p, 58) eq 1]; // Vincenzo Librandi, Sep 11 2012

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&]

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

Extensions

Definition corrected (following an observation by David Broadhurst) by M. F. Hasler, Jan 15 2016