OFFSET
1,1
EXAMPLE
13 is in the sequence because 5i, 8i, 1 and 12 are solutions of x^2 = 1 (mod 13).
MATHEMATICA
h[n_] := Flatten[Table[a + b I, {a, 0, n - 1}, {b, 0, n - 1}]]; sol[n_] := Select[h[n], Mod[#^2, n] == 1 &]; Select[Range[100], Length[sol[#]] > 2 &]
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Jun 16 2013
STATUS
approved