Compute a factorization (of a polynomial) into coprime squarefree factors.
The factorization may sometimes be finer than necessary, i.e. two factors
could have the same multiplicity.
Example
/**/ Use R ::= QQ[x,y];
/**/ f := (x^2-1)^2*(y+2)^3;
/**/ indent(SqFreeFactor(f));
record[
RemainingFactor := 1,
factors := [x^2 -1, y +2],
multiplicities := [2, 3]
]