Programmation GP
Boucle while
Ajouter
rho(n)=
{
my(x=2,y=5);
while(gcd(y-x,n)==1,
x=(x^2+1)%n;
y=(y^2+1)%n; y=(y^2+1)%n
);
gcd(n,y-x);
}
et faire
\r
rho(2^64+1)
%1 = 274177