This function is provided by the package
LLLBases.
The first n-1 columns of the matrix z form a basis of the kernel of the n integers of the list s, and the dot product of the last column of z and s is the gcd g.
The method used is described in the paper:
Havas, Majewski, Matthews,
Extended GCD and Hermite Normal Form Algorithms via Lattice Basis Reduction, Experimental Mathematics 7:2 p. 125 (1998).
For an example,
i1 : s = apply(5,i->372*(random 1000000))
o1 = {330441276, 107988252, 326524488, 42328020, 219662280}
o1 : List
|
i2 : (g,z) = gcdLLL s
o2 = (372, | 2 7 -1 -44 7 |)
| -14 -18 -9 -15 -6 |
| 7 -7 13 18 -11 |
| 18 -17 -28 -1 4 |
| -10 12 -8 47 8 |
o2 : Sequence
|
i3 : matrix{s} * z
o3 = | 0 0 0 0 372 |
1 5
o3 : Matrix ZZ <--- ZZ
|