Forum: Wagstaff PRP Search
2022-02-09, 01:54
|
Replies: 12
Views: 3,877
more simplifications (2)
t3(q)={w=(2^q+1)/3;s=4;for(i=2,q,s=(s^2-2)%(6*w));if(ispowerful(s+2),print1(q","))}
forprime(x=1,1000,t3(x))
3,5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701,
or
...
|
Forum: Wagstaff PRP Search
2022-02-09, 01:28
|
Replies: 12
Views: 3,877
more wagstaff simplifications
t3(q)={w=(2^q+1)/3;s=4;for(i=1,q,s=(s^2-2)%(6*w););if(s==14||s==194,print1(q","))}
forprime(x=1,100000,t3(x))
3,5,7,11,13,17,19,23,31,43,61,79,101,127,167,191,199,313,347,701,1709,2617,3539,
|
Forum: Wagstaff PRP Search
2022-02-08, 16:08
|
Replies: 12
Views: 3,877
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))
...
|
Forum: Wagstaff PRP Search
2022-02-08, 16:04
|
Replies: 12
Views: 3,877
more fixes in wagstaff
t3(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,t3(x))
...
|
Forum: Wagstaff PRP Search
2022-02-08, 15:23
|
Replies: 12
Views: 3,877
fix
For Mersenne:
t1(q)={w=(2^q-1);S0=4;S=S0;for(i=1,q-1,S=Mod(S^2-2,w));s1=lift(Mod(S-5-9,w));s2=lift(Mod(S-5+9,w));if(s2==2,print1(q","))}
forprime(x=3,5000,t1(x))
...
|