OFFSET
1,1
COMMENTS
Only 1 (mod 4) numbers have been found.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
t[n_] := Which[Mod[n, 4] == 3, n + 1, Mod[n, 4] == 1, n - 1, True, n]; Select[1 + Range[99000], PowerMod[2, # - 1, #] == 1 && !PrimeQ[#] && Im[PowerMod[1 + 2I, t[#], #]] == 0 &]
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, May 22 2012
STATUS
approved