OFFSET
1,2
COMMENTS
Primes: 5, 58061, 341455271,...
MATHEMATICA
Select[Range[10^5], Mod[PowerMod[2, #, 4*#-1] + #^2, 4*#-1] == 0 &] (* Giovanni Resta, Oct 09 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Sep 24 2019
EXTENSIONS
a(10)-a(27) from Giovanni Resta, Oct 09 2019
STATUS
approved