login
Primes p that have Kronecker symbol (p|86) = 1.
1

%I #13 Sep 08 2022 08:45:57

%S 3,5,17,19,23,29,31,37,41,47,61,79,97,103,127,131,149,157,163,167,179,

%T 193,211,227,239,271,277,281,311,331,337,347,349,353,359,367,373,389,

%U 401,419,421,431,439,467,479,487,491,499,523,569,571,587,599,617,653

%N Primes p that have Kronecker symbol (p|86) = 1.

%C Originally incorrectly named "primes which are squares mod 86", which is sequence A106891. - _M. F. Hasler_, Jan 15 2016

%H Vincenzo Librandi, <a href="/A191051/b191051.txt">Table of n, a(n) for n = 1..1000</a>

%t Select[Prime[Range[200]], JacobiSymbol[#,86]==1&]

%o (Magma) [p: p in PrimesUpTo(653) | KroneckerSymbol(p, 86) eq 1]; // _Vincenzo Librandi_, Sep 11 2012

%o (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

%Y Cf. A191017, A191018, A191020, A191023, A191025, A191026, A191028, A191029, A191032, A191034, A191036, A191037, A191040, A191042, A191043, A191046, A191047, A191049.

%K nonn,easy

%O 1,1

%A _T. D. Noe_, May 25 2011

%E Definition corrected by _M. F. Hasler_, Jan 15 2016