mersenneforum.org - View Single Post

JCoveiro · 2022-02-09, 22:26

t3(q)={s=4;for(i=2,q,s=(s^2-2)%(6*(2^q+1)/3));s=(s-14)%(2^(q-2)-2);if(s==0,print1(q","))}
forprime(x=3,5000,t3(x))

Code:

5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701,1709,2617,3539,

2022-02-09, 22:26	#9
JCoveiro "Jorge Coveiro" Nov 2006 Moura, Portugal 2⁴×3 Posts	t3(q)={s=4;for(i=2,q,s=(s^2-2)%(6*(2^q+1)/3));s=(s-14)%(2^(q-2)-2);if(s==0,print1(q","))} forprime(x=3,5000,t3(x)) Code: 5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701,1709,2617,3539,