View Single Post
Old 2022-02-08, 16:08   #11
JCoveiro
 
"Jorge Coveiro"
Nov 2006
Moura, Portugal

24·3 Posts
Default what if?

t(q)={w=(2^q+1)/3;s=4;for(i=1,q-1,s=(s^2-2)%(w));s1=lift(Mod(s-5-9,w));s2=lift(Mod(s-5+9,w));if((s1==0)||(s2==0),print1(q","))}
forprime(x=3,5000,t(x))

Code:
3,5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701,1709,2617,3539,
JCoveiro is offline   Reply With Quote