Programmation GP
forvec
Avec forvec :
s3(n)=
{
my(m=sqrtint(n));
forvec(v=vector(3,i,[1,m]),
if (v*v~==n,
return(v)));
}
? s3(12345)
%2 = [4,77,80]
Pour un meilleur algorithme, utiliser qfsolve.