up previous next
determine (minimal) GCD free basis of a set of integers
GCDFreeBasis(L: LIST of INT): LIST of INT
This function returns a GCD free basis of a set of integers; you can
think of this as the set of all numbers (except 1) obtainable by
performing GCD and exact division operations.
Given a set N = [ N_1,...N_k ]
we seek a basis G = [ G_1,...,G_s ]
is a product of powers of the G_j
, and the G_j
coprime; the set G
is called a GCD free basis for N
. In general the set
is not uniquely defined.
/**/ GCDFreeBasis([factorial(20), factorial(10)]);
[46189, 4, 14175]