OFFSET
1,1
COMMENTS
5, 11 and all primes congruent to 1, 4, 9, 14, 16, 26, 31, 34, 36, or 49 (mod 55). - Robert Israel, Jan 15 2016
MAPLE
S55:= {seq(x^2 mod 55, x=0..27)}:
select(t -> member(t mod 55, S55), [seq(ithprime(i), i=1..1000)]); # Robert Israel, Jan 15 2016
MATHEMATICA
Join[{5, 11}, Select[Prime[Range[250]], MemberQ[{1, 4, 9, 14, 16, 26, 31, 34, 36, 49}, Mod[#, 55]]&]] (* Harvey P. Dale, Jan 17 2022 *)
PROG
(PARI) select(p->issquare(Mod(p, 55))&&isprime(p), [1..1500]) \\ It would be more efficient to select only among primes, replacing [1..1500] by primes([1, 1500]), in which case the isprime() condition can be omitted from the selection function. But we wanted to provide a universally valid characteristic function in the 1st argument of select(). - M. F. Hasler, Jan 15 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Jan 15 2016
STATUS
approved