Thread: 1+1 Selfridges PRP test View Single Post
2020-11-20, 06:32   #19
Kebbaj

"Kebbaj Reda"
May 2018
Casablanca, Morocco

7×11 Posts

Quote:
 Originally Posted by paulunderwood You could have done n%100==63 instead. Anyway thanks. With the counterexample [62611454518799, 9] I need to go back to the drawing board. I had a PM for Dr. Sardonicus, wherein he describes a method to generate counterexamples for a=4. Perhaps he will post his method here.
Thank you also to have drive me into this question of this exceptions. And also to put me on the language PARI, usually I program on MMA.
At the beginning I wanted to translate your code into MMA to better measure it, but I couldn't translate the "kronecker" function
(if someone knows how to translate the "kronecker" function into MMA?).

Occasionally, I wanted to make a calculation status in your code, to know where "n" is in the calculation, but I didn't succeed:

Code:
{
forstep(a=3,3,1,print(a);D=a^2-4;E=D^2-4;s=187655186463;
forstep(n=177655186463,1840473667452,4*E,
if(n>=s,s=s+10000000000;print([n,s]),
if((n-floor(n%)*100)==63,
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])))))))
}
I don't know why its plant?

Last fiddled with by Kebbaj on 2020-11-20 at 06:39