|
%I #4 Jan 15 2016 20:44:46
%S 2,5,7,19,31,41,47,59,67,71,97,101,103,107,109,113,131,149,157,163,
%T 173,191,193,211,227,233,257,281,283,293,307,311,317,349,359,373,379,
%U 397,419,421,431,439,443,467,479,503,521,541,547,563,577,593,599,607,617,653,659,661,683
%N Primes which are squares (mod 31).
%o (PARI) select(p->issquare(Mod(p,31))&&isprime(p),[1..1000]) \\ This is to provide a generic characteristic function as 1st arg of select(), there are other ways to produce the sequence more efficiently. - _M. F. Hasler_, Jan 15 2016
%Y Cf. A191040.
%K nonn
%O 1,1
%A _M. F. Hasler_, Jan 15 2016
|
