%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