M2:  p=3	s[p-3]==  +  2^((p+1)/2) (mod Mp)
M3:  p=5	s[p-3]==  +  2^((p+1)/2) (mod Mp)
M4:  p=7	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M5:  p=13	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M6:  p=17	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M7:  p=19	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M8:  p=31	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M9:  p=61	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M10: p=89	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M11: p=107	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M12: p=127	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M13: p=521	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M14: p=607	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M15: p=1279	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M16: p=2203	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M17: p=2281	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M18: p=3217	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M19: p=4253	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M20: p=4423	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M21: p=9689	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M22: p=9941	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M23: p=11213	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M24: p=19937	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M25: p=21701	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M26: p=23209	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M27: p=44497	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M28: p=86243	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M29: p=110503	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M30: p=132049	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M31: p=216091	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M32: p=756839	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M33: p=859433	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M34: p=1257787	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M35: p=1398269	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M36: p=2976221	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M37: p=3021377	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M38: p=6972593	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M39: p=13466917	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M40: p=20996011	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M41: p=24036583	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M42: p=25964951	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M43: p=30402457	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M44: p=32582657	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M45: p=37156667	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M46: p=42643801	s[p-3]==  -  2^((p+1)/2) (mod Mp).
M47: p=43112609	s[p-3]==  +  2^((p+1)/2) (mod Mp).
M48: p=57885161 s[p-3]==  -  2^((p+1)/2) (mod Mp).

24+ 23-