Bill Allombert on Thu, 14 Dec 2023 11:16:43 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: foursquares-2.16.1.gp threesquares_fact


On Wed, Dec 13, 2023 at 10:34:43PM -0800, Thomas D. Dean wrote:
> What purpose does F serve in qfsolve([matdiagonal([1,1,1,-n]),F])?
> 
> ? threesquares_fact(n,F)
> %66 = [15, 0, 3]
> ? abs(qfsolve(matdiagonal([1,1,1,-n]))[1..3])~
> %67 = [15, 0, 3]

F is the factorization of discriminant of matdiagonal([1,1,1,-n]) which is n.

This is to avoid factoring n several time.

Cheers,
Bill