Reduce[Mod[x*x+3,2^37-1]==Mod[0,2^37-1]]