%I #19 Sep 08 2022 08:45:57
%S 2,5,7,19,41,47,59,67,71,97,101,103,107,109,113,131,149,157,163,173,
%T 191,193,211,227,233,257,281,283,293,307,311,317,349,359,373,379,397,
%U 419,421,431,439,443,467,479,503,521,541,547,563,577,593,599,607,617
%N Primes that are squares mod 31.
%H Vincenzo Librandi, <a href="/A191024/b191024.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[150]],JacobiSymbol[#,31]==1&]
%o (Magma) [p: p in PrimesUpTo(617) | JacobiSymbol(p, 31) eq 1]; // _Vincenzo Librandi_, Sep 10 2012
%K nonn,easy
%O 1,1
%A _T. D. Noe_, May 24 2011