View Single Post
Old 2020-11-20, 11:01   #21
Kebbaj
 
Kebbaj's Avatar
 
"Kebbaj Reda"
May 2018
Casablanca, Morocco

10011012 Posts
Default

OK that works!
I find that the functionality of PARI GP looks a lot of thinks like MMA.
Code:
{
str = fileopen("C:\\pari GP\\here.log", "a");
forstep(a=3,3,1,print(a);D=a^2-4;E=D^2-4;bound=177665186463;
forstep(n=177655186463,1840473667452,4*E,if((n-floor(n/100)*100)==63,
if(n>=bound,bound=100000000+bound;print([n,bound]),
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]);
filewrite(str, ([n,r,s]))))))));fileclose(str)
}
Thanks

I will try to find for a = 7

Last fiddled with by Kebbaj on 2020-11-20 at 11:13
Kebbaj is offline   Reply With Quote