View Single Post
Old 2020-11-17, 19:07   #1
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

2·3·599 Posts
1+1 Selfridges PRP test

This script outputs no counterexamples:

Code:
{
forstep(a=3,100,1,print(a);
D=a^2-4;E=D^2-4;
forstep(n=4*E-1,50000000*E,4*E,
if(kronecker(D,n)==-1&&!ispseudoprime(n),
r=Mod(D,n)^((n-1)/2);s=Mod(E,n)^((n-1)/2);
if(r==-1&&s==1,
print([n,r,s])))))
}
Can you find one?
paulunderwood is offline   Reply With Quote