mersenneforum.org - View Single Post

JCoveiro · 2022-02-08, 16:08

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,

2022-02-08, 16:08	#11
JCoveiro "Jorge Coveiro" Nov 2006 Moura, Portugal 2⁴×3 Posts	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,