  We show a method for constructing a polynomial interpolating roots'
multiplicities of another polynomial, that does not use companion matrices.
This leads to a modification to Guersenzvaig--Szechtman square-free
decomposition algorithm that is more efficient both in theory and in practice.
