RefineGCDFreeBasis(B: LIST of INT, N: INT): LIST of INT
This function computes a refined GCD free basis by adjoining a given
integer to it. The value returned is [NewB, N2] where NewB is the
refined basis and N2 is the part of N coprime to every element of B.